arxiv: v3 [gr-qc] 29 Jun 2015

Transcription

1 QUANTITATIVE DECAY RATES FOR DISPERSIVE SOLUTIONS TO THE EINSTEIN-SCALAR FIELD SYSTEM IN SPHERICAL SYMMETRY JONATHAN LUK AND SUNG-JIN OH arxiv: v3 [gr-qc] 29 Jn 205 Abstract. In this paper, we stdy the ftre casally geodesically complete soltions of the spherically symmetric Einstein-scalar field system. Under the a priori assmption that the scalar field φ scatters locally in the scale-invariant bonded-variation BV norm, we prove that φ and its derivatives decay polynomially. Moreover, we show that the decay rates are sharp. In particlar, we obtain sharp qantitative decay for the class of global soltions with small BV norms constrcted by Christodolo. As a conseqence of or reslts, for every ftre casally geodesically complete soltion with sfficiently reglar initial data, we show the dichotomy that either the sharp power law tail holds or that the spacetime blows p at infinity in the sense that some scale invariant spacetime norms blow p.. Introdction In this paper, we stdy the qantitative long time dynamics for the spherically symmetric dispersive spacetimes satisfying the Einstein-scalar field eqations. More precisely, these are spherically symmetric soltions M, g, φ to the Einstein-scalar field system, where g is a Lorentzian metric and φ is a real valed fnction on a 3 + dimensional manifold M, sch that M, g is ftre casally geodesically complete and φ scatters locally in the scale-invariant bonded-variation BV norm. For these spacetimes, we establish a Price-law type decay for the scalar field φ, the Christoffel symbols associated to g and all of their derivatives. To obtain the decay reslts, we do not need to assme any smallness of the initial data. Moreover, we show that the decay rates in this paper are sharp. The spherically symmetric Einstein-scalar field system, being one of the simplest model of selfgravitating matter in this symmetry class, has been stdied extensively both nmerically and mathematically. In a seminal series of papers by Christodolo [4], [5], [6], [7], [8], he achieved a complete nderstanding of the singlarity strctre of spherically symmetric spacetime soltions to this system. The clmination of the reslts shows that generic spherically symmetric initial data with one asymptotically flat end give rise to a spacetime whose global geometry is either dispersive with a Penrose diagram represented by Figre or contains a black hole region BH which terminates in a spacelike crvatre singlarity S with a Penrose diagram represented by Figre 2. In particlar, in either of these generic scenarios, the spacetime possesses a complete nll infinity I + and ths obeys the weak cosmic censorship conjectre. Moreover, in either case the maximal Cachy development of the data is inextendible with a C 2 Lorentzian metric and therefore also verifies the strong cosmic censorship conjectre. We refer the readers to [5] for a comprehensive discssion on general singlarity strctres for spherically symmetric spacetimes. The remarkable resoltion of the cosmic censorship conjectres however gives very little information on the long time dynamics for these spacetimes except for the small data 2 case [6]. In particlar, not mch is known abot the asymptotic decay of the scalar field as measred by a far-away observer at nll infinity. In the dispersive case, Christodolo showed that the Bondi in the BV class, i.e., the initial data for vrφ has bonded variation. More precisely, Christodolo showed that the non-generic set of initial data has co-dimension at least two in the BV topology. 2 i.e., when the initial data is close to that of Minkowski space.

2 i + Γ I + Σ i 0 Figre. The dispersive case S i + b Γ BH Γ H + I + Σ i 0 Figre 2. The black hole case mass decays to zero along nll infinity withot an explicit decay rate. In the black hole case, he showed that the Bondi mass approaches the mass of the black hole, from which one can infer the non-qantitative decay for the scalar field along nll infinity [6]. The long time dynamics in the case where the spacetime settles to a black hole was sbseqently stdied in the seminal work 3 of Dafermos-Rodnianski [2]. They proved a qantitative decay rate for the scalar field and its derivatives in the spacetime inclding along nll infinity I + and the event horizon H +. The proof is based on the local conservation of energy, which is sbcritical, together with techniqes exploiting the conformal geometry of the spacetime and the celebrated red-shift effect along the event horizon. The reslt in particlar jstified, in a nonlinear setting, the heristics of Price [9]. It trns ot that the qantitative decay rates, when combined with the reslts of [], also have interesting conseqences for the strong cosmic censorship conjectre in the context of the spherically symmetric Einstein-Maxwell-scalar field system. In this paper, we stdy the other generic scenerio, i.e., spherically symmetric dispersive spacetime soltions to the Einstein-scalar field system. Unlike in the black hole case, the monotonic Hawking mass is spercritical and provides no control over the dynamics of the soltion. We ths do not expect to be able to obtain qantitative decay rates for large soltions withot imposing extra assmptions. Instead, we assme a priori the non-qantitative decay of a critical qantity - the 3 In fact, they considered the more general Einstein-Maxwell-scalar field eqations. 2

3 BV norm 4 - bt only locally in a region where the area of the orbit of the symmetry grop SO3 remains niformly bonded. Under this assmption of local BV scattering, we show that the scalar field and all its derivatives decay with a qantitative rate, reminescent of the Price law decay rates in the black hole case. We refer the readers to the statement of the main theorems in Section 3 for the precise rates that we obtain. We prove, in particlar, a qantitative decay rate for the scalar field along nll infinity. Or reslts apply in particlar to the class of soltions arising from initial data with small BV norm. Christodolo [6] showed that these spacetimes are ftre casally geodesically complete. One can easily dedce from [6] that in fact these spacetimes satisfy the BV scattering assmption and therefore the soltions obey the qantitative decay estimates of or main theorem see Theorem 3.5. On the other hand, or reslts do not reqire any smallness assmptions on the initial data. We conjectre that indeed or class of spacetimes contains those arising from large data: Conjectre. There exists initial data of arbitrarily large BV norm whose maximal global development scatters locally in the BV norm. In addition to the pper bonds that we obtain in or main theorem, we also constrct examples where we prove lower bonds for the soltions with the same rates as the pper bonds. In particlar, there exists a class of initial data with compactly spported scalar field whose ftre development satrates the decay estimates in the main theorem. This shows that the decay rates are sharp. We note that the decay rate is also consistent with the nmerical stdy of Bizón-Chmaj- Rostworowski [2]. As a corollary of the main reslt on decay, we show the following dichotomy: either the qantitative decay rates are satisfied or the soltion blows p at infinity. The latter are soltions sch that some scale-invariant spacetime norms become infinite see precise definition in Definition 3.2. The decay reslt in this paper easily implies that the locally BV scattering soltions that we consider are stable against small, reglar, spherically symmetric pertrbations. More ambitiosly, one may conjectre that Conjectre 2. Spherically symmetric locally BV scattering dispersive soltions to the Einsteinscalar field eqations are stable against non-spherically symmetric pertrbations. Conjectre 2, if tre, will generalize the monmental theorem on the nonlinear stability of Minkowski spacetime of Christodolo-Klainerman [9] see also a simpler proof in [6]. For nonlinear wave eqations satisfying the nll condition, it is known [], [22] that large soltions decaying sfficiently fast are globally stable against small pertrbations. On the other hand, or main theorem shows a qantitative decay rate for spherically symmetric locally BV scattering dispersive spacetimes. Conjectre 2 can therefore be viewed as an attempt to generalize the reslts in [], [22] to the Einstein-scalar field system. We will address both Conjectres and 2 in ftre works. Acknowledgements. The athors wold like to thank Mihalis Dafermos and Igor Rodnianski for valable discssions. We also thank Jonathan Kommemi for providing the Penrose diagrams. J. Lk is spported by the NSF Postdoctoral Fellowship DMS S.-J. Oh is a Miller Research Fellow, and thanks the Miller Institte at UC Berkeley for the spport... Otline of the paper. We otline the remainder of the paper. In Section 2, we discss the set-p of the problem and in particlar define the class of soltions considered in the main theorem, i.e., the locally BV scattering soltions. In Section 3, we state the main theorems in the paper Theorems 3. and 3.2, their conseqences and additional theorems on the optimality of the decay 4 Soltions of bonded variation have been first stdied by Christodolo [6] and plays an important role in the proof of the cosmic censorship conjectres [8]. 3

4 rates. In the same section, we otline the main ideas of the proof. In Sections 4 and 5, we explain the main analytic featres of the eqations and the geometry of the class of spacetimes that we consider. Sections 6 and 7 consist of the main content of this paper. In Section 6, we prove the decay estimates for φ, v rφ and rφ, i.e., the first main theorem Theorem 3.. In Section 7, sing in particlar the reslts in Section 6, we derive the decay bonds for the second derivatives for rφ and the metric components, i.e., the second main theorem Theorem 3.2. In the remaining sections of the paper, we trn to other theorems stated in Section 3. In Section 8, we give a proof of the dichotomy allded to above, i.e., either the qantitative decay holds or the spacetime blows p at infinity. In Section 9, we sketch a proof of a refinement of the conclsions of the main theorems in the small data case. Finally, in Section 0, we prove optimality of the decay rates asserted by the main theorems. 2. Set-p In this section, we define the set-p, formlate the eqations in a doble nll coordinate system and explain the characteristic initial vale problem. This will allow s to state the main theorem in the next section. 2.. Spherically Symmetric Einstein-Scalar-Field System SSESF. We begin with a brief discssion on the derivation of the Spherically Symmetric Einstein-Scalar-Field System SSESF from the 3 + -dimensional Einstein-scalar-field system. Soltions to the Einstein-scalar field eqations can be represented by a triplet M, g µν, φ, where M, g µν is a 3 + -dimensional Lorentzian manifold and φ a real-valed fnction on M. The spacetime metric g µν and the scalar field φ satisfy the Einstein-scalar-field system: R µν 2. 2 g µνr =2T µν, µ µ φ =0. where R µν is the Ricci crvatre of g µν, R is the scalar crvatre, and µ is the covariant derivative given by the Levi-Civita connection on M, g. The energy-momentm tensor T µν is given by the scalar field φ, i.e. 2.2 T µν = µ φ ν φ 2 g µν λ φ λ φ. Assme that the soltion M, g µν, φ is spherically symmetric, i.e., the grop SO3 of three dimensional rotations acts smoothly and isometrically on M, g, where each orbit is either a point or is isometric to S 2 with a rond metric. The scalar field φ is reqired to be constant on each of the orbits. These assmptions are propagated by 2.; hence, if M, g µν, φ is a Cachy development, then it sffices to assme spherical symmetry only on the initial data. The qotient M/SO3 gives rise to a +-dimensional Lorentzian manifold with bondary, which we denote by Q, g ab. The bondary Γ consists of fixed points of the grop action. We define the area radis fnction r on Q to be Area of symmetry sphere r :=. 4π and r = 0 at Γ. Note that each component of Γ is a timelike geodesic. We assme that Γ is non-empty and connected, and moreover that there exists a global doble nll coordinates, v, i.e. a coordinate system, v covering Q in which the metric takes the form 2.3 g ab dx a dx b = Ω 2 d dv 4

5 for some Ω > 0. We remark that both assmptions are easily jstified if M, g is a Cachy development of a spacelike hypersrface homeomorphic to R 3. The metric g µν of M is characterized by Ω and r and takes the form 2.4 g µν dx µ dx ν = Ω 2 d dv + r 2 ds 2 S 2 where ds 2 S is the standard line element on the nit sphere S 2. Therefore, we may reformlate the 2 spherically symmetric Einstein-scalar-field system SSESF in terms of the triplet φ, r, Ω as r v r = r v r 4 Ω2, SSESF r 2 v log Ω = r v r + 4 Ω2 r 2 φ v φ, r v φ = r v φ v r φ, 2Ω r Ω = r 2 + r φ 2, 2Ω v r v Ω = vr 2 + r v φ 2, with the bondary condition r = 0 along Γ Basic assmptions, notations and conventions. In this sbsection, we introdce the basic assmptions on the base manifold Q, as well as some notations and conventions that will be sed in the rest of the paper. Definition of Q and M. Denote by R + the +-dimensional Minkowski space, with the standard doble nll coordinates, v. Let Q be a +-dimensional Lorentzian manifold which is conformally embedded into R + with ds 2 Q = Ω2 d dv. Given a non-negative fnction r on Q, we define the set Γ := {, v Q : r, v = 0}, called the axis of symmetry. We also define M, g µν to be the +3-dimensional Lorentzian manifold with M = Q S 2 and g µν given by 2.3; this is to be thoght of as the fll spacetime before the symmetry redction. We refer to 2. for the fll interpretation. Assmptions on the conformal geometry of Q. We assme that Γ Q is a connected set, which is the image of a ftre-directed timelike crve emanating from the point,. We also assme that C Q, where C = {, v R + : =, v < }. Frthermore, Q is assmed to be the domain of dependence of Γ and C to the ftre, in the sense that every casal crve in Q has its past endpoint on either Γ or C. Notations for the conformal geometry of Q. Denote by C resp. C v the constant resp. v crve in Q. Note that these are nll crves in Q. Given 0, v 0 Q, we define the domain of dependence of the line segment C 0 {v v 0 }, denoted by D 0, v 0, to be the set of all points p Q sch that all past-directed casal crves passing p intersects Γ C 0 {v v 0 }, pls the line segment C 0 {v v 0 } itself. Also, we define the ftre nll infinity I + to be the set of ideal points, + sch that sp C r =. Integration over nll crves. Whenever we integrate over a sbset of C resp. C v, we will se the standard line element dv d for integrals over, i.e., 2 f = f, v d C v { 2 } 2 f = f, v dv, C {v v v 2 } v 5

6 respectively. Fnctions of bonded variation. Unless otherwise specified, fnctions of bonded variation BV considered in this paper will be assmed to be right-continos. By convention, v f dv or v f will refer to the distribtional derivative of f, which is a finite signed measre, and v f dv or v f will denote the total variation measre. Unless otherwise specified, these measres will be evalated on intervals of the form v, v 2 ]. Ths, according to or conventions, 2 v v fv dv =fv 2 fv, 2 v v fv dv =T.V. v,v 2 ][f]. New variables. We introdce the following notation for the directional derivatives of r: λ := r v, r ν :=, The Hawking mass m, v is defined by the relation 2.5 2m r = a r a r = 4Ω 2 r v r. For a soltion to SSESF, the qantity m possesses sefl monotonicity properties see Lemma 4.3, which will be one of the key ingredients of or analysis. We define the mass ratio to be µ := 2m r. We also define the Bondi mass on C by M := lim v m, v, provided the limit exists. The Bondi mass M i := M = lim v m, v on the initial crve C is called the initial Bondi mass Refomlation in terms of the Hawking mass. The Hawking mass as defined in 2.5 trns ot to obey sefl monotonicity See 4.2. We therefore reformlate SSESF in terms of m and eliminate Ω. Notice that by 2.4 and 2.5, the metric is determined by r and m. We say that φ, r, m on Q is a soltion to SSESF if the following eqations hold: µ λ = µr λν, 2.6 µ v ν = µr λν, 2.7 { 2ν m = µr 2 φ 2, 2λ v m = µr 2 v φ 2, 2.8 v rφ = µλν µr φ, and moreover, the following bondary conditions hold: r = 0 and m = 0 along Γ. 6

7 We remark that sing 2.6, the wave eqation 2.8 for φ may be rewritten in either of the following two eqivalent forms: v rφ = λφ, v rφ = v νφ Choice of coordinates. Note that Q is rled by the family of nll crves C. Since a nll crve C with cannot intersect C, its past endpoint mst be on Γ. Therefore, or assmptions so far impose the following conditions on the doble nll coordinates, v on Q: is constant on each ftre-directed nll crve emanating from Γ and v is constant on each conjgate nll crve. However, these conditions are insfficient to give a niqe choice of a coordinate system, as the system SSESF and assmptions so far are invariant nder the change of coordinates U, v V v, U = V =, for any strictly increasing fnctions U and V. To remove this ambigity, we fix the choice of the coordinate system, once and for all, as follows. We first fix v on C relating it with the fnction r. Specifically, we will reqire that v = 2r + on C, which in particlar implies that 2.9 λ, v = 2. Next, in order to fix, we prescribe sch that Γ = {, v : = v}. To do so, for every otgoing nll crve C in Q, follow the incoming nll crve to the past starting from C Γ ntil the point p where it intersects the initial crve C. We then define the -coordinate vale for C to be the v-coordinate vale for p. Under sch coordinate choice, D 0, v 0 may be expressed as D 0, v 0 = {, v Q : [ 0, v 0 ], v [, v 0 ]}. Moreover, if r and φ are sfficiently reglar fnctions on Q, then or coordinate choice leads to lim v + λ + ν, v = lim v λ + ν, v = 0 and lim v + v + rφ, v = lim v v + rφ, v = 0. These conditions will be incorporated into precise formlations of soltions to SSESF with limited reglarity in the following sbsection Characteristic initial vale problem. In this paper, we stdy the characteristic initial vale problem for SSESF with data prescribed on C, nder qite general assmptions on the reglarity. In this sbsection, we give precise formlations of initial data and soltions to SSESF to be considered in this paper. We begin with a discssion on the constraint imposed by SSESF more precisely, on initial data for φ, r, m. In fact, the constraint is very simple, thanks to the fact that they are prescribed on a characteristic i.e., nll crve C. Once we prescribe φ on C, the coordinate condition 2.9 dictates the initial vales of r, and the initial vales of m are then determined by the constraint 2.7 along the v direction, as well as the bondary condition m, = 0. In other words, initial data for φ, r, m possess only one degree of freedom, namely the prescription of a single real-valed fnction φ, v, or eqivalently, v rφ, v. Following Christodolo [6], we say that an initial data set for φ, r, m is of bonded variation BV if v rφ, is a right-continos BV fnction on [, with finite total variation on,. We also define the notion of soltion of bonded variation to SSESF as follows. Definition 2. Bonded variation soltions to SSESF. A soltion φ, r, m to SSESF is called a soltion of bonded variation on Q if on every compact domain of dependence D 0, v 0, the following conditions hold: sp D0,v 0 ν < and sp D0,v 0 λ <. 7

8 2 λ is BV on each C D 0, v 0 niformly in, and ν is BV on each C v D 0, v 0 niformly in v. 3 For each a with a, a Γ, lim ν + λa, a + ɛ = 0. ɛ 0+ 4 φ is an absoltely continos fnction on each C D 0, v 0 with total variation bonded niformly in, and also an absoltely continos fnction on each C v D 0, v 0 with total variation bonded niformly in v. 5 For each a with a, a Γ, lim sp ɛ 0 0<δ ɛ lim sp ɛ 0 0<δ ɛ T.V. {a δ} a δ,a [φ] = 0, T.V. a,a+δ {a+δ} [φ] = 0, lim sp ɛ 0 0<δ ɛ lim sp ɛ 0 0<δ ɛ T.V. a ɛ,a δ {a δ} [φ] = 0, T.V. {a+δ} a+δ,a+ɛ [φ] = 0. 6 v rφ is BV on each C D 0, v 0 niformly in, and rφ is BV on each C v D 0, v 0 niformly in v. 7 For each a with a, a Γ, lim ɛ 0+ v rφ + rφ a, a + ɛ = 0. We also consider more reglar data and soltions, as follows. We say that an initial data set for φ, r, m is C if v rφ, is C on [, with sp C 2 vrφ <. In the following definition, we define the corresponding notion of a C soltion to SSESF. Definition 2.2 C soltions to SSESF. A soltion φ, r, m to SSESF is called a C soltion on Q if the following conditions hold on every compact domain of dependence D 0, v 0 : sp D0,v 0 ν < and sp D0,v 0 λ <. 2 λ, ν are C on D 0, v 0. 3 For each a with a, a Γ, lim ν + λa, a + ɛ = lim ν + λa ɛ, a = 0. ɛ 0+ ɛ 0+ 4 v rφ and rφ are C on D 0, v 0. 5 For each a with a, a Γ, lim v rφ + rφ a, a + ɛ = lim v rφ + v rφ a ɛ, a = 0. ɛ 0+ ɛ 0+ Remark 2.3. By [6, Theorem 6.3], a BV initial data set leads to a niqe BV soltion to SSESF on {, v : + δ, v } for some δ > 0. If the initial data set is frthermore C, then it is not difficlt to see that the corresponding soltion is also C persistence of reglarity. In fact, this statement follows from the argments in Section 7 of this paper; see, in particlar, the proof of Lemma Local scattering in BV and asymptotic flatness. We are now ready to formlate the precise notion of locally BV scattering soltions to SSESF, which is the class of soltions that we consider in this paper. In particlar, for this class of soltions, we make a priori assmptions on its global geometry. Definition 2.4 Local scattering in BV. We say that a BV soltion φ, r, m to SSESF is locally scattering in the bonded variation norm BV, or a locally BV scattering soltion, if the following conditions hold: 8

9 Ftre completeness of radial nll geodesics: Every incoming nll geodesic in Q has its ftre endpoint on Γ, and every otgoing nll geodesic in Q is infinite towards the ftre in the affine parameter. Moreover, there exists a global system of nll coordinates, v and Q is given by 2.0 Q = {, v : [,, v [, }. 2 Vanishing final Bondi mass: The final Bondi mass vanishes, i.e., 2. M f := lim M = 0. 3 Scattering in BV in a compact r-region: There exists R > 0 sch that for Q cpt defined to be the region {, v Q : r, v R}, we have 2.2 vrφ 2 0, v log λ 0 C Q cpt C Q cpt as. Several remarks concerning Definition 2.4 are in order. Remark 2.5. In fact, the condition 2.0 is a conseqence of ftre completeness of radial nll geodesics and the preceding assmptions. To see this, first recall or assmption that C = {, v : =, v [, }. Hence from or choice of the coordinate and ftre completeness of incoming radial nll geodesics, it follows that the range of mst be [,. Frthermore, for each [,, the range of v on C is [, by ftre completeness of otgoing radial nll geodesics and Definition 2.. More precisely, ftre completeness of C implies that it can be contined past {} [, v 0 ] as long as 0 bonded on {} [, v 0 ] for every finite v 0. Ω2 dv <, and Definition 2. implies 5 that Ω 2 = 4λν µ indeed remains Remark 2.6. For more reglar e.g., C asymptotically flat soltions, the conditions and 2 in Definition 2.4 may be replaced by a single eqivalent condition, namely reqiring the fll spacetime M, g to be ftre casally geodesically complete as a +3-dimensional Lorentzian manifold. In particlar, 2 follows from the deep work [4] of Christodolo, in which it was proved that if M f > 0 then the space-time necessarily contains a black-hole and ths is not ftre casally geodesically complete. Remark 2.7. As we will see in the proof, there exists a niversal ɛ 0 sch that 3 in Definition 2.4 can be replaced by the weaker reqirement that there exists R > 0 and U > 0 sch that vrφ 2 ɛ 0, v log λ ɛ 0 C Q cpt C Q cpt for U. To simplify the exposition, we will omit the proof of this improvement. Remark 2.8. For a sfficiently reglar, asymptotically flat soltion to SSESF, Condition in Definition 2.4 is eqivalent to reqiring that the conformal compactification of Q is depicted by a Penrose diagram as in Figre in the introdction. For more discssion on Penrose diagrams, we refer the reader to [2, Appendix C] and [5]. In fact, this eqivalence follows easily from the classification of all possible Penrose diagrams for the system SSESF given in the latter reference. We also define the precise notion of asymptotic flatness for initial data with BV or C reglarity. As we shall see soon in the main theorems, the rate of decay for the initial data, measred in r, is directly related to the rate of decay of the corresponding soltion, in both and r. 5 We refer to the proof of Proposition 5.3 below for details of estimating ν in terms of assmptions on φ, vrφ µ and λ. 9

10 Definition 2.9 Asymptotic flatness of order ω in BV or C. For ω >, we make the following definition. We say that an initial data set is asymptotically flat of order ω in BV if v rφ, BV[, and there exists I > 0 sch that 2.3 sp + r ω v rφ I <. C 2 We say that an initial data set is asymptotically flat of order ω in C if v rφ, C [, and there exist I 2 > 0 sch that 2.4 sp + r ω v rφ + sp + r ω + vrφ 2 I 2 <. C C Remark 2.0. The initial Bondi mass M i := lim v m, v can be easily bonded by CI 2 ; see Lemma 4.5. Remark 2.. Observe that both conditions imply that rφ, v tends to a finite limit as v ; in particlar, lim v φ, v = 0. This serves to fix the gage freedom φ, r, m φ + c, r, m for soltions to SSESF. 3. Main reslts 3.. Main theorems. With the definitions of locally BV scattering soltions and asymptotically flat initial data, we now have the necessary means to state the main theorems of this paper. Roghly speaking, these theorems say that locally BV scattering soltions with asymptotically flat initial data exhibits qantitative decay rates, which can be read off from the rate ω in Definition 2.9. The first theorem is for initial data and soltions in BV. Theorem 3. Main Theorem in BV. Let φ, r, m be a locally BV scattering soltion to SSESF with asymptotically flat initial data of order ω in BV. Then for ω := min{ω, 3}, there exists a constant A > 0 sch that φ A min{ ω, r ω }, v rφ A min{ ω, r ω }, rφ A ω. The second theorem is for initial data and soltions in C. Theorem 3.2 Main Theorem in C. Let φ, r, m be a locally BV scattering soltion to SSESF with asymptotically flat initial data of order ω in C. Then, in addition to the bonds , there exists a constant A 2 > 0 sch that for ω := min{ω, 3}. 2 vrφ A 2 min{ ω+, r ω+ }, 2 rφ A 2 ω+, v λ A 2 min{ 3, r 3 }, ν A 2 3. Some remarks regarding the main theorems are in order. Remark 3.3. Notice that in Theorem 3.2, the decay rates for v λ and ν are independent of the order ω of asymptotic flatness of the initial data. This is becase the scalar field terms enter the eqations for v log λ and v log ν qadratically see eqations 4.6 and 4.7 and ths as 0

11 long as ω >, their contribtions to the decay rates of v λ and ν are lower order compared to the term involving the Hawking mass. Remark 3.4. By Remark 2.3, a C initial data set gives rise to a C soltion. Hence Remark 2.6 applies, and the conditions 2 of Definition 2.4 may be replaced by a single eqivalent condition of ftre casal geodesic completeness of M, g in the case of Theorem 3.2. Remark 3.5. In general, the constants A and A 2 depend not only on the size of the initial data e.g., I, I 2, bt rather on the fll profile of the soltion. Nevertheless, for the special case of small initial total variation of v rφ, A and A 2 do depend only on the size of the initial data; see 3.3. Remark 3.6. If the initial data also verify higher derivative estimates, then the techniqes in proving Theorems 3. and 3.2 also allow s to derive decay bonds for higher order derivatives. The proof of the higher derivative decay estimates is in fact easier than the proofs of the first and second derivative decay bonds since we have already obtained sfficiently strong control of the scalar field and the geometry of the spacetime. We will omit the details. Remark 3.7. The decay rates that we obtain in these variables imply as immediate corollaries decay rates for v φ, φ, etc. See Corollaries 6.9 and 7.3. Remark 3.8. The decay rates in the main theorems are measred with respect to the doble nll coordinates, v normalized at the initial crve and the axis Γ as in 2.4. To measre the decay rate along nll infinity, one can alternatively normalize the coordinate 6 by reqiring r = 2 at ftre nll infinity. As we will show in 5.3, for the class of spacetimes considered in this paper, the decay rates with respect to this new system of nll coordinates are the same p to a constant mltiplicative factor. Remark 3.9. In view of Remark 2.7, the assmption of local BV scattering can be replaced by the bondedness of the sbcritical qantities vrφ 2 p C, C Q cpt v log λ p C, C Q cpt for p >. This is becase for fixed ɛ 0, one can choose R to be sfficiently small depending on C and apply Hölder s ineqality to ensre that vrφ 2 ɛ 0, v log λ ɛ 0. C Q cpt C Q cpt Remark 3.0. We also notice that the proof of or main theorem can be carried ot in an identical manner for locally BV scattering soltions arising from asymptotically flat Cachy data. More precisely, we can consider initial data given on a Cachy hypersrface v = f, with C f C and satisfying the constraint eqation together with the following bonds on the initial data: + r ω φ + v rφ + rφ + λ 2 + ν + Ĩ, 2 and + r ω + vrφ 2 + rφ 2 + v log λ + log ν Ĩ2. Then, if we assme in addition that the spacetime is locally BV scattering to the ftre, the conclsions of Theorems 3. and 3.2 hold. 6 In particlar, this normalization is sed in [2] for the black hole case. By changing the nll coordinates, we can ths more easily compare the decay rates in or setting and that in [2].

12 Remark 3.. Or main theorems can also be viewed as reslts on pgrading qalitative decay to qantitative decay estimates. Sch problems have been widely stdied in the linear setting withot the assmption on spherical symmetry on nontrapping asymptotically flat Lorentzian manifolds [3], [2], [7], as well as for the obstacle problem on Minkowski space [8], [20]. In the nonlinear setting, we mention the work of Christodolo-Tahvildar Zadeh [0], who stdied the energy critical 2-dimensional spherically symmetric wave map system and proved asymptotic decay for the soltion and its derivatives BV scattering and the blow-p at infinity scenerio. The condition of local BV scattering in the main theorems follows if one rles ot, a priori, a blow-p at infinity scenario. More precisely, we say that a soltion blows p at infinity if some scale-invariant spacetime norms are infinite as follows: Definition 3.2. Let φ, r, m be a BV soltion to SSESF sch that the condition of Definition 2.4 ftre completeness of radial nll geodesics holds. We say that the soltion blows p at infinity if at least one of the following holds: sp λ =, where λ Γ := lim λ, v. 2 3 Γ v + v λ φ λ v φdvd =, φ v ν rφ v φ ν rφdvd =. Remark 3.3. We do not prove in the paper the existence or non-existence of soltions that blow p at infinity. We remark that this is analogos to the blow-p at infinity scenarios which have recently been constrcted in some simpler semilinear, critical wave eqations [4]. It follows from or main theorem that if a soltion does not blow p at infinity, it obeys qantitative decay estimates. More precisely, we have Theorem 3.4 Dichotomy between blow-p at infinity and BV scattering. Let φ, r, m be a BV soltion to SSESF sch that the condition of Definition 2.4 ftre completeness of radial nll geodesics holds. Assme frthermore that the initial data for φ, r, m obey the condition 7 lim v φ, v = 0 and 3.8 vrφ 2 dv + sp v rφ <. C C Then either the soltion blows p at infinity; or 2 the soltion is globally BV scattering, in the sense that the conditions 2 and 3 of Definition 2.4 hold with R =. This theorem is established in Section 8. It then follows from or main theorems Theorems 3. and 3.2 that if a BV soltion does not blow p at infinity and possesses asymptotically flat initial data, then it obeys qantitative decay estimates Refinement in the small data in BV case. By a theorem of Christodolo [6], the maximal development of data with small BV norms does not blow p at infinity. The previos theorem applies, and ths the corresponding soltion is globally BV scattering, in the sense described in Theorem 3.4. Moreover, a closer inspection of the proof of the main theorems reveals that the following stronger conclsion holds in this case. 7 By Remark 2., note that this is the only condition on limv φ, v which is consistent with asymptotic flatness. 2

13 Theorem 3.5 Sharp decay for data with small BV norm. There exists a niversal ɛ 0 > 0 sch that for 0 < ɛ ɛ 0, the following statements hold. If the initial data set is asymptotically flat of order ω in BV and C 2 vrφ < ɛ, then the maximal development φ, r, m is globally BV scattering, in the sense that Definition 2.4 holds with arbitrarily large R > 0. Moreover, it satisfies the estimates with A C I I + ɛ. Here and similarly in 2, we se the convention that C I depends on I in a nondecreasing fashion. 8 2 If, in addition, the initial data set is asymptotically flat of order ω in C, then the maximal development also satisfies with A 2 C I2 I 2 + ɛ. The point of this theorem is that we only need to know that the initial total variation to be small in order to conclde pointwise decay rates; in particlar, I, I 2 can be arbitrarily large. In this sense, Theorem 3.5 generalizes both the small BV global well-posedness theorem [6, Theorem 6.2] and the earlier small data scattering theorem [3] for data that are small in a weighted C norm. A proof of this theorem will be sketched in Section Optimality of the decay rates. Or main theorems show pper bonds for the decay rates of the scalar field φ and its derivatives both towards nll infinity i.e., in r and along nll infinity i.e., in. For ω = ω < 3, if the decay rate of the initial data towards nll infinity satisfies also a lower bond, then we can show that both the r and decay rates in Theorem 3. are satrated. More precisely, Theorem 3.6 Sharpness of t ω tail for < ω < 3. Let < ω < 3. Sppose, in addition to the assmptions of Theorem 3., that there exists V sch that the initial data set satisfies the lower bond r ω v rφ, v L > 0, for v V. Then there exists a constant L ω > 0 sch that for sfficiently large. v rφ, v L ω min{r ω, ω }, rφ, v L ω ω, Remark 3.7. One can also infer the sharpness of the decay of φ from that of its derivatives. We will omit the details. This theorem will be proved in 0.. In fact, the proof of this theorem is similar to the proof of the pper bonds in the first main theorem Theorem 3.. We will show that after restricting to sfficiently large, the initial lower bond propagates and the nonlinear terms only give lower order contribtions. Notice also that the analogos statement is false for ω 3, since the nonlinear terms may dominate the contribtion of the initial data. For ω 3, we can show that the decay rates in Theorem 3. are sharp in the following sense: Theorem 3.8 Sharpness of t 3 tail. For arbitrarily small ɛ > 0, there exists a locally BV scattering soltion φ, r, m to SSESF which satisfies the following properties: 8 In particlar, for I sfficiently small, we have the estimate A CI + ɛ for some absolte constant C. 3

14 v rφ, v is smooth, compactly spported in the v-variable and has total variation less than ɛ, i.e., C 2 vrφ < ɛ. 2 There exists a constant L 3 > 0 sch that for sfficiently large. v rφ, v L 3 min{r 3, 3 }, rφ, v L 3 3, To prove Theorem 3.8, we will first establish a sfficient condition for the desired lower bonds in terms of non-vanishing of a single real nmber L, which is compted from information at the nll infinity. This reslt Lemma 0. is proved sing the decay rates proved in the main theorems, and we believe it might be of independent interest. In 0.3, we will complete the proof of Theorem 3.8 by constrcting an initial data set for which L can be bonded away from zero. This can be achieved by showing that the soltion is close to that of a corresponding linear problem and controlling the error terms after taking ɛ > 0 to be sfficiently small and sing Theorem Strategy of the proof of the main theorems. Roghly speaking, the proof of decay of φ and its derivatives can be split into three steps. In the first two steps, we control the incoming part 9 of the derivatives of the scalar field and metric components, i.e., v rφ, 2 vrφ and v λ. To this end, we split the spacetime into the exterior region Q ext := {, v Q : v 3} and the interior region Q int := {, v Q : v 3}. In the first step, we control the incoming part of the soltion in the exterior region. In this region, we have r v,, ths the negative r weights in the eqations give the reqired decay of φ and its derivatives. We then prove bonds in the interior region in the second step. Here, we exploit certain non-qantitative smallness in the spacetimes qantities as given by the assmption of local BV scattering to propagate the decay estimates from the exterior region to the interior region all the way p to the axis. Finally, in the third step, we control the otgoing part of the soltion, i.e., rφ, 2 rφ and ν, by showing that the decay bonds that we have proved along the axis can be propagated in the otgoing direction. We remind the readers that the above sketch is only a heristic argment and is not tre if taken literally. In particlar, in order to carry ot this procede, we need to first show that the local BV scattering assmption provides some control over the spacetime geometry. As we will show below, the estimates are derived in slightly different fashions for the first and the second derivatives of rφ. We note in particlar that carrying ot this general scheme relies heavily on the analytic strctre of the Einstein-scalar field eqations, inclding the montonicity properties as well as the nll strctre of the renormalized eqations Estimates for first derivatives of rφ. To obtain decay bonds for the first derivatives of rφ, we will rely on the wave eqation v rφ = 2mλν µr 2 φ. Notice that when we solve for the incoming radiation v rφ sing this as a transport eqation in, the right hand side does not depend explicitly on the otgoing radiation rφ. Instead, the right hand side consists of terms that are either lower order in terms of derivatives or satisfy a certain monotonicity property. 9 We call these variables incoming becase they obey a transport eqation in the direction. 4

15 In particlar, this eqations shows that as long as φ can be controlled, we can estimate v rφ by integrating along the incoming direction. On the other hand, we can also control φ once a bond on v rφ is known by integrating along the otgoing v direction. To achieve the desired decay rates for φ, v rφ and rφ, we follow the three steps otlined above: Bonds 0 for v rφ and φ in v 3: In the exterior region, we have r, v, it is therefore sfficient to prove the decay in r. First, we prove that sp C + rφ is bonded. This is achieved in a compact region by continity of the soltion and in the region of large r by integrating v rφ in the otgoing direction from the compact region. Since v rφ can in trn be controlled by φ, we get the desired bond. To improve over this bond we define B U := sp sp ω φ + r ω φ [,U] C and show via the wave eqation that r ω v rφ C + ɛ B U, where ɛ 0 as. This gives the optimal decay rate for v rφ in the exterior region p to an arbitrarily small loss, which can be estimated once B U can be controlled. 2 Bonds for v rφ and φ in v 3: For the decay of the first derivatives, the interior region {v 3} is frther divided into the intermediate region {r R} and the compact region {r R}. In these two regions, the r weight in the eqation is not sfficient to give the sharp decay rate. Instead, we start from the decay rate v rφ obtained in the first step in the exterior region and propagate this decay estimate inwards. To achieve this, we need to show that 2mλν is small when is sfficiently large. µr 2 2a r R and v 3: In the intermediate region where we still have a lower bond on r, the reqired smallness is given by the qalitative information that the Hawking mass approaches 0. Ths, from some large time onwards, 2mλν becomes sfficiently small and we can µr 2 integrate the wave eqation directly to obtain the desired decay bonds. 2b r R and v 3: In this region, we se the local BV scattering assmption to show that 2mλν {r R} 0 as. This smallness allows s to propagate the decay estimates µr 2 from the crve r = R to the region r < R. At this point, we can also recover the control for B U and close the estimates in step. This allows s to derive all the optimal decay rates for φ and v rφ 3 Bonds for rφ: To achieve the bonds for rφ, first note that along the axis we have rφ = v rφ. Ths, by the previos derived control for v rφ, we also have the decay of rφ along the axis. We then consider the wave eqation as a transport eqation in the otgoing direction for rφ to obtain the sharp decay for rφ in the whole spacetime Estimates for second derivatives of rφ. As for the first derivatives, we control the second derivatives by first integrating the eqation in the exterior region p to a crve v = 3. We then propagate the decay bonds from the exterior region to the interior region sing the estimates already derived for the first derivative of φ, as well as the local BV scattering assmption. However, at this level of derivatives, some new difficlties arise as we now describe. 0 The estimates in this region are similar to the corresponding bonds for the black hole case in [2]. There, it was observed that the qantity vrφ, which Dafermos-Rodnianski called an almost Riemann invariant, verifies an eqation sch that the right hand side has sefl weights in r and give the desired decay rates. In particlar, since we are simply sing compactness, the constants in Theorem 3. depend not only on the size of the initial data. 5

16 Renormalization and the nll strctre. The assmption of local BV scattering implies that 3.9 v φ + vrφ 2 0 C {r R} as. When combined with Christodolo s BV theory, this also implies that as v, we have 3.0 φ + rφ 2 0. C v {r R} Notice that on C resp. C v, we only control the integral of vrφ 2 and v φ resp. rφ 2 and φ. Sppose when integrating along the incoming direction to control vrφ 2 and v λ, we need to estimate terms of the form φ v φ. C v {r R} We can apply the BV theory to show that for v sfficiently large, φ ɛ. On the other hand, one can show that C v {r R} sp v φ C sp vrφ 2 C v {r R} J C v Q cpt which can be controlled by the qantity that we are estimating. However, in eqation 4.2 for 2 vrφ derived by differentiating 2.8, there are terms of the form v φ v φ sch that neither of the factors can be controlled a priori in L by the local BV scattering assmption. In other words, the eqation does not obey any nll condition. To deal with this problem, we follow [6] and introdce the renormalized variables vrφ 2 v λφ, rφ 2 νφ, v log λ λ µ µ r + vφ λ v rφ ν rφ, log ν ν µ µ r + φ λ v rφ ν rφ which have the property that the nonlinear terms arising in the eqations for these variables in fact have a nll strctre. In particlar, we can apply the above heristic procedre to obtain decay estimates in the compact region r R. Non-renormalized variables and decay towards nll infinity. While the renormalization allows s to apply the BV theory in the interior region, it does not give the optimal r decay rates in the exterior region. For example, the renormalized qantity µ λ v log λ µ r + vφ λ v rφ ν rφ decays only as r 2 towards nll infinity de to the contribtion of µ µ λ r, which is weaker than the desired r 3 decay for v log λ. Therefore, in order to obtain the optimal estimates everywhere in the spacetime, we need to se the variables 2 vrφ, 2 rφ, v λ and ν together with their renormalized versions. 6

17 Copling of the incoming and otgoing parts. Finally, an additional challenge is that nlike the estimates for the first derivatives of the scalar field, the bonds for the incoming part of the soltion vrφ 2 and v λ are copled to that for the otgoing part rφ 2 and ν. Likewise, to control rφ, 2 we need estimates for vrφ 2 and v λ. For example, in the eqation for v log λ µ λ µ r + vφ λ v rφ ν rφ, there is a term involving rφ 2 on the right hand side. In particlar, in order to obtain the desired decay for v λ, we need to at the same time prove the decay for rφ. 2 Strategy for obtaining the decay estimates. With the above difficlties in mind, we can now give a very rogh sketch of the strategy of the proof. Bonds for vrφ 2 and v λ for large r: As in the case for the first derivatives, we first prove the optimal r decay for vrφ 2 and v λ in the exterior region. To this end, we integrate the eqations satisfied by the non-renormalized variables. We note that the error terms can all be bonded sing the local BV scattering assmption and the decay estimates already proved for the first derivatives. 2 Bonds for all second derivatives: Steps 2 and 3 for the decay bonds for the first derivatives are now copled. Define B 2 U := sp sp ω vrφ 2 + ω rφ 2 + ω v λ + ω ν. [,U] C We then show that B 2 U can control the error terms arising from integrating the renormalized eqations in the sense that we can obtain an ineqality of the form weighted renormalized variables C 2 + ɛ 2 B 2 U, where ɛ 2 0 as 2. We then prove that the renormalized variables in fact control all the weighted second derivatives in B 2. After choosing 2 to be sfficiently large, we show that B 2 U is bonded independent of U and ths all the second derivatives have ω decay. 3 Optimal bonds in terms of decay: While we have obtained ω decay for the second derivatives, the decay rates are not the sharp rates claimed in the main theorem. To finally obtained the desired bonds, we integrate the eqations of the non-renormalized variables and se the preliminary estimates obtained in and 2 above. Here, we make se of the fact that the estimates obtained in step 2 above are sfficiently strong both in terms of reglarity and decay to control the error terms in the non-renormalized eqations. 4. Analytic properties of SSESF In this section, we discss the analytic properties of SSESF. These inclde scaling, monotonicity and the nll strctre of the system. All these featres will play crcial roles in the analysis. 4.. Scaling. For a > 0, SSESF is invariant nder the scaling of the coordinate system together with the scaling of the fnctions a, This in particlar implies that the BV norms v av r ar, m am, Ω Ω, φ φ. 2 vrφ, v dv and 7 v λ, v dv

18 are scale invariant. Ths the a priori assmptions 2.2 are taken with respect to localized versions of scale invariant norms Monotonicity properties. We first begin with basic monotonicity properties of r. Lemma 4. Monotonicity of r. Let φ, r, m be a BV soltion to SSESF. Then we have ν < 0 in Q, and λ > 0 when µ > 0, λ = 0 when µ = 0, λ < 0 when µ < 0. Proof. This was proved in [6, Propositions. and.2]; we reprodce the proof for the reader s convenience. Note the eqation v r 2 = 2 Ω2. which easily follows from SSESF. As r 2 = 2r r = 0 on Γ and r > 0 on Q, we easily see that ν < 0. Then from the definition of µ, the second conclsion also follows. According to the sign of λ, a general Penrose diagram Q is divided into three sbregions as follows: T := {, v Q : λ < 0}, A := {, v Q : λ = 0}, R := {, v Q : λ > 0}. These are called the trapped region, apparent horizon, and reglar region, respectively. The next lemma, which we borrow from [6], shows that the soltions to SSESF considered in this paper consist only of the reglar region R. Therefore, extensive discssion of T and A will be sppressed. Lemma 4.2 [6, Proposition.4]. Let φ, r, m be a BV soltion to SSESF. Then the casal past of Γ in Q is contained in R. In particlar, Q = R if φ, r, m satisfies the condition in Definition 2.4 ftre completeness of radial nll geodesics. Next, we trn to monotonicity properties of the Hawking mass m, which will play an important role in or paper. The following lemma is an obvios conseqence of 2.7. Lemma 4.3 Monotonicity of m. For a BV soltion φ, r, m to SSESF, we have v m 0, m 0 in R. By the monotonicity v m 0, the limit M := lim v m, v exists possibly + at this point for each. This is called the Bondi mass at retarded time. The following statement is an easy corollary of the preceding lemma. Corollary 4.4 Monotonicity of the Bondi mass. Let φ, r, m be a BV soltion to SSESF, and sppose that C R for [, 2 ]. Then the Bondi mass M is a non-increasing fnction on [, 2 ]. The following lemma shows that M i < for initial data sets considered in this paper. Lemma 4.5. Sppose that v rφ, is asymptotically flat or order ω > in the sense of Definition 2.9. Then we have 4. M i := lim v m, v CI2. 8

19 This is an easy conseqence of 2.7 and Lemma 4.; we omit its proof. By the preceding corollary, it follows that M < for each. We conclde this sbsection with additional monotonicity properties of soltions to SSESF, sefl for controlling the geometry of locally BV scattering soltions to SSESF. Lemma 4.6. Let φ, r, m be a BV soltion to SSESF. For, v R, we have λ λ, v, v, µ µ λ = v ν 0. Proof. The lemma follows from the formla λ log µ = ν r φ 2 and Nll strctre of the evoltion eqations. In this sbsection, we follow [6] and demonstrate that the evoltion eqations verify a form of nll strctre. In particlar, the nll strctre occrs in the eqations for the second derivatives of the scalar field and the metric. However, it is not apparent if we simply take the derivatives of the eqations 2.6 and 2.8. Instead, we rewrite the eqations in renormalized variables for which the nll strctre is manifest. We will perform this renormalization separately for the wave eqations for φ and for the eqations for λ and ν. - The wave eqation for φ. Taking v of the eqation 2.8, we obtain 2 vrφ = v λ φ = λ v φ + v λφ, or eqivalently, after sbstitting in the first eqation in 2.6, vrφ = 2mλν µr 2 vφ + ν µ vφ 2 φ + 2mν µr 2 vλφ 4m µr 3 λ2 νφ. Some terms on the right hand side, sch as µ ν v φ 2 φ, do not exhibit nll strctre and are dangeros near Γ. To tackle this, we rewrite Ths, from the first eqation, we derive v λφ = [ v λφ] v λ φ. 4.3 [ 2 vrφ v λφ] = λ v φ v λ φ. By switching and v, we obtain the following analogos eqations in the conjgate direction. 4.4 v 2 rφ = 2mλν µr 2 φ + λ µ φ 2 φ + 2mλ µr 2 νφ 4m µr 3 λν2 φ. 4.5 v [ 2 rφ νφ] = v ν φ ν v φ. - The eqations for λ and ν. From 2.6, we have log λ = µ µr ν, v log ν = 9 µ µr λ.