On the energy of a non-singular black hole solution satisfying the weak energy condition

The energy-momentum localization for a new four-dimensional and spherically symmetric, charged black hole solution that through a coupling of general relativity with non-linear electrodynamics is everywhere non-singular while it satisfies the weak energy condition is investigated. The Einstein and M\{o} ller energy-momentum complexes have been employed in order to calculate the energy distribution and the momenta for the aforesaid solution. It is found that the energy distribution depends explicitly on the mass and the charge of the black hole, on two parameters arising from the space-time geometry considered, and on the radial coordinate. Further, in both prescriptions all the momenta vanish.In addition, a comparison of the results obtained by the two energy-momentum complexes is made, whereby some limiting and particular cases are pointed out.


Introduction
Even though the problem of energy-momentum localization has triggered a lot of interesting research work, it still remains not fully answered for more than a century. The only step one can hope to make forward is to find a more powerful tool for accessing this issue of general relativity.
The use of different tools for energy-momentum localization, like superenergy tensors [1]- [4], quasi-local definitions [5]- [9] and the energy-momentum complexes [10]- [16] has led to the development of several interesting works. In particular, the energy-momentum complexes of Einstein [10], Landau-Lifshitz [12], Papapetrou [13], Bergmann-Thomson [14] and Weinberg [16] are pseudotensorial quantities and coordinate dependent. They can be used in Cartesian and quasi-Cartesian coordinates, more precisely in Schwarzschild Cartesian coordinates and in Kerr-Schild Cartesian coordinates, and have yielded so far many physically meaningful results [17]- [23]. The Møller energy-momentum complex [15] allows the calculation of energy and momenta in any coordinate system, including Schwarzschild Cartesian coordinates and Kerr-Schild Cartesian coordinates, and has also provided physically interesting results for many space-time geometries, in particularly for (3 + 1), (2 + 1) and (1 + 1) space-times [24]- [44]. It is worth noting that different pseudotensors yield the same energy value for any metric of the Kerr-Schild class and also for solutions more general than those of the Kerr-Schild class (for reviews and references, see the works [48], [49] and [50]). Moreover, the localization of energy-momentum has been studied in the context of teleparallel theory of gravitation whereby many similar results have been obtained [51]- [62].
The Einstein, Landau-Lifshitz, Papapetrou, Bergmann-Thomson, Weinberg and Møller prescriptions are in agreement with the definition of the quasi-local mass given by Penrose [45] and developed by Tod [46] for some gravitating systems. Recently, the new concept in the localization of energy is that of quasi-local energy-momentum associated with a closed 2-surface. In fact, it has been demonstrated that "the quasi-local quantities could provide a more detailed characterization of the states of the gravitational field" [47]. The energy localization is also connected with the quasi-local energy given by Wang and Yau [63], [64]. The rehabilitation of energy-momentum complexes concerns the searching for a common quasi-local energy value. Recently, an important discovery has been made, namely, by considering pseudotensors and quasi-local approaches in the context of the Hamiltonian formulation and with the choice of a 4D isometric Minkowski reference geometry on the boundary, it is found that for any closed 2-surface there is a common value for the quasi-local energy for all expressions that are in agreement (to linear order) with the Freud superpotential or, put simply, all the quasi-local expressions in a large class yield the same energy-momentum [65], [66]. This paper is organized as follows. In Section 2 we briefly present the non-singular black hole solution that satisfies the weak energy condition. In Section 3 the two prescriptions of Einstein and Møller used for the calculation of the energy distribution and momenta are introduced. In Section 4 we show the results obtained for the aforementioned non-singular black hole solution by applying the pseudotensorial definitions given in Section 3. Finally, in the Discussion we summarize the results and we focus on commenting about some limiting and particular cases. Throughout we have used geometrized units (c = G = 1) and for the signature the choice has been (+,−,−,−). In the case of the Einstein prescription we have used the Schwarzschild Cartesian coordinates {t, x, y, z} and for the Møller prescription the Schwarzschild coordinates {t, r, θ, φ}, respectively. Further, Greek indices take values from 0 to 3, while Latin indices run from 1 to 3.

The Non-Sinular Black Hole Solution Satisfying the Weak Energy Condition
In this section the new spherically symmetric and charged non-singular black hole solution that satisfies the weak energy condition developed by L. Balart and E. C. Vagenas [67] is introduced. This black hole solution has been found in the context of general relativity coupled to non linear electrodynamics via a term L(F ) in the action, with the Lorentz-invariant scalar F = 1 4 F µν F µν , and F µν is the electromagnetic tensor that reduces to the Maxwellian field tensor in the weak field case. Here, F µν is restricted only to the electric field. The black hole metric has been constructed by employing the Dagum distribution [68], a continuous probability distribution, given by where x > 0 , while the parameters a, b, p ∈ R + . For p = 1 a and b = 1, σ(x) can be written as Here, with the transformation x → q 2 M r the black hole metric function reads with γ ∈ R + , while M and q are the mass and the charge of the black hole, respectively. For β ≥ 3 (a+1) this black hole metric is non-singular everywhere, while for β ≤ 3 (a+1) it satisfies the weak energy condition T µν u µ u ν ≥ 0 for every timelike vector u µ , according to which the local energy density must be positive-definite. This general inequality is shown in [69] to be equivalent to two inequalities on the first and second derivative of the mass function w.r.t. r. In order to have a non-singular black hole solution that also satisfies the weak energy condition, one imposes the condition β = 3 (a+1) , so that the obtained new spherically symmetric, static and charged non-singular black hole metric has a line element of the form with , and the metric function now reads The corresponding electric field is For small values of the radial coordinate r the family of metrics derived from (5) behaves as a de Sitter black hole while asymptotically it behaves as the Schwarzschild black hole. It is important to notice that a particuar case is obtained for a = 1 (hence β = 3/2) and γ = 1 6 , so that the metric function (5) reads and the black hole solution asymptotically behaves as the Reissner-Nordström solution. For all the other values of the parameters a and γ the non-singular black hole space-time geometries with the metric function (5) asymptotically do not have the behaviour of the Reissner-Nordström solution. It must also be pointed out that only for the values a = 1 and γ = 1 6 one gets a nonlinear electrodynamics model leading to Maxwell's electrodynamics in the weak field approximation.

Einstein and Møller Prescriptions
The Einstein energy-momentum complex [10] defined for a (3 + 1) dimensional space-time is given by The superpotentials h µλ ν are expressed as and satisfy the necessary antisymmetric property In the Einstein prescription the local conservation law holds: The energy and momentum can be calculated with where θ 0 0 and θ 0 i represent the energy and momentum density components, respectively.
Applying Gauss' theorem, the energy-momentum becomes with n i the outward unit normal vector over the surface dS. In eq. (14) P 0 represents the energy. The expression for the Møller energy-momentum complex [15] is given by with the Møller superpotentials M µλ The Møller superpotentials M µλ ν are also antisymmetric Møller's energy-momentum complex also satifies the local conservation law In (15) J 0 0 is the energy density and J 0 i represents the momentum density components.
For the Møller prescription, the energy and momentum distributions are given by and the energy distribution can be evaluated with Using Gauss' theorem one obtains

Energy-Momentum Distribution of the Non-Singular Black Hole Solution Satisfying the Weak Energy Condition
In order to perform the calculations using the Einstein prescription, the metric given by the line element (4) is transformed into Schwarzschild Cartesian coordinates applying the coordinate transformation x = r sin θ cos ϕ, y = r sin θ sin ϕ, z = r cos θ. Then, the line element takes the following form: The calculations for µ = 0, 1, 2, 3 and i = 1, 2, 3 for the components of the superpotential h 0i µ (needed in (14)) in quasi-Cartesian coordinates are while the non-vanishing components of the superpotential are With the aid of the line element (22), the expression for the energymomentum distribution (14) and the expressions (24)-(26) for the superpotentials, the energy distribution in the Einstein prescription for the new charged non-singular black hole solution satisfying the weak energy condition is obtained as while, due to (23), we find that all the momentum components vanish: while all the other components vanish. Inserting the above expression (29) of the Møller superpotential into (21), we obtain the energy distribution in the Møller prescription: from which we infer that the energy in the Møller prescription can be given in terms of the energy in the Einstein prescription.
Furthermore, as expected from the vanishing of the spatial components of the Møller superpotential, all the momentum components are found to  be zero: In Fig. 3 the energy distribution in the Møller prescription for four different values of the parameter a and M = 1, q = 0.1, and γ = 1 is plotted, while in Fig. 4 the behaviour of the Møller energy distribution near the origin is presented for four different values of the parameter a and M = 1,  A comparison of the energy distributions in the Einstein and Møller prescriptions is given in Fig. 5 for a = 0.5 and M = 1, q = 0.1, and γ = 1. Additionally, in Fig. 6 an analogous comparison of the energy distributions in the Einstein and Møller prescriptions near the origin is presented for the same values of the parameters a, M , q, and γ.

Discussion
The aim of this work is the study of the energy-momentum distribution for a new spherically symmetric and charged, non-singular black hole solution satisfying the weak energy condition.
To this purpose, the Einstein and Møller energy-momentum complexes have been applied and it has been found that all the momenta vanish in both pseudotensorial prescriptions. Additionally, the energy distributions obtained have well-defined expressions showing a dependence on the mass M , the charge q, the two parameters γ and a, introduced in Section 2, and on the radial coordinate r. Furthermore, the Møller energy is expressed in terms of the Einstein energy. Both energies acquire the same value M (ADM mass) for r → ∞ or for q = 0.
We have also studied the limiting behavior of the energy distribution for r → 0 and r → ∞, and for the particular cases q = 0 and a = 1 and γ = 1 6 , respectively. In Table 1 we summarize the physically meaningful results for these limiting and particular cases.
At this point some remarks are need to clarify the results. From Table 1 we conclude that the energy distribution in both prescriptions vanishes near the origin r → 0, while for either r → ∞ or q = 0 both energies are equal to the ADM mass M , in agreement to the result obtained by Virbhadra for the energy distribution of the Schwarzschild black hole solution [50]. For the particular case a = 1 and γ = 1 6 we get an expression for the energy distribution that depends on the mass M and the charge q of the black hole, as well as on the radial coordinate r for a model of general relativity coupled to nonlinear electrodynamics, the latter corresponding to Maxwell's theory in the weak field approximation. Furthermore, the non-singular black hole solution with β = 3 2 , a = 1 and γ =   In Fig. 7 the energy distributions in the Einstein prescription E E and in the Møller prescription E M are presented, as a function of r, for the particular case a = 1 and γ = 1 6 , while in Fig. 8 the same functions for the same values of the parameters are presented near the origin. Looking closer at the behaviour of the energy distribution near zero, it is seen from Fig. 2 that the Einstein energy tends to zero from positive values, as expected from expression (27), and it is an increasing function of r. The Møller energy, on the other hand, tends also to zero but, according to Fig. 4, close enough to zero, i.e. for r 0.04, it acquires negative values before reaching zero, a result which is supported by (30). For values of r greater than 0.04, the Møller energy is positive and increasing.
In Fig. 5 and Fig. 6, we consider the particular value a = 0.5 as an example in order to compare the two energies obtained and we conclude that the Einstein energy E E is everywhere greater than the Møller energy. In fact, one can advocate that the positive energy region can be used for the effect of a convergent gravitational lens [70]- [72].
The  (30)). In fact, the energy becomes negative for r a < 2γ( q 2 M ) a . As an example, in Fig. 9 we see that the energy distribution becomes negative for 0 < r < 0.0275 if we choose a = 0.7 and γ = 1. For larger values of r, the energy distribution becomes positive and then it increases.
As pointed out in the Introduction, a modern viewpoint regarding the energy-momentum localization is represented by a quasi-local energy-momentum associated with a closed 2-surface. The pseudotensors have this property, mostly since Penrose developed his definition of quasilocal mass [45]. Further, according to [47], since a generally accepted expression for the energy-momentum has not been found until today, some important criteria for the quasi-local energy-momentum expressions must be satisfied, basically (a) the expressions should yield the standard values at spatial infinity, and (b) they should not give negative values for the energy.
In the case of the black hole solution that satisfies the weak energy condition and is under study in the present work, both the Einstein and Møller energymomentum complexes do yield the standard values at infinity. In fact, in the Einstein prescription the energy takes only positive values, while in the case of the Møller prescription the energy takes negative values when r satisfies the inequality r a < 2γ( q 2 M ) a . We believe that this apparent weakness of the Møller prescription can be justified by the properties of the particular metric considered. Indeed, a similar behaviour of the Møller energy-momentum complex was found in [73] and [74]. This strange behaviour is attributed to the particularities of these black hole solutions originating in the coupling of the gravitational field to non-linear electrodynamics.
Although the energy distribution in the Møller prescription does not exhibit the desired behaviour for every value of r, we can conclude that the two prescriptions used in the present work still constitute instructive tools for the energy-momentum localization. Given the problem that arises from the Møller prescription, it is challenging to employ other pseudo-tensorial prescriptions available, as well as the teleparallel equivalent theory of general relativity, for further investigation of the energy-momentum localization in the context of the four-dimensional spherically symmetric and charged, non-singular black hole solution satisfying the weak energy condition. We consider facing this challenge in a future work.