Stress Tensor on Null Boundaries

Using the Brown-York prescription for the definition of quasilocal gravitational energy-momentum tensor on a boundary and also complete canonical structure on a null boundary which has been found recently \cite{Aghapour:2018icu}, we propose a similar stress tensor on the null boundary. Then we exploit this stress tensor to compute the quasi-local energy and angular momentum for some well-known gravitational solutions. We have found that in addition to reference spacetime method for regularizing total energy, in the case of null boundary we can add a possible counterterm so avoiding embedding difficulties.


Introduction
In general relativity (GR), defining a local energy for gravitational field is problematic (for a review on the subject see [2]). The problem is rooted in the general covariance principle and the fact that the first derivative of metric can always vanish in a properly chosen local coordinate system. But in order to have a notion of local energy density with correct physical dimension, such quantity could be defined in terms of the metric and its first derivative.
Also, a related notion to energy density in field theories is the action functional. In fact, the standard covariant action in GR, i.e., the Hilbert-Einstein (HE) action, contains second order derivatives of the metric, which is another consequence of the general covariance principle.
As it is well-known, this action does not have a well-posed variational principle and needs to be complemented with additional terms defined on the boundaries of the spacetime. Brown and York [3] have pointed out that having an action with a well-posed variational principle, by using Hamilton-Jacobi analysis, one can define the quasilocal energy of the system.
A careful variation of the HE action in addition to suggesting the proper boundary term (boundary action), to ensure a well-posed variational principle, also provides dynamical degrees of the theory as well as its canonical structure (see e.g., [4]). An important point, in finding the complete canonical structure by such a procedure, is the necessity of the condition that one should not suppress any degree of freedom beforehand by imposing restrictions. The proper boundary action complementing the HE action, when the boundary is timelike or spacelike, is the well-known Gibbons-Hawking-York (GHY) term [5,6]. Applying variational principle in GR for null boundaries has been a subject of investigation in recent years [7][8][9][10].
In these papers, the authors are either interested in finding the proper boundary terms on null boundaries, and therefore ignored those terms that are fixed by boundary conditions, or imposed some restrictions on variations in such a way that the resulting canonical structure is not complete. One of such restrictions is that the variations keep the character of the boundary unchanged. But one of the metric degrees of freedom is responsible for variations which alter the boundary character. This can be easily seen as follows: If the boundary is specified by φ = const. for some scalar field φ, then the normal to the boundary is proportional to ∂ a φ. For a null boundary, we have ∂ a φ∂ a φ = 0, while for a general metric variation we find: ∂ a φ∂ a φ → ∂ a φ∂ a φ − δg ab ∂ a φ∂ b φ . (1) Therefore, unless we set variations of metric in direction of ∂ a φ to zero, the boundary doesn't remain null. Recently we have shown that such variations appear in canonical structure [1].
In addition, it has been shown that, in order to preserve such variations, a general doublefoliation framework is needed. The conjugate momentum to these variations is a scalar Ξ (for definition see section 3 or [1]). In this article we will show that this scalar provides the quasilocal energy density on a null boundary For the definition of quasilocal stress tensor on null boundary we follow the same method presented by Brown and York in the case of quasilocal stress tensor on timelike boundary which means that the derivative of the action with respect to the metric variations which are tangential to the boundary [3,11].
The content of the paper is organized as follows. In the next section we first briefly review the original Brown-York construction of quasi-local gravitational energy-momentum tensor.
In section 3, we follow the same analysis for the case where instead of timelike boundary we have a null one. Finally, in section 4 we calculate the energy and angular momentum in various space times using the proposed stress tensor.

The Brown-York Tensor
In this section we review the Brown-York [3] definition of gravitational stress-tensor on the boundary. To illustrate the main idea beyond Brown-York definition of quasi-local energy momentum, it is useful to start with an example in classical mechanics.

The Hamilton-Jacobi Method
Let L(q(t),q(t), t) be the classical Lagrangian for a particle. The action functional is I[q(t)] = t 2 t 1 L(q(t),q(t), t) dt, for initial and final configuration , i.e. (q 1 , t 1 ) and (q 2 , t 2 ), respectively. A general variation of this action for a given history is: Extremizing the action and imposing the boundary condition by fixing q and t at the end points, provides the equations of motion that are given by the first term in the integrand which are known as the Euler-Lagrange equations. Moreover, the Hamilton-Jacobi principal function S(q 1 , t 1 ; q 2 , t 2 ) is defined as the value of the action for a solution q(t) of the equation of motion from (q 1 , t 1 ) to (q 2 , t 2 ).
Thus, according to (2), the derivative of the principal function S with respect to q, i.e. δS δq , gives the canonical momenta p = ∂L ∂q while the derivative with respect to t, i.e. δS δt , yields the minus of energy: H = ( ∂L ∂qq − L). An important premonition is that in order to ensure that the above procedure works well, the variational principle must be well-posed beforehand. Indeed, a Lagrangian with no wellposed variational principle does not lead to the correct relations for momenta and energy of the system. For example, the action S 1 = 1 2 mẋ 2 dt leads to the equations of motion and correct momentum and energy for a free particle. However, by considering another action as S 2 = S 1 − 1 2 m d dt (xẋ)dt = − 1 2 mxẍdt , one could see that, although the difference with the first action is a total derivative and the equations of motion are unchanged, by varying this action one finds: which means that in order to get the equation of motion we need to determine both x andẋ at both end points. This leads to inconsistency with second order equation of motion, because according to that we need to just fix the position at the ends. Therefore, for this action the variational principle is not well-posed 2 . As a consequence, the derivative of the principal function does not lead to the correct definition for momentum and energy. However, we must note that all total derivatives do not spoil the variational principle, e.g. changing the action by S → S + t 2 t 1 dh dt dt for arbitrary function h(q(t), t) is allowed. Hence, the action and the principal function are not unique; as a result, momenta and energy are also not exclusive.
However, this arbitrariness can be fixed by choosing the zero point of energy, for example set the arbitrary function h so that it yields to zero energy for free particle at rest in the above example.

Stress tensor on timelike boundary
Having learned enough from the above simple mechanical example, lets begin with HE action in d dimension: in which R is the Ricci scalar. By varying the action, if we consider the boundary segments to be either time-like or space-like, one gets: where the sum is over each boundary segment B i and every corner C i at intersection of two neighboring segments of boundary. Moreover, K is the extrinsic curvature of each segment and ϑ is the angle or the boost parameter between segments , depending on their character.
In addition, δh ab is the variation of metric on each boundary while δq ab is the variation on each co-dimension two joints. For the details of the calculations and specially the method by which one can take into account the contribution of joints, see Refs. [8,10,11] or [1]. We Now consider a region in space-time as in fig 1 which is generated by evolution of a spacelike surface Σ from Σ 1 to Σ 2 and a timelike boundary T . The intersection of Σ leafs with T are co-dimension two surfaces S from S 1 to S 2 . Now in this region of space-time we calculate the variation of the action (6), and impose the equations of motions , i.e. G ab = 0.
Then we may find: Here, the symbol , h ab is the metric on each spacelike boundary and γ ab is the induced metric on the timelike boundary T . P ab and Π ab are respectively the gravitational momenta of Σ and T : where K ab and χ ab are their corresponding extrinsic curvatures. In the original work of Brown and York [3], they assumed the boundaries to be orthogonal, thus, the contribution of the joints was missing. However, after the work of Hayward [13], who emphasized the importance of this term for the variational principle to be well-posed, this term appeared in later works [11,14,15].
If we ignore these joints and consider the boundaries to be orthogonal, then in an analogous way to the mechanical example one may define the gravitational canonical momentum as the derivative of the principal functions with respect to the induced metric on spacelike segments of the boundary and gravitational energy-momentum-stress tensor as the derivative with respect to the induced metric on timelike segment. In fact, the Π ab has the same expression as for ADM canonical momentum. Thus, the energy-momentum-stress tensor (or just the stress tensor for abbreviation) will be defined as: To define conserved quantities, now we choose a spacelike co-dimension two surface S in T with unit timelike normal u a . Then, the metric γ ab is further decomposed as γ ab = q ab −u a u b .
Here, u a defines local the flow of time in T . Moreover, according to Brown-York [3] for an isometry in the boundary generated by the Killing vector ξ a , the conserved charge associated to this symmetry is defined by: In fact, there is some points regarding to the expression (9) and (10). The first one is that, as we pointed out before, these expressions are not unique. One can append a subtraction term S 0 to the boundary action without affecting variational problem when S 0 depends on fixed boundary data, S 0 = S 0 (h ab , γ ab ), which leads to ambiguities in the definitions of energy and momenta. According to Brown and York interpretation these ambiguities are consequence of freedom to choose the zero point of energy and redefine system momenta with canonical transformation [3,11]. On the other hand, Eq. (10) leads to infinities when calculated for large spheres in general systems in spacetime. The Brown-York proposal is to choose subtraction term S 0 such that the modified action S − S 0 leads to zero energy for flat spacetime. Therefore, the zero point is chosen to be the flat spacetime while the scheme for other geometries is to embed their boundary in flat spacetime. Thus, the modified expression for the stress tensor will be As another remark, if we want to consider the effect of joints or non-orthogonal boundaries, we must care about the components of δγ ab contained in δq ab appeared in the last term of (7).
Also, we must be careful about the question that with respect to which observer we desire to calculate the quasi-local quantities. Note that when the boundaries are non-orthogonal, the Eulerian observers orthogonal to Σ constant are different from the observers orthogonal to S constant in the boundary T . In these cases, a further decomposition of the induced metric γ ab with the assistance of the vector u a = N ∇ a t is required, where t is a foliation of T . By this decomposition we get [11]: where N and V a are lapse and shifts of decomposition. The components of δS δγ ab can be calculated as which are known respectively as the quasilocal energy density, tangential momentum density and spatial stress. The details of calculations for different observers can be found in [11,15].
In these cases it has been shown that a double foliation of space-time is a natural setup for calculations as described in appendices of [11]. There, the relation between quantities, as measured by different observers, has been obtained. The relations between two sets of normals to S as depicted in Fig.2 is:n The quasilocal energy density associated with the two surfaces S as seen by the observers orthogonal to Σ constant is [11]: where k ab = −q c a q d b ∇ cnd and l ab = −q c a q d b ∇ cūd are defined so that k and l are their trace, respectively. In the case for which two observers are at rest with respect to each other i.e. v = 0 or, in other words, the boundaries are orthogonal, the quasilocal energy density becomes: where k 0 is the corresponding extrinsic curvature as embedded in flat space. So, the total quasilocal energy on S becomes [3,11]: Evaluating the above integral provides an expression for E as a function of r, E = E(r).
However, we must note that because of the Eq. (17), in general the form of this function varies for different observers. For example, in the case of Schwarzschild black hole the function for static, radially infalling and boosted observers has been obtained in [15].
As the last point, let us note that for the above relations , and for different observers, it is always assumed that these observers do not move in null geodesics. In fact, null observers can not be contained in the above setup because the normal vectors (or four-velocities) are always normalized to unity. As the name indicates, due to the fact that the norm vanishes for null observers, normalization does not make sense in this case. In fact, here, the normal vectors or four-velocities are both normal and tangent to the hypersurface, the induced metric becomes degenerate, and the usual extrinsic curvature does not make sense.
In the next section we will use standard treatment for null hyper-surfaces. The resulting stress tensor can provide the quasi-local quantities as measured by null observers.

Asymptotic AdS spacetime and counterterm method
Finding proper subtraction term by embedding in reference spacetime is a difficult task. In fact, it is not possible to embed a boundary with an arbitrary intrinsic metric in the reference spacetime. In the case of asymptotically AdS spacetimes, there is an attractive proposal without necessity of embedding in reference spacetime as proposed in [16]. This approach is inspired from AdS/CFT duality [17,18] by interpreting divergences which appear in stress tensor when the boundary is moved to infinity as dual to standard ultraviolet divergences of quantum field theory ; then, they can be removed by adding local counterterms to the boundary action. For instance for AdS 4 the counterterm Lagrangian in the boundary proposed to be: where R is the scalar curvature of induced metric on the boundary, and ℓ is the AdS radius.
Adding this term to the usual GHY term in the timelike boundary, then variation of action yields to a regularized stress tensor such that energy becomes finite at the r → ∞ limit. The main advantage of this method, beside the fact that counterterms are covariant and do not spoil the variational principle, is the lack of embedding difficulties in BY method. However, these counterterms are known just for asymptotically AdS spaces. For asymptotic flat space, since there is no length scale ℓ, finding such counterterms is problematic. Also, the limit ℓ → ∞ in Eq. (20) does not lead to a unique covariant expression [19].

Hamilton-Jacobi analysis on null boundary
In this section we are going to repeat the calculation of previous section when instead of a timelike boundary we have a null one. The corresponding spacetime region is illustrated in Fig.3. In order to calculate the variation of the action in this region we need to express

The set up
A segment N of the spacetime boundary, characterized by φ 0 = 0, is called a null hypersurface if ∇ a φ 0 ∇ a φ 0 = 0. This feature of null boundary indicates that the normal vector to the null surface is also tangent to it. This property is the origin of some difficulties when dealing with such hypersurfaces; because, as a consequence, the induced metric becomes degenerate and therefore constructing a projector to the null surface just from its normal is not possible.
One standard remedy to this problem is to introduce an auxiliary null vector k a which lays out of the hypersurface and therefore ℓ a k a = 0, when ℓ a is the null normal to the boundary, e.g., ℓ a ∝ ∇ a φ 0 on the boundary. For more details about the geometry of null hypersurfaces we refer the interested reader to [20,21].
By defining ℓ a as the normal to the null boundary, we introduce the auxiliary null form k a and take the normalization of these null forms to be everywhere as: ℓ a ℓ a = 0 , k a k a = 0 and ℓ a k a = −1 .
With the aid of ℓ a and k a , we can define the projector as: This projector is not in fact a projector on null surface; instead, it essentially projects spacetime vectors onto the co-dimension two surface S, to which ℓ a and k a are orthogonal.
A systematic approach to define these co-dimension two surfaces is to use a doublefoliation by two scalar fields (φ 0 , φ 1 ). The intersection of level surfaces of φ 0 and φ 1 are the co-dimension two surfaces S. In this foliation ℓ a and k a can be expanded generally as: Three out of four coefficients in the above expansions are determined by normalization conditions (21) and one remains free, due to the rescaling or boost gauge freedom (ℓ a → αℓ a , k a → 1 α κ a ). On the boundary, we set the coefficient B = 0, so that the boundary is a level surface of φ 0 and we have ℓ a B = A ∇ a φ 0 , and it is null because ∇ a φ 0 ∇ a φ 0 ∝ ℓ a ℓ a = 0. Note that ℓ a and k a are form fields defined in whole space-time and everywhere we have ℓ 2 = k 2 = 0, whereas ∇ a φ 0 ∇ a φ 0 = 0 is satisfied just on the boundary.
In this double foliation framework, the spacetime metric becomes: By covariant differentiation of vectors ℓ a and k a and projecting them in different directions, using q a b , ℓ a and k a , we can define the following geometric objects from ∇ a ℓ b and ∇ a ℓ b : These relations are generalizations of the relation ∇ a n b = −K ab + n a a b to current case where decomposition has been done with two null vectors. The definitions are as follows: where Θ ab and Ξ ab are extrinsic curvatures of S while ω a , η a andη a are twists. In addition, a a andā a are tangent accelerations of ℓ a and k a to S, respectively. Moreover, κ andκ are in-affinity parameters 4 .
As a side remark, let us point out that this general double foliation described above is different from double null foliation as considered by various authors [22,23] . Although, the double null foliation is useful for initial value problem as shown by Sachs [24], one could show that for variational principle it suffers from having a partial gauge fixing condition which is a disadvantage [25]. This is because requiring that the level sets of a coordinate φ being null, fixes one of the metric components: The double null folation is a special case of the above set up if we set B = C = 0 in whole spacetime. If fact, we will see in the next section that the variation of such components is important for finding complete canonical structure, and more importantly the canonical momenta of such variations is the quasilocal energy density of system.

Variation of Hilbert-Einstein action
The main point in calculating the variations on a null surface is that the operators δ for variations and the covariant derivative ∇ are defined in spacetime, whereas the condition for boundary, in the above setup B = 0, is valid only on the null boundary. Thus, we must first apply these operators and then impose the condition B = 0. For example consider variation of the vector ℓ a given by Eq. (23). By variation of this vector on the boundary, one finds δℓ a = δA∇ a φ 0 + δB∇ a φ 1 . On the other hand, if we set B to zero first and then vary the 4 The quantity κ is called surface gravity in the case of a black hole horizon null surface. equation, we find that δℓ a = δA∇ a φ 0 . This varaition is not valid, because variation of one degrees of freedom in metric has been killed in this relation, i.e. δB ∝ ℓ a ℓ b δg ab = 0. This variation is responsible for taking out the boundary from being null. To avoid losing any degree of freedom, in the following we will not set B = 0 until the end of calculations.
We consider a hypersurface to be a leaf of one of the foliations, let it be φ 0 = const.
Variation of HE action on such hypersurface in the double foliation determined by eqs. (23) to (25) leads to [1]: where H is the determinant of H ij defined above. The first line as a total variation suggests appropriate boundary term that must be subtracted in order to have a well-posed variational principle. The total variation term in the boundary integral, i.e.
has a nice geometric interpretation. It can be rewritten as Bk a ) is the normal to the hyperfurface, and K a bc = −(k a ∇ b ℓ c + l a ∇ b k c ) is defined such that K a is its trace on two last indices 5 .
On a null boundary now we set B = 0; thus, we find the action with a well-posed variational principle to be Let us remind that, as mentioned in previous subsection, there is a scaling boost gauge symmetry in null hypersurface description by two vectors ℓ a and k a . We can use this gauge freedom to set D = 1 in this gauge and by putting B = 0 we have √ H = A. Therefore, the action becomes: This is what found in [8,9]. In [7] the authors also have set the lapse A = 1, so the last term vanishes. Using the above action, the variation of the principal function on the null boundary becomes: Similar expression for canonical structure on the null boundary has been found in [9]. There, the authors considered only the variations that keeps the boundary to remain null, thus, the last term of the first line was missing in their analysis. Note that in this expression δq ab , δβ 1a and (−2δB) are all variations of the metric components tangential to the boundary as easily can be seen from the decomposition (25). This is an important point because, as in the non-null case, the variational principle tells us that for Dirichlet boundary condition we only need to fix the tangential metric components. Here also the number of degrees of freedom matches with the non-null case; for example, in four dimensions δh ab has six components while (δq ab , δβ 1a , −2B) altogether have 3 + 2 + 1 = 6 components.

The stress tensor
Having found the variation of action on the null boundary, we can use BY prescription to find the stress tensor on the boundary. In doing so we must differentiate the principal function in space-time region, illustrated in Fig 3, with respect to the metric component tangential to the null boundary segments N . Note that differentiating with respect to metric components in Σ 1 and Σ 2 gives the canonical momenta P ab similar to the previous section. According to the Eq.(25), the metric components tangential to N are (δq ab , δβ 1a , −2B). Note that, in contrast to previous timelike boundary presented in the last section, there is no induced metric on null surface; thus, we must differentiate with respect to each component of metric separately according to: ).
On the other hand, the variational principle does not fix the boundary action completely; therefore, we can subtract any functional S 0 of fixed boundary data from the action (33).
where ω a 0 and Ξ 0 are defined as: According to BY prescription, eqs. (36) to (38), may be used to define quasilocal energy density ǫ, tangential momentum density j a and spatial stress s ab , respectively.
The subtraction term S 0 is related to the zero point of energy and must be well chosen so that the stress tensor leads to finite energy for different systems and conventionally zero energy for Minkowski space time. The expression (36) is similar to the Eq. (18) for timelike boundary. Ξ is defined by Ξ = −q ab ∇ a k b while k reads as k = −q ab ∇ a n b . The similarity is in the sense that for null boundary the only vector pointing out of the boundary is k a . For nonnull boundaries, studied in the previous section, the normal to boundary was also pointing out of it; however, note that for null hypersurfaces the normal is within the boundary.

Quasi-local quantities
In this section we examine the proposed stress tensor by evaluating conserved charges of various space-times and then we compare it with the known results. Here, as to the timelike boundary, we propose the Quasi-local quantity to be: where ℓ a has the role of time flow on the boundary while ξ a is a Killing vector which generates an isometry of the boundary. If ℓ a is the generator of time translation symmetry, then the total energy becomes: where the quasi-local energy density ǫ is defined by (36). When there is a rotational Killing vector ζ a , then its corresponding angular momentum is: In the following we will calculate the above quantities for some well-known gravitational solutions.

Minkowski space
The simplest example for investigation is Minkowski spacetime. In retarded-spherical coordinates the metric is as follows: where dΩ 2 = γ AB dσ A dσ B is the metric on a unit sphere. It is evident that u = const. is a null surface. Comparing the above line element with Eq.(25) and using the foliation relations (23) and (24) yields ℓ a = ∇ a u, k a = ∇ a r + 1 2 ∇ a u and q AB = r 2 γ AB . From the definitions (29) one easily finds: As a result Ξ = − 1 r , and the integral dθdφ √ q Ξ = − dθdφ r sin 2 θ will be infinite as r → ∞. Subtracting Ξ 0 leads to zero energy for flat space as expected.

Schwarzschild black hole
Our next simple example is Schwarzschild black hole. In retarded Eddington-Finkelstein coordinates the metric is given by: In this case we have ℓ a = ∇ a u and k a = 1 2f (r) ∇ a u + ∇ a r. Here, by calculating the quasilocal energy density, we get: By replacing f (r) = 1 − 2M r , we obtain ǫ = M 4πr 2 . Thus the total energy becomes: It is interesting that the quasilocal energy, as calculated for the null observers in N , is independent of the distance r for which the energy is calculated. At first sight, it may seems strange because the usual BY expression for energy is just equal to the ADM mass at r → ∞. However, notice that as mentioned previously, in general the function E(r) is observer dependent. The in-dependency of energy to the distance also has been observed for boosted foliation of Schwarzschild black hole in [11].

Slow-rotating black hole
Here we examine conserved quantities for slow rotating Kerr black holes. Again we write the metric in the retarded Eddington-Finkelstein coordinates: Comparing with metric (25) we arrive at: Up to the first power in J, one finds the same expression for energy as the Schwarzschild case. The angular momentum quantity is related to ω a , using its definition we find: where ξ = ∂ φ is rotational Killing symmetry. Therefore, the total angular momentum becomes:

Asymptotic flat spacetime and the Bondi mass
Here, we want to study the quasi-local quantities for gravitational theories in which the metric has an asymptotic flat space behaviour. In order to do that, we suppose the metric in the Bondi coordinates. In this gauge, the most general four-dimensional metric takes the form: where ∂ r det( q AB r 2 ) = 0. By comparison the above equation with Eq.(25), we find: The expressions for Θ AB and Ξ AB in the metric (25) have been calculated explicitly in [1] and are: where D is covariant derivative on two sphere S, compatible with the metric q AB 6 . Using standard asymptotic expansions [26]: we can easily find the following leading terms for Ξ: By embedding in flat space we have the reference term Ξ 0 = − 1 r . The term D A W A is a total derivative on compact two sphere S which vanishws by integration. Thus, finally the total energy becomes: which is the expression known as Bondi mass.

A possible counterterm for asymptotic flat spacetime
We have seen that for asymptotic AdS space there is counterterm method for regularizing quasi-local quantities that has some advantages. The first benefit is that the dependency on the reference spacetime and mathematical difficulty of embedding is avoided. Furthermore, the counterterms have direct interpretation in the dual field theory and has important role for building the dictionary of AdS/CFT. Then the natural question is that: why there is no similar counterterms for flat space? For flat space there is not a length scale counterpart to the AdS length scale; thus, we can not have terms like: ℓ 2n−1 R n .
In addition, the counterterms should not spoil the variational principle so that the extrinsic curvature terms are forbidden. Therefore, on the time-like boundary we don't have a viable candidate as counterterm in asymptotic flat spacetime. However, as we have seen in this paper, null boundaries are special. Consider the following term on the null boundary: with arbitrary numerical coefficient α. The first point is that adding such term is compatible with Dirichlet boundary conditions. Also, we must note that although the above term vanishes on the null boundary, its variation is none zero. Therefore, it has a contribution to the stress tensor and energy though it is zero for on-shell action. Fortunately, we can set the coefficient α = − 1 2 so that its contribution to energy make the total finite. By adding this term to the boundary action, the quasilocal energy density becomes: For asymptotic flat spacetime studied previously, using (53), Θ is easily computed and its leading terms is as: Therefore the expression (57) leads to correct total energy without needing to embedding and reference spacetime.

Conclusion and outlook
In this work we have extended the Brown-York prescription to the case of null boundaries.
We have seen that the variation of one metric components is responsible for taking away the boundary from being null. In fact, it was shown that an exact derivation with this variation gives the expression for quasilocal energy density. In addition, the expressions for total energy and angular momentum were examined for some known spacetimes. Moreover, a special property of the calculated energy for null observers was found to be the in-dependency to radial distance. Furthermore, it was shown that in the case of null boundaries, there is a possible counterterm which is consistent with variational principle and can be added to boundary action so that the total energy becomes finite, without necessity of embedding in reference spacetime.
Regarding the similarity with the standard AdS/CFT dictionary, one may interpret the introduced stress tensor as expectation value of stress tensor in dual field theory to flat space.
Because, according to Penrose diagrams for asymptotic flat spacetimes, the null hypersurfaces I + and I − are regarded as the boundaries of spacetime; thus, one application of the stress tensor found in this paper may be flat holography.
Another application of the formalism presented in this article and [1] is to revisit the gravity in light-front. The usual investigation of gravity in light-front uses the double-null foliation of spacetime [27,28]. However, as noted above, the double-null foliation leads to partial gauge fixing of the metric. It is in contrast to field theory formulation in the lightfront coordinate which no gauge fixing is required. As we have seen a general double foliation can preserve all degrees of freedom. Thus, it is motivating to revisit gravity in light-front.