Hörmander’s method for the characteristic Cauchy problem and conformal scattering for a nonlinear wave equation

The purpose of this note is to prove the existence of a conformal scattering operator for the cubic defocusing wave equation on a non-stationary background. The proof essentially relies on solving the characteristic initial value problem by the method developed by Hörmander. This method consists in slowing down the propagation speed of the waves to transform a characteristic initial value problem into a standard Cauchy problem.


The result
Consider a globally hyperbolic smooth manifold (M,ĝ), whose metricĝ satisfies the vacuum Einstein equations. Let be a Cauchy hypersurface in M and T be a futureoriented timelike vector normal to . Consider the Cauchy problem for the defocusing cubic equation∇ α∇ αφ =φ 3 B Jérémie Joudioux jeremie.joudioux@aei.mpg.de http://www.jeremiejoudioux.eu 1 Albert-Einstein-Institut, Max-Planck-Institut für Gravitationsphysik, Am Mühlenberg 2, 14476 Potsdam, Germany where∇ is the Levi-Civita connection associated withĝ. Since (M,ĝ) is globally hyperbolic, there exists a time foliation on the manifold, and the global existence of solutions and the existence of scattering operator can be discussed for this problem. Nonetheless, since the metric is a priori not static, the standard approach to prove the global existence, for instance, through Strichartz estimates, and the existence of the scattering operator, which would require a proof of the local energy decay, might fail. The cubic wave equation (1) enjoys nonetheless good properties regarding conformal transformations. Indeed, if one considers a metric g conformal to the original metricĝ one obtains the following identity where ∇ α is the Levi-Civita connection associated with g. Hence, ifφ is a solution to the problem (1), then −1φ is a solution to the corresponding geometric equation arising from the conformal metric g. This property can be exploited to prove the global existence of solutions to the problem (1) by considering the local existence for the conformal problem, and construct a scattering operator by relying on that operation. Penrose introduced in the 1960s a class of space-times admitting a conformal compactification. These space-times, known as asymptotically simple space-times, model isolated bodies. The compactification is obtained by completing the manifold with two disconnected hypersurfaces, null for the conformal metric when the cosmological constant vanishes, which represent the extremities of future and past null geodesics. These hypersurfaces are denoted by I + and I − , respectively. The strategy for the construction of the scattering operator exploiting the existence of the conformal compactification is the following. For convenience, the construction is done for smooth data with compact support.
Consider smooth compactly supported initial data (φ 0 ,φ 1 ) on for the Cauchy problem (1). These data are appropriately transformed into data (φ 0 , φ 1 ) for the conformal cubic wave equation: Since the problem is now set on a space compact in time, the global existence problem for the Cauchy problem (1) turns into a much easier local existence problem (2) for the rescaled equation. The existence result of Cagnac-Choquet-Bruhat [3] can be used to address the existence in this context. The radiation profile of the functionφ can be obtained by considering the trace of −1φ . We define the mapping, generalising the inverse wave operators: Energy estimates can be proven for the cubic wave operator and we prove in fact The existence of the wave operators defined on H 1 (I ± ), inverses of the trace operators T ± −1 , is obtained by solving the characteristic Cauchy problem for the conformal cubic wave equation with data on I ± . The generalisation of the standard scattering operator is defined as the operator In this paper, we prove the existence of a bi-Lipschitz conformal scattering operator, which extends the notion of scattering operator to space-times which are non-static: Theorem There exists a locally bi-Lipschitz operator S : H 1 (I − ) → H 1 (I + ), which associates with past radiation profiles of solutions to the problem (1) the corresponding future radiation profiles.

Conformal techniques, scattering, asymptotic behaviour
The well-posedness of the Cauchy problem for this equation, and in this geometric setting, has been addressed in [3]. The existence of a scattering operator on Minkowski space-time has been considered, on flat space-times in [15], for general semi-linear wave equations, and by conformal methods on flat space-time in [2]. The existence of a scattering operator, under more constraining assumptions, has been considered by the author in [14], and this work extends this previous work to the generic cubic wave equation. In particular, we remove an artificial assumption on the decay of the coefficient of the nonlinearity. The purpose of this assumption was to compensate for the blow-up of the Sobolev constant associated with the Sobolev embeddings of H 1 into L 6 in dimension 3. Since the metric is not stationary, this geometrical setting is a priori not amenable to standard analytic techniques to prove the existence of a scattering operator. To construct this operator, we exploit the conformal invariance of the equation, in conjunction with the conformal compactification of the space-time. The rescaled solution −1φ = φ can be extended up toM by means of Eq. (2). The trace of −1φ on the boundaries is the radiation profiles, in this conformal setting, discussed by Friedlander [9]. Hence, the trace operators T ± , associating with initial data for the Cauchy problem for the conformal wave equation (2) to the traces of the corresponding solutions to (2) on the boundaries I ± , generalise the standard inverse wave operators of the classical scattering theory, see [18]. The actual existence of the scattering operator is performed by inverting the trace operators, that is to say, within the appropriate function spaces, solving the characteristic Cauchy problem with data on the boundaries of M.
Conformal methods to study global problems for partial differential equations in relativity go back to Penrose and Sachs and the peeling of higher spin fields. It was further studied in the context of the Cauchy problem in the mid-80s, early '90s (see, for instance, [3]). Friedlander [9] pioneered scattering theory in relativity. The reader who wishes to have a more exhaustive bibliography should refer to [13,21,27]. In the past years, this approach by conformal techniques to understand the asymptotic behaviour of solutions of field equations in general relativity has been studied in various contexts: Mason and Nicolas [18] obtained the first analytic result for linear fields, followed by a peeling result on the Schwarzschild background, for scalar wave [19], spin 1/2 and spin 1 fields [20]. This has been later extended to nonlinear waves on Kerr black holes [25]. As mentioned before, the conformal scattering construction was extended to a nonlinear wave equation [14]. Similar constructions based on local energy estimates [23] were obtained on Reissner-Nordström black holes [22] for the Maxwell equations. More recently, the existence of a conformal scattering operator has been proven for Yang-Mills fields on the de Sitter background [30].

Description of the work
We now review the key points of the result. As already mentioned, this work is an extension of previous work [14]. We are in particular building on the estimates performed in the asymptotic spatial region i 0 .
The main improvement of this paper is the possibility to handle a more general nonlinearity. In the previous paper, we assumed that the nonlinearity was multiplied by a decaying function b. This restriction has its origin in the treatment of the characteristic initial value problem on I ± to construct the inverse of the trace operators T ± . The technique relies on a fixed point argument in the energy space obtained by foliating the interior of the light cones at infinity I ± . Because of the power nonlinearity, Sobolev embeddings are required. Nonetheless, the volume of the leaves goes to zero, and a well-known consequence is the blow-up of the related Sobolev constant, associated with the Sobolev embedding from H 1 into L 6 . The nonlinearity is multiplied by a function decaying sufficiently fast to compensate for this blow-up of the Sobolev constant.
To circumvent that issue, the strategy to approach the characteristic initial value problem is changed. In particular, this strategy avoids the use of the Sobolev embeddings on a shrinking foliation and relies on the work by Hörmander [12] for the wave equation. The principle is the following. One assumes the existence of a solution for the Cauchy problem for the considered geometric equation. 1 One considers an initial characteristic surface C. A time function being chosen to split the metric where h t is the Riemannian metric on the leaves of the foliation induced by the time function ( t ), the propagation speed of the metric is slowed down The initial characteristic hypersurface C becomes spacelike for g λ , and solutions to the Cauchy problem for the geometric equation associated with g λ can be solved. In the case of the cubic wave equation, this result has been obtained by Choquet-Bruhat-Cagnac [3]. Given initial data, a family of solutions is then obtained, depending on λ. Standard compactness arguments can then be used to prove the existence of an accumulation point as λ goes to one which satisfies the characteristic initial value problem. This method by Hörmander has been extended to metrics of weak regularity in [24] for the linear wave equation. This paper clarifies some points of [24]; in particular, in the way some energy estimates are performed and in the use of trace theorems for intermediate derivatives.
This technique by Hörmander can only be used in a neighbourhood of timelike infinity, where the Penrose compactification leads indeed to a compact space. Near spacelike infinity i 0 , the metric coincides with the Schwarzschild metric, and the chosen conformal factor does not lead to a compactification in the neighbourhood of i 0 . There, the existence of solutions to the characteristic initial value problem can be proven by relying on an explicit foliation. Global solutions to the Cauchy problem in the past of I + (resp. the future of I − ) are obtained by appropriately gluing the different solutions to the characteristic initial value problem in the neighbourhood of i ± and the neighbourhood of i 0 .

Organisation of the paper
Section 2 contains the precise geometric framework, in Sect. 2.2, and reminders on the cubic wave equation, in Sect. 2.3, in particular, the global existence of solutions to the defocusing wave equations, see Proposition 2.2 and some estimates for the characteristic initial value problem, see Proposition 2.4. Section 3 solves the characteristic initial value problem for the cubic wave equation, by Hörmander's method, see Theorem 3.2. Section 4 addresses the proof of the uniqueness of solutions to the Cauchy problem, see Proposition 4.5. The next section, Sect. 5, addresses the problem of solving the characteristic initial value problem, see Proposition 5.4. The global characteristic initial value problem is solved in Sect. 6. Finally, the trace operators, as well as the scattering operators, are obtained in Sect. 7, see Theorem 7.1, and its corollary, Corollary 7.2. "Appendix A" contains reminders on a trace theorem, see in particular Theorem A.1, and its use when energy estimates can be proven.

Conventions, notations
All along the paper, the metrics have signature (− + ++), and we use the Einstein summation conventions. The expression a b, where a and b are two functions defined on M means that there exists a constant C > 0 depending on the geometry such that a ≤ Cb on M.
If a b and b a, then one writes a ≈ b.

Regular asymptotically simple space-time
Penrose introduced the concept of asymptotically simple space-times [26,Chapter 9.6] as idealised space-times modelling a spatially localised gravitational source. These space-times do not need to satisfy Einstein equation. Nonetheless, they can be conceived as resulting from a small perturbation of Minkowski space-time (when the cosmological constant vanishes), in which case they are a strong form of asymptotically flat space-times (see [8,Chapter 2.3] for further discussion). They can be obtained by combining the result of gluing by Corvino-Schoen and Chruściel-Delay [5][6][7], and the result by Andersson-Chruściel [1] on the existence of space-time with smooth I + arising from regular enough data on a hyperboloid. Nonetheless, in the context of stability theorems of Minkowski space-times, the conformal boundary exists but it not smooth [4]. More importantly, smoothness at conformal infinity has been discarded as physically relevant for standard physical systems (see [4]). One important point of Hörmander's method is precisely that the regularity of the initial characteristic surface could be low, see [24], though we are not working in this frame. Though it is not necessary, we assume that the Einstein equations in vacuum with vanishing cosmological constant are satisfied. In particular, the scalar curvature of the physical metric vanishes. (1) The interior ofM is M.
The framework of this paper imposes the use of energy estimates to prove the existence of a solution of the Cauchy problem set up on the characteristic cones at infinity. These energy estimates usually rely on the use of Stokes' theorem, requiring that the metric is regular enough at the tip of the cone. This is why an extra regularity assumption is made at tips of the boundary of the manifold.
Assumption 1 Let (M,ĝ) be regular asymptotically simple space-time whose conformal compactification is denoted by (M, g = 2ĝ ). We assume that there exist a neighbourhood U 1 of i 0 in M and a system of coordinates (t, r , θ, φ) on U 1 ∩M such that, in U 1 ∩M, the metric is the Schwarzschild metric Another important assumption which is made on the structure of the asymptotically simple space-time and which is the most important restriction on the geometry: Assumption 2 Assume that there exists a neighbourhood U 2 of i ± in M, and an embedding of U 2 into a manifoldM, such that the metric g extends to a C ∞ -metric ontō M.

Conformal changes for the cubic wave equation
Consider the conformally related metrics g andĝ = −2 g. Then, it is a straightforward calculation thatφ is a solution to the cubic wave equation on the manifoldM if, and only if, onM, the function φ = −1φ satisfies the equation

Estimates for the defocusing wave equation
It should be noticed that we are working here with the defocusing wave equation. Given a time foliation ( t ) onM, with unit future-oriented normal T , we define the energy E(t) at time t where D is the induced connection on t . E(t) is approximately conserved along the evolution, in the sense that Hence, for the defocusing wave equation, it is possible to prove global existence for large data in The key tool is a Sobolev embedding from H 1 ( t ) into L 6 ( t ) in dimension 3. As long as these embeddings can be performed uniformly (for instance, when all the leaves are diffeomorphic to R 3 , or the 3-sphere), then the global existence for all data is ensured.
The strategy to perform the estimates is the following: • whenever we are working with a finite time interval, we use the standard energy estimates for the linear wave equations, similar to [28, Chapter 1, §3], where, for the local existence, H s -norms are propagated; • we otherwise use the conservation of the energy with the nonlinear term.

The Cauchy problem for the cubic wave equation
The section contains some preliminary known results about the global existence of solutions to the Cauchy problem for the cubic wave equation admits a global solution in C 1 (R, L 2 (X )) ∩ C 0 (R, H 1 (X )). Furthermore, the a priori estimate holds, for all t ∈ R:

Remark 2.3
We could add a function b in the equation as follows: Cagnac and Choquet-Bruhat in [3] addressed the local well-posedness for this equation under the following assumptions on the function b: the function b is bounded on M, once differentiable function on M, admits a C 1 -extension to M and, for a given futureoriented unit timelike vector field T a on the unphysical space-time M, there exists a constant c such that |T a ∇ a b| b.
The conformal scattering operator could be constructed in that situation in a similar fashion. This extension has very little meaning in the context of that work, and its purpose, and is therefore ignored.

A priori estimates for the characteristic initial value problem
In previous work, we have proven various a priori estimates. Since these estimates are the basis of Hörmander's method for solving the characteristic Cauchy problem, we recall in this section the various energy estimates which were previously proved. It is important to note that, in this previous work, as noted in [14], the a priori energy estimates are proved without the decay assumptions at the tip of the vertex on the function b, as considered in Remark 2.3.
The first result is a local a priori energy estimate in the neighbourhood of the tip of the vertex. The geometric and analytic settings are the following: let p be a point in M and U an open neighbourhood of p in M, with compact closure. Assume that one can define the past light cone C − ( p) globally in U . Let t be a time function defining a foliation { t } t=0···1 of the past of C − ( p). The unit future-oriented normal to t is denoted by T a . On C − ( p), let l be a generator of future-oriented null geodesics ending up at p. Let n be a null vector field transverse to C − ( p) (see Fig. 1 for a representation with respect to a given past light cone) such that The family (l, n) is completed into a family (l, n, e 1 , e 2 ) so that the family (l, e 1 , e 2 ) is tangent to the cone C − ( p) and (e 1 , e 2 ) is orthogonal to (l, n) . The derivatives with respect to the vectors e 1 , e 2 are denoted by ∇ S 2 .
We consider the standard Sobolev spaces on t , H 1 ( t ) and L 2 ( t ). The set of smooth functions on the cone C − ( p) in the future of 0 with compact support away from the tip p is endowed with the norm: where n dμ[ĝ] is the contraction of the ambient four-dimensional volume form dμ [ĝ] with the null vector n transverse to the light cone C − ( p). The completion for this norm of the set of smooth functions with compact support away from p is denoted by . Using standard energy estimates, one gets the following proposition: The following inequalities hold: and Furthermore, for all t, one has , one immediately gets the same inequality without the L 4 norms.

Characteristic Cauchy problem à la Hörmander
The purpose of this section is to establish an a priori well-posedness result of the characteristic Cauchy problem for the wave equation on a curved space-time based on the result by Hörmander [12], extended in [24] and partially used in [14] to prove the existence and uniqueness of a weak solution to the characteristic Cauchy problem in H 1 (M). More regularity is required when establishing the Lipschitz continuity of the wave operator. The method used by Hörmander is based on a reduction of the propagation speed of the considered wave equation so that one can resort to the standard existence result of the wave. In the following, one restricts oneself to a light cone, but the result can be extended in the same way to arbitrary weakly lightlike hypersurfaces.
The geometric setting is the following: consider a point p in M. One denotes by C − ( p) the past light cone from p. Consider a time function defined in the interior of C − ( p). The induced foliation is denoted by t . The part of C − ( p) between 0 and T is denoted C T . One assumes that the slice 0 is in the past of p and the slice containing p is denoted by T . The unit future-oriented normal vector to the time slices is denoted by T a . One considers a 3+1 splitting of the metric in the following form: where h t is a Riemannian metric on t and the lapse N is given by: The functional setting is given by: • on the time slice t , one defines the energy of a function φ: • on the light cone C − ( p), the following H 1 norm is considered: where l is a non-vanishing generator of the null directions of C − ( p) and T a dμ[g] is the contraction of the space-time volume form with the vector field T . The space is defined as being the completion of the space of smooth functions whose compact support does not contain p. One also define in similar way H 1 (C T ).

Remark 3.1
One could have defined, like Hörmander, H 1 (C − ( p)) by transporting the H 1 structure of a timelike slice on C − ( p). It happens that the two definitions coincide [14].
Consider the characteristic Cauchy problem Note that, for the discussion of this section, the scalar curvature term is excluded from the discussion, but can be treated similarly, exploiting the fact the curvature is bounded.
The purpose of this section is to prove that there exists a strong solution to the characteristic Cauchy problem up to the hypersurface 0 in the past of p.
Using Proposition 2.2, one can then prove the following theorem:

Theorem 3.2 The characteristic Cauchy problem (2.1) admits a global strong solution down to
The following a priori estimates furthermore hold: Proof Assuming that there exists a solution in C 0 (R, H 1 ( τ )), the a priori estimates are a consequence of in [14,Proposition 4.15]. The proof of the theorem will require at a point the use of Sobolev embeddings. This requires to work with an extension of the foliation ( τ ) to a cylinder. One then considers a smooth isometric embedding of (J − ( p) ∩ J + ( 0 ), g) into a compact cylinder (U , g), see Fig. 1. The foliation ( τ ) is extended on the cylinder U in a spacelike foliation of U for the extension of the metric g. We denote by (˜ τ ) this extension; all the leaves of this foliation are now topological 3-spheres endowed with a Riemannian metric, and as a consequence, all the Sobolev spaces H k (˜ τ ) considered on the leaf˜ τ are equivalent since the leaf is compact. Furthermore, the cone C T is extended as a weakly spacelike hypersurfaceC, which by construction is a Cauchy surface. The spacelike foliation ( τ ) is extended in the past ofC, for values of τ in a compact interval I = [a, b], with T ≤ b.
We consider the initial datum θ on C T in H 1 (C T ). This function θ is extended as a functionθ in H 1 (C) (see, for instance, [29, Chapter V, Theorems 5 and 5']) such that where the constant C depends solely on the geometry of the boundary of C T . In particular, this constant blows up as T goes to 0. The proof of the well-posedness relies on the method introduced by Hörmander [12] which consists in slowing down the propagation speed of the waves. A detailed proof of Hörmander's paper is given in [24] when the metric is only Lipschitz. The following proof follows step by step this paper (and especially the scheme given in the proof of Theorem 4.3), modifying it when necessary.
The reduction of the propagation speed is realised by introducing a parameter λ in 1 2 , 1 and a family of metrics g λ : g λ = −λ 2 N 2 dt 2 + h˜ t so that the hypersurfaceC becomes spacelike for the metric g λ . According to Proposition 2.2, for any λ in 1 2 , 1 , the Cauchy problem where c λ is a constant which depends continuously on the scalar curvature of g λ . Since the interval under consideration is compact and since, as a consequence of the previous remark, c λ depends continuously on λ, one can replace c λ by its supremum over 1 2 , 1 . Furthermore,C is compact and, as consequence, all the Sobolev spaces associated with a smooth metric are equivalent. This includes H 1 (C). The family Consider now a sequence (λ n ) converging to 1 and one denotes by (φ n ) the associated sequence. One denotes by U the volume delimited by 0 and C − ( p).

Remark 3.3
Before studying the convergence process, let us remind the reader that, using a priori estimates of Proposition 2.4, the energy on the initial time slice˜ 0 controls all the energies on the time slices˜ τ for τ ∈ I . As a consequence, a convergence stated for H 1 ( 0 ) actually holds in L ∞ (I , H 1 (˜ τ )).
One now proceeds by extracting successively sequences of (λ n ) as follows: Hence, using Kakutani's theorem, there exists a sub-sequence of (φ n ) converging weakly in H 1 (U ) towards a function φ in H 1 (U ). Furthermore, using Rellich-Kondrachov's theorem, since H 1 (U ) is compactly embedded in H σ (U ) for σ < 1, a diagonal extraction process gives: We have also used that φ n lies for all n in C 0 (I , L 2 (˜ τ )) which is closed in L ∞ (I , L 2 (˜ τ )). and, as a consequence, (2.10) (6) Furthermore, since g λ is smooth, λ φ λ also converges in the sense of distributions towards φ. Using the convergence (2.10), (φ 3 n ) converges towards φ 3 in the sense of distributions. This means that the function φ satisfies, in the sense of distributions, φ = φ 3 .
(7) It remains to prove that φ satisfies the initial conditions. This is a direct consequence of Eq. (2.6): The next step of the proof consists in proving that the solution is in fact a strong solution of the equation: (1) Since we are working on a compact space, the trace operator is in fact compact from H 1 (U ) in L 2 (C). Up to an extraction, (φ n ) converges then strongly in L 2 (C). As a consequence, φ| C T is equal to φ 0 in L 2 (C − ( p)). (2) One already has, as a consequence of the a priori estimates, (3) One then considers the operator where α is a smooth function on U , L 1 is a first-order purely spatial operator, and L 2 is a second-order purely spatial operator. It is clear from this decomposition that if φ lies in L ∞ (I , H 1 (˜ τ )) and ∂ τ φ in L ∞ (I , L 2 (˜ τ )), then L 0 φ lies L ∞ (I , H −1 (˜ τ )) since H 1 (˜ τ ) and L 2 (˜ τ ) embed themselves continuously and uniformly over the foliation in H −1 (˜ τ ).
(b) The a priori estimates and Eq. (2.11) give: As a consequence (this fact is proved in [24, p. 535]), we have: In this context, on can prove the a priori estimates over the foliation˜ τ : The energy is then locally Lipschitz continuous. As a consequence, since φ is in C 0 (I , w − H 1 (˜ τ )), φ is also in C 0 (I , H 1 (˜ τ )).

Estimates for the characteristic Cauchy problem
In the previous section, it has been proved that the Cauchy problem admits a local solution to the characteristic Cauchy problem, relying only on the existence of a priori estimates for the solutions. The uniqueness has, so far, not be proved yet, as well as the continuity in the initial data. This can be achieved in various ways. The path we chose to follow is based on the work of Baez-Segal-Zhou [2] and relies on an a priori existence result of solutions to the Cauchy problem, which has been handled in a specific way separately in Sect. 3. The estimates are established by relying on a reduction to the Cauchy problem. The setting is both the framework introduced by Hörmander and Baez-Segal-Zhou. The geometric setting and notations are the same as in the previous section, see Fig. 1, and the proof of Theorem 3.2.

Remark 4.1
The setting introduced by Hörmander considers only space-times of the form X × R, where X is a compact manifold, a priori without boundary. The result can nonetheless be immediately extended to the case when the manifold X has a boundary. This can be seen using the following remark: the energy space remains H 1 (X ) (and notḢ 1 (X )). Assume that the manifold with boundary X is embedded in a bigger compact manifold without boundary (which can always be done). If the boundary of X is at least C 1 , the space H 1 (X ) extends continuously into H 1 (X ). As a consequence, all energy estimates involving X can be brought back onto X . In this specific case, the boundary of T is the intersection of T with J − ( p) and is, as a consequence, smooth.
One introduces first the following operator: which associates with initial data (φ, ψ) the trace over C T of the solution of the linear wave equation with initial data (φ, ψ) on T . Following Hörmander [12], this operator has the following properties: Theorem 4.2 (Hörmander) The linear operator T is one-one, onto and bi-continuous, that is to say that there exists a constant C depending only on the geometry of the manifold such that, for all (φ, ψ) in H 1 ( 0 ) × L 2 ( 0 ): We introduce the extension operator E C and E 0 , as considered in [29, Chapter V, Theorems 5 and 5']:

Lemma 4.3 Let δ 0 be an initial data set in H 1 (C T ). Let H be a function defined over the foliation
. Then, the solution δ of the Cauchy problem can be extended to the past of˜ T by solving the Cauchy problem on˜ 0 with initial data E 0 T −1 (δ 0 ) for the equation: where the function H is extended by 0 outside J − (C T ).
For such an equation, the energy estimates are simple to obtain, since the foliation of reference to establish them have non-vanishing volume. For the sake of consistency, these estimates are nonetheless proved here. To that end, we introduce the stress-energy tensor The energy of a function φ on any weakly spacelike hypersurface S is defined as where is the Hodge dual with respect to the metric g. In particular, one checks that

Lemma 4.4 Let δ be as in Lemma 4.3. There exists an increasing function C
Proof Since the techniques under considerations are standard, the proof is only sketched. The error term associated with the considered stress-energy tensor is Applying Stokes theorem between 0 and s , one gets immediately the following inequality: One deals with the nonlinear term using Sobolev embeddings over the compact foliation τ , whose volume of leaves does not go to 0: One finally gets and using Grönwall's inequality closes the estimates.
Let now be θ and ξ be two functions in H 1 (C T ) and consider the two characteristic Cauchy problems: The function δ defined as being the difference between u and v satisfies the wave equation Using Lemma 4.3, we extend δ to the whole cylinder by considering δ to the problem with initial data E 0 T −1 (θ − ξ) . Using the continuity of trace operator, the continuity of the extension operator, and Lemma 4.4, one proves:

Proposition 4.5 There exists an increasing function C such that the following inequalities hold:
Proof The proof of these energy inequalities is a direct consequence of Proposition 4.4. It suffices to notice that, for H 2 = φ 2 + φψ + ψ 2 , using the triangular inequality, Using a priori estimates such as the one proved in [14, Proposition 6.2], the later terms can be bounded by either Finally, an important consequence of the energy estimate (4.5) and Theorem 3.2 is the following proposition: Proposition 4. 6 The characteristic Cauchy problem: admits at most one solution up to time T in C 0 ([0, T ], H 1 ( t )) ∩ C 1 ([0, T ], L 2 ( t )).

Solving the characteristic Cauchy problem in the neighbourhood of
i 0

Preliminary result
The purpose of this section is to explain how the characteristic Cauchy problem can be solved in the neighbourhood of the spacelike infinity i 0 of the asymptotically simple manifold. The neighbourhood of i 0 is assumed to be isometric to the Schwarzschild space-time to agree with the work of the Corvino-Schoen and Chruściel-Delay. The Schwarzschild metric is given, in the standard spherical coordinateŝ Performing the change of coordinates the metric takes the form A conformal rescaling with a conformal factor = 1 R is finally performed to define the unphysical metric: Consider the domain + u 0 = {u ≤ u 0 }, for some u 0 to be chosen later. It has been proved [14,18] that the following lemma holds: There exists u 0 < 0, |u 0 | large enough, such that the following decay estimates in the coordinate (u, r , θ, ψ) hold: Furthermore, the vector field The parameter will be chosen later when we perform the energy estimates. We define, in + u 0 = {t > 0, u < u 0 }, the following hypersurfaces, for u 0 given in R (Fig. 2): The volume form associated with g in the coordinates (R, u, ω S 2 ) is then: (4.1) If φ is function defined on + u 0 , one defines the following energies, in the domain Furthermore, if one introduces the parameter τ : the following energy estimates holds: • Energy decay: • A priori estimates:

Local existence for the Characteristic Cauchy problem on I +
The purpose of this section is to prove that the Cauchy problem: where H 1 (S u 0 ) is defined as being the completion of the space of traces of smooth functions on S u 0 for the norm T lin ab being the stress-energy tensor associated with the linear wave equation. In a previous work [14], the following result has been proved: • using the fact that the leaves of the foliation H s endowed with its induced metric are uniformly equivalent to a cylinder, it is possible to establish energy estimates for the difference of solutions of the nonlinear equation; • it is possible to establish an existence result for solutions of the nonlinear equation for small data, using either a Picard iteration (as in [13]) or Hörmander slowing down process.

Proposition 5.3
Let φ and ψ be two solutions of the characteristic Cauchy problem and denote by δ = φ − ψ their difference. One denotes the energy of δ by where T lin ab is the stress-energy tensor for δ associated with the linear wave equation. Then, δ satisfies the following energy estimates on the foliation H s ; In particular, the Cauchy problem Proof Let φ and ψ be two global solutions of (4.3). Their difference δ = φ − ψ satisfies the equation: Using the Stokes theorem between I + u 0 , S u 0 and H s , one gets, since δ vanishes on S u 0 and I + u 0 : As already said, Sobolev estimates can be performed uniformly on the foliation (H s ). As a consequence, using Hölder's inequality, one gets: where T a is the approximate Morawetz vector field The constant depends only on the L ∞ -bound of the scalar curvature. Finally, using Grönwall lemma, one gets that As a consequence, the mapping This defines a mapping One proves, using standard energy estimates that, for s ≤ , (4.5) In particular, φ belongs to C 0 [0, ], H 1 (H s ) . Let R be a positive real number. Assume that Then, there exists an small enough, depending on R, such that, for all ψ such that Furthermore, let ψ,ψ be two functions in C 0 [0, ], H 1 (H s ) , and consider φ,φ the corresponding solution to the problem (4.6). Their difference φ −φ satisfies One easily proves, as in Proposition 5.3, that where C is an increasing function of R. Hence, R fixed, for small enough, the mapping φ → ψ is a contraction of the ball of radius  The existence of global solutions to the Cauchy problem is proved using a gluing process: the Cauchy problem with data in H 1 (I + ) is constructed up to a spacelike hypersurface and, then, considering the traces of the solution on two given hypersurfaces in the neighbourhood of i 0 solved up to the initial time slice 0 . The main issue arising in this process is that the constants arising in the energy estimates depend on the L ∞ -bounds of the metric and its inverse on the bounded neighbourhood of i 0 . Since the compactification we are working with has the particularity to have an asymptotic end at i 0 , these constants are not bounded on the whole future of 0 . This problem is avoided as follows (see Fig. 3): • Let u 0 be in R such as in Proposition 5.2.
• Let be a given spacelike hypersurface such that: is transverse to I + ; -is in the past of S u 0 and in the future of 0 ; in particular, coincides with 0 far from i 0 . These choices of u 0 and are made once for all. in the neighbourhood of i 0 . Furthermore, ifφ is another solution with initial dataθ, the following energy estimates holds: there exists a increasing function f such that Proof The proof of the theorem relies on the consecutive use of Theorem 3.2, Proposition 4.6 and of Proposition 5.4. Let θ be in H 1 (I + ). Using Theorem 3.2, one proves that the Cauchy problem admits a solution φ up to , as defined earlier. Then, using Proposition 5.4, there exists a > 0 such that the Cauchy problem given by Consider the time slice H . This time slice can be extended intoM into a spacelike Cauchy slice denoted by H. The timelike vector field T is extended as a timelike vector to H. One now considers the functions (ψ 0 , ψ 1 ) defined piecewise by: and By construction, since H is a spacelike Cauchy surface on the unphysical spacetime, the result from Cagnac-Choquet-Bruhat [3] can be applied immediately to prove the existence of a unique solution to the Cauchy problem with data (ψ 0 , ψ 1 ) on H up to the Cauchy surface 0 . Furthermore, this solution belongs to Consequently, there exists a global solution of the Cauchy problem:

Existence and regularity of the scattering operator
The purpose of this section is to prove the regularity and the existence of conformal scattering operator for the nonlinear wave equation In the context of asymptotically simple space-times with regular i ± , such as the one considered in [14,18], it has been proved by Mason and Nicolas, for the scalar wave equation, that the trace operators can be obtained from the inverse wave operators˜ ± as follows: let F ± be the null geodesic flows identifying the hypersurface 0 with I ± ; then, T ± are given by where F is the pullback by the null geodesic flows F ± . The same result will hold in the context of the nonlinear equation. The conformal scattering result which is stated here should consequently be considered as a standard scattering result for a nonstationary metric for a nonlinear wave equation. The existence of scattering operators was obtained on the flat background for sub-critical wave equations by [10,11,15,32,33].
One considers the operators defined by where u is the unique solution of the Cauchy problem on the conformally compactified space-time One finally introduces the conformal scattering operator defined as The operators T ± : H 1 ( 0 ) × L 2 ( 0 ) → H 1 (I ± ) are well-defined, invertible and locally bi-Lipschitz, that is to say that T ± and (T ± ) −1 are Lipschitz on any ball in H 1 ( 0 ) × L 2 ( 0 ) and H 1 (I ± ), respectively.
Proof The existence and continuity of the operators guarantee the existence of solution to the Cauchy problem for the conformal wave equation, see Proposition 2.2. The existence of their inverses and their continuity are obtained thanks to the Lipschitzness guaranteed by the energy estimate of Proposition 6.1.
The final product of this paper, the existence of a scattering operator, is achieved by composing the operator (T − ) −1 with T + : Corollary 7.2 There exists a locally bi-Lipschitz invertible operator S : . This operator generalises the notion of classical scattering operators (defined when the metric is static) to non-stationary metrics.
Proof As mentioned before, the construction of S is a straightforward consequence of Theorem 7.1. The fact that the operator coincides with the classical notion of scattering operator is stated in [18].
Let a, b be two elements in R ∪ {±∞} and m be an positive integer. One considers the following space: The following trace theorem then holds, see [16, theorems 2.3 and 3.1, chapter 1,]: W (a, b). Then, for all j in {0, . . . , m − 1}, the derivatives of φ satisfy: Remark A. 2 The theory developed by Lions and Magenes applies to general Hilbert spaces, independently of the geometric context (in particular, if one works on a space with boundary or on a compact manifold).
One applies this result to the following situation: let M be a compact manifold and consider The interpolation spaces between both are given by [31,  Let t 0 be in I . The purpose is to prove that f is continuous in t 0 .
Let (w k ) be a sequence of smooth function converging towards w in H s . Let be a positive real and consider: Since, using Cauchy-Schwarz inequality, and since (w k ) converges towards w in H s , and since φ is in L ∞ (I , H s (M)), there exists a K such that: One now considers the interpolation operators between L 2 (M) and H s (M) (for arbitrary s): This is a family of unbounded, self-adjoint operators for which: Since w K is smooth, it admits a pre-image by F −1 s • F σ . As a consequence, one has: Cauchy-Schwarz inequality implies then: Since φ is in C 0 (I , H σ (M)), there exists an open neighbourhood U of t 0 in I such that, for all t in U ⊂ I : As a consequence, for all t in U , that is to say that f is continuous in t 0 .
One now considers a model case fitting the context of the wave equation. Let (I × M, −N 2 dt 2 + g t ) be a Lorentzian manifold, I being a compact interval. One defines the energy of the function as being the function of t and φ in H 1 (M) and ψ in L 2 (M) by: The final lemma required to achieve the wanted regularity for the solution of the characteristic Cauchy problem for the wave equation is the following: Lemma A. 6 One assumes that, for all (φ, ψ) in H 1 (M) × L 2 (M), the function: is continuously differentiable in I . Let φ be a function such that: • φ is in C 0 (I , H 1 (M) − w); as a consequence, φ lies in L ∞ (I , H 1 (M)) • ∂ t φ is in C 0 (I , L 2 (M) − w); as a consequence, φ lies in L ∞ (I , L 2 (M)).
The function: is furthermore assumed to be continuous in t.
Proof The proof of this fact can be found in [16, Chapter 3, Section 8.4, p. 279]. For the sake of self-consistency, the proof is quoted here, with some adaptations to our framework. Let t be in I and consider (t n ) n a sequence of elements of I converging towards t and define the quantity: Since φ is in C 0 (I , H 1 (M) − w) and ∂ t φ is in C 0 (I , L 2 (M) − w), the last two terms converge towards −2E(t, φ(t), ∂ t φ(t)). Since the energy is assumed to be continuous, the first tow terms converge towards 2E(t, φ(t), ∂ t φ(t)).
The sequence (ξ n ) converges then towards 0 when n grows. As a consequence, φ and ∂ t φ are in C 0 (I , H 1 (M)) and C 0 (I , L 2 (M)), respectively.