Consistent thermodynamics and topological classes for the four-dimensional Lorentzian charged Taub-NUT spacetimes

In this paper, we derive the consistent thermodynamics of the four-dimensional Lorentzian Reissner-Nordström-NUT (RN-NUT), Kerr–Newman-NUT (KN-NUT), and RN-NUT-AdS spacetimes in the framework of the (ψ-N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi -\mathscr {N}$$\end{document})-pair formalism, and then investigate their topological numbers by using the uniformly modified form of the generalized off-shell Helmholtz free energy. We find that these solutions can be included into one of three categories of those well-known black hole solutions, which implies that these spacetimes should be viewed as generic black holes from the perspective of the topological thermodynamic defects. In addition, we demonstrate that although the existence of the NUT charge parameter seems to have no impact on the topological number of the charged asymptotically locally flat spacetimes, it has a remarkable effect on the topological number of the charged asymptotically locally AdS spacetime.

In this paper, within the framework of the (ψ − N )pair formalism, we will first utilize the generalized Komar super-potential [38] to derive the consistent thermodynamics of the four-dimensional Lorentzian Reissner-Nordström-NUT (RN-NUT), Kerr-Newman-NUT (KN-NUT), and RN-NUT-AdS spacetimes, and to investigate their topological number via the uniformly modified form of the generalized off-shell Helmholtz free energy.We find that these spacetimes should also be viewed as generic black holes from the thermodynamic topological perspective.
The remaining part of this paper is organized as follows.In Sec.II, we give a brief review of the novel thermodynamic topological approach for Taub-NUT-type spacetimes.In Sec.III, we first derive the consistent formulation of thermodynamic properties of the four-dimensional Lorentzian RN-NUT spacetime and then investigate its topological number.In Sec.IV, we turn to discuss the case of the Lorentzian KN-NUT spacetime.In Sec.V, we then extend to discuss the more general Lorentzian RN-NUT-AdS 4 spacetime.Finally, we present our conclusions in Sec.VI.

II. THERMODYNAMIC TOPOLOGICAL APPROACH FOR TAUB-NUT-TYPE SPACETIMES
In accordance with the thermodynamic topological approach suggested in Ref. [37], it is possible to introduce the modified form of the generalized off-shell Helmholtz free energy for a Taub-NUT-type black hole thermodynamical system with the mass M, the entropy S, and the Misner potential ψ, as well as the gravitational Misner charge N , where τ is an extra variable that can be treated as the inverse temperature of the cavity surrounding the Taub-NUT-type black hole.Only when τ = T −1 , the modified form of the generalized Helmholtz free energy (1) is on-shell and reduces to the common Helmholtz free energy: F = M − T S − ψN of the Taub-NUT-type black holes [39][40][41].According to Ref. [24], a core vector φ is defined as where the two parameters satisfy the ranges: 0 < r h < +∞, 0 ≤ Θ ≤ π, respectively.The component φ Θ is divergent at Θ = 0 and Θ = π, indicating that the direction of the vector is outward there.One can define the topological current by using Duan's φmapping topological current theory [42][43][44] as follows: where Since it is easy to prove that the above current (3) is conserved, and one can quickly obtain ∂ µ j µ = 0 and then indicate that the topological current is a δ -function of the field configuration [43,44] where the three dimensional Jacobian J µ (φ /x) obeys: It is easy to demonstrate that j µ equals to zero only when φ a (x i ) = 0, and one can easily obtain the topological number W as follows: where β i is the positive Hopf index that counts the number of the loops of the vector φ a in the φ -space when x µ are around the zero point z i , while η i = sign(J 0 (φ /x) z i ) = ±1 is the Brouwer degree, and w i is the winding number for the i-th zero point of φ that is contained in the domain Σ.In addition, if two different closed curves Σ 1 and Σ 2 enclose the same zero point of φ , the corresponding winding number must equal.On the other hand, if there is no zero point of φ in the enclosed region, one must have W = 0. Note that the local winding number w i can be used to characterize the local thermodynamic stability, with positive and negative values corresponding to thermodynamically stable and unstable black holes, respectively, and the global topological number W represents the difference between the number of thermodynamically stable black holes and the number of thermodynamically unstable black holes of a classical black hole solution at a fixed temperature [24].Therefore, not only can one distinguish between different black hole phases (thermodynamically stable or unstable) of the same black hole solution at a specific temperature based upon the local winding number, but also can classify the black hole solutions based upon the global topological number.Furthermore, according to this classification, black holes with the same global topological number (even if they are of different geometric types) have similar thermodynamical topology properties.

III. RN-NUT SPACETIME
As the simplest charged case, we will investigate the fourdimensional Lorentzian RN-NUT solution [45][46][47], and adopt the following line element in which the Misner strings are symmetrically distributed along the polar axis: where in which m, n and q are the mass, the NUT charge, and the electric charge parameters, respectively.The event horizon radius: r h = m+ m 2 + n 2 − q 2 is the largest root of the equation: f (r h ) = 0.In addition, the electromagnetic gauge potential one-form is given by with which a convenient gauge choice is made so that its temporal component vanishes at infinity.The metric ( 6) with the Abelian gauge potential ( 7) is an exact solution to the field equations derived from the Lagrangian density: For latter convenience, we introduce a generalized Komar super-potential (see Eq. (5.20) of Ref. [38]) as follows: associated with a Killing vector ξ and the Faraday-Maxwell field strength tensor: It can be shown that ∇ b Ξ ab = 0.

A. Consistent thermodynamics
As we are only focused on the purely electrically charged RN-NUT solutions, so we will first rederive the consistent thermodynamics of the four-dimensional Lorentzian RN-NUT spacetime within the framework of the (ψ − N )-pair formalism.
The Bekenstein-Hawking entropy is one quarter of the area of the event horizon the Gibbons-Hawking temperature is proportional to the surface gravity κ on the event horizon in which a prime represents the partial derivative with respective to its variable.
Secondly, the total electric charge distribution over a twodimensional sphere with a finite radius r is given by the Gauss' integral which clearly shows it is radial-dependently distributed and should be on the Misner string singularities [48].So the electric charge on the event horizon is The corresponding electrostatic potential at the event horizon simply reads where χ = ∂ t is the timelike Killing vector normal to the event horizon.
In the Lorentzian RN-NUT spacetime, there are additional Killing horizons (north/south pole axes) connected to the Misner strings with the associated Misner potential being In the language of exterior differential forms, the Hodge dual two-form corresponding to the generalized Komar superpotential ( 8) can be obtained as sin θ dθ ∧ dφ , (14) for the timelike Killing vector χ = ∂ t .
Using the above generalized Komar superpotential twoform (14) to replace the ordinary Komar one and following the same pattern of the (ψ − N )-pair formalism as did in Refs.[40,49] (namely, one deliberately separates the integral into three parts: the spatial infinity, the horizon, and two Misner string tubes), one can derive the integral Bekenstein-Smarrlike mass formula and then verify that the differential first law can also be satisfied: In the derivation of the Smarr-like formula, one can define the conserved mass as which exactly coincides with the Komar mass calculated via the usual Komar integral at infinity, since the well fall-off asymptotic behavior of the Maxwell field.In the computation, we has used the determinant √ −g = r 2 + n 2 sin θ .On the hand hand, we have instead The last thermodynamic quantity of the Misner charge can be easily determined via the above Bekenstein-Smarr-like mass formula as which is non-globally conserved.Alternately, it can also be evaluated via the Misner tubes integral , which reproduces the above expression ( 19) after using the identity m = r 2 h − n 2 + q 2 /(2r h ).It can be further identified that all the above thermodynamic quantities are related to the Gibbs free energy of the fourdimensional Lorentzian RN-NUT spacetime [50] whose expression can be obtained via a Wick-rotated back procedure from the Euclidean action of the Euclidean RN-NUT spacetime: where h is the determinant of the induced metric h i j , K is the trace of the extrinsic curvature tensor defined on the boundary with this induced metric, and K 0 is the subtracted one of the massless uncharged Taub-NUT solution as the reference background.The computation of the Euclidean action integral yields the following expression for the Gibbs free energy where β = 1/T is the interval of the time coordinate.
By the way, it should be noted that the above results are completely consistent with those given in Ref. [51] without any "derivation".

B. Topological number
Next, we will investigate the topological number of the four-dimensional Lorentzian RN-NUT spacetime.We note that the Helmholtz free energy simply reads Replacing T with 1/τ in Eq. ( 23) and using m = r 2 h − n 2 + q 2 /(2r h ), then the generalized off-shell Helmholtz free energy is Adopting the definition of Eq. ( 2), the components of the vector φ can be easily calculated as follows: By solving the equation: φ r h = 0, one can arrive at a curve on the r h − τ plane.For the four-dimensional RN-NUT spacetime, one can obtain We point out that Eq. ( 26) consistently reduces to the one obtained in the case of the four-dimensional RN black hole [24] when the NUT charge parameter n is turned off.Note that the generation point satisfies the constraint conditions: Taking q/r 0 = 1 and n/r 0 = 1 for the four-dimensional Lorentzian RN-NUT spacetime, we plot in Fig. 1 and Fig. 2, respectively, for the zero points of the component φ r h , and for the unit vector field n on a portion of the Θ − r h plane with τ = 7r 0 in which r 0 is an arbitrary length scale set by the size of a cavity enclosing the RN-NUT spacetime.From Fig. 1, one generation point can be found at τ/r 0 = τ c /r 0 = 5.13.It is clear that the RN-NUT spacetime behaves like the RN black hole, showing that the NUT charge parameter appears to have no impact on the thermodynamic topological classification for the static charged asymptotically locally flat spacetime.Consequently, it would be fascinating to learn more about the relationship between geometric topology and thermodynamic topology: for example, it would be very interesting to investigate the topological number of ultraspinning black holes [52][53][54][55][56][57][58][59] and their usual counterparts.In Fig. 2, the zero points are located at (r h /r 0 , Θ) = (0.23, π/2), and (0.86, π/2), respectively.Thus, one can read the winding numbers w i for the blue contours C i : w 1 = 1, w 2 = −1, which are similar to those of the RN black hole [24].In terms of the topological global properties, one can easily obtain the topological number W = 0 for the RN-NUT spacetime from Fig. 2, which is also the same one as that of the RN black hole.As a result, based upon the viewpoint of the thermodynamic topological numbers, the Lorentzian RN-NUT spacetime should be welcomed into the black hole family.Furthermore, it can be indicated that, while the RN-NUT spacetime and RN black hole are evidently distinguished in geometric topology, they belong to the same class in terms of thermodynamic topology.

IV. KN-NUT SPACETIME
In this section, we will focus on the case of a rotating charged Taub-NUT spacetime by considering the fourdimensional KN-NUT solution [60][61][62][63], whose line element with the Misner strings symmetrically distributed along the rotation axis is written in the Boyer-Lindquist coordinates as: where in which m, n, a and q are the mass, the NUT charge, the rotation and the electric parameters, respectively.The event horizon radius is r h = m + m 2 + n 2 − a 2 − q 2 .In addition, the electromagnetic gauge potential one-form is given by in a gauge that its temporal component vanishes at infinity.

A. Consistent thermodynamics
Now, we investigate the consistent thermodynamics of the four-dimensional Lorentzian KN-NUT spacetime within the framework of the (ψ − N )-pair formalism.The Bekenstein-Hawking entropy is taken as one quarter of the event horizon area: The Hawking temperature is proportional to the surface gravity κ on the event horizon The angular velocity at the event horizon and the Misner potential are, respectively, The electric charge Q on the event horizon can be computed as and its corresponding electrostatic potential at the event horizon is where χ = ∂ t + Ω∂ ϕ is the Killing vector normal to the event horizon.
As for the conserved mass, one can compute it just like the non-rotating case and get for the timelike Killing vector ∂ t .One can note that it exactly coincides with the Komar mass evaluated via the usual Komar integral at infinity.
One can anticipate that both the first law and the Bekenstein-Smarr mass formula for the Lorentzian KN-NUT spacetime should read from which one can first solve J h in terms of N from the integral Smarr-like formula (37), and then solve N = N (r h , q, n, a) from the differential first law.After abandoning an integration constant, one can finally get the expressions for the gravitational Misner charge and the angular momentum as follows: We would like to point out that the above two expressions are identical to those of Ñ and J given by Eqs.(3.15) and (3.16) in Ref. [41] in the case the magnetic charge parameter is turned off, namely, the asymptotic magnetic charge vanishes.The first law and the Bekenstein-Smarr mass formula precisely correspond to the magnetic version of the full cohomogeneity first law [41] when the magnetic charge parameter is set to zero.On the other hand, one can utilize the generalized Komar superpotential (8) rather than the usual Komar one with respect to the Killing vector χ = ∂ t + Ω∂ ϕ and follow the same paradigm of the (ψ −N )-pair formalism as did in Ref. [64] to derive the integral Bekenstein-Smarr-like mass formula (37) and then check that the above thermodynamic quantities simultaneously satisfy the differential first law as well.In this way, it is facilitated to use the GRTensor II package to perform the algebraic manipulation to get the above involved expressions of J h and N .Here, we will not repeat the "derivation" but just provide a simple and equivalent way to evaluate the horizon angular momentum J h by using our generalized Koamr superpotential (8): which reproduces the above expression after using m = (r 2 h − n 2 + a 2 + q 2 /(2r h ).
Adopting the same procedure as did in Sec.III A, the calculation of the Euclidean action integral (21) of the KN-NUT spacetime yields the Gibbs free energy which coincides with the result of Eq. (3.1) given in Ref. [41] in the case when the magnetic charge parameter g = 0 is turned off.Furthermore, it can be also identified as

B. Topological number
In order to obtain the thermodynamic topological number of the KN-NUT spacetime, we need to get the expression of the generalized off-shell Helmholtz free energy in advance.The Helmholtz free energy is given by It is a simple matter to obtain the generalized off-shell Helmholtz free energy as Then, the components of the vector φ are given by where Therefore, by solving the equation: φ r h = 0, one can obtain as the zero point of the vector field φ r h .We also point out that Eq. ( 46) consistently reduces to the one obtained in the case of the four-dimensional Kerr-Newman black hole [27] when the NUT charge parameter n vanishes.
Taking n/r 0 = 1 and q/r 0 = 1 as well as n/r 0 = 1 for the KN-NUT spacetime, we plot the zero points of the component φ r h in Fig. 3, and the unit vector field n on a portion of the Θ − r h plane in Fig. 4 with τ/r 0 = 50, respectively.In Fig. 3, one generation point can be found at τ/r 0 = τ c /r 0 = 41.49.At τ = τ 1 , there are one thermodynamically unstable KN-NUT spacetime and one thermodynamically stable KN-NUT spacetime, just like the Kerr-Newman black hole [27].In Fig. 4, one can observe that the zero points are located at (r h /r 0 , Θ) = (1.70,π/2), and (3.41, π/2), respectively.Based upon the local property of the zero points, we can obtain the topological number: W = 1 − 1 = 0 for the KN-NUT spacetime, which is the same one as that of the Kerr-Newman black hole [27].Therefore, the four-dimensional KN-NUT spacetime should be present in the large family of black holes.Additionally, it can be concluded that even though the KN-NUT spacetime and Kerr-Newman black hole have undoubtedly distinct geometric topologies, they are the same type from the viewpoint of the thermodynamic topology, just like the RN-NUT spacetime and RN black hole, which have been addressed in Sec.III and Ref. [24], respectively.V. RN-NUT-ADS 4 SPACETIME In this section, we turn to explore the Lorentzian charged Taub-NUT spacetime with an negative cosmological constant, namely, the Lorentzian RN-NUT-AdS 4 spacetime, whose metric and Abelian gauge potential are still given by Eqs. ( 6)-( 7), but now f (r) = r 2 − 2mr − n 2 + q 2 + r 4 + 6n 2 r 2 − 3n 4 /l 2 , in which the AdS radius l is related to the thermodynamic pressure P = 3/ 8πl 2 of the four-dimensional AdS black hole [65][66][67].One can show that the generalized Komar superpotential now obey an identity: ∇ b Ξ ab = −6ξ a /l 2 in the present case.

A. Consistent thermodynamics
We now investigate the thermodynamical properties within the (ψ − N )-pair formalism of the four-dimensional Lorentzian RN-NUT-AdS spacetime.Since we are extending the results already appeared in the subsection III A, so we will mainly collect the needed expressions and just outline the different aspects.For the event horizon r h , which is the location of the largest root of the radial function: f (r h ) = 0, the Bekenstein-Hawking entropy is while the Hawking temperature has a different expression On the event horizon, the electric charge and and its corresponding electrostatic potential The conformal mass can be evaluated as which is associated with the timelike Killing vector: χ = ∂ t .The Misner potential is Using the metric determinant: √ −g = r 2 + n 2 sin θ , one can define the thermodynamic volume which is conjugate to the pressure: P = 3/ 8πl 2 .Within the framework of the extended phase space, one can substitute the above thermodynamical quantities into the Bekenstein-Smarr mass formula and use the identity: m = (r 2 h − n 2 + q 2 )/(2r h ) + r 4 h + 6n 2 r 2 h − 3n 4 /(2l 2 r h ) to acquire the expression of the gravitational Misner charge: Then one can verify that they also completely satisfy the first law: We would like to point out that the mass formulae presented here exactly correspond to the magnetic-type first law and Smarr-like mass formula of the unconstrained ψ −N pair formalism of the consistent thermodynamics of the dyonic RN-NUT-AdS 4 spacetimes [49,50] when the asymptotic magnetic charge is turned off.In particular, the above expression for the Misner charge (54) coincides with that of N (2) explicitly given by Eq. ( 57) in Ref. [49] after setting the magnetic charge parameter to zero.
One can also follow the same steps as did in Refs.[40,49] to derive the above Smarr-like formula.To do so, in addition to use the generalized Komar superpotential two-form (14), one must also introduce a dual Killing co-potential ω to cancel the divergence at infinity.We shall not repeat this algebraic excise here.Instead, there is another simple way to regulate the divergence by making a subtraction from the massless pure NUT-charged background.
One can obtain the Gibbs free energy [50] which coincides with those computed via the Euclidean action integral, namely G = I/β .In order to obtain this result, one can calculate the Euclidean action [50,68] for the Euclidean spacetime where K and R(h) are the extrinsic curvature and Ricci scalar of the boundary metric h µν , respectively.In order to remove the divergence, the action includes, in addition to the ordinary Einstein-Hilbert term, the Gibbons-Hawking boundary term and the corresponding AdS boundary counterterms [69][70][71][72][73].
Note that the Gibbs free energy (56) should also be identified with B. Topological number In the following, we will investigate the topological number of the four-dimensional Lorentzian RN-NUT-AdS spacetime.The Helmholtz free energy simply reads Replacing T with 1/τ and substituting l 2 = 3/(8πP), then the generalized off-shell Helmholtz free energy is given by (60) Thus, the components of the vector φ are obtained as follows: from which one can obtain the zero point of the vector field φ r h as which consistently reduces to the one obtained in the fourdimensional RN-AdS 4 black hole case [24] when the NUT charge parameter n is turned off.We also point out that the annihilation point satisfies the constraint conditions: Taking the pressure Pr 2 0 = 0.2 and the NUT charge parameter n/r 0 = 1 as well as the electric charge parameter q/r 0 = 1 for the RN-NUT-AdS 4 spacetime, we plot the zero points of φ r h in the r h − τ plane in Figs. 5 and the unit vector field n on a portion of the Θ − r h plane with τ = r 0 in Fig. 6, respectively.In Fig. 5, one annihilation point can be found at τ/r 0 = τ c /r 0 = 1.10.From Fig. 6, one can find that the zero points are located at (r h /r 0 , Θ) = (0.74, π/2), and (1.84, π/2), respectively.According to the conclusions in Ref. [15], it can be inferred that the second-order phase transition occurs in the RN-NUT-AdS 4 spacetime system.Thus, it is very interesting to explore the phase transitions of the RN-NUT-AdS 4 spacetime, such as the Hawking-Page phase transitions [74] and the P-V criticality [75] to check the correctness of the above conjecture.In addition, for the RN-NUT-AdS 4 spacetime, we observe that the topological number is: W = 0, and is different from that of the RN-AdS 4 black hole, which has: W = 1 [24], because the third zero point of the RN-AdS 4 black hole solution vanishes once the NUT charge parameter is introduced.Therefore, it indicates that the NUT charge parameter has a remarkable effect on the topological number for the static charged asymptotically local AdS spacetime.As a result, at least according to the viewpoint of the thermodynamic topological approach, the Lorentzian RN-NUT-AdS 4 spacetime should be included into a member of the black hole family.The red arrows represent the unit vector field n on a portion of the r h − Θ plane with q/r 0 = 1, n/r 0 = 1, Pr 2 0 = 0.2 and τ/r 0 = 1 for the RN-NUT-AdS 4 spacetime.The zero points (ZPs) marked with black dots are at (r h /r 0 , Θ) = (0.74, π/2), (1.84, π/2) for ZP 1 and ZP 2 , respectively.The blue contours C i are closed loops surrounding the zero points.

VI. CONCLUSIONS
Our results found in the present paper are now summarized in the following Table I.
In this paper, we first derive the consistent thermodynamics of the four-dimensional Lorentzian charged RN-NUT, KN-NUT, and RN-NUT-AdS spacetimes within the framework of  first law [41,49,50] of the dyonic NUT-charged spactimes when the magnetic charge parameter vanishes.Then we investigate their topological numbers by using the uniformly modified form of the generalized off-shell Helmholtz free energy.We found that the RN-NUT spacetime has: W = 0, which is the same one as that of the RN black hole [24].We showed that the KN-NUT spacetime has: W = 0, which is identical to that of the Kerr-Newman black hole [27].In addition, we also indicated that the RN-NUT-AdS 4 spacetime has: W = 0, which is different from that of the RN-AdS 4 black hole (W = 1) [24].Therefore, one can conclude that although the existence of the NUT charge parameter seems to have no impact on the topological number of the charged asymptotically locally flat spacetimes, it has an important effect on the topological number of the charged asymptotically locally AdS spacetime.Furthermore, it can be demonstrated that the four-dimensional RN-NUT, KN-NUT and RN-NUT-AdS spacetimes should be treated as generic black holes from the standpoint of the thermodynamic topological approach.There are two promising further topics that can be pursued in the future.As mentioned above, one intriguing topic is to explore the phase transitions of the RN-NUT-AdS 4 spacetime.Another one is to extend the present work to the more general dyonic cases [76] and higher-even dimensional cases [77,78].

FIG. 1 .
FIG.1.Zero points of the vector φ r h shown in the r h − τ plane with n/r 0 = 1 and q/r 0 = 1.The generation point for the RN-NUT spacetime is represented by the black dot with τ c .At τ = τ 1 , there are two RN-NUT spacetimes.

FIG. 3 .
FIG.3.Zero points of the vector φ r h shown in the r h − τ plane with n/r 0 = 1, q/r 0 = 1 and a/r 0 = 1.The generation point for the KN-NUT spacetime is represented by the black dot with τ c .At τ = τ 1 , there are two KN-NUT spacetimes.Obviously, the topological number is: W = 1 − 1 = 0.

Solutions
− N )-pair formalism, which exactly correspond to the magnetic version of the full cohomogeneity (unconstrained)

TABLE I .
The topological number W , numbers of generation and annihilation points for the four-dimensional Lorentzian charged Taub-NUT spacetimes.