Perturbations of the Gravitational Energy in the TEGR: Quasinormal Modes of the Schwarzschild Black Hole

We calculate the gravitational energy spectrum of the perturbations of a Schwarzschild black hole described by quasinormal modes, in the framework of the teleparallel equivalent of general relativity (TEGR). We obtain a general formula for the gravitational energy enclosed by a large surface of constant radius $r$, in the region $m\,<<\,r\,<<\infty$, where $m$ is the mass of the black hole. Considering the usual asymptotic expression for the perturbed metric components, we arrive at finite values for the energy spectrum. The perturbed energy depends on the two integers $n$ and $l$ that describe the quasinormal modes. In this sense, the energy perturbations are discretised. We also obtain a simple expression for the decrease of the flux of gravitational radiation of the perturbations.


Introduction
The response of a black hole or neutron star to external, nonradial perturbations is described by quasinormal modes. A comprehensive review of the physics related to this interesting phenomena is found in refs. [1,2]. These modes are damped oscillations of the space-time geometry that may be used to characterize the intrinsic properties of the physical system. Investigations on quasinormal modes are carried out both analytically and numerically. The modes are characterized by a spectrum of discrete, complex valued frequencies. The real part of the frequency is related to the oscillation frequency, and the imaginary part yields the rate at which each mode is damped as a result of emission of radiation. Thus, quasinormal modes may be important to gravitational waves astrophysics. In addition, the study of these modes is, to some extent, a testing ground for ideas in quantum gravity. Chandrasekhar once stated [3] that one relevant way of investigating a physical system is by perturbing it, and analysing the response of the system. This is precisely the role of QNM in the context of black holes.
In the analysis of Einstein's equations for a perturbed Schwarzschild space-time one finds a number of equations that are eventually reduced to two one-dimensional wave equations for the perturbed metric components (for axial and polar perturbations) [1]. However, the nature of the potential in these equations precludes exact solutions in terms of known functions [4]. In particular, the equations for the metric perturbations admit solutions provided the frequencies are discrete and under the imposition of special boundary conditions. The solutions of these equations are the quasinormal modes.
In this article, we calculate the energy spectrum of the quasinormal modes in the framework of the teleparallel equivalent of general relativity (TEGR). A proper definition for the energy-momentum of the gravitational field is of ultimate importance for a comprehensive understanding of Einstein's general relativity, and yet there is no general agreement regarding an acceptable expression. There is not an unique approach to the problem (pseudotensors, quasilocal expressions, for instance), and even within each approach there is no preferred definition for the energy-momentum. In the context of the field equations of the TEGR, however, one finds a suitable and consistent framework for the definitions of the gravitational energy-momentum and 4angular momentum [5,6], that satisfy the algebra of the Poincaré group. This definition of the gravitational energy-momentum has been applied to several configurations of the gravitational field, and leads to consistent results. It satisfies important requirements that any gravitational energy-momentum definition must satisfy [6].
We will consider the perturbed Schwarzschild space-time and calculate the gravitational energy enclosed by a surface of constant radius r. The energy contained within this region of the unperturbed space-time is known and may be easily evaluated out of the definition that arises in the TEGR, as well as by means of several quasilocal definitions for the gravitational energy. Thus, by subtracting the unperturbed energy from the total energy that includes the perturbations, we obtain the energy of the perturbations only. In the present analysis, we will consider both axial and polar perturbations. We are restricting the analysis to gravitational perturbations (we do not address scalar and electromagnetic perturbations). We find that the perturbed energy (i) oscillates and at the same time is damped, and (ii) depends on the integers n (the overtone number) and l (the angular momentum number). To our knowledge, a similar analysis regarding the spectrum of the energy perturbations of the Schwarzschild space-time has not been presented so far. We conjecture that the perturbed energy in arbitrary (nonspherical) volumes of the Schwarzschild space-time is also characterized by integers.

Review of the gravitational energy-momentum definition in the TEGR
We assume that the space-time geometry is defined by the tetrad field e a µ only. In this case, the only possible nontrivial definition for the torsion tensor is given by T aµν = ∂ µ e aν − ∂ ν e aµ . In the TEGR, it is possible to rewrite Einstein's equations in terms of e a µ and T aµν . The Lagrangian density of the theory is defined by L M stands for the Lagrangian density for the matter fields. The Lagrangian density L is invariant under the global SO(3,1) group. Invariance under the local SO(3,1) group is verified as long as we take into account the total divergence that arises in the identity where R(e) is the scalar Riemannian curvature. However, the field equations derived from eq. (1) are invariant under local SO(3,1) transformations, and are equivalent to Einstein's equations. They read where δL M /δe aµ = eT aµ . The definition of the gravitational energy-momentum may be established in the framework of the Lagrangian formulation defined by (1), according to the procedure of ref. [5] (we will make c = 1 = G). Equation (3) may be rewritten as where T λµ = e a λ T aµ and t λµ is defined by In view of the antisymmetry property Σ aµν = −Σ aνµ , it follows that ∂ λ e e a µ (t λµ + T λµ ) = 0 .
The equation above yields the continuity (or balance) equation, Therefore we identify t λµ as the gravitational energy-momentum tensor [5], as the total energy-momentum contained within a volume V , as the gravitational energy-momentum flux [7,8], and as the energy-momentum flux of matter [8]. In view of (4), eq. (8) may be written as where Π aj = −4ke Σ a0j . A summary of all issues discussed above may found in ref. [9] Equation (11) is the definition for the gravitational energy-momentum presented in ref. [6], obtained in the framework of the vacuum field equations in Hamiltonian form. It is invariant under coordinate transformations of the three-dimensional space and under time reparametrizations. Note that (6) is a true energy-momentum conservation equation. We also remark that for finite volumes of integration, the free index of the integral of the energymomentum 4-vector density must be a Lorentz index, such as a in the left hand side of definition (11), because under a global SO(3,1) transformation both sides of (11) transform consistently. In contrast, the integral on the right hand side of a hypothetical quantity such as P µ = d 3 x eV µ (with spacetime index µ) would be coordinate dependent (the integral of a space-time vector density eV µ is neither invariant nor covariant under coordinate transformations), and therefore ill-defined for finite, three-dimensional volumes of integration.
In the ordinary formulation of arbitrary field theories, energy, momentum, angular momentum and the centre of mass moment are frame dependent field quantities, that transform under the global SO(3,1) group. In particular, energy transforms as the zero component of the energy-momentum four-vector. These features of special relativity must also hold in general relativity, since the latter yields the former in the limit of weak (or vanishing) gravitational fields. As an example, consider the total energy of a black hole, represented by the mass parameter m. As seen by a distant observer, the total energy of a static Schwarzschild black hole is given by E = mc 2 . However, at great distances the black hole may be considered as a particle of mass m, and if it moves with constant velocity v, then its total energy as seen by the same distant observer is E = γmc 2 , where γ = (1−v 2 /c 2 ) −1/2 . Likewise, the gravitational momentum, angular momentum and the centre of mass moment are naturally frame dependent field quantities, whose values vary from frame to frame, and are different for different observers in arbitrary space-times.
Before closing this section, we note that the existence of a scalar density (on a three-dimensional spacelike hypersurface) like ∂ j Π aj = −4k∂ j (e Σ a0j ) in eq. (11) is natural in the framework of theories constructed out of third order tensors such as the torsion tensor T aµν . Naively, we observe that the contraction of two space-time indices in a third order tensor yields a vector of the type φ µ , and therefore ∂ µ (eφ µ ) is a well defined scalar density that under integration yields a well behaved surface integral, in similarity to eq. (11). The existence of these well defined scalar densities is natural in the TEGR, as in the identity below eq. (2), but of course, after a number of manipulations and rearrangements, these scalar densities could also be established in Einstein-Cartan type theories.

3
Axial perturbations of the Schwarzschild black hole The analysis of the quasinormal modes in the Schwarzschild space-time consists in solving the differential equations for the perturbed metric components, with appropriate boundary conditions. These equations were first written down by Regge and Wheeler [10]. They obtained the general form of the simplest nonspherical perturbations for the Schwarzschild black hole, namely, the axial and polar perturbations. Regge and Wheeler showed that the equations that describe the axial perturbations may be separated if the perturbed metric tensor h µν is expanded in tensorial spherical harmonics.
The general form of h µν is obtained and further simplified by taking into account the gauge symmetry (coordinate invariance) of the field equations. The simplest form of the nonvanishing axial perturbations are given by [10] where P l (θ) are the Legendre polynomials. The functions h 0 (t, r) and h 1 (t, r) satisfy the equations [1, 10] where f (r) = 1 − 2m/r. Equation (15) is a consequence of (13) and (14). The quantity ∂h 0 /∂t in (13) may be substituted in (14). Thus, (13) and (14) may be combined into one equation. Defining Ψ(t, r) = (1/r)f (r)h 1 (t, r), the resulting equation reads In terms of the tortoise coordinate x = r + 2m ln(r/2m − 1), eq. (16) may be rewriten as V RW (x) is the Regge-Wheeler potential. It is given by The time dependence of h µν is assumed to be of the type e −iωt , where ω is a complex frequency. Thus, it follows that Ψ(t, r) = e −iωt ψ(r). As a consequence, the function ψ(r) satisfies the time independent equation In view of the asymptotic behaviour of the potential V RW (x), the expressions of ψ(x) in the asymptotic limits x → ∞ (r → ∞) and x → −∞ (r → 2m) are given by ψ(x) ∼ = Ae iωx and ψ(x) ∼ = Ae −iωx , respectively, where A << 1 in order to ensure the perturbative character of the solution. Expressions with positive exponential correspond to waves that are purely outgoing at infinity, and the ones with negative exponential, to waves that are purely ingoing at the horizon. For weakly excited states, the quasi-normal frequencies ω may be obtained from the third order WKB approximation. In the case of the gravitational perturbation, the fundamental state is given by n = 0 and l = 2 [11,12]. These semi-analytic approximations yield values very close to the numerical values [13]. The stable solutions decay with time. Thus, if ω r and ω i represent the real and imaginary parts, respectively, of the complex valued frequency, then we must have ω i < 0.
The possible values of ω for highly damped modes (large values of n) are independent of l, and may be obtained by means of numerical investigations. They are given by [1] where n is an integer. This expression has also been obtained by means of analytic procedure in ref. [14]. However, the frequency to be employed in the construction of the Figures below, in subsections 3.1 and 4.1, does not correspond to highly damped modes. The problems regarding the divergent behaviour of h 1 (t, r) in the limit r → ∞ have been discussed in refs. [1] (section 3.1.2) and [15]. However, the asymptotic behaviour of h 1 (t, r) poses no problem to the result of our analysis, since we will integrate (11) over a surface of constant (finite) radius r << ∞.
We need the asymptotic behaviour of h 0 (t, r). Given that ∂ 0 h 0 = iω h 0 , and which follows from (13) in the limit r → ∞, we find The m/r correction to the expressions above will be negligible in the development of our analysis. By means of straightforward calculations, we obtain from (22) an expression that holds in the limit r → ∞ and will be useful in the following subsection. We find

The gravitational energy of the axial perturbations
In this subsection, we will address the axially perturbed Schwarzschild spacetime and calculate the gravitational energy contained within a simple three dimensional volume, namely, a spherical volume of radius r. In order to simplify the calculations, we will consider a large sphere so that the field quantities are considered in the limit m/r << 1. The purpose is to find the expression of the energy perturbations in this limit. To our knowledge, no expression for the gravitational energy perturbations has been obtained so far in the present physical context. In the course of the calculations we find that the perturbed gravitational energy depends on the square of the field quantities h 0 and h 1 that appear in (12). Therefore in order to calculate the first order contribution to the perturbed energy we must keep the terms quadratic in h 03 and h 13 .
We start with the metric tensor [10] Definition (11) for the gravitational energy-momentum is frame dependent. Thus, we choose a configuration of tetrad fields that has a clear physical interpretation. Tetrad fields are interpreted as reference frames adapted to preferred fields of observers in space-time. This interpretation is possible by identifying the e (0) µ components of the frame with the four-velocities u µ of the observers, e (0) µ = u µ [16,17]. In the present analysis we will establish a set of tetrad fields adapted to static observers in space-time. Therefore, we require e (0) i = 0. This condition fixes 3 components of the frame. The other three components are fixed by choosing a orientation of the frame in the three-dimensional space. Thus, e (0) µ is parallel to the worldline of the observers, and e (k) µ are the three unit vectors orthogonal to the timelike direction. We fix e (k) µ such that e (1) µ , e (2) µ and e (3) µ in cartesian coordinates are unit vectors along the x, y and z directions. The tetrad field in (t, r, θ, φ) coordinates that satisfies these conditions is given by with the following definitions: g 03 (−g 00 ) 1/2 , E 2 = 1 + g 11 g 2 03 + g 00 g 2 13 r 2 sin 2 θ , The frame above satisfies e (0) i = 0. It is possible to show that if we neglect g 03 and g 13 , the frame components e (k) µ (t, x, y, z) in the limit m/r << 1 are given by e (1) µ ∼ = (0, 1, 0, 0) , e (2) µ ∼ = (0, 0, 1, 0) , and e (3) µ ∼ = (0, 0, 0, 1). We denote by g and e the determinants of g µν and e a µ , respectively. We find g = −e 2 = −r 4 sin 2 θ E 2 . Thus, we have The gravitational energy is given by the a = (0) component of (11). Transforming the volume integral into a surface integral and considering the definition Π aj = −4k Σ a0j , we find that The integration will be carried out on a surface S of constant radius r.
The non-vanishing components of the torsion tensor that are relevant to the evaluation of the gravitational energy are The evaluation of Σ (0)01 requires a large number of algebraic manipulations, but otherwise is simple. The full expression of eΣ (0)01 is given by By keeping terms of order (h 0 ) 2 and (h 1 ) 2 only in the integrand of (30), we obtain the approximate expression that is relevant to the present analysis. The leading term in the energy expression is m, the total Schwarzschild energy. It is easy to verify that which is a well known result. We recall that B is given by (26). The gravitational energy of the perturbations will be denoted by δP (0) , and is obtained out of the second order terms in (30), assuming −g 00 ∼ = 1 and g 11 ∼ = 1. Therefore, for fixed values of r such that m/r << 1 we obtain, after some simplifications, Recalling that k = 1/16π (we are assuming c = 1 = G), we find Now we make use of (22) and (23), and we arrive at where we have used π 0 dθ sin θ∂ θ (P l (cos θ)) 2 = 2l(l + 1) 2l + 1 .
In the absence of the perturbation, i.e., when h 0 = 0 = h 1 , δP (0) vanishes. The approximate solution (34) may be compared with the numerical solution presented in Figure 1. The numerical solution, represented by the continuous line, is obtained directly from the exact expression for the total gravitational energy minus the black hole energy, requiring the function ψ to be a solution of eq. (19). The boundary conditions for the numerical solution of eq. (19) are ψ(r 1 ) = A exp −iωx 1 and ψ(r 2 ) = A exp +iωx 2 , where x 1 = r 1 + 2m ln(r 1 /2m − 1), x 2 = r 2 + 2m ln(r 2 /2m − 1), r 1 = 2m + m/4 and r 2 = 200m. In Figure 1 we display the initial instant of time of the perturbation. We see that the larger is the value of the coordinate r, more intense is the energy of the perturbation, probably because the farther one is from the black hole, the more energy is necessary to perturb it. Note that r is not a dynamical coordinate, it is just a parameter of the coordinate system that fixes the radius of integration. The evolution of the gravitational energy perturbation is stable, i.e., the latter tends to zero when t → ∞, as we see in Figure 2.  Considering the real part of the solution given by eq. (34), we have We may calculate the variation of the perturbation energy in a halflife interval T /2, for weakly damped solutions, but neglecting the damping effects, considering only ω = 2ω r . By making use of we obtain Since the black hole looses energy in the course of time, as we see in Figure  (2), it is important to analyse the rate at which the energy decay occurs. This issue may be described by the energy flux of the radiated energy. The present formalism for the gravitational energy-momentum allows a definition for the gravitational energy-momentum flux [8], which has been applied to the radiation of gravitational energy in the Bondi space-time. The gravitational energy flux is given by the a = (0) component of Φ a g defined by (9), or simply by In the present context, we find Again we define ω = 2ω r and find that after a period T = 2π/ω the decrease in the flux of gravitational radiation is Φ (0)

Polar perturbations of the Schwarzschild black hole
In this section, we repeat the analysis developed in the previous section for the axial perturbations. In the context of polar perturbations, the calculations are more intricate, in spite of the similarities with the calculations carried out for the axial perturbations. We will consider the perturbed metric tensor as given by g µν =ḡ µν + h µν . The line element for the unperturbed metricḡ µν tensor is and the polar perturbations h µν are described by [10] h µν = e −iωt P l (cos θ) where f = 1 − 2m/r and H 0 , H 1 , H 2 , K are functions of t, r. Therefore, the perturbed metric tensor is written as where we define H 0 = H 0 e −iωt P l (cos θ) , Einstein's equations reduce to seven non-trivial equations: one algebraic equation, three first order differential equations for the metric tensor components, and three second order equations. The algebraic equation is eventually reduced to [10] The three first order differential equations are [18] dK dr The three remaining equations eventually yield the algebraic equation where λ ≡ 1 2 (l − 1)(l + 2). By rewriting the field quantities according to [18] and and defining the tortoise coordinate x as we obtain the relations α(r) = λ(λ + 1)r 2 + 3λmr + 6m 2 r 2 (λr + 3m) , It follows that Θ = dψ dx (56) and where is the Zerilli potential.
Although the algebraic forms of the Zerilli and Regge-Wheeler potentials are quite distinct, they share strong similarities [1]. For very large values of the coordinate r, we have V Ze → 0. Thus, in the limit r >> m we obtain a simple solution for the function ψ, given by where A is a constant with dimension of length. In order to ensure the perturbative character of the solution, the amplitude A is required to satisfy A << 1, in natural units. In the limit r >> m, the solution ψ grows with the radial distance. Thus, since the approximations are made assuming |ψ 2 | << |ψ|, we must ensure that |A 2 e 2iωx | << |Ae iωx |. The frequencies of the polar perturbations are the same as for the axial perturbations [19], provided we consider the same values of (n, l) when comparing the frequencies.
Considering as usual ω = ω r + iω i , we have It is important to remark that the stability of the solutions demand that the function Ψ = e −iωt ψ decay sufficiently fast with time. Consequently, the stable solutions are those for which ω i < 0. The solutions obtained in the following approximations are valid as long as Therefore, the approximate analytic solutions will be valid provided the condition −∞ << x << 1 |ω i | ln A is satisfied. Just like in the axial case, we must find the asymptotic limits for the perturbative functions. Thus, in the limit r >> m, equations (52,53,54,55) become As a consequence, equations (49,50) are approximated by and Making use of eq. (48) in the limit r >> m, together with the equations (65,66), we obtain H ≈ −ω 2 rψ .
The expressions above given by eqs. (65,66,67) will be used in the following subsection, where we will evaluate the gravitational energy of the polar perturbations.

The gravitational energy of the polar perturbations
In similarity to the analysis in subsection 3.1, where we addressed the gravitational energy of the axial perturbations, here we also construct a set of tetrad fields adapted to static observers (with respect to the asymptotic flat space-time limit of the black hole), that yields the metric tensor (43). A set of tetrad fields that satisfy the necessary conditions is given by where and whose determinant is e = D 2 ACr 2 sin θ.
With the help of these quantities, we find The expression above leads to the gravitational energy given by eq. (11). An approximate analytic expression for the quantity above may be obtained by considering terms up to the second order in H 0 , H 1 and K. The unperturbed expression of the momentum given by eq. (72) is simply Π (0)1 g = 4k sin θr(f 1/2 − 1). The integration of this quantity on a surface at spacelike infinity yields the expected expression for the total gravitational energy, E = m. Thus, the quantity that yields the gravitational energy of the perturbations is δΠ (0)1 = Π (0)1 − Π (0)1 g . The approximate analytic expression of the latter is The first line in (73) contains the first order terms δΠ (0)1 f irst , whereas the other terms correspond to the second order terms δΠ (0)1 sec . The total energy within a spherical surface of radius r is obtained by summing these two quantities, δΠ (0)1 = δΠ (0)1 f irst +δΠ (0)1 second and integrating over the spherical surface. The first order terms depend on the variable θ in the form P l (cos θ). The integration of this quantity vanishes, i.e., π 0 d(cosθ)P l (cos θ) = 0. Therefore, only the second order terms yield a non-trivial result when integrated over the whole spherical shell.
Using the results given by eqs. (65,66,67), we integrate (73) over a spherical surface of radius r, and obtain the gravitational energy of the perturbation. We find where Ψ(t, r) = e −iωt ψ(r). In view of the orthogonality property of the Legendre Polinomials, π 0 P l (cos θ) 2 sin θdθ = 2 2l + 1 , we obtain In similarity to the axial case, the final result given by eq. (75) has real and imaginary parts, both in the exponential and in the frequency ω. As before, we take the real part of the whole expression. In the case of highly damped oscillations, we have ω r << ω i , and therefore we may use ω ≈ ω i .
The solution given by eq. (75) may be compared with the one obtained directly from the numerical integration of eq. (72). In Figure 3, the continuous line represents the numerical solution, whereas eq. (75) is represented by the dashed line. In the numerical solution, the function ψ was obtained by considering the full expression of the Zerilli potential, and the boundary conditions are exactly the same as considered in the case of axial perturbations. As before, the numerical solution was obtained considering r 1 = 2m + m/4 and r 2 = 200m in the boundary conditions for eq. (57). We believe that the errors in the approximate analytic solution, compared to the numerical solution, are due to the approximation made in eq. (57), where we neglected the Zerilli potential in the limit r >> m. In Figure 3, we are considering the initial instant of time of the perturbation. The time evolution of the energy of the perturbation is stable, since the gravitational energy of the perturbation is damped in the course of time, as wee see in Figure 4.

Final remarks
In this paper we obtained the expression of the gravitational energy of the perturbations due the quasinormal modes in the Schwarzschild space-time. The expressions given by (34,75) hold in the limit m << r << ∞, and depend on the integers n and l. The variation in space and in time of these quantities are displayed in Figures (1-4). These are the major results of the present analysis. The dependence on the coordinate r in expressions (34,75) is due to the very nature of the QNM. The functions h 0 and h 1 diverge at spacelike infinity, which is a feature of the QNM. However, the perturbation is assumed to be localized, and does not make sense in the limit r → ∞. In our understanding, the axial and polar quasi-normal perturbations represent ripples in the space-time geometry, as they do represent in the metric formulation of general relativity. However, the physical description and results displayed by Figures (1,2,3,4) can only be obtained in the TEGR, since the  the gravitational energy-momentum P a given by eqs. (8) and (11) cannot be established in the standard formulation of general relativity. One conclusion is that the gravitational energy of the black hole oscillates. However, since the energy of the perturbations is concentrated far from the event horizon, as we see in Figures (1,3), we conclude that the mass of the black hole does not vary in time, i.e., the mass of the black hole (restricted to interior of the event horizon) neither increases nor decreases with time. Figures (2,4) show that energy of the perturbations rapidly decays for large instants of time. Here we make a distinction between the gravitational energy of the black hole (the zeroth component of the gravitational energymomentum 4-vector P a ), and the mass of the black hole (an invariant of the Poincaré group, P a P a = −m 2 ). Note, in addition, that m is the mass of the black hole established in a frame where the black hole is at rest. In view of these considerations, we may conclude that the energy of the perturbations is a non-local effect, since it takes place sufficiently far from the event horizon. We may conjecture that the farther one is from the black hole, i.e., the farther an external perturbation is imparted to the black hole, the more energy is required to perturb it.
The dependence of expressions (34,75) on the integers n and l implies a discretization of the energy perturbations. This discretization must be further investigated. However, it is known that the very discretization of the QNM leads to a discretization of the area of the black hole's event horizon [20]. A possible interplay between the discretizations of the energy perturbations and of the area of the black hole is an issue of relevant interest.
The flux of gravitational radiation emitted in the process of damping of the quasinormal modes of a Schwarzschild black hole has been investigated in refs. [21,22,23], by means of the Landau-Lifshitz pseudotensor [24], assuming that at large distances from the black hole the space-time perturbation is represented by a plane gravitational wave. The analysis in the latter references was made assuming n = 0 , l = 2. In the present investigation, the expression for the flux of gravitational radiation given by eq. (9) is invariant under coordinate transformations, in contrast to the pseudotensors, that are coordinate dependent expressions. It would be interesting to address a specific physical configuration like the collision of two black holes, and compare the energy radiated after the collision (i) by means of eq. (39), and (ii) through the approach based on the latter references, in the context of the Landau-Lifshitz pseudotensor. This issue will be addressed elsewhere.
The amplitude A in the function ψ arises as a constant of integration in the solutions of eqs. (19) and (57). Let us consider the tortoise coordinate x = r+2M ln r 2M − 1 . If we make x = x+ 1 iω ln A, then it is easy to see that ψ = e iωx = Ae iωx . Therefore, there is an arbitrariness in the establishment of the amplitude A. Since the profiles of the gravitational energy perturbations in Figures (1-4) is very much similar, with a suitable choice of amplitudes we may locally enhance the similarity between the two perturbations, as we see in Figure (5). However, note that ψ is dimensionless in the axial perturbation, but has dimension (of distance) in the polar perturbation. But still, the similarities between the two perturbations is very clear. The procedure presented in this article may also be applied to extended formulations of gravity, namely, the f(T) type theories. In some of these approaches there are proposals for the gravitational energy, as for instance in refs. [25], [26], [27], [28], [29]. There are various motivations to address such extended theories, as for instance the analysis of the dark energy and dark matter problems, the problem of singularities of black holes, and the investigation of the degrees of freedom carried by gravitational waves. The analysis of the gravitational energy perturbations in these extended frameworks is, of course, of relevant interest.