the notion of  conformal infinity  introduced by penrose almost forty years ago is one of the most fruitful concepts within einstein s theory of gravitation . most of the modern developments in the theory are based on or at least influenced in one way or another by the conformal properties of einstein s equations in general or , in particular , by the structure of null infinity : the study of radiating solutions of the field equations and the question of fall - off conditions for them ; the global structure of space - times ; the structure of singularities ; conserved quantities ; multipole moments ; asymptotic quantization procedures ; -spaces and the twistor programme ; the null hypersurface formulation of general relativity ; the conformal field equations and their importance for the numerical evolution of space - times . obviously , there exists a vast amount of material related to the subject of conformal infinity which can not be covered adequately within this review article . we will discuss here those issues of  conformal infinity  which are relevant for numerical applications . on the one hand , this is a restriction to a subtopic which is reasonably narrow and yet broad enough to encompass the central ideas and new developments . on the other hand , questions concerning the numerical treatment of gravitational radiation and , in particular the problems which arise from the attempt to numerically model infinitely extended systems , suggest that the conformal methods can be useful not only for rigorous arguments but also for numerical purposes . indeed , we will show that the conformal picture has matured enough to provide an approach to applications in numerical relativity which relies on a very sound theoretical basis not only with respect to the physical appropriateness , but also with respect to the mathematical well - posedness of the problems considered . what is even more remarkable is the fact that the numerical implementations of the conformal picture via the conformal field equations are numerically well - defined in the sense that there are no spurious instabilities in the codes ( which so often are the stumbling blocks for the traditional approaches via the adm equations ) , so that the computed solutions demonstrably converge to the order of the discretization scheme . the conformal approach based on the hyperboloidal initial - value problem allows us to compute ( semi-)global space - times including their asymptotic structure which in turn enables us to rigorously , i.e. without any further approximation beyond the discretization , determine the radiation coming out from the system under consideration . the work which has been devoted to this approach up to now clearly shows its power . the results obtained have been checked against exact results ( exact solutions or known theorems ) , and there is no doubt that the results are correct . the geometric concepts which have been devised by penrose now turn out to be very useful in practical applications in the sense that they provide the solution to all conceptual problems posed by the notion of  gravitational radiation  in connection with numerical computations . already , it is obvious that with this tool one can achieve results which have not been feasible by any other numerical method . furthermore , the analysis at space - like infinity ( see section 3.5 ) has finally put within reach the feasibility of computing the global structure of asymptotically flat space - times from asymptotically flat cauchy data . this article will begin with a survey of the general background , i.e. the physical motivations behind the idea of asymptotically flat space - times , and a short account of the historical development which led to our current understanding . next , the necessary mathematical ideas will be introduced and the regular conformal field equations will be discussed . the last section is concerned with current issues in the numerical implementation of the conformal methods . we will see that the conformal field equations provide a very powerful method to study global problems in numerical relativity such as gravitational wave propagation and detection as well as the emergence of singularities and their horizons . isolated systems provide valuable insights into the workings of a physical theory by restricting one s attention to small subsystems . they serve as models for systems in the real world allowing us to deduce statements about their behaviour , and to attribute to them various physical properties such as mass , momentum , emitted radiation , etc . therefore , it is desirable that a theory should allow within its mathematical framework the characterization of such systems . in general relativity , this is a difficult problem . the reason is a familiar one : the metric which , in other theories , provides a background structure on which the physical fields act , is itself a dynamical object in general relativity . in this section we discuss some of the issues which lead us to focus on asymptotically flat space - times as models for isolated systems in relativity and hence as realistic gravitationally radiating systems . the conformal structure of space - times has found a wide range of interesting applications in general relativity with various motivations . of particular importance for us is the emergence of conformal geometric ideas in connection with isolated systems ( see also the related discussions in [ 66 , 56 ] ) . as an illustration , consider a gravitating system ( e.g. a binary system or a star ) somewhere in our universe , evolving according to its own gravitational interaction , and possibly reacting to gravitational radiation impinging on it from the outside . we are interested in detecting and evaluating these waves because they provide us with important information about the physics governing the system . for several obvious reasons , it is desirable not only to have a description of such situations within the theoretical framework of the theory but , furthermore , to have the ability to simulate them numerically . two problems arise : first , we need to idealize the physical situation in an appropriate way , since it is hopeless to try to analyze the behaviour of the system in its interaction with the rest of the universe . we are mainly interested in the behaviour of the system and not so much in other processes taking place at large distances from the system . since we would like to ignore those regions , we need a way to isolate the system from their influence . we might want to do this by cutting away the uninteresting parts of the universe along a time - like cylinder t enclosing the system . thereby , we effectively replace the outer part by data on t. the evolution of our system is determined by those data and initial data on some space - like hypersurface s. but now we are faced with the problem of interpreting the data . it is well known that initial data are obtained from some free data by solving elliptic equations . it is very difficult to give a physical meaning to initial data obtained in this way , and it is even more difficult , if not impossible , to specify a system , i.e. to determine initial data , exclusively from ( local ) physical properties of the constituents of the system like energy - momentum , spin , material properties , etc . in a similar spirit , the data on the time - like boundary t are complicated and only to a rather limited extent do they lend themselves to physical interpretation . for instance , it is not known how to extract from those data any piece which would unambiguously correspond to the radiation emitted by the system . how can we be sure that we capture essentially the same behaviour independently of how we define t ? thus , we are led to consider a different kind of  isolation procedure  . alone in the universe  in the sense that we assume it being embedded in a space - time manifold which is asymptotically flat . how to formulate this is a priori rather vague . somehow we want to express the fact that the space - time  looks like  minkowski space - time  at large distances  from the source . certainly , fall - off conditions for the curvature have to be imposed as one recedes from the source and these conditions should be compatible with the einstein equations . this means that there should exist solutions of the einstein equations which exhibit these fall - off properties . we would then , on some initial space - like hypersurface s , prescribe initial data which should , on the one hand , satisfy the asymptotic conditions . on the other hand , the initial data should approximate in an appropriate sense the initial conditions which give rise to the real behaviour of the system . our hope is that the evolution of these data provides a reasonable approximation of the real behaviour . as before , the asymptotic conditions which in a sense replace the influence of the rest of the universe on the system should not depend on the particular system under consideration . furthermore , we would hope that the conditions are such that there is a well defined way to allow for radiation to be easily extracted . it turns out that all these desiderata are in fact realized in the final formulation . these considerations lead us to focus on space - times which are asymptotically flat in the appropriate sense . even if we choose the idealization of our system as an asymptotically flat space - time manifold , we are still facing the task of adequately simulating the situation numerically . this is a formidable task , even when we ignore complications arising from difficult matter equations . the simulation of gravitational waves in an otherwise empty space - time coming in from infinity , interacting with themselves , and going out to infinity is a challenging problem . the reason is obvious : asymptotically flat space - times necessarily have infinite extent while computing resources are finite . the conventional way to overcome this apparent contradiction is the introduction of an artificial boundary  far away from the interesting regions  . during the simulation this boundary evolves in time thus defining a time - like hypersurface in space - time . there one imposes conditions which , it is hoped , approximate the asymptotic conditions . however , introducing the artificial boundary is nothing but the reintroduction of the time - like cylinder t on the numerical level with all its shortcomings . instead of having a  clean  system which is asymptotically flat and allows well defined asymptotic quantities to be precisely determined , one is now dealing again with data on a time - like boundary whose meaning is unclear . even if the numerical initial data have been arranged so that the asymptotic conditions are well approximated initially by the boundary conditions on t , furthermore , the numerical treatment of an initial - boundary value problem is much more complicated than an initial value problem because of instabilities which can easily be generated at the boundary . what is needed , therefore , is a definition of asymptotically flat space - times which allows to overcome both the problem of  where infinity is  and the problem of simulating an infinite system with finite resources . the key observation in this context is that  infinity  is far away with respect to the space - time metric . but , what if we replaced these metre sticks by ones which grow in length the farther out we go ? then it might be possible that only a finite number of them suffices to cover an infinite range , provided the growth rate is just right . this somewhat naive picture can be made much more precise : instead of using the physical space - time metric g to measure distance and time , we use a different metric g = g , which is  scaled down  with a scale factor . if  can be arranged to approach zero at an appropriate rate , then this might result in  bringing infinity in to a finite region  with respect to the unphysical metric g. we can imagine attaching points to the space - time which are finite with respect to g but which are at infinity with respect to g . in this way we can construct a boundary consisting of all the end points of the succession of finitely many rescaled metre sticks arranged in all possible directions . this construction works for minkowski space and so it is reasonable to define asymptotically flat space - times as those for which the scaling - down of the metric is possible . we arrived at this idea by considering the metric structure only  up to arbitrary scaling  , this is the conformal structure of the space - time manifold in question . by considering the space - time only from the point of view of its conformal structure we obtain a picture of the space - time which is essentially finite but which leaves its causal properties , and hence the properties of wave propagation unchanged . this is exactly what is needed for a rigorous treatment of radiation emitted by the system and also for the numerical simulation of such situations . the way we have presented the emergence of the conformal structure as the essence of asymptotically flat space - times is not how it happened historically . since it is rather instructive to see how various approaches finally came together in the conformal picture we will present in the following section a short overview of the history of the subject . the conformal treatment of infinity has grown out of many roots and it has a history which goes back to the middle of the past century . the natural question arising was how to characterize wave - like solutions of einstein s vacuum field equations ( allowing for a possible cosmological constant  ) : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${g_{ab } } + \lambda { g_{ab } } = 0.$$\end{document } 
 the existence of gravitational waves was predicted by einstein   shortly after he had found the general theory of relativity . however , due to the mathematical and physical complexity of the theory of gravitation there was confusion for a long time about whether the field equations really rigorously do admit solutions with a wave - like character . for instance , rosen   came to the conclusion that there were no such solutions because in the class of plane symmetric waves every non - flat solution of the field equations became singular on a two - dimensional submanifold . this result was , however , due to the lack of understanding of the different kinds of singularities which can occur in a covariant theory and the singularity appearing in the plane wave space - times later turned out to be a coordinate singularity . thus , one of the early problems in the research area of gravitational waves was the invariant characterization of radiation . in 1957 , pirani   started the investigation with the suggestion that the algebraic properties of the riemann ( more specifically the weyl ) tensor should be considered as indicating the presence of radiation . in particular , under the assumption that gravitational radiation can in fact be characterized by the curvature tensor and that it propagates with the local speed of light , he proposed the definition that gravitational radiation is present if the curvature tensor has petrov types { 211 } or { 31}. he arrived at this characterization by the observation that a gravitational wave - front would manifest itself as a discontinuity of the riemann tensor across null hypersurfaces , these being the characteristics for the einstein equations . , pirani did not enforce the field equations , but towards the end of the paper he proposed to look at the equations arbcd = 0 , which follow from the bianchi identity for vacuum space - times . this is the first hint at the importance of the bianchi identity for the study of gravitational waves . in   lichnerowicz the next important step in the development of the subject was trautman s study of the question of boundary conditions for the gravitational field equations . he wanted to obtain conditions which were general enough to allow for gravitational radiation of an isolated system of matter , but still strong enough to guarantee uniqueness ( appropriately defined ) of the solution given he gave an asymptotic fall - off condition for the metric coefficients with respect to a certain class of coordinate systems . it was obtained by analogy to the situation with the scalar wave equation and maxwell theory , where the fields can be required to satisfy sommerfeld s  ausstrahlungsbedingung  . in those cases it is known that there exist unique solutions for given initial data , while in the gravitational case this was not known at the time and , in fact , remained unknown until only quite recently . he defined an energy - momentum four - vector at infinity which is well defined as a consequence of the boundary condition . it is obtained as the limit of an integral of the energy - momentum pseudo - tensor over a space - like manifold with boundary as it stretches out to infinity . by application of stokes theorem , the three - dimensional integral can be converted to a surface integral over the boundary , the sphere at infinity , of certain components of the so - called  superpotentials  for the energy - momentum pseudo - tensor . nowadays these are recognized as certain special cases of the nester - witten two - form   and , hence , trautmans energy - momentum integral coincides with the bondi - sachs expression or the adm expression depending on how the limit to infinity is taken . he considered the difference of energy - momentum between two space - like hypersurfaces and concluded that it must be due to radiation crossing the time - like cylinder which together with the two hypersurfaces bounds a four - dimensional volume . if the limit would have taken out to infinity along null directions , then this result would coincide with the bondi - sachs mass - loss formula . finally , trautman observed that the definitions of pure radiation fields given by pirani and lichnerowicz are obeyed not exactly but only asymptotically by gravitational fields satisfying his boundary condition . hence , he concluded that such solutions approach radiation fields in the limit of large distances to the source . based on the ideas of pirani and trautman and guided by his own investigations   of the structure of retarded linearized gravitational fields of particles , sachs   proposed an invariant condition for outgoing gravitational waves . the intuitive idea was that at large distances from the source , the gravitational field , i.e. the riemann tensor , of outgoing radiation should have approximately the same algebraic structure as does the riemann tensor for a plane wave . sachs analyzed these deviations in detail and obtained rather pleasing qualitative insights into the behaviour of the curvature in the asymptotic regime . in contrast to the earlier work , sachs used more advanced geometrical methods . based on his experience in the study of algebraically special metrics [ 34 , 124 , 86 ] he first analyzed the geometry of congruences of null curves . after the introduction of an appropriate null tetrad he used the bianchi identity for the curvature tensor in a form which follows from the vacuum field equations to obtain the characteristic fall - off behaviour of the curvature components which has been termed the  peeling property  . the riemann tensor of a vacuum space - time has this property , if for any given null geodesic with an affine parameter r which extends to infinity , the curvature falls off along the curve in such a way that to order 1/r it is null ( petrov type n or { 4 } ) with a quadruple pnd ( principal null direction ) along the curve . to order 1/r it has type iii ( { 31 } ) with the triple pnd pointing along the curve , and to order 1/r it has type ii ( { 211 } ) with the double pnd oriented along the curve . to order 1/r it is algebraically general ( type i or { 1111 } ) but one of the pnd s lies in the direction of the geodesic . to order 1/r the curvature symbolically one can express this behaviour in the form 
 ( 1)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c = \frac{{{c_4}}}{r } + \frac{{{c_3}}}{{{r^2 } } } + \frac{{{c_2}}}{{{r^3 } } } + \frac{{{c_1}}}{{{r^4 } } } + o({r^ { - 5}}),$$\end{document } 
 where c on the left hand side denotes the weyl tensor , and each cr on the right hand side stands for a tensor which has an r - fold pnd along the null geodesic and which is independent of r. the important point is that the part of the curvature which has no relation to the null direction of the outgoing geodesic goes as 1/r . sachs postulated the outgoing radiation condition to mean that a bounded source field is free of mixed ( i.e. a non - linear superposition of in- and outgoing ) radiation at large distances if and only if the field has the peeling property . later it was realized [ 26 , 102 ] that this condition does not exclude ingoing radiation . instead it is possible to have an ingoing wave profile provided that it falls off sufficiently fast as a function of an advanced time parameter . the study of gravitational waves and the related questions was the main area of research of the group around h. bondi and f. pirani at king s college , london , in the years between 1955 and 1967 . in a series of papers [ 92 , 93 , 24 , 120 , 94 , 125 , 25 , 127 ] they analyzed several problems related to gravitational waves of increasing complexity . the most important of their results certainly is the work on axisymmetric radiating systems by bondi , van der burg and metzner . they used a different approach to the problem of outgoing gravitational waves . instead of looking at null geodesics they focused on null hypersurfaces , and instead of analyzing the algebraic structure of the curvature using the bianchi identity they considered the full vacuum field equations . their work was concerned with axisymmetric systems but shortly afterwards sachs   removed this additional assumption . this is a scalar function u whose level surfaces are null hypersurfaces opening up towards the future . based on the assumption that such a function exists , one can introduce an adapted coordinate system , so - called bondi coordinates , by labeling the generators of the null hypersurfaces with coordinates on the two - sphere and introducing the luminosity distance r ( essentially the square root of the area of outgoing wave fronts ) along the null generators . the metric , when written in this kind of coordinate system , contains only six free functions . asymptotic conditions were imposed to the effect that one should be able to follow the null geodesics outwards into the future for arbitrarily large values of r. then the metric was required to approach the flat metric in the limit of infinite distance . additionally , it was assumed that the metric functions and other quantities of interest ( in particular the curvature ) were analytic functions of 1/r . the field equations in bondi coordinates have a rather nice hierarchical structure which is symptomatic for the use of null coordinates and which allows for a simpler formal analysis compared to the related cauchy problem . in essence their procedure amounts to the formulation of a certain characteristic initial value problem ( see [ 128 , 129 ] ) and the identification of the free data . it turns out that the freely specifiable data are two functions , essentially components of the metric , on an initial null hypersurface u = u0 and two similar functions at these latter functions are bondi s news functions whose non - vanishing is taken to indicate the presence of gravitational radiation . the most important consequence is the demonstration that outgoing gravitational waves carry away energy from the source and hence diminish its mass . this is the consequence of the bondi - sachs mass loss formula which relates the rate of the mass decrease to the integral over the absolute value of the news . another consequence of the analysis was the peeling property : for space - times which satisfy the vacuum field equations and the bondi - sachs boundary conditions , the curvature necessarily has the asymptotic behaviour ( 1 ) as predicted by sachs direct analysis of the riemann tensor using the vacuum bianchi identity . thus , the bondi - sachs conditions imply the covariant outgoing radiation condition of sachs and also the boundary condition proposed by trautman . the group of coordinate transformations which preserve the form of the metric and the boundary conditions was determined . this infinite dimensional group which became known as the bms group is isomorphic to the semi - direct product of the homogeneous lorentz group with the abelian group of so - called super - translations . the emergence of this group came as a surprise because one would have expected the poincar group as the asymptotic symmetry group but one obtained a strictly larger group . however , the structure of the bms group is sufficiently similar to the poincar group . in particular , it contains a unique abelian normal subgroup of four dimensions , which can be identified with the translation group . this result forms the basis of further investigations into the nature of energy - momentum in general relativity . the bms group makes no reference to the metric which was used to derive it . therefore , it can be interpreted as the invariance group of some universal structure which comes with every space - time satisfying the bondi - sachs boundary conditions . the bms group has been the subject of numerous further investigations since then . for some of them we refer to [ 95 , 96 , 97 , 101 , 111 , 118 , 126 ] . at about the same time it combined the spinor methods which had been developed earlier by penrose   with the ( null-)tetrad calculus used hitherto . in particular , they constructed a coordinate system which was very similar to the ones used by bondi et al . and the only difference was their use of an affine parameter instead of luminosity distance along the generators of the null hypersurfaces of constant retarded time . based on these coordinates and an adapted null frame they showed that the single assumption 0 = o(1/r ) ( and the technical assumption of the uniformity of the angular derivatives ) already implied the peeling property as stated by sachs . the use of 0 as the quantity whose properties are specified on a null hypersurface was in accordance with a general study of characteristic initial value problems for spinor equations and in particular for general relativity undertaken by penrose . an important point in this work was the realization that the bianchi identity could be regarded as a field equation for the weyl tensor . it might be useful here to point out that it is a misconception to consider the bianchi identity as simply a tautology and to ignore it as contributing no further information , as it is done even today . it is an important piece of the structure on a riemannian or lorentzian manifold which relates the ( derivatives of the ) ricci and weyl tensors . if the ricci tensor is restricted by the einstein equations to equal the energy - momentum tensor , then the bianchi identity provides a differential equation for the weyl tensor . its structure is very similar to the familiar zero rest - mass equation for a particle with spin 2 . in fact , in a sense one can consider this equation as the essence of the gravitational theory . newman and unti   carried the calculations started in   further and managed to solve the full vacuum field equations asymptotically for large distances . the condition of asymptotic flatness was imposed not on the metric but directly on the weyl tensor in the form suggested by newman and penrose , namely that the component 0 of the weyl tensor should have the asymptotic behaviour 0 = 00/r+o(1/r ) . from this assumption alone ( and some technical requirements similar to the ones mentioned above ) they obtained the correct peeling behaviour of the curvature , the form of the metric up to the order of 1/r , in particular its flatness at large distances , and also the bondi - sachs mass loss formula . later , the procedure developed by newman and unti to integrate the vacuum field equations asymptotically has been analyzed by dixon , who showed that it can be carried out consistently to all orders of 1/r . it is remarkable how much progress could be made within such a short time ( only about four years ) . the trigger seems to have been the use of the structure of the light cones in one form or another in order to directly describe the properties of the radiation field : the introduction of the retarded time function , the use of an adapted null - tetrad , and the idea to  follow the field along null directions  . schcking had emphasized the conformal invariance of the massless free fields , a fact which had been established much earlier by bateman   and cunningham   for the wave equation and the free maxwell field , and by mclennan   for general spin . this had led to the idea that conformal invariance might play a role also in general relativity and , in particular , in the asymptotic behaviour of the gravitational radiation field ( see   for a personal account of the development of these ideas ) . finally , penrose   outlined a completely different point of view on the subject , arrived at by taking the conformal structure of space - time as fundamental . he showed that if one regarded the metric of minkowski space - time to be specified only up to conformal rescalings g  g for some arbitrary function  , then one could treat points at infinity on the same basis as finite points . minkowski space - time could be completed to a highly symmetrical conformal manifold by adding a  null - cone at infinity  . the well known zero rest - mass fields which transform covariantly under conformal rescalings of the metric are well defined on this space and the condition that they be finite on the null - cone at infinity translates into reasonable fall - off conditions for the fields on minkowski space . on the infinite null - cone one could prescribe characteristic data for the fields which correspond to the strength of their radiation field . he suggested that asymptotically flat space - times should share at least some of these properties . this point of view proved successful . in a further paper by penrose , the basic qualitative picture we have today is developed . roughly speaking ( see the next section for a detailed account ) , the general idea is to attach boundary points to the  physical  space - time manifold which idealize the end - points at infinity reached by infinitely extended null geodesics . the unphysical manifold is equipped with a metric which is conformal to the physical metric with a conformal factor  which vanishes on the boundary . in addition , the structure of the conformal boundary is uniquely determined by the physical space - time . the metric of minkowski space - time in polar coordinates is 
 ( 2)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde g = d{t^2 } - d{r^2 } - { r^2}d{\sigma ^2},$$\end{document } 
 where d is the metric of the unit sphere . to perform the conformal rescaling we introduce null coordinates u = t - r and v = t + r. this puts the minkowski line - element into the form 
 ( 3)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde g = dudv - \frac{1}{4}{(v - u)^2}d{\sigma ^2},$$\end{document } 
 where d is the metric of the unit sphere . the coordinates u and v each range over the complete real line , subject only to the condition v - u  0 . this infinite range is compactified by transforming with an appropriate function , e.g. 
 ( 4)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$u = \tan u,\:\;\;\;\;\;v = \tan v,$$\end{document } 
 thus introducing new null coordinates u and v , in terms of which the metric takes the form 
 ( 5)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde g = \frac{1}{{4{{\cos } ^2}u{{\cos } ^2}v}}\ { 4du\;dv - { \sin ^2}(v - u)d{\sigma ^2}\ } .$$\end{document } 
 the coordinates u , v both range over the open interval ( -/2 , /2 ) with the restriction v - u  0 . obviously , the minkowski metric is not defined at points with u = /2 or v = /2 . now we define a different metric g = g , conformally related to g by the conformal factor  = 2 cosu cosv . thus , 
 ( 6)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g = 4du\;dv - { \sin ^2}(v - u)d{\sigma ^2}.$$\end{document } 
 this metric is perfectly regular at the points mentioned above and , in fact , g is the metric of the einstein cylinder { \mathcal e}. this can be verified by defining an appropriate time and radius coordinate . with t = u + v , r = v - u we have 
 ( 7)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g = d{t^2 } - d{r^2 } - { \sin ^2}rd{\sigma ^2}.$$\end{document } thus , we may consider minkowski space to be conformally embedded into the einstein cylinder . the minkowski metric determines the structure of the boundary , namely the two 3-dimensional null - hypersurfaces  and  which represent ( future and past )  null - infinity  . they are given by the conditions v = /2,|u| < /2 (  ) and u = -/2,|v| < /2 (  ) . they represent  future and past time - like infinity  , the start respectively the end point of time - like geodesics , while i is a point with u = -/2 , v = /2 . it is the start and end point of all space - like geodesics , hence it represents  spacelike infinity  . figure 1the embedding of minkowski space into the einstein cylinder 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal e}$$\end{document}. the embedding of minkowski space into the einstein cylinder 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal e}$$\end{document}. the conformal boundary of minkowski space - time consists of the pieces  , i , and i. these are fixed by the minkowski metric . in contrast to this , the conformal manifold into which minkowski space - time is embedded ( here the einstein cylinder ) is not fixed by the metric . obviously , had we chosen a different conformal factor  =  with some arbitrary positive function  we would not have obtained the metric of the einstein cylinder but a different one . we see from this that although the conformal boundary is unique , the conformal extension beyond the boundary is not . first of all it simplifies the discussion of problems at infinity which would involve complicated limit procedures when viewed with respect to the physical metric . transforming to the unphysical metric attaches a boundary to the manifold so that issues which arise at infinity with respect to the physical metric can be analyzed by local differential geometric arguments in the neighbourhood of the boundary . this is particularly useful for discussing solutions of conformally invariant field equations on space - time . the basic idea is the following : consider a space - time which allows us to attach a conformal boundary , thus defining an unphysical manifold conformally related to the given space - time . suppose we are also given a solution of a conformally invariant equation on this unphysical manifold . because of the conformal invariance of the equation , there exists a rescaling of that unphysical field with a power of the conformal factor , which produces a solution of the equation on the physical manifold . then the physical solution will have a characteristic asymptotic behaviour which is entirely governed by the conformal weight of the field , i.e. by the power of the conformal factor used for the rescal - ing . thus , the regularity requirement of the unphysical field translates into a characteristic asymptotic fall - off or growth behaviour of the physical field , depending on its conformal weight . penrose used this idea to show that solutions of the zero rest - mass equations for arbitrary spin on a space - time , which can be compactified by a conformal rescaling , exhibit the peeling property in close analogy to the gravitational case as discovered by sachs . take as an example the spin-2 zero rest - mass equation for a tensor kabcd with the algebraic properties of the weyl tensor 
 ( 8)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^a}{k_{abcd } } = 0.$$\end{document } 
 this is the equation for linear perturbations of the gravitational fields propagating on a fixed background . it is conformally covariant , in the sense that it remains unchanged provided the field is rescaled as kabcd  kabcd , i.e. , it is a conformal density with weight -1 . using the geometric technique of conformal compactification , penrose was able to establish the peeling property also for general ( non - linear ) gravitational fields . furthermore , he showed that the group of transformations of the conformal boundary leaving the essential structure invariant was exactly the bms group . this geometric point of view suggested that the asymptotic behaviour of the gravitational field of an isolated radiating gravitational system can be described entirely in terms of its conformal structure . the support for this suggestion was overwhelming from an aesthetical point of view , but a rigorous support for this claim was provided essentially only from the examination of the formal expansion type solutions of bondi - sachs and newman - unti and the analysis of explicit stationary solutions of the field equations . the geometric point of view outlined above is the foundation on which many modern developments within general relativity are based . let us now discuss the notion of asymptotically flat space - times and some of their properties in more detail . we have seen in section 2.1 that the question of how to define isolated systems in general relativity has led to the mathematical idealization of asymptotically flat space - times . they are defined by the requirement that they allow the attachment of a smooth conformal boundary . the precise definition is : 
 definition 1 
 a smooth ( time- and space - orientable ) space - time ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , gab ) is called asymptotically simple , if there exists another smooth lorentz manifold ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , gab ) such that 
 ( i)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } is an open submanifold of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } with smooth boundary \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \tilde { \mathcal m } = i;$$\end{document}(ii)there exists a smooth scalar field  on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , such that gab = gab 
 on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , and so that  = 0 , d  0 on ;(iii)every null geodesic in 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } 
 acquires a future and a past endpoint on . 
 an asymptotically simple space - time is called asymptotically flat , if in addition rab = 0 in a neighbourhood of . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } is an open submanifold of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } with smooth boundary \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \tilde { \mathcal m } = i;$$\end{document } there exists a smooth scalar field  on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , such that gab = gab 
 on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , and so that  = 0 , d  0 on  ; every null geodesic in 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } 
 acquires a future and a past endpoint on . thus , asymptotically flat space - times are a subclass of asymptotically simple space - times , namely those for which the einstein vacuum equations hold near . examples for asymptotically simple space - times which are not asymptotically flat include the de sitter and anti - de sitter space - times , both solutions of the einstein equations with non - vanishing cosmological constant . we will concentrate here on asymptotically flat space - times . according to the first condition , the space - time ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , gab ) , which we call the physical space - time can be considered as part of a larger space - time ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , gab ) , the unphysical space - time . as a submanifold of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } the physical space - time can be given a boundary which is required to be smooth . the unphysical metric gab is well - defined on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } and , in particular , on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , while the physical metric gab is only defined on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } and can not be extended in a well - defined sense to the boundary of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } or even beyond . the metrics generate the same conformal structure , they are conformally equivalent in the sense that on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } they define the same null - cone structure . note that although the extended manifold \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } and its metric are called unphysical , there is nothing unphysical about this construction . as we shall see below , the boundary of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } is uniquely determined by the conformal structure of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } and , therefore , it is just as physical as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. the extension beyond the boundary , given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } is not unique , as we have already seen in section 2.2 , but this is of no consequence for the physics in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } because the extension is causally disconnected from \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. the second condition fixes the behaviour of the scaling factor on  as being  of the order 1/r  as one approaches  from within \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. the condition ( iii ) is a completeness condition to ensure that the entire boundary is included . in some cases of interest , this condition is not satisfied . in the schwarzschild space - time , for instance , there are null - geodesics which circle around the singularity , unable to escape to infinity . this problem has led to a weakening of definition 1 to weakly asymptotically simple space - times ( see e.g. ) . such space - times are essentially required to be isometric to an asymptotically simple space - time in a neighbourhood of the boundary . a different completeness condition has been proposed by geroch and horowitz . in the following discussion of the analytic and geometric issues , weakly asymptotically simple space - times will not play a role so that we can assume our space - times to be asymptotically simple . of course , for applications weakly asymptotically simple space - times are important because they provide interesting examples of space - times with black holes . we defined asymptotically flat space - times by the requirement that the einstein vacuum equation holds near the boundary , i.e. , that asymptotically the physical space - time is empty . there are ways to relax this condition by imposing strong enough fall - off conditions on the energy - momentum tensor without violating any of the consequences . since we are concerned here mainly with the asymptotic region , we are not really interested in including any matter fields . therefore , we will assume henceforth that the physical space - time is a vacuum space - time . this does not mean that the following discussion is only valid for vacuum space - times , it simply allows us to make simpler statements . the conformal factor  used to construct the boundary  is , to a large extent , arbitrary . this raises the important question about the uniqueness of the conformal boundary as a point set and as a differential manifold . if this uniqueness were not present , then the notion of  points at infinity it could then happen that two curves which approach the same point in one conformal boundary for a space - time reach two different points in another conformal completion . or , similarly , that two conformal extensions which arise from two different conformal factors were not smoothly related . in fact , it can be shown that between two smooth extensions there always exists a diffeomorphism which is the identity on the physical space - time , so that the two extensions are indistinguishable from the point of view of their topological and differential structure . it also follows from schmidt s so called b - boundary construction [ 131 , 132 , 134 ] . from the condition that the vacuum einstein equation holds , one can derive several important consequences for asymptotically flat space - times : 
  is a smooth null hypersurface in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. is shear - free. has two connected components , each with topology s  .the conformal weyl tensor vanishes on . 
 the first part of statement ( a ) follows from the fact that  is given by the equation  = 0 . since  has a non - vanishing gradient on  , regularity follows . furthermore , from the einstein vacuum equations one has  = 0 on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. hence , equation ( 111 ) implies on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\omega^{2}\lambda-\frac{1}{4}\omega\square \omega+\frac{1}{2}\nabla^{a}\omega\nabla_{a}\omega=0.$$\end{document } 
 this equation can be extended smoothly to the boundary of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , yielding there the condition nan = 0 for the co - normal na = -a of . hence , the gradient of the conformal factor is null , and  is a null hypersurface . is a smooth null hypersurface in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}.  has two connected components , each with topology s  . the conformal weyl tensor vanishes on . as such it is generated by null geodesics . the statement ( b ) asserts that the congruence formed by the generators of  has vanishing shear . to show this we look at equation ( 110 ) and find from ab = 0 that 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\omega\phi_{ab}+\nabla_{a}\nabla_{b}\omega-\frac{1}{4}g_{ab}\square \omega=0,$$\end{document } 
 whence , on  we get ( writing mab for the degenerate induced metric on  ) 
 ( 9)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$2\displaystyle \nabla_{a}n_{b}=\mathcal{l}_{n}m_{ab}=-\frac{1}{2}m_{ab}\square \omega,$$\end{document } 
 whence the lie - derivative of mab along the generators is proportional to mab , which is the shear - free condition for null geodesic congruences with tangent vector n ( see [ 75 , 118 ] ) . to prove statement ( c ) we observe that since  is null , either the future or the past light cone of each of its points has a non - vanishing intersection with \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. this shows that there are two components of  , namely  on which null geodesics attain a future endpoint , and  where they attain a past endpoint . these are the only connected components because there is a continuous map from the bundle of null - directions over \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } to  , assigning to each null direction at each point p of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } the future ( past ) endpoint of the light ray emanating from p in the given direction . if  were not connected then neither would be the bundle of null - directions of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , which is a contradiction ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } being connected ) . to show that the topology of  is s   requires a more sophisticated argument which has been given by penrose   ( a different proof has been provided by geroch ) . it has been pointed out by newman   that these arguments are only partially correct . he rigorously analyzed the global structure of asymptotically simple space - times and he found that , in fact , there are more general topologies allowed for . however , his analysis was based on methods of differential topology not taking the field equations into account . indeed , we will find later in theorem 6 that the space - time which evolves from data close enough to minkowski data will have a  with topology s  . the proof of statement ( d ) depends in an essential way on the topological structure of . we refer again to . the vanishing of the weyl curvature on  is the final justification for the definition of asymptotically flat space - times : anishing ricci curvature implies the vanishing of the weyl tensor and hence of the entire riemann tensor on . the physical space - time becomes flat at infinity . but there is another important property which follows from the vanishing of the weyl tensor on . consider the weyl tensor cbcd of the unphysical metric gab which agrees on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } with the weyl tensor cbcd of the physical metric gab because of the conformal invariance ( 107 ) . on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , cbcd satisfies the vacuum bianchi identity 
 ( 10)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde \nabla _ a}{\tilde c^a}_{\;bcd } = 0.$$\end{document } 
 this equation looks superficially like the zero rest - mass equation ( 8) for spin-2 fields . however , the conformal transformation property of ( 10 ) is different from the zero rest - mass case . the equation is not conformally invariant since the conformal rescaling of a vacuum metric generates ricci curvature in the unphysical space - time by equation ( 108 ) , which then feeds back into the weyl curvature via the bianchi identity ( cf . equation ( 112 ) ) . however , we can define the field 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$k_{\;\;bcd}^a = { \omega ^ { - 1}}c_{\;\;bcd}^a$$\end{document } 
 on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. as it stands , kbcd is not defined on . but the vanishing of the weyl tensor there and the smoothness assumption allow the extension of kbcd to the boundary ( and even beyond ) as a smooth field on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. it follows from equation ( 10 ) that this field satisfies the zero rest - mass equation 
 ( 11)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a}k_{bcd}^a = 0$$\end{document } 
 on the unphysical space - time \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } with respect to the unphysical metric . therefore , the rescaled weyl tensor kbcd is a genuine spin-2 field with the natural conformal behaviour . in fact , this is the field which most directly describes the gravitational effects , in particular its values on the boundary are closely related to the gravitational radiation which escapes from the system under consideration . it propagates on the conformal space - time in a conformally covariant way according to equation ( 11 ) which looks superficially like the equation ( 8) for a ( linear ) spin-2 zero rest - mass field . however , there are highly non - linear couplings between the connection given by a and the curvature given by kabcd . in the physical space - time , where the conformal factor is unity , the field kabcd coincides with the weyl tensor which is the source of tidal forces acting on test particles moving in space - time . for these reasons ( 11 ) and the regularity on  follows a specific fall - off behaviour of the field kbcd and hence of the weyl tensor which is exactly the peeling property obtained by sachs . it arises here from a reasoning similar to the one presented towards the end of section 2.2 . it is a direct consequence of the geometric assumption that the conformal completion be possible and of the conformal invariance of equation ( 11 ) . this equation for the rescaled weyl tensor is an important sub - structure of the einstein equation because it is conformally invariant in contrast to the einstein equation itself . in a sense it is the most important part also in the system of conformal field equations which we consider in the next section . the possibility of conformal compactification restricts the lowest order structure of the gravitational field on the boundary . this means that all asymptotically flat manifolds are the same in that order , so that the conformal boundary and its structure are universal features among asymptotically flat space - times . this is nicely illustrated by the weyl tensor which necessarily vanishes on the conformal boundary , but the values of the rescaled weyl tensor kabcd = cabcd are not fixed there . in summary , our qualitative picture of asymptotically flat space - times is as follows : such space - times are characterized by the property that they can be conformally compactified . more importantly , these points together form a three - dimensional manifold which is smoothly embedded into a larger extended space - time . smoothness of the resulting manifold with boundary translates into asymptotic fall - off conditions for the physical metric and the fields derived from it . the boundary emerges here as a geometric concept and not as an artificial construct put in by hand . boundary condition at infinity  as we and the early works sometimes did . in this subsection we continue our discussion of the prototype of an asymptotically flat space - time , namely minkowski space - time . the motivation for doing so is partly to get more acquainted with the idea of conformal compactification and partly to show why this concept is the correct one for the description of radiation processes . for the following discussion we refer to the rescaled metric 
 ( 12)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g = 4du\,dv - { \sin ^2}(v - u)d{\sigma ^2},$$\end{document } 
 which is the metric of the einstein cylinder . the long side of the triangle consists of all the points in the centre , r = 0 ( i.e. u = v ) . the other two sides of the triangle correspond to null - infinity . the points i are points in the centre with u = v = /2 , while i is a point with u = -/2 , v = /2 . the lines meeting at i are lines of constant t , while the lines emanating from i and converging into i are lines of constant r. in the four - dimensional space - time the lines of constant t correspond to three - dimensional space - like hypersurfaces which are asymptotically euclidean . in the case of minkowski space - time , the metric can be extended in a regular way to three points representing future and past time - like infinity and space - like infinity , but this is not generally so . already in the case of the schwarzschild metric , which is also an asymptotically flat metric , there are , strictly speaking , no such points because any attempt to extend the metric yields a singularity . however , it is common language to refer to this behaviour by saying that  the points i and/or i are singular  . the reason for this is related to the presence of mass . for any space - time which has a non - vanishing adm - mass , the point i is necessarily singular while the singularity of the time - like infinities is , in general , related to the fall - off properties of the energy - momentum tensor in time - like directions . in the case of the schwarzschild solution ( like in any stationary black - hole solution ) it is the presence of the static ( stationary ) black hole which is responsible for the singularity of i. let us now assume that there is a particle which moves along the central world - line r = 0 , emitting radiation . for the sake of simplicity we assume that it emits electro - magnetic radiation which travels along the outgoing null - cones to null - infinity . the null - cones are symbolized in figure 2 by the straight lines going off the particle s world - line at 45. we are now interested in the behaviour of the signal along various space - like hypersurfaces . in figure 3 the generic features discussed below will be the same for any asymptotically flat space - time as long as we stay away from the corners of the diagram . the vertical dashed line is the world - line of the particle which defines the time axis . we have displayed two asymptotically euclidean space - like hyper - surfaces intercepting the time axis at two different points and reaching out to space - like infinity . furthermore , there are two hypersurfaces which intersect the time axis in the same two points as the asymptotically flat ones . they reach null - infinity , intersecting in a two - dimensional space - like surface . this geometric statement about the behaviour of the hypersurfaces in the unphysical space - time translates back to the physical space - time as a statement about asymptotic fall - off conditions of the induced ( physical ) metric on the hypersurfaces , namely that asymptotically the metric has constant negative curvature . this is , in particular , a property of the space - like hyperboloids in minkowski space . thus , such hypersurfaces are called hyperboloidal hypersurfaces . an important point to keep in mind is that the conformal space - time does not  stop  at  but that it can be extended smoothly beyond . the extension is not uniquely determined as we have already discussed in connection with the embedding of minkowski space into the einstein cylinder ( cf . thus , the extension plays no role for the concept of null - infinity but it can be very helpful for technical reasons , in particular when numerical issues are discussed . we now imagine that the central particle radiates electro - magnetic waves of uniform frequency , i.e. proportional to sin( ) , where  is the particle s proper time . this gives rise to a retarded electro - magnetic field on the entire space - time which has the form   sin(u ) , where u is a retarded time coordinate on minkowski space with u =  on the central world - line . we ignore the fall - off of the field because it is irrelevant for our present purposes . let us now look at the waves on the various hypersurfaces . in the physical space - time , the hypersurfaces extend to infinity and we can follow the waves only up to an arbitrary but finite distance along the hypersurfaces . these are surfaces of constant minkowski time which implies that the signal is again a pure sine . note , however , that this is only true for these special hypersurfaces . even in minkowski space - time we could choose space - like hypersurfaces which are not surfaces of constant minkowski time but which nonetheless are asymptotically euclidean . on such surfaces the wave would look completely different . on the reason for the decrease in frequency is the fact that these surfaces tend to become more  characteristic  as they extend to infinity , thus approaching surfaces of constant phase of the retarded field . the final diagram shows the signal obtained by an idealized observer which moves along the piece of  between the two intersection points with the hyperboloidal surfaces . the signal is recorded with respect to the retarded time u which , in the present case , is a so called bondi parameter ( see section 4.3 ) . therefore , the observer measures a signal at a single frequency for a certain interval of this time parameter . a different bondi time would result in a signal during a different time interval but with a single , appropriately scaled , frequency . using an arbitrary time parameter would destroy the feature that only one frequency is present in the signal . this is , in fact , the only information that can be transmitted from the emitter to the receiver under the given circumstances . the wave - forms of the signal as they appear in the conformal space - time , i.e. with respect to a coordinate system which covers a neighbourhood of  , are shown in figure 5 . in the specific case of minkowski space - time the signal on the asymptotically euclidean surfaces shows the  piling up  of the waves as they approach spacelike infinity . since the field and the surfaces are both smooth across  , the signal can continue on across null - infinity without even noticing its presence . the points where  is crossed are indicated in the diagram as two little crosses . the values of the field at these points are the same as the boundary values of the signal in the third diagram . here the signal on the same region of  as in figure 4 is displayed , but with respect to the coordinate u = arctan(u ) which is not a bondi parameter . what these diagrams teach us is the following : it has been convenient in relativity to decompose space - time into space and time by slicing it with a family of space - like hypersurfaces . in most of the work on existence theorems of the einstein equations it has been convenient to choose them to be cauchy surfaces and thus asymptotically euclidean . also , in most numerical treatments of einstein s equations the same method is used to evolve space - times from one space - like hypersurface to the next ( see section 4 ) . here the hypersurfaces used are finite because the numerical grids are necessarily finite . in the approaches based on the standard einstein equations it makes no difference whether the grid is based on a finite portion of an asymptotically euclidean or a hyperboloidal hypersurface . the fact that the space - time should be asymptotically flat has to be conveyed entirely by a suitable boundary condition which has to be imposed at the boundary of the finite portion of the hypersurface ( i.e. at the little crosses in figure 3 ) . however , this implies that the accuracy of the wave - form templates obtained with such approaches depends to a large extent on the quality of that boundary condition . so far there exists no suitable boundary condition which would be physically reasonable and lead to stable codes . in the conformal approach one has the option to  include infinity  by using the conformal field equations ( see section 3 ) . the diagrams show that the hyperboloidal surfaces are very well suited to deal with the radiation problems . they provide a foliation of the conformal space - time on which one can base the evolution with the conformal field equations . the solution obtained will be smooth near  and we  only  need to locate  on each hypersurface to read off the value of the radiation data ( as indicated in the second diagram of figure 5 ) . as we indicated already in the introduction , the amount of results and developments related to the conformal structure of space - times and , in particular , to conformal infinity is overwhelming . we take the opportunity to refer to several other developments which are not treated here in detail . there exist several reviews of the subject from different points of view , e.g. by geroch , by penrose , by schmidt , by newman and tod , by ashtekar [ 6 , 7 ] , and by friedrich [ 53 , 56 ] . a large part of the literature on null - infinity is concerned with  conserved quantities  . there it can be defined in terms of limits of integrals , called linkages , over spheres which approach a cut of null - infinity  , where the integrals are taken over certain vectorfields in the physical space - time which suitably approximate the infinitesimal generators of asymptotic symmetries . penrose , who had earlier   reexpressed the original bondi - sachs expressions in terms of genuine geometric quantities at  , has also derived them from his quasi - local mass proposal . helicity lowering  of the rescaled weyl tensor at  using a two - index asymptotic twistor . other approaches ( see   for a review ) start from a hamiltonian or lagrangian formulation of the theory and derive the energy - momentum expressions via noether theorems or the moment - map of symplectic geometry ( see e.g. [ 9 , 12 ] ) . these formulations also provide a framework for  asymptotic quantization  , a scheme which is geared towards a scattering - matrix description for gravity . the universal structure of  provides the necessary background structure for the definition of a phase - space of the radiative modes of the gravitational field and its subsequent quantization . while the energy - momentum expressions all coincide , there is still disagreement about the various angular - momentum expressions ( see e.g. the review article by winicour ) . this difficulty is caused by the group structure of the bms group which does not allow to single out a unique lorentz subgroup ( it is obtained only as a factor group ) . hopefully these discrepancies will be resolved once the structure of the gravitational fields at i is completely understood . all the  conserved quantities  are associated with a ( space - like ) cut of null - infinity which is used for evaluation of the surface integrals , and an infinitesimal generator of the asymptotic symmetry group used in defining the integrand . the prime example is again the bondi - sachs energy - momentum , which obeys the famous bondi - sachs mass - loss formula which relates the values of the energy - momentum at two given cuts with a negative definite  flux integral  over the part of  between the two cuts . furthermore , there exist the somewhat mysterious newman - penrose constants , five complex quantities which are also defined by surface integrals over a cut of . in contrast to the previous conserved quantities , the np constants are absolutely conserved in the sense that they do not depend on the particular cut which is used for the evaluation of the integrals . in space - times which have a regular point i , the np constants turn out to be the value of the gravitational field at i. if i is singular , then the np constants are still well - defined , although now they should probably be considered as the value of the gravitational field at an ideal point i. other interpretations relate them to certain combinations of multipole moments of the gravitational field [ 102 , 118 ] . people have tried to give an interpretation of the np constants in terms of a lagrangian or symplectic framework [ 71 , 70 , 122 ] , but these results are still somewhat unsatisfactory . very recently , friedrich and knnr   were able to connect the np constants defined at null - infinity to initial data on a space - like asymptotically euclidean ( time - symmetric ) hypersurface . finally , we want to mention the recent formulation of general relativity as a theory of null hypersurfaces , see . this theory has its roots in the observation that one can reconstruct the points of minkowski space - time from structures defined on null - infinity . the future light - cone emanating from an arbitrary point in minkowski space - time is a shear - free null hypersurface intersecting  in a cut . the shear - free property of the light - cone translates into the fact that the cut itself is given as a solution of a certain differential equation , the  good cut equation  on . conversely , it was realized that in flat space the solution space of the good cut equation is isometric to minkowski space - time ( in particular , it carries a flat metric ) . attempts to generalize this property led to newman s -space construction   which associates with each ( complexified ) asymptotically flat and ( anti-)self - dual space - time a certain complex four - dimensional manifold which carries a ricci - flat metric . it is obtained as the solution space of the complex good cut equation . trying to avoid the unphysical complexification has finally led to the above mentioned null surface formulation of general relativity . at this point construction   which associates with each anti self - dual vacuum space - time a certain three - dimensional complex manifold . the interpretation of these manifolds at the time was that they should provide the one - particle states of the gravitational field in a future quantum theory of gravity . the non - linear gravitons themselves have led to remarkable developments in pure mathematics ( see e.g. the contributions in ) . it is clear that penrose s method of conformally compactifying space - time does provide a convenient and elegant framework for discussing questions of asymp - totics . what is not obvious , however , is how to answer the question of existence of space - times which actually do possess the asymptotic structure suggested by the conformal picture . on the one hand , a very specific geometric property of the space - time is required , namely the possibility to attach a smooth conformal boundary , and on the other hand , the einstein equation , a differential equation on the metric of the space - time , has to be satisfied . it is not clear whether these two different requirements are indeed compatible . to be more precise , one would like to know how many asymptotically flat solutions of the field equations actually exist , and whether this class of solutions contains the physically interesting ones which correspond to radiative isolated gravitating systems . the only possible avenue to answering questions of existence consists of setting up appropriate initial value problems and proving existence theorems for solutions of the einstein equations subject to the appropriate boundary conditions . the conventional cauchy problem which treats the einstein equation as a second order partial differential equation for the metric field is already rather complicated by itself ( see e.g. the review by choquet - bruhat and york ) . but to obtain statements of the type mentioned above is further complicated by the fact that in order to discuss the asymptotic fall - off properties of solutions one would need to establish global ( long time and large distance ) existence together with detailed estimates about the fall - off behaviour of the solution . the geometric characterization of the asymptotic conditions in terms of the conformal structure suggests to discuss also the existence problem in terms of the conformal structure . the general idea is as follows : suppose we are given an asymptotically flat manifold which we consider as being embedded into an appropriate conformally related unphysical space - time . the einstein equations for the physical metric imply conditions for the unphysical metric and the conformal factor relating these two metrics . it turns out that one can write down equations which are regular on the entire unphysical manifold , even at those points which are at infinity with respect to the physical metric . existence of solutions of these  regular conformal field equations  on the conformal manifold then translate back to ( semi-)global results for asymptotically flat solutions of the field equations in physical space . this approach towards the existence problem has been the programme followed by friedrich since the late 1970 s . in this section we will discuss the conformal field equations which have been derived by friedrich , and the various subproblems which have been successfully treated by using the conformal field equations . before deriving the equations we need to define the arena where the discussion is taking place . definition 2 
 a conformal space - time is a triple ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , gab ,  ) such that 
 ( i)(\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , gab ) is a ( time- and space - orientable ) lorentz manifold;(ii) is a smooth scalar field on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } 
 such that the set 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m } = \ { p \in { \mathcal m}:\omega ( p ) > 0\ } $ $ \end{document } 
 is non - empty and connected;(iii)the gravitational field kbcd = cbcd extends smoothly to all of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. 
 two conformal space - times ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , gab ,  ) and ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , ab , ^ ) are equivalent if 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } 
 and 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\hat { \mathcal m}$$\end{document } 
 are diffeomorphic and if , after identification of 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } 
 and 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } 
 with a suitable diffeomorphism , there exists a strictly positive scalar field  on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } 
 such that ^ =  and ^gab = gab . ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , gab ) is a ( time- and space - orientable ) lorentz manifold ;  is a smooth scalar field on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } 
 such that the set 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m } = \ { p \in { \mathcal m}:\omega ( p ) > 0\ } $ $ \end{document } 
 is non - empty and connected ; the gravitational field kbcd = cbcd extends smoothly to all of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. from this definition follows that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } is an open sub - manifold of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } on which a metric gab = gab is defined , which is invariant in the sense that two equivalent conformal space - times define the same metric gab . the space - time ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , gab ) allows the attachment of a conformal boundary which is given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$i = \ { p \in { \mathcal m}:\omega ( p ) = 0,\;d\omega   \ne 0\ } $ $ \end{document}. the above definition of conformal space - times admits much more general situations than those arising from asymptotically flat space - times ; this generality is sometimes needed for numerical purposes . under the conditions of definition 2 , it follows that the weyl tensor vanishes on  because the gravitational field ( i.e. the rescaled weyl tensor ) is smooth on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. note that we make no assumptions about the topology of . if each null geodesic which starts from the inside of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } has a future and a past endpoint on  , then \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } is asymptotically simple in the sense of definition 1 . if , in addition , the metric gab is a vacuum metric then  has the implied topology s  . note also that it is quite possible to have situations , where gab is a vacuum metric and where the topology of  is not s   , but e.g. t  . then , necessarily , there must exist null geodesics which do not reach . in the special case when  is empty , the conformal factor  is strictly positive , i.e. \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m } = { \mathcal m}$$\end{document } , and the conformal space - time is isometric to the physical space - time ( choosing  =  ) . our goal is to express the vacuum equations in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } in terms of geometric quantities on the unphysical space - time . when expressed in terms of unphysical quantities it reads ( see the formulae of appendix 7 ) 
 ( 13)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$0=\displaystyle \tilde{g}_{ab}=g_{ab}-\frac{2}{\omega}(\nabla_{a}\nabla_{b}\omega - g_{ab}\square \omega)-\frac{3}{\omega^{2}}g_{ab}\nabla_{c}\omega\nabla^{c}\omega.$$\end{document } 
 this equation can be interpreted as the einstein equation for the metric gab with a source term which is determined by the conformal factor . if we assume  to be known , then it is a second order equation for gab , which is formally singular on  , where  vanishes . therefore , it is very hard to make any progress towards the existence problem using this equation . to remedy this situation , friedrich [ 43 , 44 , 45 ] suggested to consider a different system of equations on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } which can be derived from the geometric structure on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , the conformal transformation properties of the curvature and the vacuum einstein equation on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. it consists of equations for a connection a , its curvature and certain other fields obtained from the curvature and the conformal factor . let us assume that a is a connection on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } which is compatible with the metric gab so that 
 ( 14)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ c}{g_{ab } } = 0$$\end{document } 
 holds . we will write down equations for the following unknowns : 
  the connection a, the conformal factor  , a one - form a and a scalar function s, a symmetric trace - free tensorfield ab , and a completely trace - free tensorfield kbcd which has the symmetries of the weyl tensor . we introduce the zero - quantity 
 ( 15)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$z = ( t_{bc}^a,\:{{\mathcal q}_{abcd}},\:{{\mathcal p}_{abc}},\:{{\mathcal b}_{bcd}},\:{{\mathcal s}_a},\:{{\mathcal s}_{ab}},\:{{\mathcal d}_a}),$$\end{document } 
 where tbca is the torsion tensor of a and the other components of z are defined in terms of the unknowns by 
 ( 16)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal q}_{abcd } } \equiv { r_{abcd } } - \omega { k_{abcd } } + 2{g_{c[a}}{\phi _ { b]d } } - 2{g_{d[a}}{\phi _ { b]c } } - 4{g_{c[a}}{g_{b]d}}\lambda , $ $ \end{document } 
 ( 17)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal p}_{abc } } \equiv 2{\nabla _ { [ c}}{\phi _ { a]b } } + 2{g_{b[c}}{\nabla _ { a]}}\lambda   - { k_{cab}}^d{\sigma _ d},$$\end{document } 
 ( 18)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal b_{bcd } } \equiv { \nabla _ a}{k^a}_{bcd},$$\end{document } 
 ( 19)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_a } \equiv { \nabla _ a}\omega   - { \sigma _ a},$$\end{document } 
 ( 20)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_{ab } } \equiv { \nabla _ a}{\sigma _ b } - { g_{ab}}s + \omega { \phi _ { ab}},$$\end{document } 
 ( 21)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal d}_a } \equiv { \nabla _ a}s + { \phi _ { ab}}{\sigma ^b } - \omega { \nabla _ a}\lambda   - 2\lambda { \sigma _ a}.$$\end{document } 
 in addition , we consider the scalar field 
 ( 22)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t } \equiv 2\omega s - 2{\omega ^2}\lambda   - { \sigma _ a}{\sigma ^a}$$\end{document } 
 on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. the equations z = 0 are the regular conformal vacuum field equations . in contrast to equation ( 13 ) this system is regular1 on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , even on  because there are no terms containing .  the conformal factor  , a one - form a and a scalar function s ,  a symmetric trace - free tensorfield ab , and  a completely trace - free tensorfield kbcd which has the symmetries of the weyl tensor . consider the equation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal b}_{bcd } } = { \nabla _ a}{k^a}_{bcd } = 0$$\end{document}. this subsystem lies at the heart of the full system of conformal field equations because it feeds back into all the other parts . it was pointed out in section 2.2 that the importance of the bianchi identity had been realized by sachs . however , it was first used in connection with uniqueness and existence proofs only by friedrich [ 45 , 44 ] . its importance lies in the fact that it splits naturally into a symmetric hyperbolic system of evolution equations2 and constraint equations . energy estimates for the symmetric hyperbolic system naturally involve integrals over a certain component of the bel - robinson tensor , a well known tensor in general relativity which has certain positivity properties . the usefulness of the conformal field equations is documented in 
 theorem 1 
 suppose that a 
 is compatible with gab 
 and that z = 0 on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document}. if 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t } = 0$$\end{document } 
 at one point of 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , then 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t } = 0$$\end{document } 
 everywhere and , furthermore , the metric gab 
 is a vacuum metric on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document}. 
 proof : from the vanishing of the torsion tensor it follows that a is the levi - civita connection for the metric gab . then , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal q}_{abcd } } = 0$$\end{document } is the decomposition of the riemann tensor into its irreducible parts which implies that the weyl tensor cbcd = kbcd , that ab is the trace - free part of the ricci tensor , and that  = 24r . the equation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_a } = 0$$\end{document } defines a in terms of  , and the trace of the equation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_ab } = 0$$\end{document } defines s = 1/4. the trace - free part of that equation is the statement that ab = 0 , which follows from the conformal transformation property ( 110 ) of the trace - free ricci tensor . with these identifications the equations \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal b}_{abc } } = 0$$\end{document } resp . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal p}_{abc } } = 0$$\end{document } do not yield any further information because they are identically satisfied as a consequence of the bianchi identity on ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , g ) , resp . ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } , g ) . finally , we consider the field \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t}$$\end{document}. taking its derivative and using \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_{ab } } = 0$$\end{document } and 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal d}_a } = 0$$\end{document } , we obtain \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a}{\mathcal t } = 0$$\end{document}. hence , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t}$$\end{document } vanishes everywhere if it vanishes at one point . it follows from the transformation ( 111 ) of the scalar curvature under conformal rescalings that \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t } = 0$$\end{document } implies  = 0 . it is easy to see that the conformal field equations are invariant under the conformal rescalings of the metric specified in definition 2 and the implied transformation of the unknowns . the conformal invariance of the system implies that the information it contains depends only on the equivalence class of the conformal space - time . the reason for the vanishing of the gradient of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t}$$\end{document } is essentially this : if we impose the equation ab = 0 for the trace - free part of the ricci tensor of a manifold , then by use of the contracted bianchi identity we obtain a = 0 . expressing this in terms of unphysical quantities leads to the reasoning in theorem 1 . the special case  = 1 reduces to the standard vacuum einstein equations , because then we have kabcd = cabcd and a = 0 . then \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal s}_{ab } } = 0$$\end{document } implies ab = 0 and s = 0 , while \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal t } = 0$$\end{document } forces  = 0 . the other equations are identically satisfied . given a smooth solution of the conformal field equations on a conformal manifold , theorem 1 implies that on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde { \mathcal m}$$\end{document } we obtain a solution of the vacuum einstein equation . in particular , since the weyl tensor of gab vanishes on  due to the smoothness of the gravitational field , this implies that the weyl tensor has the peeling property in the physical space - time . therefore , if existence of suitable solutions of the conformal field equations on a conformal manifold can be established , one has automatically shown existence of asymptotically flat solutions of the einstein equations . the main advantage of this approach is the fact that the conformal compactification supports the translation of global problems into local ones . note that the use of the conformal field equations is not limited to vacuum space - times . it is possible to include matter fields into the conformal field equations provided the equations for the matter have well - defined and compatible conformal transformation properties . this will be the case for most of the interesting fundamental field equations ( maxwell , yang - mills , scalar wave [ 77 , 78 ] etc . ) we have set up the conformal field equations as a system of equations for the conformal geometry of a conformal lorentz manifold . as such they are invariant under general diffeomorphisms and , as we have seen , under conformal rescalings of the metric . in this form they are not yet very useful for treating questions of existence of solutions or even for numerical purposes . for existence results and numerical evolution the geometric equations have to be transformed into partial differential equations for tensor components which can then be used to set up well - posed initial value problems for hyperbolic systems of evolution equations . this process , sometimes referred to as  hyperbolic reduction  consists of several steps . by imposing suitable gauge conditions one can specify a coordinate system , a linear reference frame and a conformal factor . then the equations can be written as equations for the components of the geometric quantities with respect to the chosen frame in the chosen conformal gauge and as functions of the chosen coordinates . in the next step , one needs to extract from the equations a subsystem of propagation equations which is hyperbolic so that it has a well - posed initial value problem . finally , one has to make sure that solutions of the reduced system give rise to solutions of the full system . this step may involve the verification that the gauge conditions imposed are compatible with the propagation equations , or that other equations ( constraints ) not included in the reduced system are preserved under the propagation . the first two steps , choice of gauge and extraction of the reduced system , are very much related . furthermore , the gauge conditions should be such that they can be imposed locally without loss of generality . the gauge freedom present in the conformal field equations can easily be determined . the freedom to choose the coordinates amounts to four scalar functions while the linear reference frame , which we take to be orthogonal , once the geometric equations have been transformed into equations for components , the next step is to extract the reduced system . these are equations for the components of the geometric quantities defined above as well as for gauge - dependent quantities : the components of the frame with respect to the coordinate basis , the components of the connection with respect to the given frame and the conformal factor . there are several well - known choices for coordinates ( harmonic , gau , bondi , etc . ) , as well as for frames ( fermi - walker transport , newman - penrose , etc . ) . these are usually  hard - wired  into the equations and one has no further control on the properties of the gauge . gau coordinates for instance have the tendency to become singular when the geodesic congruence which is used for their definition starts to self - intersect . similarly , bondi coordinates are attached to null - hypersurfaces which have the tendency to self - intersect thus destroying the coordinate system . in the context of existence proofs and the numerical evolution of the equations it is of considerable interest to have additional flexibility in order to prevent the coordinates or the frame from becoming singular . the goal is to  fix the gauge  in as flexible a manner as possible and to obtain reduced equations which still have useful properties . a scheme to obtain the reduced equations in symmetric hyperbolic form while still allowing for arbitrary gauges has been devised by friedrich   ( see also   for various examples ) . cartan s structure equations which express the torsion and curvature tensors in terms of tetrad and connection coefficients are two - form equations : they are skew on two indices and the information contained in the equations is not enough to fix the tetrad and the connection by specifying the torsion and the curvature . normally , this is achieved by reducing the number of variables , in this case the number of tetrad components and connection coefficients . however , one can just as well add appropriate further equations to have enough equations for all unknowns . the additional equations should be chosen so that the ensuing system has  nice  properties . consider , in flat space with coordinates ( x ) = ( t , x , x , x ) , a one- form  which we require to be closed : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ \mu } { \omega _ \nu } - { \partial _ \nu } { \omega _ \mu } = 0.$$\end{document } 
 from this equation we can extract three evolution equations , namely 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \partial _ t}{\omega _ 1 } - { \partial _ 1}{\omega _ 0 } = 0,\\ 
 { \partial _ t}{\omega _ 2 } - { \partial _ 2}{\omega _ 0 } = 0,\\ 
 { \partial _ t}{\omega _ 3 } - { \partial _ 3}{\omega _ 0 } = 0 . \end{array}$$\end{document } 
 obviously , these three equations are not sufficient to reconstruct  from appropriate initial data . one possibility to proceed from here is to specify one component of  freely and then obtain equations for the other three . however , it is easily seen that only by specifying 0 we can achieve a pure evolution system . otherwise , we get mixtures of evolution and constraint equations . so we may note that proceeding in this way leads to a restriction of possibilities as to which components should be specified freely and , in general , it also entails that derivatives of the specified component appear . another possible procedure is to enlarge the system by adding an equation for the time derivative of 0 . doing this covariantly implies that we should add an equation in the form of a divergence 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial ^\mu } { \omega _ \mu } = f,$$\end{document } 
 where f is an arbitrary function . this results in the system 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}{l } } 
   { { \partial _ t}{\omega _ 0 } - { \partial _ 1}{\omega _ 1 } - { \partial _ 2}{\omega _ 2 } - { \partial _ 3}{\omega _ 3 } = f , } \\ { \;\;\;\;\;\;\;\;\;{\partial _ t}{\omega _ 1 } - { \partial _ 1}{\omega _ 0 } = 0 , } \\ 
   { \;\;\;\;\;\;\;\;\;{\partial _ t}{\omega _ 2 } - { \partial _ 2}{\omega _ 0 } = 0 , } \\ 
   { \;\;\;\;\;\;\;\;\,{\partial _ t}{\omega _ 3 } - { \partial _ 3}{\omega _ 0 } = 0 , } 
 \end{array}$$\end{document } 
 which is symmetric hyperbolic for any choice of f. note also that f appears as a source term and only in undifferentiated form . clearly , our influence on the component 0 is now very indirect via the solution of the system , while before we could specify it directly . in a similar way note , however , that this way of fixing a gauge is not at all specific to these equations . since it depends essentially only on the form of cartan s structure equations it is applicable in all cases where these are part of the first order system . the cartan equations can be regarded as exterior equations for the one - forms a dual to a tetrad ea and the connection one - forms \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\omega _ { { a^\nu } } ^\mu   = \omega _ c^\mu { \nabla _ a}e_\nu ^c$$\end{document}. similar to the system above , the equations involve only the exterior derivative of the one - forms and so we expect that we should add equations in divergence form , namely 
 ( 23)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \nabla ^a}\omega _ a^\mu   = { f^\mu } , \\ 
 { \nabla ^a}\omega _ { { a^\nu } } ^\mu   = { f^\mu } _ \nu , 
 \end{array}$$\end{document } 
 with arbitrary gauge source functions f for fixing coordinates and fv for choosing a tetrad . note that a = ax implies f = x;. in a given gauge ( i.e. , coordinates and frame field are specified ) the gauge sources can be determined from 
 ( 24)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f^{\mu}=\square   x^{\mu},$$\end{document } 
 ( 25)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${f^\mu } _ \nu   = { \nabla ^a}(\omega _ c^\mu { \nabla _ a}e_\nu ^c).$$\end{document } 
 in fact , these equations are exactly the same equations as ( 23 ) except that they are written in a more invariant form . now it is obvious that the gauge sources contain information about the coordinates and the frame used . what needs to be shown is that any specification of the gauge sources fixes a gauge . in fact , suppose we are given functions ^f and ^fv on  then there exist ( locally ) coordinates ^x and a frame ^ea so that in that coordinate system the gauge sources are just the prescribed functions ^f and ^fv . this follows from the equations 
 ( 26)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\square \hat{x}^{\mu}=\hat{f}^{\mu}(\hat{x}^{\nu } ) , $ $ \end{document } 
 ( 27)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^a}(\hat \omega _ c^\mu { \nabla _ a}\hat e_\nu ^c ) = { \hat f^\mu } _ { \;\;\nu } ( { \hat x^\rho } ) .$$\end{document } 
 these are semi - linear wave equations which determine a unique solution from suitably given initial data close to the initial surface . note that on the right hand side of ( 26 ) there is a function of the ^x , and not a source term . the equations can be solved in steps . once the coordinates ^x have been determined from ( 26 ) , the right hand side of ( 27 ) finally , we need to discuss the gauge freedom in the choice of the conformal factor . in many discussions of asymptotic structure the conformal factor is chosen in such a way that null - infinity is divergence - free , in addition to the vanishing of its shear , which is a consequence of the asymptotic vacuum equations . that means that infinitesimal area elements remain unchanged in size as they are parallelly transported along the generators of . since they also remain unchanged in form due to the vanishing shear of  , they remain invariant and hence they can be used to define a unique metric on the space of generators of . this choice simplifies many calculations on  , still leaving the conformal factor quite arbitrary away from . yet , in numerical applications this choice of the conformal factor may be too rigid and so one needs a flexible method for fixing the conformal factor . it turns out that one can introduce a gauge source function for the conformal gauge as well . consider the change of the scalar curvature under the conformal rescaling gab  ^gab = gab , ^ =  , : it transforms according to 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda   \mapsto \hat \lambda   = { \theta ^ { - 2}}(\lambda   + \frac{{\square \theta } } { { 4\theta } } ) \:.$$\end{document } 
 reading this transformation law as an equation for  we obtain 
 ( 28)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\square \theta   = 4\theta ( { \theta ^2}\hat \lambda   - \lambda ) .$$\end{document } 
 it follows from this equation that we may regard the scalar curvature as a gauge source function for the conformal factor : for , suppose we specify the function ^ arbitrarily on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal m}$$\end{document } , then equation ( 28 ) is a non - linear wave equation for  which can be solved given suitable initial data . this determines a unique  , hence a unique ^ and ^gab such that the scalar curvature of the rescaled metric ^gab has scalar curvature ^. note that these considerations are local . however , the problem of identifying and fixing a gauge globally is very difficult but also very important because only when the gauges are globally known one can really compare two different space - times . having established that the gauge sources do in fact , locally , fix a unique gauge we can now split the system of conformal field equations into evolution equations and constraints . we introduce an arbitrary time - like unit vectorfield t which has a priori no relation to the tetrad field used for framing . we split all the tensorial quantities into the parts which are parallel and orthogonal to that vector field using the projector ha = ab - tat . we introduce the covariant derivatives of the vectorfield t by 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\chi ^a}_b = { h^c}_a{\nabla _ c}{t^a},\;\;\;\;\;{\chi ^a } = { t^c}{\nabla _ c}{t^a}.$$\end{document } 
 they account for half ( 9 + 3 ) of the four - dimensional connection coefficients . the other half is captured by defining a covariant derivative a which has the property that it annihilates both t and ta and agrees with a when acting on tensors orthogonal to t ( see equations ( 40 ) ) . note that we have not required that t be the time - like member of the frame , nor have we assumed that it be hypersurface orthogonal . in the latter case , ab is the extrinsic curvature of the family of hypersurfaces orthogonal to t and hence it is symmetric . furthermore , the derivative a agrees with the levi - civita connection of the metric hab induced on the leaves by the metric gab . we write the equations in terms of the derivative a and the  time derivative   which is defined in a way similar to a ( see equation ( 40 ) ) , because in this form it is quite easy to see the symmetric hyperbolicity of the equations . as they stand , the equations ( 90 , 91 , 92 , 93 , 94 , 95 , 96 , 97 , 98 , 99 , 100 , 101 , 102 , 103 ) form a symmetric hyperbolic system of evolution equations for the collection of 65 unknowns 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\{c^{\mu},c_{a}^{\mu},\chi^{a},\gamma_{{a}^{b}_{c}},\lambda^{a}{\;_{b}},\phi_{ab},\phi_{a},e_{ab},b_{ab},\omega,\sigma_{a},\sigma , s\}.$$\end{document } 
 this property is present irrespective of the particular gauge . for any choice of the gauge source functions f , fa , fab and  , the system is symmetric hyperbolic . the fact that the gauge sources appear only in undifferentiated form implies that one can specify them not only as functions of the space - time coordinates but also as functions of the unknown fields . in this way , one can feed information about the current status of the evolution back into the system in order to influence the future development . other ways of specifying the coordinate gauge , including the familiar choice of a lapse function and a shift vector , are not as flexible because then not only these functions themselves appear in the equations , but also their derivatives . specifying them as functions of the unknown fields alters the principal part of the system and , hence , the propagation properties of the solution . this may not only corrupt the character of the system but it may also be disastrous for the numerical applications because an uncontrolled change of the local propagation speeds implies that the stability of a numerical scheme can break down due to violation of the cfl condition ( see   for a more detailed discussion of these issues ) . however , due to the intuitive meaning of lapse and shift they are used ( almost exclusively ) in numerical codes . apart from various possibilities to specify the gauges which result in different systems with different numbers of unknowns , one can also set up the equations using spinorial methods . this was the method of choice in almost all of friedrich s work ( see e.g.   and also ) . the ensuing system of equations is analogous to those obtained here using the tetrad formalism . the main advantages of using spinors is the fact that the reduction process automatically leads to a symmetric hyperbolic system , that the variables are components of irreducible spinors which allows for the elimination of redundancies , and that variables and equations become complex and hence easier to handle . another possibility is to ignore the tetrad formalism altogether ( or , more correctly , to choose as a basis for the tangent spaces the natural coordinate frame ) . this also results in a symmetric hyperbolic system of equations ( see [ 55 , 81 ] ) , in which the gauge dependent variables are not the frame components with their corresponding connection coefficients but the components of the spatial metric together with the usual christoffel symbols and the extrinsic curvature . the fact that the reduced equations form a symmetric hyperbolic system leads , via standard theorems , to the existence of smooth solutions which evolve uniquely from suitable smooth data given on an initial surface . we have the 
 theorem 2 
 [ friedrich ] for functions f , fa , fab ,  on  
 and data given on some initial surface let u be the solution of the reduced equations . if u satisfies the conformal field equations ( 16 , 17 , 18 , 19 , 20 , 21 ) on the initial surface then , in fact , it satisfies them on the entire domain of dependence of the initial surface in the space - time defined by u. the proof of this theorem relies on the existence of a  subsidiary system  of equations for the zero - quantity z ( see equation ( 15 ) ) , whose vanishing indicates the validity of the conformal field equations . thus , one has uniqueness of the solutions so that z vanishes in the domain of dependence of the initial surface if it vanishes on the surface . it can be shown that solutions obtained from different gauge source functions are in the same conformal class , so they lead to the same physical space - time . from the physical point of view the most interesting scenario is the following one : a gravitating material system ( like e.g. several extended bodies ) evolves from a given initial state , possibly interacting with incoming gravitational radiation and emitting outgoing gravitational radiation until it reaches a final state . initial data are prescribed on the blue parts , i.e. on a hyperboloidal hypersurface and the part of  
 which is in its future . note that the two cones  
 and  
 are separated to indicate the non - trivial transition between them . initial data are prescribed on the blue parts , i.e. on a hyperboloidal hypersurface and the part of  
 which is in its future . note that the two cones  
 and  
 are separated to indicate the non - trivial transition between them . in more mathematical terms , this requires the solution of an initial value problem : we provide appropriate initial data , describing the initial configuration of the matter and the geometry , on a hyperboloidal hypersurface 0 , and appropriate boundary data , describing the incoming gravitational radiation , on the piece of  which is in the future of 0 . then we have to show that there is a unique solution of the conformal field equations coupled to the matter equations which exists for some time . if the situation is  close enough  to a newtonian situation , i.e. , the gravitational waves are weak and the matter itself is rather  tame , then one would expect that there is a solution , i.e. a space - time , which is regular on arbitrary hyperboloidal hypersurfaces intersecting . in general , however , we can not expect to have a regular point i representing time - like infinity . the reason is not so much the incorporation of matter into the conformal field equations but a more fundamental one . space - like infinity i is a singularity for the conformal structure of any space - time , which has a non - vanishing adm - mass . without the proper understanding of i there will be no way to bridge the gap between past and future null - infinity because i provides the link between the incoming and the outgoing radiation fields . the results obtained so far are concerned only with the pure radiation problem , i.e. the vacuum case . in   christodoulou and klainerman prove the global non - linear stability of minkowski space , i.e. the existence of global solutions of the einstein vacuum equations for  small enough  cauchy data which satisfy certain fall - off conditions at space - like infinity . however , they do not recover the peeling property for the weyl tensor but a weaker fall - off , which implies that in this class of solutions the conformal compactification would not be as smooth as it was expected to be . this raises the question whether their results are sharp , i.e. , whether there are solutions in this class which indeed have their fall - off behaviour . in that case , one would probably have to strengthen the fall - off conditions of the initial data at space - like infinity in order to establish the correct peeling of the weyl tensor . then an interesting question arises as to what the physical meaning of these stronger fall - off conditions is . an indication that maybe more restrictive conditions are needed is provided by the analysis of the initial data on hyperboloidal hypersurfaces ( see below ) . the first result   obtained with the conformal field equations is concerned with the asymptotic characteristic initial value problem ( see figure 7 ) in the analytic case . figure 7the geometry of the asymptotic characteristic initial value problem : characteristic data are given on the blue parts , i.e. an ingoing null surface and the part of  which is in its future . the geometry of the asymptotic characteristic initial value problem : characteristic data are given on the blue parts , i.e. an ingoing null surface and the part of  which is in its future . note that the ingoing surface may develop self - intersections and caustics . in this kind of initial value problem , one specifies data on an ingoing null hypersurface \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n}$$\end{document } and that part of  which is in the future of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n}$$\end{document}. the data which have to be prescribed are essentially the so - called null data on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n}$$\end{document } and  , i.e. those parts of the rescaled weyl tensor which are entirely intrinsic to the respective null hypersurfaces . in the case of theorem 3 
 [ knnr ] for given smooth null data on an ingoing null hypersurface 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n}$$\end{document } 
 and a smooth radiation field on the part 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal i}$$\end{document } 
 of  
 which is to the future of the intersection s of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n}$$\end{document } 
 with  
 and certain data on s , there exists a smooth solution of einstein s vacuum equations in the future of 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n } \cup { \mathcal i}$$\end{document } 
 which implies the given data on 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal n } \cup { \mathcal i}$$\end{document}. the result is in complete agreement with sachs earlier analysis of the asymptotic characteristic initial value problem based on formal expansion methods . another case these are vacuum solutions characterized by the fact that they are smoothly extensible through past time - like infinity , i.e. by the regularity of the point i. this case has been treated in [ 49 , 51 ] . a solution of this kind is uniquely characterized by its radiation field , i.e. the intrinsic components of the rescaled weyl tensor on . in the analytic case , a formal expansion of the solution at i can be derived , and growth conditions on the coefficients can be given to ensure convergence of the formal expansion near i. furthermore , there exists a surprising relation between this type of solutions and static solutions , summarized in 
 theorem 4 
 [ friedrich ] with each asymptotically flat static solution of einstein s vacuum field equations can be associated another solution of these equations which has a smooth conformal boundary  
 and for which the point i 
 is regular . this result establishes the existence of a large class of purely radiative solutions . for applications , however , the most important type of initial value problem so far , in the sense that the asymptotic behaviour can be controlled , has been the hyperboloidal initial value problem where data are prescribed on a hyperboloidal hypersurface . this is a space - like hypersurface whose induced physical metric behaves asymptotically like a surface of constant negative curvature ( see section 2.4 ) . in the conformal picture , a hyperboloidal hypersurface is characterized simply by the geometric fact that it intersects  transversely in a two - dimensional space - like surface . prototypes of such hypersurfaces are the space - like hyperboloids in minkowski space - time . in the minkowski picture they can be seen to become asymptotic to null cones which suggests that they reach null - infinity . however , the picture is deceiving : the conformal structure is such that the hyperboloids always remain space - like , the null - cones and the hyperboloids never become tangent . the intersection is a two - dimensional surface s , a  cut  of . the data implied by the conformal fields on such a hypersurface are called hyperboloidal initial data . the first result obtained for the hyperboloidal initial value problem states that if the space - time admits a hypersurface which extends smoothly across  with certain smooth data given on it , then the smoothness of  will be guaranteed at least for some time into the future . this is contained in 
 theorem 5 
 [ friedrich ] smooth hyperboloidal initial data on a hyperboloidal hypersurface  determine a unique solution of einstein s vacuum field equations which admits a smooth conformal boundary at null - infinity in the future of . there exists also a stability result which states that there are solutions which behave exactly like minkowski space near future time - like infinity : 
 theorem 6 
 [ friedrich ] if the hyperboloidal initial data are in a sense sufficiently close to minkowskian hyperboloidal data , then there exists a conformal extension of the corresponding solution which contains a point i 
 such that  
 is the past null cone of that point . it should be emphasized that this result implies that the physical metric of the corresponding solution is regular for all future times . thus , the theorem constitutes a ( semi-)global existence result for the einstein vacuum equations . that such data exist follows already from theorem 4 because one can construct hyper - boloidal hypersurfaces together with data on them in any of the radiative solutions whose existence is guaranteed by that theorem . however , one can also construct such data sets in a similar way to the construction of cauchy data on an asymptotically euclidean hypersurface by solving the constraint equations implied on the cauchy surface . let  be a hyperboloidal hypersurface in an asymptotically flat vacuum space - time which extends out to  , touching it in a two - surface  which is topologically a two - sphere . the assumptions on  are equivalent to the fact that on     s is a smooth riemannian manifold with boundary , which carries a smooth metric hab and a smooth function  , obtained by restriction of the unphysical metric and the conformal factor . the conformal factor is a defining function for the boundary  ( i.e. it vanishes only on  with non - vanishing gradient ) , and together with the metric it satisfies 
 ( 29)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde h_{ab } } = { \omega ^ { - 2}}{h_{ab}}$$\end{document } 
 on  , where hab is the metric induced on  by the physical metric . furthermore , let ab be the extrinsic curvature of  in the physical space - time . together with hab it satisfies the vacuum constraint equations , 
 ( 30)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde r - { \tilde \chi _ { ab}}{\tilde \chi ^{ab } } + { ( { \tilde \chi ^a}_{\;\;a})^2 } = 0,$$\end{document } 
 ( 31)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde \partial _ a}{\tilde \chi ^a}_{\;\;b } - { \tilde \partial _ b}\tilde \chi _ { \;\;a}^a = 0,$$\end{document } 
 where a is the levi - civita connection of hab and r is its scalar curvature . we discuss only the simplest case here , referring to the literature for the more general results . assume that the extrinsic curvature is a pure trace term , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde \chi _ { ab } } = \frac{1}{3}c{\tilde h_{ab}}.$$\end{document } 
 the momentum constraint ( 31 ) implies that c is constant while the hyperboloidal character of  implies that c  0 . with these simplifications and a rescaling of hab with a constant factor , the hamiltonian constraint ( 30 ) takes the form 
 ( 32)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde r =   - 6.$$\end{document } 
 a further consequence of the condition ( 29 ) is the vanishing of the magnetic part bab of the weyl tensor . for any defining function  of the boundary , the conformal factor has the form  = . expressing equation ( 32 ) in terms of the unphysical quantities hab and  yields the single second - order equation 
 ( 33)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$8{\omega ^2}\delta \phi   - 8\omega { \partial ^a}\omega { \partial _ a}\phi   - [ { \omega ^2}r + 4\omega \delta \omega   - 6{\partial ^a}\omega { \partial _ a}\omega ] \phi   = 6{\phi ^5}.$$\end{document } 
 this equation is a special case of the lichnrowicz equation and is sometimes also referred to as the yamabe equation . for a given metric hab and boundary defining function it is a second - order , non - linear equation for the function . note that the principal part of the equation degenerates on the boundary . therefore , on the boundary , the yamabe equation degenerates to the relation 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\phi { \partial ^a}\omega { \partial _ a}\omega   = { \phi ^5}.$$\end{document } 
 note also that  = 0 is a solution of ( 33 ) which , however , is not useful for our purposes because it would correspond to a conformal factor with vanishing first derivative on . therefore , we require that  be non - vanishing on the boundary , i.e. bounded from below by a strictly positive constant . then the relation above determines the boundary values of  in terms of the function  . taking derivatives of equation ( 33 ) , one finds that also the normal derivative of  is fixed on the boundary in terms of the second derivative of . a given metric hab does not fix a unique pair ( hab ,  ) . therefore , equation ( 33 ) has the property that , for fixed  , rescaling the metric hab with an arbitrary smooth non - vanishing function  on   s according to hab  hab results in a rescaling of the solution  of ( 33 ) according to    and , hence , a change in the conformal factor   . now we define the trace - free part sab = ab - 1/3habc of the projection ab of the trace - free part of the unphysical ricci tensor onto  and consider the equations 
 ( 1)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { s_{ab } } =   - { \omega ^ { - 1}}({\partial _ b}\omega   - \frac{1}{3}{h_{ab}}\delta \omega ) \:,\\ 
 { e_{ab } } =   - { \omega ^ { - 1}}({r_{ab } } - \frac{1}{3}{h_{ab}}r - { s_{ab}})\ : , 
 \end{array}$$\end{document } 
 which follow from the equations ( 20 ) and ( 16 ) , respectively . together with the fields hab ,  , s = 1/3 they provide initial data for all the quantities appearing in the evolution equations under the given assumptions . as they stand , these expressions are formally singular at the boundary and one needs to worry about the possibility of a smooth extension of the field to . this question was answered in , where the following theorem was proved : 
 theorem 7 
 suppose (  , h ) is a three - dimensional , orientable , compact , smooth riemannian manifold with boundary . then there exists a unique solution  of ( 33 ) , and the following conditions are equivalent : 
 the function  as well as the tensor fields sab 
 and eab 
 determined on the interior  from h and  =  extend smoothly to all of .the conformal weyl tensor computed from the data vanishes on .the conformal class of h is such that the extrinsic curvature of  with respect to its embedding in (  , h ) is pure trace . the function  as well as the tensor fields sab 
 and eab 
 determined on the interior  from h and  =  extend smoothly to all of . the conformal weyl tensor computed from the data vanishes on . the conformal class of h is such that the extrinsic curvature of  with respect to its embedding in (  , h ) is pure trace . condition ( 3 ) is a weak restriction of the conformal class of the metric h on  , since it is only effective on the boundary . it is equivalent to the fact that in the space - time which evolves from the hyperboloidal data , null - infinity  is shear - free . interestingly , the theorem only requires  to be orientable and does not restrict the topology of  any further . this theorem gives the answer in a highly simplified case because the freedom in the extrinsic curvature has been suppressed . but there are also several other , less restrictive , treatments in the literature . in [ 2 , 3 ] the assumption ( 29 ) is dropped allowing for an extrinsic curvature which is almost general apart from the fact that the mean curvature is required to be constant . in also this requirement is dropped ( but , in contrast to the other works , there is no discussion of smoothness of the implied conformal initial data ) , and in   the existence of hyperboloidal initial data is discussed for situations with a non - vanishing cosmological constant . the theorem states that one can construct the essential initial data for the evolution once equation ( 33 ) has been solved . the data are given by expressions which are formally singular at the boundary because of the division by the conformal factor . this is of no consequences for the analytical considerations if condition ( 3 ) in the theorem is satisfied . however , even then it is a problem for the numerical treatments because one has to perform a limit process to get to the values of the fields on the boundary . it is clear from equations ( 79 , 80 , 81 , 82 , 83 , 84 , 85 , 86 , 87 , 88 , 89 ) that the conformal constraints are regular as well . some of the equations are rather simple but the overall dependencies and interrelations between the equations are very complicated . at the moment there exists no clear analytical method ( or even strategy ) for solving this system . an interesting feature appears in connection with condition ( 3 ) of the theorem and analogous conditions in the more general cases . the necessity of having to impose this condition seems to indicate that the development of hyperboloidal data is not smooth but in general at most \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\mathcal c}^2}$$\end{document}. if the condition were not imposed then logarithms appear in an expansion of the solution of the yamabe equation near the boundary , and it is rather likely that these logarithmic terms will be carried along with the time evolution , so that the developing null - infinity looses differentiability . thus , the conformal boundary is not smooth enough and , consequently , the weyl tensor need not vanish on  which , in addition , is not necessarily shear - free . one can show   that generically hyperboloidal data fall into the class of  poly - homogeneous  functions which are ( roughly ) characterized by the fact that they allow for asymptotic expansions including logarithmic terms . this behaviour is in accordance with other work   on the smoothness on  , in particular with the bondi - sachs type expansions which were restricted by the condition of analyticity ( i.e. no appearance of logarithmic terms ) . solutions of the hyperboloidal initial value problem provide pieces of space - times which are semi - global in the sense that their future ( or past ) development is determined . however , the domain of dependence of a hyperboloidal initial surface does not include space - like infinity and one may wonder whether this fact is the reason for the apparent generic non - smoothness of null - infinity . is it not conceivable that the possibility of making a connection between  and  across i to build up a global space - time automatically excludes the non - smooth data ? if we let the hyperboloidal initial surface approach space - like infinity it might well be that condition ( 3 ) imposes additional conditions on asymptotically flat cauchy data at spatial infinity . these conditions would make sure that the development of such cauchy data is an asymptotically flat space - time , in particular that it has a smooth conformal extension at null - infinity . these questions give some indications about the importance of gaining a detailed understanding of the structure of gravitational fields near space - like infinity . one of the difficulties in obtaining more information about the structure at space - like infinity is the lack of examples which are general enough . they possess a part of a smooth null - infinity which , however , is incomplete . this is a general problem because the existence of a complete null - infinity with non - vanishing radiation restricts the possible isometry group of a space - time to be at most one - dimensional with space - like orbits . some of the boost - rotation symmetric space - times even have a regular i , thus they have a vanishing adm - mass . other examples exist of space - times which are solutions of the einstein - maxwell   or einstein - yang - mills   equations . however , they were constructed in a way which enforces the field to coincide with the schwarzschild or the reissner - nordstrom solutions near i. so they are not general enough to draw any conclusions about the generic behaviour of asymptotically flat space - times near i. as indicated above , the problem at i is one of the most urgent and important ones related to the conformal properties of isolated systems . however , it is also one of the most complicated ones and a thorough discussion of all its aspects is not possible here . we will try to explain in rough terms what the new developments are , but we have to refer to   for the rigorous statements and all the details . there have been several approaches over the years towards a treatment of space - like infinity . geroch   gave a geometric characterization along the same lines as for null - infinity based on the conformal structure of cauchy surfaces . he used his construction to define multipole moments for static space - times [ 63 , 64 ] , later to be generalized to stationary space - times by hansen . it was shown by beig and simon [ 19 , 136 ] that the multipole moments uniquely determine a stationary space - time and vice versa . different geometric characterizations of spatial infinity in terms of the four - dimensional geometry were given by sommers , ashtekar and hansen [ 10 , 6 ] , and by ashtekar and romano . the difficulties in all approaches which try to characterize the structure of gravitational fields at space - like infinity in terms of the four - dimensional geometry arise from the lack of general results about the evolution of data near spatial infinity . since there are no radiating solutions which are general enough at spatial infinity to provide hints , one is limited more or less to one s intuition . reasonable  asymptotic conditions on the gravitational field at i and from them derive certain nice properties of space - times which satisfy these conditions . but there is no guarantee that there are indeed solutions of the einstein equations which exhibit the claimed asymptotic behaviour . in a sense , all these characterisations are implicit definitions of certain classes of space - times ( namely those which satisfy the imposed asymptotic conditions ) . what is needed is an analysis at space - like infinity which is not only guided by the geometry but which also takes the field equations into account ( see e.g. [ 18 , 17 ] for such attempts using formal power series ) . recently , friedrich   has given such an analysis of space - like infinity which is based exclusively on the initial data , the field equations and the conformal structure of the space - time . first , in order to simplify the analysis , an assumption on the initial data ( metric and extrinsic curvature ) on an asymptotically euclidean hypersurface  is made . since the focus is on the behaviour of the fields near space - like infinity , the topology of  is taken to be . it is assumed that the data are time - symmetric ( ab = 0 ) and that on  a ( negative definite ) metric hab with vanishing scalar curvature is given . let  be the conformal completion of (  , hab ) which is topologically s , obtained by attaching a point i to  , and assume furthermore that there exists a smooth positive function  on  with (i ) = 0 , d(i ) = 0 and hess (i ) negative definite . furthermore , the metric hab = hab extends to a smooth metric on s. thus , the three - dimensional conformal structure defined by hab is required to be smoothly extensible to the point i. from these assumptions follows that the conformal factor near i is determined by two smooth functions u and w , where u is characterized by the geometry near i while w collects global information because \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$w(i ) = { \textstyle{1 \over 2}}{m_{{\rm{adm}}}}$$\end{document } , while u(i ) = 1 . with this information the rescaled weyl tensor , the most important piece of initial data for the conformal field equations , near i is found to consist of two parts , a  massive  and a  mass - less  part . under suitable conditions , the mass - less part , determined entirely by the local geometry near i , can be extended in a regular way to i , while the massive part always diverges at i as 1/|x| in a normal coordinate system ( x , x , x ) at i unless the adm - mass vanishes . in order to analyze the singular behaviour of the initial data in more detail , the point i is blown up to a spherical set i essentially by replacing it with the sphere of unit vectors at i. roughly speaking , this process yields a covering space \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal c}$$\end{document } of ( a suitable neighbourhood of i in )  projecting down to  which has the following properties : the pre - image i of i is an entire sphere while any other point on  has exactly two pre - image points . there exists a coordinate r on \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\mathcal c}$$\end{document } which vanishes on i and which is such that on each pair of pre - image points it takes values r and -r , respectively . the actual blowup procedure involves a rather involved bundle construction which also takes into account the tensorial ( respectively spinorial ) nature of the quantities in question . the next important step is the realization that , in order to take full advantage of the conformal structure of space - time , it is not enough to simply allow for metrics which are conformally equivalent to the physical metric but that one should also allow for more general connections . instead of using a connection which is compatible with a metric in the conformal class , one may use a connection a which is compatible with the conformal structure , i.e. which satisfies the condition 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ c}{\tilde g_{ab } } = 2{\lambda _ c}{\tilde g_{ab}}$$\end{document } 
 for some one - form a . if a is exact , then one can find a metric in the conformal class for which a is the levi - civita connection . this generalization is motivated by the use of conformal geodesics as indicated below , and its effect is to free up the conformal factor , which we call  to distinguish it from the conformal factor  given on the initial surface  , from the connection ( recall that two connections which are compatible with metrics in the same conformal class differ only by terms which are linear in the first derivative of the conformal factor relating the metrics ) . as a consequence , the conformal field equations , when expressed in terms of a generalized connection , do not any longer contain an equation for the conformal factor . it appears , instead , as a gauge source function for the choice of conformal metric . additionally , a free one - form appears which characterizes the freedom in the choice of the conformal connection . to fix this freedom , they are given in terms of a system of ordinary differential equations ( ode s ) for their tangent vector together with a one - form along them . in coordinates their crucial property is that they are defined entirely by the four - dimensional conformal structure with no relation to any specific metric in that conformal structure . a time - like congruence of such curves is used to set up a  gau coordinate system in a neighbourhood of i and to define a conformal frame , a set of four vectorfields which are orthonormal for some metric in the conformal class . this metric in turn defines a conformal factor  which rescales it to the physical metric . the one - form determined by the conformal geodesics defines a conformal connection a , thus fixing the freedom in the connection . in this way if the physical space - time is a vacuum solution of the einstein equations then one can say more about the behaviour of the conformal factor  along the conformal geodesics : it is a quadratic function of the natural parameter  along the curves , vanishing at exactly two points if the initial conditions for the curves are chosen appropriately . the vanishing of  indicates the intersection of the curves with . the intersection points are separated by a finite distance in the parameter . now one fixes an initial surface  with data as described above , and the conformal geodesics are used to set up the coordinate system and the gauge as above . when the blow - up procedure is performed for  , a new finite representation of space - like infinity figure 8the geometry near space - like infinity : the  point  i 
 has been blown up to a cylinder which is attached to  
 and . the physical space - time is the exterior part of this  stovepipe  . the  spheres   
 and  
 are shown in blue and light green , respectively . the geometry near space - like infinity : the  point  i 
 has been blown up to a cylinder which is attached to  
 and . the physical space - time is the exterior part of this  stovepipe  . the  spheres   
 and  
 are shown in blue and light green , respectively . the point i on the initial surface has been replaced by a sphere i which is carried along the conformal geodesics to form a finite cylinder . the surfaces are the surfaces on which the conformal factor  vanishes . they touch the cylinder in the two spheres i , respectively . the conformal factor  vanishes with non - vanishing gradient on i and on  while on the spheres there is for the first time a clean separation of the issues which go on at space - like infinity : the spheres i are the places where   touches i  while the finite cylinder i serves two purposes . on the one hand , it represents the endpoints of space - like geodesics approaching from different directions , while , on the other hand , it serves as the link between past and future null - infinity . the part  outside  the cylinder where r is positive between the two null surfaces  corresponds to the physical space - time , while the part with r negative is not causally related to the physical space - time but constitutes a smooth extension . for easy reference . the conformal field equations , when expressed in the conformal gau gauge of this generalized conformal framework , yield a system of equations which has similar properties as the earlier version : it is a system of equations for a frame , the connection coefficients with respect to the frame , and the curvature , split up into the ricci and the weyl parts ; they allow the extraction of a reduced system which is symmetric hyperbolic and propagates the constraint equations . the bianchi identities , which form the only sub - system consisting of partial differential equations , again play a key role in the system . due to the use of the conformal gau gauge , all other equations are simply transport equations along the conformal geodesics . the reduced system is written in symbolic form as 
 ( 35)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\rm{a}}^\tau } \frac{\partial } { { \partial \tau } } { \rm{u } } + { { \rm{a}}^r}\frac{\partial } { { \partial r}}{\rm{u + } } { { \rm{a}}^\theta } \frac{\partial } { { \partial \theta } } { \rm{u + } } { { \rm{a}}^\phi } \frac{\partial } { { \partial \phi } } { \rm{u = bu,}}$$\end{document } 
 with symmetric matrices a , a , a , a , and c which depend on the unknown u and the coordinates (  , r ,  ,  ) . this system as well as initial data for it , first defined only on the original space - time , can be extended in a regular way to an extended neighbourhood of space - like infinity which allows for the setup of a regular initial value problem at space - like infinity . its properties are most interesting : when restricted to i , the initial data coincide necessarily with minkowski data , which together with the vanishing of  implies that on the entire cylinder i the coefficient matrix a vanishes . thus , the system ( 35 ) degenerates into an interior symmetric hyperbolic system on i. therefore , the finite cylinder i is a total characteristic of the system . the two null - infinities  are also characteristics , and at the intersections i between them and i the system degenerates : the coefficient matrix a which is positive definite on i and  looses rank on i. the fact that i is a total characteristic implies that one can determine all fields on i from data given on i. i is not a boundary on which one could specify in- or outgoing fields . this is no surprise , because the system ( 35 ) yields an entirely structural transport process which picks up data delivered from  via i and moves them to  via i. it is also consistent with the standard cauchy problem where it is known that one can not specify any data  at infinity  . the degeneracy of the equations at i means that one has to take special precautions to make sure that the transitions from and to  are smooth . in fact , not all data  fit through the pipe  : friedrich has derived restrictions on the initial data of solutions of the finite initial value problem which are necessary for regularity through i. they are conditions on the conformal class of the initial data , stating that the cotton tensor and all its symmetrized and trace - removed derivatives should vanish at the point i in the initial surface . if this is not the case , then the solution of the intrinsic system will develop logarithmic singularities which will then probably spread across null - infinity , destroying its smoothness . so here is another concrete indication that initial data have to be restricted albeit in a rather mild way in order for the smooth picture of asymptotic flatness to remain intact . it is not known whether this condition is also sufficient nor what its physical implications are . note that the conditions on the cotton tensor are entirely local - at - infinity . the setting described in the above paragraphs certainly provides the means to analyze the consequences of the conformal einstein evolution near spacelike infinity and to understand the properties of gravitational fields in that region . the finite picture allows the discussion of the relation between various concepts which are defined independently at null and space - like infinity . as one application of this kind friedrich and knnr   have related the newman - penrose constants which are defined by a surface integral over a cut of  to initial data on . the cut of  is pushed down towards i where it is picked up by the transport equations of system ( 35 ) . in a similar way , one can relate the bondi- and adm - masses of a space - time . readers who are interested in obtaining information about the current status of existence theorems in general relativity are referred to rendall s recently updated living reviews article . although we have restricted ourselves to the case of a vanishing cosmological constant this does not mean that other cases have not been analyzed . friedrich   has shown the existence of asymptotically simple de sitter type solutions of the einstein - maxwell - yang - mills equations and stability of de sitter space . more recently , with the general version of his conformal field equations , he was able to demonstrate the existence of asymptotically simple anti - de sitter type solutions . finally , it should also be mentioned that recently some work has gone into generalizing the bondi - sachs calculations to poly - homogeneous space - times [ 30 , 143 , 144 ] . in this section we discuss some of the issues which have to be confronted in a numerical implementation of the einstein equations and , in particular , of the conformal field equations to solve the hyperboloidal initial value problem . the material presented in this section is based to a large extent on work by hbner [ 76 , 77 , 78 , 79 , 81 , 82 ] and frauendiener [ 38 , 39 , 40 , 41 , 42 ] . these are also the references which should be consulted for a more detailed discussion of the various numerical schemes . the presentation here will be rather cursory in order to keep it reasonably short . at the moment they are based on three different formulations of the einstein equations : the standard cauchy problem using the einstein equations in either adm - like form or in one of the current hyperbolic formulations , the characteristic initial value problem for the standard einstein equations in bondi- or newman - penrose form , and the cauchy problem for the conformal field equations . let us discuss the main properties of each of these approaches with respect to their treatment of infinity . in this connection it is useful to look at conformal diagrams . in all three cases a source , indicated by a shaded brown region , evolves between the time - like past and future infinities i and thereby emits gravitational radiation which escapes towards future null infinity . the three methods differ in how they introduce an evolution process in order to simulate the physical system . the goal is to obtain the radiation data which are received by an observer on  as accurately as possible ( see   for a discussion of the issues involved in the reception process ) . in figure 9 the geometry is sketched which is implied by solving the cauchy problem for the standard einstein equations in any of the various forms which exist today : adm formulation , old   and improved , or any one of the hyperbolic formulations like , e.g. in [ 1 , 61 ] . in this approach the physical space - time is decomposed into space and time by the introduction of a time coordinate . the hypersurfaces of constant physical time are cauchy surfaces and , hence , asymptotically euclidean . the blue line corresponds to a ( normally time - like ) hypersurface which has to be introduced because for the standard einstein equations , formulated in the physical space - time , null - infinity is infinitely far away . the numerical grids are distributed along the constant time hypersurfaces , and with each time - step they evolve from one hypersurface to the next . they have outermost grid - points which , in the course of the evolution , move along that time - like hypersurface . the green lines are surfaces of constant time . the blue line indicates the outer boundary . the green lines are surfaces of constant time . the blue line indicates the outer boundary . thus , inevitably , we are facing the task of solving an initial boundary value problem ( ibvp ) for the einstein equations . this , in fact , is one of the big problems in scientific computing as witnessed by the following quotation : 
 the difficulties caused by boundary conditions in scientific computing would be hard to overemphasize . boundary conditions can easily make the difference between a successful and an unsuccessful computation , or between a fast and a slow one . yet in many important cases , there is little agreement about what the proper treatment of the boundary should be . one of the sources of difficulty is that although some numerical boundary conditions come from discretizing the boundary conditions for the continuous problem , other  artificial  or the reason is that the number of boundary conditions required by a finite difference formula depends on its stencil , not on the equation being modeled . thus even a complete mathematical understanding of the initial boundary value problem to be solved ( which is often lacking ) is in general not enough to ensure a successful choice of numerical boundary conditions . open  ,  radiation  ,  absorbing  ,  non - reflecting  , or boundary conditions , which are numerical artifacts designed to limit a computational domain in a place where the mathematical problem has no boundary at all . 
 despite these remarks , perhaps the most basic and useful advice one can offer concerning numerical boundary conditions is this : before worrying about discretization , make sure to understand the mathematical problem . if the ibvp is ill - posed , no amount of numerical cleverness will lead to a successful computation . this principle may seem too obvious to deserve mention , but it is surprising how often it is forgotten . 
 while these remarks apply to any numerical problem , they are particularly relevant to numerical relativity . the difficulties caused by boundary conditions in scientific computing would be hard to overemphasize . boundary conditions can easily make the difference between a successful and an unsuccessful computation , or between a fast and a slow one . yet in many important cases , there is little agreement about what the proper treatment of the boundary should be . one of the sources of difficulty is that although some numerical boundary conditions come from discretizing the boundary conditions for the continuous problem , other  artificial  or the reason is that the number of boundary conditions required by a finite difference formula depends on its stencil , not on the equation being modeled . thus even a complete mathematical understanding of the initial boundary value problem to be solved ( which is often lacking ) is in general not enough to ensure a successful choice of numerical boundary conditions . open  ,  radiation  ,  absorbing  ,  non - reflecting  , or boundary conditions , which are numerical artifacts designed to limit a computational domain in a place where the mathematical problem has no boundary at all . despite these remarks , perhaps the most basic and useful advice one can offer concerning numerical boundary conditions is this : before worrying about discretization , make sure to understand the mathematical problem . if the ibvp is ill - posed , no amount of numerical cleverness will lead to a successful computation . this principle may seem too obvious to deserve mention , but it is surprising how often it is forgotten . until only quite recently no result on this situation has changed with the work in   where it is proved that the ibvp for the einstein equations in the physical space - time is well - posed for a certain class of boundary conditions . in rough terms , these boundary conditions can be understood as follows : for an arbitrarily specified time - like unit vector tangent to the boundary we can define a null - tetrad by taking the two real null - vectors as lying in the two - dimensional space spanned by that vector and the unit vector normal to the boundary . this fixes the null - tetrad up to boosts in that plane and rotations in the plane orthogonal to it . with respect to that tetrad we can define the ( complex ) weyl tensor components 0 and 4 . then , the boundary conditions have the form 
 ( 36)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$q =   - { \psi _ 4 } + \alpha { \psi _ 0 } + \beta { \bar \psi _ 0}.$$\end{document } 
 here  ,  and q are complex functions on the boundary . q can be freely specified , thus representing the free data , while  and  are restricted by some condition , which implies in particular that || + ||  1 . they can be regarded as specifying the reflection and transmission properties of the boundary . in the frame 4 can be interpreted as the part of the weyl tensor which propagates orthogonally across the boundary in the outward or inward direction , respectively . the boundary conditions ( 36 ) essentially say that we can specify the value of the ingoing part freely and couple some parts of the outgoing field back to it . consider the boundary condition with q = 0 and  =  = 0 viz . this looks like a completely absorbing boundary condition or , in other words , like the condition of no incoming radiation . recall that the weyl tensor components are defined with respect to a time - like unit vector . but there is no distinguished time - like unit - vector available so that they are all equivalent . this means that the boundary conditions , and hence the presence or absence of incoming radiation , are observer dependent . clearly , the condition of no incoming radiation is not compatible with the einstein equations anyway , because outgoing waves will always be accompanied by backscattered waves generated by the non - linear nature of the equations . the impossibility of a local completely absorbing boundary condition for the scalar wave equation in flat space has been emphasized already in numerical analysis . the deep reason behind this fact seems to be the lorentz invariance of the equations . while the conditions in ( 36 ) give rise to a well - posed ibvp , they have not found their way into the numerical work yet . returning to figure 9 we see that there are several fundamental difficulties with this ( currently standard ) approach . recall that the goal is to obtain radiative information near . to achieve this goal we need to arrange for two things : we need to make sure that we simulate an asymptotically flat space - time and we need to have a stable code in order to get meaningful results . in this approach , both these requirements can only be realized by an appropriate boundary condition . the parts of space - time above that line are influenced by the boundary condition . it is obvious from the causal dependencies that the radiative information which can be accessed using this approach is to a large extent affected by the boundary condition . a bad choice of the boundary condition , in particular one which is not compatible with the einstein evolution , or one which is not adequate for asymptotic flatness , can ruin any information gained about the radiation . . it must be 
 physically reasonable : in some way it should express the asymptotic flatness of the space - time manifold , mathematically reasonable : the ibvp should be well - posed with this boundary condition , numerically reasonable : codes which implement the boundary condition should be stable . 
 devising a boundary condition which has these properties there have been several proposals reaching from specifying minkowski or schwarzschild data to matching with linearized solutions of the einstein equations . although these conditions are physically reasonable to some extent , they are not compatible with the einstein evolution and therefore they do not possess the other two properties . such inappropriate boundary conditions can be expected to perform well for small amplitudes and for short times when the non - linear effects can be neglected . physically reasonable : in some way it should express the asymptotic flatness of the space - time manifold , mathematically reasonable : the ibvp should be well - posed with this boundary condition , numerically reasonable : codes which implement the boundary condition should be stable . we can also see from the diagram that the choice of asymptotically euclidean hypersurfaces is unfortunate for the following reason . a natural thing to do in order to increase the accuracy of the results is to push the boundary further outward towards infinity . however , pushing along surfaces of constant time brings us closer to i and not to the interesting parts in the neighbourhood of . this implies that the numerical codes not only increase in size but in addition that they have to run longer ( in the sense of elapsed physical time ) until the interesting regions are covered . one idea   to overcome the problems encountered in the standard cauchy approach is the so - called cauchy - characteristic matching procedure . in this approach one introduces a time - like hypersurface , the interface , along which one tries to match a cauchy problem with a characteristic initial value problem . the cauchy part provides the interior boundary condition for the characteristic part while the latter provides the outer boundary condition for the former . thus , the space - time is foliated by hypersurfaces which are partly space - like and partly null , where the causal character changes in a continuous but non - differentiable way across a two - dimensional surface ( see figure 10 ) . the blue line indicates the interface between the ( inner ) cauchy part and the ( outer ) characteristic part . the blue line indicates the interface between the ( inner ) cauchy part and the ( outer ) characteristic part . the characteristic part of the hypersurfaces reaches out all the way to . the advantage of this procedure is that the problem of finding the correct boundary condition has been eliminated by the introduction of the characteristic part of the scheme . the characteristic part has been implemented numerically in a very successful way ( see the living reviews article   for a recent review on that topic ) . it is based on the bondi - sachs equations for the gravitational field on null - hypersurfaces . using a  compactified coordinate along the generators of the null - hypersurfaces it is even possible to  bring null - infinity to finite places  . however , this is different from the conformal compactification because the metric is not altered . this necessarily leads to equations which  notice where infinity is  because they degenerate there . because of this degeneracy the solutions can be singular unless very specific boundary values are prescribed . the critical place is , of course , the interface where two completely different codes have to be matched . this is not merely a change in the numerical procedures but the entire setup changes across the interface : the geometry , the independent and dependent variables etc . test calculations have been performed successfully for space - times with symmetries ( cylindrical , spherical ) and/or model equations like non - linear wave equations . currently , a combined code for doing the pure gravitational problem in three dimensions is being developed . finally , we want to discuss the approach based on the conformal compactification using the conformal field equations ( see figure 11 ) . in this case , the arena is not the physical space - time but some other unphysical manifold which is conformally related via some conformal factor . the geometry of the conformal approach . the physical space - time is the part of the conformal manifold on which  is positive . introducing an appropriate time - coordinate in the conformal manifold leads to a foliation by space - like hypersurfaces which also cover the physical manifold . those hypersurfaces which intersect  transversely are hyperboloidal hypersurfaces in the physical space - time . it is important to note that they are submanifolds of the conformal manifold so that they do not stop at  but continue smoothly across  which is just another submanifold of the conformal manifold , albeit a special one . hyperboloidal initial data ( see 3.4 ) are given on one such hypersurface and are subsequently evolved with the symmetric hyperbolic ( reduced ) system of evolution equations . in contrast to the behaviour discussed above in connection with the characteristic evolution , the conformal field equations do not know where  is . they know that there might be a  ( because the conformal factor is one of the variables in the system ) but they have to be told its location on the initial data surface , and from there on this information is carried along by the evolution . first of all , one might have thought that the question of the correct boundary condition for the ibvp for the einstein equations has now merely been shifted to the even more complicated problem of finding a boundary condition for the conformal field equations . this is true but irrelevant for the following reason : suppose we give hyperboloidal initial data on that part of the initial surface which lies inside the physical region . the structure of the characteristics of the conformal field equations is such that ( provided the gauges are chosen appropriately ) the outermost characteristic is the light cone . hence  is a characteristic for the evolution equations , and this implies that the extension of the initial data into the unphysical region can not influence the interior . similarly , it is true that also the conformal field equations need a boundary condition . but since this boundary has been shifted into the unphysical region there is no need for it to be physically reasonable . it is enough to have a boundary condition which is numerically reasonable , i.e. which leads to stable codes . the information generated at the boundary travels at most with the velocity of light and so it can not swap into the physical region . numerically , things are somewhat different but one can expect that any influence from the outside will die away with the order of accuracy of the discretization scheme used . a further advantage of the conformal approach is the possibility to study global properties of the space - times under consideration . not only do the hyperboloidal hypersurfaces extend up to and beyond null - infinity but it is also possible to study the long - time behaviour of gravitational fields . if the initial data are small enough so that future time - like infinity i is a regular point ( see theorem 6 ) then one can determine in a finite ( conformal and computational ) time the entire future of the initial surface and therefore a ( semi-)global space - time up to ( and even beyond ) i. this has been successfully demonstrated by hbner   in the case of spherically symmetric calculations and more recently also in higher dimensions . the calculation of a global space - time up to i has the effect of shrinking the region on the computational domain which corresponds to the physical space - time because  seems to move through the grid . this means that more and more grid - points are lost for the computation of the physical part . then the available gauge freedom can be used to gain complete control over the movement of  through the grid . the conformal field equations lead to a first order system of approximately sixty equations depending on the specific formulation and gauge conditions . in the version derived in appendix 6 however , on the one hand one has to compare this with current proposals for hyperbolic systems for the einstein equations , which advocate a number of equations of the same order ( although there are not even equations for a conformal factor ) . on the other hand , the variables which appear in the conformal field equations are of a very geometric nature , and for further investigations of the space - times ( finding the geodesics , computing forces etc . ) , these variables would have to be computed anyway . but we can be somewhat more specific : the conformal field equations have a total number of equations which is roughly 60 . the complexity ( in particular , the memory ) scales roughly with the number of gridpoints and for a 3d code that is n , where n is the number of grid points per dimension . the upshot of this is that the memory requirements are not dictated by the number of equations because this is a fixed factor but by the dimensionality of the code , and this affects both codes in the same way . a final word concerning the inclusion of matter into the conformal approach : it is clear that this is more complicated than it is for the standard einstein equations . the reason is that the conformal field equations also contain the ( rescaled ) weyl curvature which couples to the energy - momentum tensor via its derivatives . this means that one needs equations not only for those matter fields which appear in the energy - momentum tensor but also for their derivatives . furthermore , the fact that under conformal rescalings the trace - free part of the energy - momentum tensor behaves differently compared to its trace , causes additional difficulties . it was mentioned in section 3.1 that for  = 1 the conformal field equations reduce to the standard einstein equations together with the vacuum bianchi identity . this suggests that one should specify the conformal factor initially to be equal to unity in a region which contains the sources ( assuming the sources have spatially compact support ) and then decreases until it vanishes on . furthermore , the time - derivative of  should be adjusted such that it vanishes in a neighbourhood of the source . in fact , it is possible to do so , but it is not yet known , whether and , if so , for how long one can keep up this condition of constant conformal factor around the source . the proposed procedure does not eliminate the problem that additional equations for the matter variables are needed , but it might reduce the complexity somewhat . we will now discuss some issues related to the construction of hyperboloidal initial data and the implementation of the evolution equations . in order to start the evolution one has to determine the initial data . so far , the numerical methods used for their construction are based entirely on theorem 7 . it is assumed that the extrinsic curvature of the hyperboloidal initial surface is pure trace . then the yamabe equation is solved for a conformal factor which chooses from a given conformal class of riemannian metrics the one which has constant negative curvature . this is a difficult problem if one tries to prove existence of solutions of this equation because the loss of ellipticity on the boundary means that one can not simply appeal to known theorems . the yamabe equation is solved in a more or less standard way by using a richardson iteration scheme to reduce the solution of the non - linear equation to a series of inversions of a linear operator . the degeneracy of the equation on the boundary forces the solution and its derivative by regularity to have certain well defined boundary values . this and the global nature of the elliptic equations suggests to use pseudo - spectral ( or collocation ) methods . such methods are well suited for problems for which it is known beforehand that the solution will be sufficiently regular . then , the solution can be expanded into a series of certain basis functions which are globally analytic . we refer readers interested in these methods to a recent review paper by bonazzola et al . . the more difficult problem in the determination of the initial data is the fact that the curvature components are obtained by successive division by . since  vanishes on the boundary one has to compute a value which is of the form 0/0 . although the theorem tells us that this is well defined ( provided that some boundary conditions are satisfied ) , it still poses a numerical problem , because a straightforward implementation of lhpital s rule looses accuracy . the first one   is again based on pseudo - spectral methods . essentially , in this approach the multiplication with the function  is expressed as a linear operation between the expansion coefficients of  and the other factor . the pseudo - inversion of this linear operator ( in the sense of finding its moore - penrose inverse ) then corresponds to division by . although this methods seems to work well , it is not generally applicable because it depends heavily on the geometric setup . a more general method has been devised by hbner . here , the problem of dividing by  has been reformulated into a problem for solving a second - order pde for the quotient . this pde is similar to the yamabe equation in the sense that it also degenerates on the boundary where  vanishes , but it is linear . the initial data constructed in one of the ways described above must be evolved with the propagation equations extracted from the conformal field equations . as we remarked upon above , there are several ways to write the evolution equations , and therefore there are also several ways to implement them numerically . one is based on the adm - like formulation of the conformal field equations , while the other uses a tetrad formalism . in the adm - like formulation developed by hbner , one introduces the coefficients of the intrinsic metric on the space - like slices and their extrinsic curvature as variables , while in the other formulation by frauendiener , a tetrad and connection coefficients with respect to that tetrad are the basic variables . not much difference concerning efficiency and usefulness between these formulations is to be expected , provided the numerical schemes are the same . one minor advantage of the tetrad formulation might be the additional freedom in the orientation of the tetrad which can be used on  to make the computational tetrad agree with the preferred tetrad on  ( see below ) . one of the schemes which have been used is the higher dimensional lax - wendroff scheme , also called the rotated richtmyer scheme , a discretization scheme with second order accuracy . it has been employed alone   or together with strang splitting   to treat the principal part of the equations differently from the source part . since a second order scheme requires much more computing resources compared to higher order methods to achieve the same accuracy , hbner started to use the method of lines   with a fourth order scheme to compute the spatial derivatives and fourth order runge - kutta for the evolution in time . he reports   that the fourth order method is very much superior to the second order scheme in terms of efficiency . ( the feasibility of the method of lines in relativity has been studied by hungerbhler   using pseudo - spectral methods for the spatial derivatives and a combination of adams - bashforth and adams - moulton schemes for the time evolution . ) the conformal field equations and the propagation equations derived from them are quasi - linear . this implies that the characteristics of the system depend on the current solution and this , in turn , means that one has to be able to change the time - step t between successive time - slices in order to keep a stable evolution scheme . this is necessary because the schemes are explicit schemes and , therefore , subject to the courant - friedrichs - lewy condition which states that the numerical domain of dependence of a point should always include the analytical domain of dependence . this requirement already excludes the popular leapfrog scheme which is nevertheless used sometimes also for evolving the einstein equations . a general criterion for computing the maximal time - step allowed in each iteration in arbitrary dimension has been derived in . another important point in the development of evolution codes is the numerical treatment of the boundaries . as explained already above , it is one of the advantages in the conformal approach that the outer boundary is not as influential as it is in the conventional approach using the standard cauchy problem . it was also pointed out that in this case it is enough to impose a boundary condition which results in a numerically stable code because the outer boundary is located in the unphysical region and , therefore , can not influence the physical space - time . the proper way to treat the boundary is to prescribe conditions which are compatible with the full conformal field equations , in particular with their restriction to the boundary manifold . this has not been done so far . since the outer boundary is not important for the physical effects , other ways of dealing with the boundary one way is to forget about the restrictions of the conformal field equations to the boundary and to analyze the possible boundary conditions for the propagation equations . to first order then a sufficiently general boundary condition will be obtained by specifying the ingoing fields in terms of the outgoing ones . although this boundary treatment is not compatible with the restriction of the conformal field equations to the boundary , it is compatible with the evolution equations . this means that the evolution can remain stable although the solution will not satisfy the constraints in the domain of influence of the boundary , which , however , is always in the unphysical part of space - time . another , very clever method for dealing with the boundary has been found by hbner . he realized that it is sufficient to solve the conformal field equations in the physical space - time only , and not necessary to solve them in the unphysical region as long as the characteristics remain such that the information created in the unphysical part of the computational domain can not reach the physical part . consequently , in his treatment the grid is divided into three zones : the inner zone , the outer zone , and a transition zone . the inner zone covers the physical space - time ( flagged by a positive conformal factor ) and some part of the adjacent unphysical region . on this part of the grid the conformal field equations are solved . in the outer zone , which is located in a neighbourhood of the grid boundary , one solves an advection equation which propagates outwards , off the grid . in the transition zone the effect is that the boundary condition which has to be imposed on the grid boundaries is very simple and that the noise which is generated in the transition region is propagated away from the physical region outwards towards the grid boundary . our next point is concerned with the extraction of the radiative information from the numerically generated data . this is the part of the entire numerical process where the superiority of the conformal approach becomes apparent . first of all , one needs to find  on the current time - slice . since  is the surface on which the conformal factor  vanishes and since  is explicitly known during the evolution the location of  is a simple task . the next problem is concerned with the orientation of the tetrad on . the asymptotic quantities are defined with respect to a specific geometrically characterized tetrad , a bondi frame . therefore , one needs to find the transformation from one to the other at each point of . without going into too much details ( see [ 39 , 42 , 41 ] ) we remark that most asymptotic quantities , in particular the radiation field , are of a local character so they can be read off without constructing a bondi frame . this is rather fortunate because there are global issues involved in the transformation from the computational tetrad to the bondi frame . these have implications for the determination of global quantities like the bondi energy - momentum four - vector , but they have no effect on the radiation field , which is defined as that ( complex ) component of kabcd which is entirely intrinsic to  : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\psi _ 4 } = { k_{abcd}}{n^a}{\bar m^b}{n^c}{\bar m^d}.$$\end{document } 
 here , n is a null - vector tangent to the generators of  , i.e. , na  a , and m is any complex space - like null - vector which is orthogonal to n. it is useful to require the space - like vector to be tangent to the intersection of  with the current time - slice . augmenting these two vectors by a further real null - vector l yields a null - tetrad which is fixed up to rotations in the ( two - dimensional ) tangent space of that intersection and boosts in the plane orthogonal to it . the behaviour of 4 under these transformations is that of a ghp - weighted quantity [ 67 , 117 ] with boost weight 2 and spin weight 2 . however , 4 really depends only on the null - vector n. for , suppose we perform a null - rotation around n , then m transforms into m + n for some complex valued function  on . but 4 is invariant under this transformation . so in order to find 4 it is only necessary to transform from the given computational tetrad to the tetrad specified above which is rather straightforward . in fact , the computation of 4 involves only the combination of certain components of the gravitational field with powers of a. the final step in the correct determination of the radiation is to find the correct time parameter . suppose we follow a specific null generator of  crossing through successive time - slices . on each slice then we obtain the radiation emitted by the source into the direction specified by the generator as a function of our computational time parameter . since the time coordinate is rather arbitrary , this means that the wave form determined so far has no physical meaning . the problem is already present in maxwell s theory : suppose we have an emitter which sends out a pure sine wave . a detector far away from the source can not determine the absolute frequency of the signal because the relative velocity of emitter and receiver might be non - zero but the detector should also find a pure sine signal . however , this will be true only if the detector records the signal as a function of proper time . any other time parameter along the detector s world - line will not produce a pure sine . what one needs to do in the general case is to select among all parameters along the generator a specific , geometrically distinguished one , namely a bondi parameter . a generator and such a parameter along it can be understood as a certain limit of freely falling observers with proper time clocks as they move towards infinity . bondi parameters are obtained as solutions of an ordinary , linear , second order differential equation , which is conformally invariant . the computation of the bondi energy - momentum is a global procedure , i.e. , it depends on properties of the entire cut of  with the current time - slice . first , one needs to obtain the asymptotic translation group ( see e.g. ) on each cut . this provides four functions on the cut which are then , in a second step , integrated against the  mass aspect  which is another function obtained from the  coulomb  part 2 of the gravitational field , and the  news function  which is a combination of components of the ricci tensor and connection coefficients . the first step , the determination of the translation group , is the global step because it involves solving a second order elliptic equation on the cut . all these procedures for finding the relevant data on  have been worked out analytically and they have also been tested ( at least in part ) numerically . the tests have been performed under the assumption that null - infinity admits toroidal cuts , which has the advantage that one can actually compare the numerical results with analytical expressions because a whole class of exact solutions [ 135 , 80 ] is known to exist . admittedly , such space - times are rather unphysical , but since most of the extraction procedures are local there is no doubt that they will also work in more realistic cases . let us now discuss some of the achievements of the approach based on the conformal field equations in some more detail . they range from investigations of gravitational fields coupled to a scalar field in spherical symmetry to pure gravitational interactions studied in two space dimensions and , recently , in the most general case of three space dimensions . the numerical treatment of the conformal field equations was started by hbner   who analyzed the asymptotic structure of spherically symmetric space - times in which a scalar field propagates under the influence of the gravitational attraction due to its own energy density . this is a system which has been investigated rigorously in detail earlier by christodoulou ( see e.g.   and the references therein ) . so the numerical results can be judged against the very detailed information found by analytical work . in all cases considered , the numerical results agreed with the analytical ones . as a specific example , let us look at figure 12 ( taken from ) where the  upper part  of a space - time with a singularity is shown . it is obtained from the numerical evolution of initial data which are supercritical in the sense that the initial energy density ( specified by a parameter a ) is so large that the scalar field collapses down to a black hole . the dashed line is  , while the thin line is the locus of vanishing divergence of outgoing light rays , i.e. an apparent horizon . the dashed line is  , while the thin line is the locus of vanishing divergence of outgoing light rays , i.e. an apparent horizon . this is indicated by the appearance of trapped surfaces and the subsequent formation of a singularity . the boundary of the region where the trapped surfaces exist is indicated by the thin line in the figure . it is the apparent horizon on which the divergence out of the outgoing light rays vanishes . another part of the investigation was concerned with the radiation at infinity . in figure 13 ( also from ) the scalar radiation field at null - infinity as a function of proper time of an observer on  is shown . decay of the radiation at null - infinity . in this example , the initial data was subcritical so that the scalar field , which initially collapses , subsequently disperses again . note the long time - scale , ranging over approximately six orders of magnitude in proper time . this is a remarkable achievement because so far no other numerical method has been able to monitor the evolution of relativistic space - times for such a long period of time . the next step in the application of the conformal field equations to numerical problems has been the implementation of 2-d codes for the solution of a3-like space - times [ 38 , 39 ] . these provide the first examples of vacuum space - times with gravitational radiation . of course , they can not be taken seriously as models of isolated systems because the topology of their  is not the physically distinguished s  . however , they provide important test cases for the codes and in particular for methods to extract radiation . since exact solutions with this kind of global structure are known [ 135 , 80 ] one can again compare the numerical results with their exact counterparts . the radiation field 4 and the bondi mass for a particular case figure 14the radiation field 4 
 and the bondi mass for a radiating a3-like space - time . the radiation field 4 
 and the bondi mass for a radiating a3-like space - time . in both diagrams note that this is the first time that a fully non - linear wave - form has been computed which agrees with an exact solution . as a final example of the conformal method in numerical relativity we consider the schwarzschild space - time which has recently been evolved with hbner s 3-d code . figure 15 is a numerical version of the kruskal diagram , i.e. a diagram for the conformal structure of the schwarzschild solution . what is clearly visible here are the two null - infinities ( blue lines ) and the horizons ( red lines ) . null - geodesic , i.e. the locus where the kruskal null - coordinates u and v ( see e.g. ) are equal . the dashed lines are  right going  null - geodesics , moving away from the left - hand . the diagram shows the cross - over where the two horizons ( and the central line ) intersect and , accordingly , we see a large part of the region iii , which is below the cross - over , the regions i and iv with their corresponding  s and some part of region ii where the future singularity is located . the non - symmetric look of the diagram is , of course , due to the fact that the coordinates used in the code have nothing to do with the kruskal coordinates with respect to which one usually sees the kruskal diagram of the extended schwarzschild solution . in all three cases mentioned here , there is a clear indication that long - time studies of gravitational fields are feasible . all three cases have been checked against exact results ( exact solutions or known theorems ) so that there is no doubt that the numerical results are correct . these contributions show beyond any reasonable doubt that the conformal field equations can not only be used for the analytical discussion of global properties of space - times but also for the numerical determination of semi - global solutions . clearly the problems with the artificial boundary have evaporated , the asymptotic region can accurately be determined and the wave - forms can reliably be computed . together with the analysis of i there is now good hope that the numerical computation of global space - times can be achieved in the near future . now we want to discuss some of the problems which still have to be dealt with . since the conformal approach has been tested numerically so far mostly in 2d cases with unphysical global structure , it is necessary to implement full 3d codes to run on general enough data . work on this is well underway and first results which confirm the expectations based on the two - dimensional codes are available . apart from the requirement of having a general 3d - code , there are other problems which need more consideration . . it would be interesting to see how the results of   translate to the conformal field equations . their implementation could result in stable codes which do not need any additional transition zones beyond  , and which are compatible with the evolution equation unlike the procedure used in . the necessity of dividing by the conformal factor to construct the initial data is annoying . while it has been possible to do this in the spherically symmetric case , no result is available for the general case . it would be very desirable to have another way of constructing the initial data , because then one could construct more easily data which evolve into space - times with multiple black holes . as a last problem in connection with the conformal approach one should mention that friedrich s work on i provides a way to evolve cauchy data specified on an asymptotically euclidean hypersurface to hyperboloidal initial data . a code which does that kind of evolution can provide the initial data for an evolution code for the hyperboloidal initial value problem . surely there will be difficult problems in the numerical treatment of the transport equations related to the total characteristic at spatial infinity . but the work on this problem is worthwhile because it would provide the final step to the ultimate goal of a global simulation of an isolated system . a problem which affects all the work in numerical relativity today is the obscure nature of the gauge conditions . currently there is not much understanding of the effects of a gauge condition on the resulting nature of the coordinates ( frame , conformal factor ) . most of the work done on these problems is related to the choice of a lapse function ; in particular several proposals have been made for selecting a time coordinate . these are mostly dictated by formal considerations like the need of making a system hyperbolic or of easy implementation . to some extent this is justified because the physics can not depend on the coordinates which are used . but it is also well known that there are  good  and  bad  coordinates . good  and  bad  means depends to a large extent on what the goal is . ideally , coordinates should be tied to the geometry so that they obtain a more invariant nature . in 1d cases one can set up a system of double null - coordinates ( and a derived system of time and space coordinates ) this provides a gauge which is good as long as the geometry is well behaved . but , unfortunately , this gauge can not be generalized in a straightforward way to higher dimensions ( some attempts have been made in ) . probably one should assume a pragmatic viewpoint towards the problem of finding appropriate coordinates in the sense that one should regard the gauge sources as knobs which have to be adjusted by trial and error . maybe it is possible , at least to some extent , to let the code do the  twiddling  automatically . this requires that one should be able to formulate criteria for a  good solution  which can be checked by the computer . furthermore , it is also necessary that a change in the gauge sources does not change the characteristics of the system , because otherwise it is easy to get into situations where the system is not hyperbolic anymore . finally , there is no doubt that the issue of the correct boundary condition for codes based on the standard einstein equations also needs more attention . although it is not clear what these conditions mean physically , chances are good that they will produce stable codes ( see e.g. the discussion of numerical boundary conditions in ) . if this is the case then one could try to check which ones work best by comparison with exact radiating solutions , like the boost - rotation symmetric solutions discussed in . another test should be a numerical comparison with the data computed by a hyperboloidal evolution code . such tests are important because so far the standard codes have only been tested against linearized solutions . this is the regime where one would expect them to work because the boundary condition is still benign . in this way , one could not only select physically reasonable boundary conditions for the standard codes , one could also check how well they perform in comparison with the conformal codes . in particular , one could see how accurate the radiation extraction can be done with those codes and whether the accuracy is good enough for ligo wave forms . then one can compute them with a safe conscience using the standard codes , provided they are more efficient than the conformal codes . in this appendix we show how to perform the reduction process for the conformal field equations to obtain the symmetric hyperbolic system of evolution equations and the constraints . we assume that we are given a time - like unit vector t with respect to which the reduction is done . any vector v may be decomposed into parts perpendicular and parallel to t , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v^a } = { \hat v^a } + v{t^a},\;\;\;\;\;{\rm{with}}\;{\hat v^a}{t_a } = 0,\;\;\;\;\;v = { v^a}{t_a},$$\end{document } 
 and similarly , for one - forms wa . we call a vector spatial with respect to t if it is orthogonal to t. in particular , the metric gab gives rise to a spatial metric hab by the decomposition 
 ( 37)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${g_{ab } } = { h_{ab } } + { t_a}{t_b}.$$\end{document } 
 the volume four - form abcd which is defined by the metric also gives rise to a decomposition as follows : 
 ( 38)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\varepsilon _ { abcd } } = 4{\varepsilon _ { [ abc}}{t_{d]}}\:\ ; \leftrightarrow \:\;\;{\varepsilon _ { abc } } = { \varepsilon _ { abcd}}{t^d}.$$\end{document } 
 the covariant derivative operator a is written analogously , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a } = { d_a } + { t_a}d,$$\end{document } 
 thus defining two derivative operators da with tda = 0 and d = ta . the covariant derivative of t itself is an important field . it gives rise to two component fields defined by 
 ( 39)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\chi _ a}^b = { d_a}{t^b},\:\;\;\;\;\;\;{\chi ^b } = d{t^b}.$$\end{document } 
 note that a is spatial in both its indices and that there is no symmetry implied between the two indices . it is useful to define two new derivative operators  and a by the following relations : 
 ( 40)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}{v_b } = { d_a}{v_b } + { t_b}{\chi _ a}^e{v_e } - { \chi _ { ab}}{t^e}{v_e},\;\;\;\;\;\;\;\partial { v_b } = d{v_b } + { t_b}{\chi ^e}{v_e } - { \chi _ b}{t^e}{v_e}.$$\end{document } 
 these operators have the property that they are compatible with the spatial metric hab and that they annihilate t and ta . if t is the unit normal of a hypersurface , i.e. , if ta is hypersurface orthogonal , then ab is symmetric , hab is the induced ( negative definite ) metric on the hypersurface , and a is its levi - civita connection . in general this is not the case and so the operator a possesses torsion . in particular , we obtain the following commutators ( acting on scalars and spatial vectors ) : 
 ( 41)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$2{\partial _ { [ a}}{\partial _ { b]}}f + 2{\chi _ { [ ab]}}\partial f = 0,$$\end{document } 
 ( 42)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$2{\partial _ { [ a}}{\partial _ { b]}}{v_c } + 2{\chi _ { [ ab]}}\partial { v_c } =   - { h_a}^r{h_b}^s{h_c}^t{r_{rst}}^d{v_d } + 2{\chi _ { [ a|c|}}{\chi _ { b]}}^e{v_e},$$\end{document } 
 ( 43)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$[\partial , \:{\partial _ a}]f + { \chi _ a}\partial f + { \chi _ a}^c{\partial _ c}f = 0,$$\end{document } 
 ( 44)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$[\partial , \:{\partial _ a}]\;{v_b } + { \chi _ a}\partial { v_b } + { \chi _ a}^c{\partial _ c}{v_b } = { \chi _ b}{\chi _ a}^e{v_e } - { \chi _ { ab}}{\chi ^e}{v_e } - { t^c}{h_a}^e{h_b}^f{r_{cef}}^d{v_d}.$$\end{document } 
 these commutators are obtained from the commutators between the derivative operators da and d by expressing them in terms of a and  on the one hand , and by the four - dimensional connection a on the other hand . this procedure yields two equations for the derivatives of t , 
 ( 45)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$2{\partial _ { [ a}}{\chi _ { b]c } } - 2{\chi _ { [ ab]}}{\chi _ c } + { h_a}^r{h_b}^s{h_c}^t{r_{rst}}^d{t_d } = 0,$$\end{document } 
 ( 46)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \chi _ { ab } } - { \partial _ a}{\chi _ b } + { \chi _ a}^c{\chi _ { cb } } + { \chi _ a}{\chi _ b } - { t^c}{t^d}{h_a}^r{h_b}^s{r_{crds}}.$$\end{document } 
 the information contained in the commutator relations and in the equations ( 45 ) and ( 46 ) is completely equivalent to the cartan equations for a which define the curvature and torsion tensors . this completes the preliminaries and we can now go on to perform the splitting of the equations . out intention is to end up with a system of equations for all the spatial parts of the fields . in order not to introduce too many different kinds of indices , all indices refer to the four - dimensional space - time , but they are all spatial , i.e. , any transvection with ta and t vanishes . if we introduce hypersurfaces with normal vector t then there exists an isomorphism between the tensor algebra on the hypersurfaces and the subalgebra of spatial four - dimensional tensors . we decompose the fields into various spatial parts and insert these decompositions into the conformal field equations defined by ( 17 , 18 , 19 , 20 , 21 ) . the fields are decomposed as follows : 
 ( 47)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${k_{abcd } } = 4{t_{[a}}{e_{b][c}}{t_{d ] } } + 2{t_{[a}}{b_{b]}}^e{\varepsilon _ { ecd } } - 2{\varepsilon _ { ab}}^e{b_{e[c}}{t_{d ] } } + { \varepsilon _ { abe}}{e^{ef}}{\varepsilon _ { fcd}},$$\end{document } 
 ( 48)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\phi _ { ab } } = { \phi _ { ab } } + 2{t_{(a}}{\phi _ { b ) } } + { t_a}{t_b}\phi , $ $ \end{document } 
 ( 49)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\sigma _ a } = { \sigma _ a } + { t_a}\sigma .$$\end{document } 
 the function  is fixed in terms of ab because ab is trace - free . inserting the decomposition of kabcd into equation ( 18 ) , decomposing the equations into various spatial parts and expressing derivatives in terms of the operators a and  yields four equations : 
 ( 50)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial ^a}{b_{ac } } + { \varepsilon _ { ead}}{e_c}^e{\chi ^{ad } } + { \varepsilon _ { ecd}}{e_a}^e{\chi ^{ad } } = 0,$$\end{document } 
 ( 51)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial ^a}{e_{ac } } - { \varepsilon _ { ead}}{b_c}^e{\chi ^{ad } } - { \varepsilon _ { ecd}}{b_a}^e{\chi ^{ad } } = 0,$$\end{document } 
 ( 52)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { e_{bc } } + { \varepsilon _ { ae(b}}{\partial ^a}{b^e}_{c ) } = \\ 
 \quad \;\;\;\,\ , - 2\chi { e_{bc } } + 2{\chi ^a}_{(b}{e_{c)a } } - { h_{bc}}{e_{ad}}{\chi ^{ad } } + { e_{a(b}}{\chi _ { c)}}^a - 2{\varepsilon _ { ea(b}}b_{c)}^e{\chi ^a } , 
 \end{array}$$\end{document } 
 ( 53)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { b_{bc } } - { \varepsilon _ { ae(b}}{\partial ^a}e_{c)}^e = \\ \quad \quad \,\ , - 2\chi { b_{bc } } + 2{\chi ^a}_{(b}{b_{c)a } } - { h_{bc}}{b_{ad}}{\chi ^{ad } } + { b_{a(b}}{\chi _ { c)}}^a + 2{\varepsilon _ { ea(b}}{e^e}_{c)}{\chi ^a}. 
 \end{array}$$\end{document } 
 here we have defined  = a . treating the other fields and equations in a similar way , we obtain equation ( 17 ) in the form of four equations : 
 ( 54)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \phi _ a } - { \partial _ a}\phi   = { \chi ^e}{\phi _ { ae } } - 2{\chi _ a}^e{\phi _ e } - { \partial _ a}\lambda   - { e_a}^d{\sigma _ d},$$\end{document } 
 ( 55)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ { [ c}}{\phi _ { a ] } } - \chi _ { [ c}^e{\phi _ { a]e } } + { \chi _ { ca}}\phi   - \frac{1}{2}{\varepsilon _ { cae}}{b^{ed}}{\sigma _ d } = 0,$$\end{document } 
 ( 56)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \phi _ { ab } } - { \partial _ a}{\phi _ b } =   - 2{\chi _ { ( a}}{\phi _ { b ) } } - { \chi _ a}^e{\phi _ { be } } + { \chi _ { ab}}\phi   + { h_{ab}}\partial \lambda   + { e_{ab}}\sigma   + { b_a}^e{\varepsilon _ { bde}}{\sigma ^d},$$\end{document } 
 ( 57)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ { [ c}}{\phi _ { a]b } } = { \chi _ { [ ac]}}{\phi _ b } - { \phi _ { [ a}}{\chi _ { c]b } } + { h_{b[a}}{\partial _ { c]}}\lambda   - \frac{1}{2}{\varepsilon _ { cae}}{b^e}_b\sigma   + \frac{1}{2}{\varepsilon _ { cae}}{e^{ef}}{\varepsilon _ { fbd}}{\sigma ^d}.$$\end{document } 
 the equation ( 19 ) for the conformal factor are rather straightforward . we obtain 
 ( 58)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \omega   = \sigma , $ $ \end{document } 
 ( 59)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}\omega   = { \sigma _ a},$$\end{document } 
 while equation ( 20 ) yields four equations : 
 ( 60)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}{\sigma _ b } + { \chi _ { ab}}\sigma   - { h_{ab}}s + \omega { \phi _ { ab } } = 0,$$\end{document } 
 ( 61)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ b}\sigma   - { \chi _ b}^e{\sigma _ e } + \omega { \phi _ b } = 0,$$\end{document } 
 ( 62)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \sigma _ b } =   - { \chi _ b}\sigma   - \omega { \phi _ b},$$\end{document } 
 ( 63)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \sigma   = { \chi ^e}{\sigma _ e } + s - \omega \phi .$$\end{document } 
 finally , the equation ( 21 ) for s gives two equations 
 ( 64)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial s =   - { \phi _ b}{\sigma ^b } - \phi \sigma   + \omega \partial \lambda   + 2\sigma \lambda , $ $ \end{document } 
 ( 65)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}s + { \phi _ { ab}}{\sigma ^b } + { \phi _ a}\sigma   - \omega { \partial _ a}\lambda   = 0.$$\end{document } 
 this completes the gauge independent part of the equations . in order to deal with the gauges we now have to introduce an arbitrary tetrad and arbitrary coordinates we extend the time - like unit vector to a complete tetrad ( t , eia ) with taeia = 0 for i = 1 , 2 , 3 . let x with  = 0,1 , 2 , 3 be four arbitrary functions which we use as coordinates . application of  and a to the coordinates yields 
 ( 66)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { x^\mu } = { c^\mu } , \:\;\;\;\;\;\;{\partial _ a}{x^\mu } = c_a^\mu .$$\end{document } 
 the four functions c and the four one - forms ca may be regarded as the 16 expansion coefficients of the tetrad vectors in terms of the coordinate basis  = /x because of the identity 
 ( 67)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial   = { c^\mu } { \partial _ \mu } \;\ ; \leftrightarrow \;\;\partial { x^\mu } = { c^\mu } , $ $ \end{document } 
 ( 68)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$e_i^a{\partial _ a } = c_i^\mu { \partial _ \mu } \;\ ; \leftrightarrow \;\;\;e_i^a{\partial _ a}{x^\mu } = e_i^ac_a^\mu .$$\end{document } in a similar spirit , we apply the derivative operators to the tetrad and obtain 
 ( 69)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial e_i^a = { \lambda ^a}_be_i^b,$$\end{document } 
 ( 70)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ c}e_i^a = \gamma _ { c\;\;b}^{\;\;a}e_i^b.$$\end{document } 
 again , transvection with t on any index of b and cb vanishes . furthermore , both ab and cab are antisymmetric in their ( last two ) indices . together with the 12 components of  and a these fields provide additional 12 components which account for the 24 connection coefficients of the four - dimensional connection a with respect to the chosen tetrad . since we will keep the tetrad fixed here , we may , however , regard them as defining tensorfields whose components happen to coincide with them in the specified tetrad . in order to extract the contents of the first of cartan s structure equations one needs to apply the commutators ( 41 ) and ( 43 ) to the coordinates to obtain 
 ( 71)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial c_a^\mu   - { \partial _ a}{c^\mu } =   - { \chi _ a}{c^\mu } - \chi _ a^bc_b^\mu , $ $ \end{document } 
 ( 72)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ b}c_a^\mu   + 2{\chi _ { [ ab]}}{c^\mu } = 0.$$\end{document } 
 similarly , the second of cartan s structure equations is exploited by applying the commutators to the tetrad vectors . equation ( 16 ) is then used to substitute for the riemann tensor in terms of the gravitational field , the trace - free part of the ricci tensor , and the scalar curvature . apart from the equations ( 45 ) and ( 46 ) , which come from acting on t , this procedure yields 
 ( 73)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { \gamma _ { abc } } - { \partial _ a}{\lambda _ { bc } } = \\ 
 \quad \;\;2{\gamma _ { a[b}}^e{\lambda _ { c]e } } - { \chi _ a}{\lambda _ { bc } } - { \chi _ a}^e{\gamma _ { ebc } } - 2{\chi _ { a[b}}{\chi _ { c ] } } + \omega { b_a}^e{\varepsilon _ { ebc } } + 2{h_{a[b}}{\phi _ { c ] } } , 
 \end{array}$$\end{document } 
 ( 74)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \partial _ a}{\gamma _ { bcd } } - { \partial _ b}{\gamma _ { acd } } + 2{\chi _ { [ ab]}}{\lambda _ { cd}}\\ 
 \quad \;\;\ , + 2{\chi _ { [ a|d}}{\chi _ { b]c } } + \omega { \varepsilon _ { ab}}^e{e_{ef}}{\varepsilon ^f}_{cd } - 2{h_{c[a}}{\phi _ { b]d } } + 2{h_{d[a}}{\phi _ { b]c } } - 4{h_{c[a}}{h_{b]d}}\lambda   = 0 . \end{array}$$\end{document } 
 now we have collected all the equations which can be extracted from the conformal field equations and cartan s structure equations . what remains to do is to separate them into constraints and evolution equations . before doing so , we notice that we do not have enough evolution equations for the tetrad components and the connection coefficients . the remedy to this situation the missing equations for the coordinates are obtained by applying the dalembert operator to the coordinates . expressing the wave operator in terms of  and a yields the additional equations 
 ( 75)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { c^\mu } + { \partial ^a}c_a^\mu   = { \chi ^a}c_a^\mu   + \chi { c^\mu } - { f^\mu } .$$\end{document } 
 in order to find the missing equations for the tetrad we need to compute the gauge source functions 
 ( 76)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { f_{0k } } = { g_{ab}}{\nabla ^c}({t^a}{\nabla _ \end{array}$$\end{document } 
 in a similar way as explained above , we may regard these functions as components of tensorfields fab and fa whose components happen to agree with them in the specified basis . thus , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${f_{ik } } = { f_{ab}}e_i^ae_k^b,\:\;\;\;\;\;\;{f_{0i } } = { f_a}e_i^a.$$\end{document } 
 computing these tensorfields from ( 76 ) gives 
 ( 77)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${f_{ab } } = \partial { \lambda _ { ab } } - { \partial ^e}{\gamma _ { eab } } - { \chi ^e}{\gamma _ { eab } } + \chi { \lambda _ { ab } } + 2{\gamma _ { ec[a}}{\gamma ^{ec}}_{b]},$$\end{document } 
 ( 78)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${f_a } =   - \partial { \chi _ a } - { \partial ^e}{\chi _ { ea } } + { \chi ^e}{\chi _ { ea } } - { \chi _ a}{\chi ^e}_e - { \chi _ { ec}}{\gamma ^{ec}}_a - { \chi ^c}{\lambda _ { ca}}.$$\end{document } 
 now we are ready to collect the constraints : 
 ( 79)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ b}c_a^\mu   + 2{\chi _ { [ ab]}}{c^\mu } = 0,$$\end{document } 
 ( 80)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \partial _ a}{\gamma _ { bcd } } - { \partial _ b}{\gamma _ { acd } } + 2{\chi _ { [ ab]}}{\lambda _ { cd}}\\ 
 \quad \;\;\ , + 2{\chi _ { [ a|d}}{\chi _ { b]c } } + \omega { \varepsilon _ { ab}}^e{e_{ef}}{\varepsilon ^f}_{cd } - 2{h_{c[a}}{\phi _ { b]d } } + 2{h_{d[a}}{\phi _ { b]c } } - 4{h_{c[a}}{h_{b]d}}\lambda   = 0 , 
 \end{array}$$\end{document } 
 ( 81)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}{\chi _ { bc } } - { \partial _ b}{\chi _ { ac } } - 2{\chi _ { [ ab]}}{\chi _ c } - { \varepsilon _ { ab}}^e{b_{ec } } + 2{h_{c[a}}{\phi _ { b ] } } = 0,$$\end{document } 
 ( 82)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ c}{\phi _ a } - { \partial _ a}{\phi _ c } - 2{\chi _ { [ c}}^e{\phi _ { a]e } } + 2{\chi _ { [ ca]}}\phi   - { \varepsilon _ { cae}}{b^{ed}}{\sigma _ d } = 0,$$\end{document } 
 ( 83)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 { \partial _ c}{\phi _ { ab } } - { \partial _ a}{\phi _ { cb } } - 2{\chi _ { [ ac]}}{\phi _ b } + 2{\phi _ { [ a}}{\chi _ { c]b}}\\ 
 \quad \quad \ ; - 2{h_{b[a}}{\partial _ { c]}}\lambda   + { \varepsilon _ { cae}}{b^e}_b\sigma   - { \varepsilon _ { cae}}{e^{e\,f}}{\varepsilon _ { fbd}}{\sigma ^d } = 0 , 
 \end{array}$$\end{document } 
 ( 84)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial ^a}{b_{ac } } + { \varepsilon _ { ead}}{e_c}^e{\chi ^{ad } } + { \varepsilon _ { ecd}}{e_a}^e{\chi ^{ad } } = 0,$$\end{document } 
 ( 85)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial ^a}{e_{ac } } - { \varepsilon _ { ead}}{b_c}^e{\chi ^{ad } } - { \varepsilon _ { ecd}}{b_a}^e{\chi ^{ad } } = 0,$$\end{document } 
 ( 86)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}\omega   - { \sigma _ a } = 0,$$\end{document } 
 ( 87)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ b } + { \chi _ { ab}}\sigma   - { h_{ab}}s + \omega { \phi _ { ab } } = 0,$$\end{document } 
 ( 88)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ b}\sigma   - { \chi _ b}^e{\sigma _ e } + \omega { \phi _ b } = 0,$$\end{document } 
 ( 89)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\partial _ a}s + { \phi _ { ab}}{\sigma ^b } + { \phi _ a}\sigma   - \omega { \partial _ a}\lambda   - 2\lambda { \sigma _ a } = 0.$$\end{document } 
 finally , we collect the evolution equations : 
 ( 90)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial c_a^\mu   - { \partial _ a}{c^\mu } =   - { \chi _ a}{c^\mu } - { \chi _ a}^bc_b^\mu , $ $ \end{document } 
 ( 91)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { c^\mu } + { \partial ^a}c_a^\mu   = { \chi ^a}c_a^\mu   + \chi { c^\mu } - { f^\mu } .$$\end{document } 
 ( 92)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { \gamma _ { abc } } - { \partial _ a}{\lambda _ { bc } } = 2{\gamma _ { a[b}}^e{\lambda _ { c]e } } - { \chi _ a}{\lambda _ { bc}}\\ 
 \quad \quad \quad \quad \quad \quad   - { \chi _ a}^e{\gamma _ { ebc } } - 2{\chi _ { a[b}}{\chi _ { c ] } } + \omega { b_a}^e{\varepsilon _ { ebc } } + 2{h_{a[b}}{\phi _ { c ] } } , 
 \end{array}$$\end{document } 
 ( 93)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \lambda _ { ab } } - { \partial ^e}{\gamma _ { eab } } = { \chi ^e}{\gamma _ { eab } } - \chi { \lambda _ { ab } } - 2{\gamma _ { ec[a}}{\gamma ^{ec}}_{b ] } + { f_{ab}},$$\end{document } 
 ( 94)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \chi _ a } + { \partial ^e}{\chi _ { ea } } = { \chi ^e}{\chi _ { ea } } - { \chi _ a}{\chi ^e}_e - { \chi _ { ec}}{\gamma ^{ec}}_a - { \chi ^c}{\lambda _ { ca } } - { f_a},$$\end{document } 
 ( 95)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \chi _ { ab } } - { \partial _ a}{\chi _ b } =   - { \chi _ a}^c{\chi _ { cb } } - { \chi _ a}{\chi _ b } - \omega { e_{ab } } + { \phi _ { ab } } + { h_{ab}}(\phi   - 2\lambda ) , $ $ \end{document } 
 ( 96)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \phi _ a } - { \partial _ a}\phi   = { \chi ^e}{\phi _ { ae } } - 2{\chi _ a}^e{\phi _ e } - { \partial _ a}\lambda   - { e_a}^d{\sigma _ d},$$\end{document } 
 ( 97)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { \phi _ { ab } } - { \partial _ { ( a}}{\phi _ { b ) } } =   - 2{\chi _ { ( a}}{\phi _ { b ) } } - \chi \:({a^e}{\phi _ { b)e}}\\ 
 \quad \quad \quad \quad \quad \;\ ; + \chi ( ab)\phi   + { h_{ab}}\partial \lambda   + { e_{ab}}\sigma   + { b_{(a}}^e{\varepsilon _ { b)de}}{\sigma ^d } , 
 \end{array}$$\end{document } 
 ( 98)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { e_{bc } } + { \varepsilon _ { ae(b}}{\partial ^a}{b^e}_{c ) } =   - 2\chi { e_{bc}}\\ 
 \quad \quad \quad \ ; + 2{\chi ^a}_{(b}{e_{c)a } } - { h_{bc}}{e_{ad}}{\chi ^{ad } } + { e_{a(b}}{\chi _ { c)}}^a - 2{\varepsilon _ { ea(b}}{b^e}_{c)}{\chi ^a } , 
 \end{array}$$\end{document } 
 ( 99)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{l } 
 \partial { b_{bc } } - { \varepsilon _ { ae(b}}{\partial ^a}{e^e}_{c ) } =   - 2\chi { b_{bc}}\\ 
 \quad \quad \;\;\;\ , + 2{\chi ^a}_{(b}{b_{c)a } } - { h_{bc}}{b_{ad}}{\chi ^{ad } } + { b_{a(b}}{\chi _ { c)}}^a + 2{\varepsilon _ { ea(b}}{e^e}_{c)}{\chi ^a } , 
 \end{array}$$\end{document } 
 ( 100)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \omega   = \sigma , $ $ \end{document } 
 ( 101)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial { \sigma _ b } =   - { \chi _ b}\sigma   - \omega { \phi _ b},$$\end{document } 
 ( 102)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial \sigma   = { \chi ^e}{\sigma _ e } + s - \omega \phi , $ $ \end{document } 
 ( 103)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\partial s =   - { \phi _ b}{\sigma ^b } - \phi \sigma   + \omega \partial \lambda   + 2\sigma \lambda .$$\end{document } 
 this is the complete system of evolution equation which can be extracted from the conformal field equations . it is important to keep in mind that with our conventions the spatial metric hab is negative . define a = a/. the levi - civita connection a of the new metric is given in terms of the levi - civita connection a of gab by its action on an arbitrary vectorfield v , 
 ( 104)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a}{v^b } = { \tilde \nabla _ a}{v^b } + c_{ac}^b{v^c } = { \tilde \nabla _ a}{v^b } + ( \delta _ a^b{\upsilon _ c } + \delta _ c^b{\upsilon _ a } - { g^{bd}}{g_{ac}}{\upsilon _ d}){v^c}.$$\end{document } 
 from the action on vectorfields we can obtain the action on tensors of arbitrary valence in the usual way . it is useful to split the riemann tensor into several pieces which transform differently under conformal rescalings . we write 
 ( 105)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${r^{ab}}_{cd } = { c^{ab}}_{cd } - 4\delta _ { [ c}^{[a}{p^{b]}}_{d]}.$$\end{document } 
 the tensor cbcd is , of course , weyl s conformal tensor , characterised by the property of having the same symmetries as the riemann tensor with all traces vanishing . the other piece , the tensor pab , can be uniquely expressed in terms of the ricci tensor 
 ( 106)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p_{ab } } =   - \frac{1}{2}\left ( { { r_{ab } } - \frac{1}{6}{g_{ab}}r } \right ) = { \phi _ { ab } } - { g_{ab}}\lambda .$$\end{document } 
 the tensor ab is proportional to the trace - free part of the ricci tensor , while \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda   = 24r = { \textstyle{1 \over 4}}{g^{ab}}{p_{ab}}$$\end{document } is a multiple of the scalar curvature . under the conformal rescaling gab  gab = gab , the different parts of the curvature transform as follows : 
 ( 107)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde c^a}_{\;\;bcd } \mapsto { c^a}_{bcd } = \tilde c_{\;\;bcd}^a,$$\end{document } 
 ( 108)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde p_{ab } } \to { p_{ab } } = { \tilde p_{ab } } - { \tilde \nabla _ a}{\upsilon _ b } + { \upsilon _ a}{\upsilon _ b } - \frac{1}{2}{\tilde g_{ab}}{\tilde g^{cd}}{\upsilon _ c}{\upsilon _ d}.$$\end{document } 
 thus , the weyl tensor is invariant under conformal rescalings . when pab is expressed entirely in terms of the transformed quantities we get the relation 
 ( 109)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde p_{ab } } = { p_{ab } } + { \nabla _ a}{\upsilon _ b } + { \upsilon _ a}{\upsilon _ b } - \frac{1}{2}{g_{ab}}{g^{cd}}{\upsilon _ c}{\upsilon _ d},$$\end{document } 
 from which we can deduce ( note that the contractions are performed with the transformed metric ) 
 ( 110)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\tilde \phi _ { ab } } = { \phi _ { ab } } + \frac{1}{\theta } \left ( { { \nabla _ a}{\nabla _ b}\theta   - \frac{1}{4}{g_{ab}}\square\theta } \right),$$\end{document } 
 ( 111)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\tilde \lambda   = { \theta ^2}\lambda   - \frac{1}{4}\theta\square \theta   + \frac{1}{2}{\nabla ^a}\theta { \nabla _ a}\theta .$$\end{document } 
 next consider the bianchi identity [erab]cd = 0 . inserting the decomposition ( 105 ) and taking appropriate traces allows us to write it as two equations , 
 ( 112)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a}{c^a}_{bcd } + 2{\nabla _ { [ c}}{p_{d]b } } = 0,$$\end{document } 
 ( 113)\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla _ a}{p^a}_b - { \nabla _ b}{p^a}_a = 0.$$\end{document }
the notion of conformal infinity has a long history within the research in einstein s theory of gravity . today , 
 conformal infinity  is related with almost all other branches of research in general relativity , from quantisation procedures to abstract mathematical issues to numerical applications . 
 this review article attempts to show how this concept gradually and inevitably evolved out of physical issues , namely the need to understand gravitational radiation and isolated systems within the theory of gravitation and how it lends itself very naturally to solve radiation problems in numerical relativity . 
 the fundamental concept of null - infinity is introduced . 
 friedrich s regular conformal field equations are presented and various initial value problems for them are discussed . 
 finally , it is shown that the conformal field equations provide a very powerful method within numerical relativity to study global problems such as gravitational wave propagation and detection .