Next Article in Journal
Gravitational Interaction in the Chimney Lattice Universe
Next Article in Special Issue
Variational Principles in Teleparallel Gravity Theories
Previous Article in Journal
Curved Momentum Space, Locality, and Generalized Space-Time
Previous Article in Special Issue
An Improved Framework for Quantum Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
José Wadih Maluf
1,*,
Sérgio Ulhoa
1,
Fernando Lessa Carneiro
1 and
Karlúcio H. C. Castello-Branco
2
1
Instituto de Física, Universidade de Brasília, Brasília, DF 70919-970, Brazil
2
Instituto de Ciências da Educação, Universidade Federal do Oeste do Pará, Santarém, PA 68040-070, Brazil
*
Author to whom correspondence should be addressed.
Universe 2021, 7(4), 100; https://doi.org/10.3390/universe7040100
Submission received: 1 March 2021 / Revised: 7 April 2021 / Accepted: 8 April 2021 / Published: 14 April 2021
(This article belongs to the Special Issue Teleparallel Gravity: Foundations and Observational Constraints)

Abstract

:
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 < < , 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 discretized. We also obtain a simple expression for the decrease of the flux of gravitational radiation of the perturbations.

1. 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 [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 then analyzing 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 no unique approach to the problem (e.g., pseudotensors and quasilocal expressions), 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 4-angular momentum [5,6], which 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 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 restrict 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 thus far. We conjecture that the perturbed energy in arbitrary (nonspherical) volumes of the Schwarzschild space-time is also characterized by integers.
Notation: space-time indices μ , ν , and SO(3,1) indices a , b , run from 0 to 3. Time and space indices are indicated according to μ = 0 , i , a = ( 0 ) , ( i ) . The tetrad field is denoted e a μ , and the torsion tensor reads T a μ ν = μ e a ν ν e a μ . The flat, Minkowski space-time metric tensor raises and lowers tetrad indices and is fixed by η a b = e a μ e b ν g μ ν = ( 1 , + 1 , + 1 , + 1 ) . The determinant of the tetrad field is represented by e = det ( e a μ ) .

2. 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 = k e ( 1 4 T a b c T a b c + 1 2 T a b c T b a c T a T a ) 1 c L M k e Σ a b c T a b c 1 c L M ,
where k = c 3 / 16 π G , T a = T b b a , T a b c = e b μ e c ν T a μ ν and
Σ a b c = 1 4 ( T a b c + T b a c T c a b ) + 1 2 ( η a c T b η a b T c ) .
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
e R ( e ) e ( 1 4 T a b c T a b c + 1 2 T a b c T b a c T a T a ) + 2 μ ( e T μ ) .
where R ( e ) is the scalar Riemannian curvature. However, the field equations derived from Equation (1) are invariant under local SO(3,1) transformations and are equivalent to Einstein’s equations. They read
e a λ e b μ ν ( e Σ b λ ν ) e ( Σ b ν a T b ν μ 1 4 e a μ T b c d Σ b c d ) = 1 4 k c e T a μ ,
where δ L M / δ e a μ = e T 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 make c = 1 = G ). Equation (3) may be rewritten as
ν ( e Σ a λ ν ) = 1 4 k e e a μ ( t λ μ + T λ μ ) ,
where T λ μ = e a λ T a μ and t λ μ is defined by
t λ μ = k ( 4 Σ b c λ T b c μ g λ μ Σ b c d T b c d ) .
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,
d d t V d 3 x e e a μ ( t 0 μ + T 0 μ ) = S d S j [ e e a μ ( t j μ + T j μ ) ] .
Therefore, we identify t λ μ as the gravitational energy–momentum tensor [5],
P a = V d 3 x e e a μ ( t 0 μ + T 0 μ ) ,
as the total energy–momentum contained within a volume V,
Φ g a = S d S j ( e e a μ t j μ ) ,
as the gravitational energy–momentum flux [7,8], and
Φ m a = S d S j ( e e a μ T j μ ) ,
as the energy–momentum flux of matter [8]. In view of (4), Equation (8) may be written as
P a = V d 3 x j Π a j = S d S j Π a j ,
where Π a j = 4 k e Σ a 0 j . A summary of all issues discussed above can be found in [9].
Equation (11) is the definition for the gravitational energy–momentum presented in [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 reparameterizations. 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 energy–momentum 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 e V μ (with space-time index μ ) would be coordinate dependent (the integral of a space-time vector density e V μ 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 center of mass moment are frame dependent field quantities, which 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 = m c 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 = γ m c 2 , where γ = ( 1 v 2 / c 2 ) 1 / 2 . Likewise, the gravitational momentum, angular momentum, and the center 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) such as j Π a j = 4 k j ( e Σ a 0 j ) in Equation (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 Equation (11). The existence of these well defined scalar densities is natural in the TEGR, as in the identity below Equation (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]
h 03 = h 0 ( t , r ) sin θ θ P l ( θ ) , h 13 = h 1 ( t , r ) sin θ θ P l ( θ ) .
where P l ( θ ) are the Legendre polynomials. The functions h 0 ( t , r ) and h 1 ( t , r ) satisfy the equations [1,10]
1 f ( r ) h 0 t [ f ( r ) h 1 ] r = 0 ,
1 f ( r ) [ 2 h 1 t 2 2 h 0 t r + 2 r h 0 t ] + 1 r 2 [ l ( l + 1 ) 2 ] h 1 = 0 ,
f ( r ) 2 [ 2 h 0 r 2 2 h 1 t r 2 r h 1 t ] + 1 r 2 [ r f r 1 2 l ( l + 1 ) ] h 0 = 0 ,
where f ( r ) = 1 2 m / 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
2 Ψ t 2 f r r f ( r Ψ ) r + 2 f 2 r 2 ( r Ψ ) r + f r 2 [ l ( l + 1 ) 2 ] Ψ = 0 .
In terms of the tortoise coordinate x = r + 2 m ln ( r / 2 m 1 ) , Equation (16) may be rewritten as
Ψ 2 t 2 Ψ 2 x 2 + V R W ( x ) Ψ ( t , x ) = 0 .
V R W ( x ) is the Regge–Wheeler potential, which is given by
V R W ( x ) = 1 2 m r l ( l + 1 ) r 2 6 m r 3 ,
and x = x ( r ) .
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
d 2 ψ d r 2 + [ ω 2 V R W ( x ) ] ψ ( x ) = 0 .
In view of the asymptotic behavior of the potential V R W ( x ) , the expressions of ψ ( x ) in the asymptotic limits x ( r ) and x ( r 2 m ) are given by ψ ( x ) A e i ω x and ψ ( x ) A e 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 quasinormal 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]
ω n = ln 3 8 π m i 4 m n + 1 2 + O ( n 1 / 2 ) ,
where n is an integer. This expression has also been obtained by means of analytic procedure in [14]. However, the frequency to be employed in the construction of the figures below, in Section 4 and Section 5, does not correspond to highly damped modes.
The problems regarding the divergent behavior of h 1 ( t , r ) in the limit r were discussed by ref. [1] (Section 4) and ref. [15]. However, the asymptotic behavior of h 1 ( t , r ) poses no problem to the result of our analysis, since we integrate (11) over a surface of constant (finite) radius r < < .
We need the asymptotic behavior of h 0 ( t , r ) . Given that 0 h 0 = i ω h 0 , and
0 h 0 = f [ f ( r ) h 1 ] r 1 h 1 ,
which follows from (13) in the limit r , we find
h 0 i f ω r Ψ r Ψ .
The m / r correction to the expressions above is 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 is useful in the following subsection. We find
1 h 0 0 h 1 2 r h 1 2 Ψ .

4. The Gravitational Energy of the Axial Perturbations

In this section, we address the axially perturbed Schwarzschild space-time and calculate the gravitational energy contained within a simple three dimensional volume, namely a spherical volume of radius r. To simplify the calculations, we 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 thus 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, 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]
d s 2 = 1 2 m r d t 2 + 1 2 m r 1 d r 2 + + r 2 d θ 2 + r 2 sin 2 θ d ϕ 2 + 2 h 0 sin θ θ P l ( cos θ ) d t d ϕ + 2 h 1 sin θ θ P l ( cos θ ) d r d ϕ .
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 establish a set of tetrad fields adapted to static observers in space-time. Therefore, we require e ( 0 ) i = 0 . This condition fixes three 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
e a μ = A 0 0 D 0 B sin θ cos ϕ r cos θ cos ϕ E r sin θ sin ϕ + F sin θ cos ϕ 0 B sin θ sin ϕ r cos θ sin ϕ E r sin θ cos ϕ + F sin θ sin ϕ 0 B cos θ r sin θ F cos θ ,
with the following definitions:
A = ( g 00 ) 1 / 2 , B = ( g 11 ) 1 / 2 , D = g 03 ( g 00 ) 1 / 2 , E 2 = 1 + g 11 g 03 2 + g 00 g 13 2 r 2 sin 2 θ , F = g 13 ( g 11 ) 1 / 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
e = r 2 sin θ 1 + f 1 g 03 2 r 2 sin 2 θ f g 13 2 r 2 sin 2 θ .
The gravitational energy is given by the a = ( 0 ) component of (11). Transforming the volume integral into a surface integral and considering the definition Π a j = 4 k Σ a 0 j , we find that
P ( 0 ) = S d S j Π ( 0 ) j = S d S j 4 k e Σ ( 0 ) 0 j .
The integration is 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
T 001 = A 1 A , T 301 = D 1 A , T 003 = A 0 D , T 103 = B 0 F , T 303 = r E 0 ( r E ) sin 2 θ D 0 D + F 0 F , T 013 = A 1 D , T 113 = B 1 F , T 212 = r ( 1 B ) , T 313 = r E ( 1 ( r E ) B ) sin 2 θ D 1 D + F 1 F , T 223 = r F .
The evaluation of Σ ( 0 ) 01 requires a large number of algebraic manipulations, but otherwise it is simple. The full expression of e Σ ( 0 ) 01 is given by
e Σ ( 0 ) 01 = r sin θ 1 + 1 f g 03 2 r 2 sin 2 θ f g 13 2 r 2 sin 2 θ × { 2 ( f 1 / 2 1 ) g 03 2 + 2 f [ f g 13 2 + 2 r 2 f 1 / 2 sin 2 θ r 2 sin 2 θ 1 + 1 + 1 f g 03 2 r 2 sin 2 θ f g 13 2 r 2 sin 2 θ ] + r f 1 / 2 g 03 ( 1 g 03 0 g 13 ) } × 4 g 03 2 + f f g 13 2 + r 2 sin 2 θ 1 .
By keeping terms of order ( h 0 ) 2 and ( h 1 ) 2 only in the integrand of (31), 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
4 k r d θ d ϕ [ sin θ ( g 00 ) 1 / 2 r ( 1 B ) ] = m ,
which is a well known result. We recall that B is given by (27).
The gravitational energy of the perturbations is denoted by δ P ( 0 ) and is obtained out of the second-order terms in (31), assuming g 00 1 and g 11 1 . Therefore, for fixed values of r such that m / r < < 1 , we obtain, after some simplifications,
δ P ( 0 ) = k d θ d ϕ sin θ r ( P l ( cos θ ) ) 2 h 0 2 3 h 1 2 r h 0 1 h 0 0 h 1 .
Recalling that k = 1 / 16 π (we are assuming c = 1 = G ), we find
δ P ( 0 ) = 1 8 d θ sin θ r ( P l ( cos θ ) ) 2 h 0 2 3 h 1 2 r h 0 1 h 0 0 h 1 .
Now, we make use of (22) and (23), and we arrive at
δ P ( 0 ) l ( l + 1 ) 2 l + 1 r Ψ 2 ,
where we use
0 π d θ sin θ θ ( P l ( cos θ ) ) 2 = 2 l ( l + 1 ) 2 l + 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 Equation (19). The boundary conditions for the numerical solution of Equation (19) are ψ ( r 1 ) = A exp i ω x 1 and ψ ( r 2 ) = A exp + i ω x 2 , where x 1 = r 1 + 2 m ln ( r 1 / 2 m 1 ) , x 2 = r 2 + 2 m ln ( r 2 / 2 m 1 ) , r 1 = 2 m + m / 4 and r 2 = 200 m .
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 shown in Figure 2.
Considering the real part of the solution given by Equation (34), we have
δ P ( 0 ) = l ( l + 1 ) 2 l + 1 r A 2 cos 2 ω r ( x t ) e 2 ω i ( x t ) .
We may calculate the variation of the perturbation energy in a half-life interval T / 2 , for weakly damped solutions, but neglecting the damping effects, considering only ω = 2 ω r . By making use of
1 T / 2 0 T / 2 d t sin ω t = 2 π ,
we obtain
δ P 1 / 2 ( 0 ) = 2 r A 2 π l ( l + 1 ) 2 l + 1 .
Since the black hole looses energy in the course of time, as shown in Figure 2, it is important to analyze 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 is applied to the radiation of gravitational energy in the Bondi space-time. The gravitational energy flux is given by the a = ( 0 ) component of Φ g a defined by (9), or simply by
Φ g ( 0 ) = d P ( 0 ) d t = d ( δ P ( 0 ) ) d t .
In the present context, we find
Φ g ( 0 ) = 2 l ( l + 1 ) 2 l + 1 A 2 r [ ω r sin [ 2 ω r ( x t ) ] + ω i cos [ 2 ω r ( x t ) ] ] e 2 ω i ( x t ) .
Again, we define ω = 2 ω r and find that after the period T = 2 π / ω the decrease in the flux of gravitational radiation is Φ g ( 0 ) ( t + T ) Φ g ( 0 ) ( t ) = Φ g ( 0 ) ( t ) ( 1 e 2 π ω i / ω r ) , and therefore
Φ g ( 0 ) ( t + T ) = Φ g ( 0 ) ( t ) e 2 π ω i / ω r .
Expressions (39) and (40) are clearly connected to the variation in time of the axial energy perturbation displayed in Figure 2.

5. 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, despite the similarities with the calculations carried out for the axial perturbations. We consider the perturbed metric tensor as given by g μ ν = g ¯ μ ν + h μ ν . The line element for the unperturbed metric g ¯ μ ν tensor is
d s 2 = f d t + f 1 d r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ,
and the polar perturbations h μ ν are described by [10]
h μ ν = e i ω t P l ( cos θ ) H 0 f H 1 0 0 H 1 H 2 f 1 0 0 0 0 r 2 K 0 0 0 0 r 2 K sin 2 θ ,
where f = 1 2 m / r and H 0 , H 1 , H 2 , and K are functions of t , r . Therefore, the perturbed metric tensor is written as
g μ ν = f ( 1 H 0 ) H 1 0 0 H 1 f 1 ( 1 + H 2 ) 0 0 0 0 r 2 ( 1 + K ) 0 0 0 0 r 2 ( 1 + K ) sin 2 θ ,
where we define
H 0 = H 0 e i ω t P l ( cos θ ) ,
H 1 = H 1 e i ω t P l ( cos θ ) ,
H 2 = H 2 e i ω t P l ( cos θ ) ,
K = K 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]
H 2 = H 0 H .
The three first-order differential equations are [18]
d K d r + r 3 m r ( r 2 m ) K 1 r H + 1 2 l ( l + 1 ) i ω r 2 H = 0 ,
d H d r + r 3 m r ( r 2 m ) K r 4 m r ( r 2 m ) H + i ω r r 2 m + 1 2 l ( l + 1 ) i ω r 2 H 1 = 0 ,
d H 1 d r + i ω r r 2 m K + i ω r r 2 m H + 2 m r ( r 2 m ) H 1 = 0 .
The three remaining equations eventually yield the algebraic equation
6 m r + 2 λ H + 2 λ 2 ω 2 r 3 r 2 m + 2 m ( r 3 m ) r ( r 2 m ) K + 2 i ω r + M l ( l + 1 ) i ω r 2 H 1 = 0 ,
where λ 1 2 ( l 1 ) ( l + 2 ) .
By rewriting the field quantities according to [18]
K = α ( r ) ψ + β ( r ) Θ ,
and
R 1 ω H 1 = γ ( r ) ψ + ρ ( r ) Θ ,
and defining the tortoise coordinate x as
d x = 1 f d r x = r + 2 m ln r 2 m 1 ,
we obtain the relations
α ( r ) = λ ( λ + 1 ) r 2 + 3 λ m r + 6 m 2 r 2 ( λ r + 3 m ) ,
β ( r ) = 1 ,
γ ( r ) = i λ r 2 + 3 λ m r + 3 m 2 ( r 2 m ) ( λ r + 3 m ) ,
ρ ( r ) = i r 2 r 2 m .
It follows that
Θ = d ψ d x
and
d 2 ψ d x 2 + ω 2 V Z e ψ = 0 ,
where
V Z e = f ( r ) 2 λ 2 ( λ + 1 ) r 3 + 6 λ 2 m r 2 + 18 λ m 2 r + 18 m 3 r 3 ( λ r + 3 m ) 2
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 Z e 0 . Thus, in the limit r > > m , we obtain a simple solution for the function ψ , given by
ψ ( x ) e i ω x ψ ( x ) = A e i ω x ,
where A is a constant with dimension of length. 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 2 i ω x | < < | A e 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
| A e i ω r x e ω i x | < < 1 A < < e ω i x = e | ω i | x .
It is important to remark that the stability of the solutions demand that the function Ψ = e i ω t ψ decay sufficiently quickly with time. Consequently, the stable solutions are those for which ω i < 0 . The solutions obtained in the following approximations are valid as long as
x < < 1 | ω i | ln A .
Therefore, the approximate analytic solutions are valid provided the condition < < x < < 1 | ω i | ln A is satisfied.
As in the axial case, we must find the asymptotic limits for the perturbative functions. Thus, in the limit r > > m , Equations (52)–(55) become
α ( r ) λ + 1 r ,
β ( r ) = 1 ,
γ ( r ) i ,
ρ ( r ) i r .
As a consequence, Equations (49) and (50) are approximated by
K Θ = d ψ d x = λ + 1 r + i ω ψ
and
H 1 ω 2 r ψ .
Making use of Equation (48) in the limit r > > m , together with Equations (65) and (66), we obtain
H ω 2 r ψ .
The expressions above given by Equations (65)–(67) are used in the following subsection, where we evaluate the gravitational energy of the polar perturbations.

The Gravitational Energy of the Polar Perturbations

Similar to the analysis in Section 4, where we address the gravitational energy of the axial perturbations, here we construct a set of tetrad fields adapted to static observers (with respect to the asymptotic flat space-time limit of the black hole), which yields the metric tensor (43). A set of tetrad fields that satisfy the necessary conditions is given by
e a μ = A B 0 0 0 C sin θ cos ϕ D r cos θ cos ϕ D r sin θ sin ϕ 0 C sin θ sin ϕ D r cos θ sin ϕ D r sin θ cos ϕ 0 C cos θ D r sin θ 0 ,
where
A = f 1 / 2 ( 1 H ) 1 / 2 , A B = H 1 , A C = [ 1 H 2 + H 1 2 ] 1 / 2 , D = ( 1 + K ) 1 / 2 ,
and whose determinant is e = D 2 A C r 2 sin θ .
The components of the torsion tensor T λ μ ν and the tensor Σ λ μ ν that are needed in the calculations below are
T 212 = 1 2 1 ( D 2 r 2 ) D C r , T 313 = 1 2 1 ( D 2 r 2 sin 2 θ ) D C r sin 2 θ , T 202 = 1 2 0 ( D 2 r 2 ) , T 303 = 1 2 0 ( D 2 r 2 sin 2 θ ) ,
and
Σ 001 = 1 2 ( g 01 g 01 g 00 g 11 ) ( g 22 T 212 + g 33 T 313 ) , Σ 101 = 1 2 ( g 01 g 01 g 00 g 11 ) ( g 22 T 202 + g 33 T 303 ) .
With the help of these quantities, we find
Π ( 0 ) 1 = 4 k sin θ r 2 2 C 1 ( D 2 ) + D 2 r C D r B r 2 2 A C 0 ( D 2 ) .
The expression above leads to the gravitational energy given by Equation (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 Equation (72) is simply Π g ( 0 ) 1 = 4 k 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 Π g ( 0 ) 1 . The approximate analytic expression of the latter is
δ Π ( 0 ) 1 4 k sin θ r 2 2 f 1 / 2 1 K + r K f 1 / 2 1 2 1 2 r f 1 / 2 H + 4 k sin θ [ r 2 f 1 / 2 3 4 H 1 H 1 2 K H + r 8 K 2 + r 2 2 f 1 / 2 H 1 0 K r 2 4 f 1 / 2 H 1 K ] .
The first line in (73) contains the first order terms δ Π f i r s t ( 0 ) 1 , whereas the other terms correspond to the second-order terms δ Π s e c ( 0 ) 1 . The total energy within a spherical surface of radius r is obtained by summing these two quantities, δ Π ( 0 ) 1 = δ Π f i r s t ( 0 ) 1 + δ Π s e c o n d ( 0 ) 1 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 ( c o s θ ) 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 Equations (65)–(67), we integrate (73) over a spherical surface of radius r and obtain the gravitational energy of the perturbation. We find
δ P ( 0 ) 1 16 r 3 ω 4 Ψ ( t , r ) 2 0 π P l ( cos θ ) 2 sin θ d θ ,
where Ψ ( t , r ) = e i ω t ψ ( r ) . In view of the orthogonality property of the Legendre polynomials,
0 π P l ( cos θ ) 2 sin θ d θ = 2 2 l + 1 ,
we obtain
δ P ( 0 ) 1 8 r 3 ω 4 2 l + 1 Ψ ( t , r ) 2 .
Similar to the axial case, the final result given by Equation (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 Equation (75) may be compared with the one obtained directly from the numerical integration of Equation (72). In Figure 3, the continuous line represents the numerical solution, whereas Equation (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 = 2 m + m / 4 and r 2 = 200 m in the boundary conditions for Equation (57). We believe that the errors in the approximate analytic solution, compared to the numerical solution, are due to the approximation made in Equation (57), where we neglected the Zerilli potential in the limit r > > m . In Figure 3, we consider 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 shown in Figure 4.

6. Final Remarks

In this paper, we obtain the expression of the gravitational energy of the perturbations due the quasinormal modes in the Schwarzschild space-time. The expressions given by (34) and (75) hold in the limit m < < r < < and depend on the integers n and l. The variation in space and time of these quantities are displayed in Figure 1, Figure 2, Figure 3 and Figure 4. These are the major results of the present analysis. The dependence on the coordinate r in expressions (34) and (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 quasinormal 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 in Figure 1, Figure 2, Figure 3 and Figure 4 can only be obtained in the TEGR, since the the gravitational energy–momentum P a given by Equations (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 shown in Figure 1 and Figure 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. Figure 2 and Figure 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 energy–momentum 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) and (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 is investigated in [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 Equation (9) is invariant under coordinate transformations, in contrast to the pseudotensors, which are coordinate dependent expressions. It would be interesting to address a specific physical configuration, such as the collision of two black holes, and compare the energy radiated after the collision: (i) by means of Equation (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 Equations (19) and (57). Let us consider the tortoise coordinate x = r + 2 M ln r 2 M 1 . If we make x = x + 1 i ω ln A , then it is easy to see that ψ = e i ω x = A e i ω x . Therefore, there is an arbitrariness in the establishment of the amplitude A. Since the profiles of the gravitational energy perturbations in Figure 1, Figure 2, Figure 3 and Figure 4 are very similar, with a suitable choice of amplitudes, we may locally enhance the similarity between the two perturbations, as shown in Figure 5. Note that ψ is dimensionless in the axial perturbation but has dimension (of distance) in the polar perturbation. However, the similarities between the two perturbations are 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 (see, e.g., [25,26,27,28,29]). There are various motivations to address such extended theories, 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.

Author Contributions

Conceptualization, The four authors contributed to the conceptualization; Formal analysis, The four authors contributes to the formal analysis; Investigation, S.U., F.L.C. and K.H.C.C.-B.; Methodology, S.U., F.L.C. and K.H.C.C.-B.; Software, F.L.C.; Supervision, J.W.M.; Writing–original draft, J.W.M.; Writing–review & editing, J.W.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work received no funding other than the regular institutional support from Universidade de Brasília and Universidade Federal do Oeste do Pará.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

There is no conflict of interest, of any kind, in the development and conclusion of this work.

References

  1. Nollert, H.-P. Quasinormal modes: The characteristic ‘sound’ of black holes and neutron stars. Class. Quantum Grav. 1999, 16, R159. [Google Scholar]
  2. Kokkotas, K.D.; Schmidt, B.G. Quasi-Normal Modes of Stars and Black Holes. Living Rev. Rel. 1999, 2, 2. [Google Scholar]
  3. Chandrasekhar, S. Colloquium at Syracuse University; Syracuse University: Syracuse, NY, USA, 1992. [Google Scholar]
  4. Schutz, B.F.; Will, C.W. Black hole normal modes—A semianalytic approach. Astrophys. J. 1985, 291, L33. [Google Scholar]
  5. Maluf, J.W. The gravitational energy-momentum tensor and the gravitational pressure. Ann. Phys. 2005, 14, 723. [Google Scholar]
  6. Maluf, J.W.; da Rocha-Neto, J.F.; Toríbio, T.M.L.; Castello-Branco, K.H. Energy and angular momentum of the gravitational field in the teleparallel geometry. Phys. Rev. D 2002, 65, 124001. [Google Scholar]
  7. Maluf, J.W.; Faria, F.F.; Castello-Branco, K.H. The gravitational energy–momentum flux. Class. Quantum Grav. 2003, 20, 4683. [Google Scholar]
  8. Maluf, J.W.; Faria, F.F. On gravitational radiation and the energy flux of matter. Ann. Phys. 2004, 13, 604. [Google Scholar]
  9. Maluf, J.W. The teleparallel equivalent of general relativity. Ann. Phys. 2013, 525, 339. [Google Scholar]
  10. Regge, T.; Wheeler, J.A. Stability of a Schwarzschild singularity. Phys. Rev. 1957, 108, 1063. [Google Scholar]
  11. Iyer, S.; Will, C.M. Black-hole normal modes: A WKB approach. I. Foundations and application of a higher-order WKB analysis of potential-barrier scattering. Phys. Rev. D 1987, 35, 3621. [Google Scholar]
  12. Iyer, S. Black-hole normal modes: A WKB approach. II. Schwarzschild black holes. Phys. Rev. D 1987, 35, 3632. [Google Scholar]
  13. Konoplya, R.A. Quasinormal behavior of the D-dimensional Schwarzschild black hole and the higher order WKB approach. Phys. Rev. D 2003, 68, 024018. [Google Scholar]
  14. Berti, E.; Cardoso, V.; Starinets, A.O. Quasinormal modes of black holes and black branes. Class. Quantum Grav. 2009, 26, 163001. [Google Scholar]
  15. Nollert, H.-P.; Schmidt, G.G. Quasinormal modes of Schwarzschild black holes: Defined and calculated via Laplace transformation. Phys. Rev. D 1992, 45, 2617. [Google Scholar]
  16. Maluf, J.W.; Faria, F.F.; Ulhoa, S.C. On reference frames in spacetime and gravitational energy in freely falling frames. Class. Quantum Grav. 2007, 24, 2743. [Google Scholar]
  17. Hehl, F.H.; Lemke, J.; Mielke, E.W. Two Lectures on Fermions and Gravity. In Geometry and Theoretical Physics; Debrus, J., Hirshfeld, A.C., Eds.; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  18. Zerilli, F.J. Effective potential for even parity Regge-Wheeler gravitational perturbation equations. Phys. Rev. Lett. 1970, 24, 737. [Google Scholar]
  19. Chandrasekhar, S.; Detweiler, S. The Quasi-Normal Modes of the Schwarzschild Black Hole. Proc. R. Soc. Lond. A Math. Phys. Sci. 1975, 344, 441. [Google Scholar]
  20. Maggiore, M. Physical Interpretation of the Spectrum of Black Hole Quasinormal Modes. Phys. Rev. Lett. 2008, 100, 141301. [Google Scholar] [PubMed] [Green Version]
  21. Cunningham, C.T.; Price, R.H.; Moncrief, V. Radiation from collapsing relativistic stars. I. Linearized odd-parity radiation. Astrophys. J. 1978, 224, 643. [Google Scholar]
  22. Nagar, A.; Rezolla, L. Gauge-invariant Non-spherical Metric Perturbations of Schwarzschild Black-Hole Spacetimes. Class. Quantum Grav. 2005, 22, R167. [Google Scholar]
  23. Nakano, H.; Ioka, K. Second Order Quasi-Normal Mode of the Schwarzschild Black Hole. Phys. Rev. D 2007, 76, 084007. [Google Scholar] [CrossRef] [Green Version]
  24. Landau, L.D.; Lifshitz, E.M. The Classical Theory of Fields, 4th ed.; Pergamon Press: Oxford, UK, 1975. [Google Scholar]
  25. Capozziello, S.; Capriolo, M.; Transirico, M. The gravitation energy–momentum pseudotensor: The cases of F(R) and F(T) gravity. Int. J. Geom. Methods Mod. Phys. 2018, 15, 1850164. [Google Scholar]
  26. Capozziello, S.; Capriolo, M.; Transirico, M. The gravitational energy-momentum pseudo-tensor of higher-order theories of gravity. Ann. Phys. 2017, 529, 1600376. [Google Scholar] [CrossRef] [Green Version]
  27. Capozziello, S.; Capriolo, M.; Caso, L. Gravitational waves in higher order teleparallel gravity. Class. Quantum Grav. 2020, 37, 235013. [Google Scholar]
  28. Capozziello, S.; Capriolo, M.; Caso, L. Weak field limit and gravitational waves in higher-order gravity. Int. J. Geom. Methods Mod. Phys. 2019, 16, 1950047. [Google Scholar]
  29. Ulhoa, S.C.; Spaniol, E.P. On the Gravitational Energy-momentum Vector in f(T) Theories. Int. J. Mod. Phys. D 2013, 22, 1350069. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Comparison between the real parts of the gravitational energy of the axial perturbation. The continuous line represents the values of the numerical integration of the energy, whereas the dashed line represents the values obtained from the approximate analytic expression (34). The parameters used are A = 10 6 , m = 1 . The frequency was required to be the fundamental frequency ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained from the third-order WKB method. The data represent the initial instant of time of the perturbation, i.e., t = 0 .
Figure 1. Comparison between the real parts of the gravitational energy of the axial perturbation. The continuous line represents the values of the numerical integration of the energy, whereas the dashed line represents the values obtained from the approximate analytic expression (34). The parameters used are A = 10 6 , m = 1 . The frequency was required to be the fundamental frequency ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained from the third-order WKB method. The data represent the initial instant of time of the perturbation, i.e., t = 0 .
Universe 07 00100 g001
Figure 2. Comparison between the real parts of the gravitational energy of the axial perturbation, as a function of time. The continuous line represents the values of the numerical integration of the energy, whereas the dashed line represents the values obtained from the approximate analytic expression (34). The parameters used are A = 10 6 , m = 1 . The frequency was required to be the fundamental frequency ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained from the third-order WKB method. The radius of integration was chosen to be r = 60 , in natural units.
Figure 2. Comparison between the real parts of the gravitational energy of the axial perturbation, as a function of time. The continuous line represents the values of the numerical integration of the energy, whereas the dashed line represents the values obtained from the approximate analytic expression (34). The parameters used are A = 10 6 , m = 1 . The frequency was required to be the fundamental frequency ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained from the third-order WKB method. The radius of integration was chosen to be r = 60 , in natural units.
Universe 07 00100 g002
Figure 3. Comparison between the real parts of the gravitational energy of the polar perturbation. The continuous line represents the numerical value of the evaluation, and the dashed line corresponds to the value resulting of the approximate analytic expression (75). The parameters used in the analysis are A = 10 6 , m = 1 , and the frequency considered is the fundamental one ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The data represent the initial instant of time of the perturbation, i.e., t = 0 .
Figure 3. Comparison between the real parts of the gravitational energy of the polar perturbation. The continuous line represents the numerical value of the evaluation, and the dashed line corresponds to the value resulting of the approximate analytic expression (75). The parameters used in the analysis are A = 10 6 , m = 1 , and the frequency considered is the fundamental one ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The data represent the initial instant of time of the perturbation, i.e., t = 0 .
Universe 07 00100 g003
Figure 4. Comparison between the real parts of the gravitational energy of the polar perturbation. The continuous line represents the numerical value of the evaluation, and the dashed line corresponds to the value resulting of the approximate analytic expression (75). The parameters used in the analysis are A = 10 6 , m = 1 , and the frequency considered is the fundamental one ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The integration is made on a surface of radius r = 60 , in natural units.
Figure 4. Comparison between the real parts of the gravitational energy of the polar perturbation. The continuous line represents the numerical value of the evaluation, and the dashed line corresponds to the value resulting of the approximate analytic expression (75). The parameters used in the analysis are A = 10 6 , m = 1 , and the frequency considered is the fundamental one ( n = 0 , l = 2 ) ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The integration is made on a surface of radius r = 60 , in natural units.
Universe 07 00100 g004
Figure 5. Comparison between the real parts of the gravitational energy of the perturbations. The dashed line represents the numerical expression of the energy for the axial perturbations and the dotted line represents the numerical energy of the polar perturbations. In natural units (whenever necessary), we have A a x i a l = 10 6 , A p o l a r = 5 1 / 2 10 6 , m = 1 , and the frequency (fundamental, ( n = 0 , l = 2 )) is ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The data represent the initial instant of time of the perturbations, i.e., t = 0 .
Figure 5. Comparison between the real parts of the gravitational energy of the perturbations. The dashed line represents the numerical expression of the energy for the axial perturbations and the dotted line represents the numerical energy of the polar perturbations. In natural units (whenever necessary), we have A a x i a l = 10 6 , A p o l a r = 5 1 / 2 10 6 , m = 1 , and the frequency (fundamental, ( n = 0 , l = 2 )) is ω = 0.373162 i 0.089217 , obtained by means of the third-order WKB method. The data represent the initial instant of time of the perturbations, i.e., t = 0 .
Universe 07 00100 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maluf, J.W.; Ulhoa, S.; Carneiro, F.L.; Castello-Branco, K.H.C. Perturbations of the Gravitational Energy in the TEGR: Quasinormal Modes of the Schwarzschild Black Hole. Universe 2021, 7, 100. https://doi.org/10.3390/universe7040100

AMA Style

Maluf JW, Ulhoa S, Carneiro FL, Castello-Branco KHC. Perturbations of the Gravitational Energy in the TEGR: Quasinormal Modes of the Schwarzschild Black Hole. Universe. 2021; 7(4):100. https://doi.org/10.3390/universe7040100

Chicago/Turabian Style

Maluf, José Wadih, Sérgio Ulhoa, Fernando Lessa Carneiro, and Karlúcio H. C. Castello-Branco. 2021. "Perturbations of the Gravitational Energy in the TEGR: Quasinormal Modes of the Schwarzschild Black Hole" Universe 7, no. 4: 100. https://doi.org/10.3390/universe7040100

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop