New Taub-NUT Black Holes with Massive Spin-2 Hair

Introduction.- In the beginning of last century, Einstein discovered General Relativity (GR), which has two important predictions. One is the black hole and the other is gravitational wave. In this century, people have directly detected gravitational waves LIGOScientific:2017vwq ; LIGOScientific:2017zic and taken images of black holesEventHorizonTelescope:2019dse ; EventHorizonTelescope:2019uob ; EventHorizonTelescope:2019jan ; EventHorizonTelescope:2019ths ; EventHorizonTelescope:2019pgp ; EventHorizonTelescope:2019ggy . These great successes validate GR as a classical theory. However, GR is not a perfect theory; it is not renormalizable and can have have an infinite possibility of higher order corrections tHooft:1974toh . It was pointed out that one can obtain a renormalizable theory by including also the quadratic curvature terms to the Einstein-Hilbert action. The price to pay is that these higher-derivative terms that make the theory renormalizable will inevitably introduce ghost modes Stelle:1976gc . Many works have been done on addressing the ghost phenomena of higher derivative gravity Horava:2009uw ; Lu:2009em ; Li:2008dq ; Lu:2010ct ; Lu:2010cg ; Lu:2011zk and it was pointed out that the ghost problem might not be severe to the theory Smilga:2013vba .

There are three quadratic curvature terms in general, Ricci scalar square, Ricci tensor square and Riemann tensor square($R^{2},R_{\mu\nu}R^{\mu\nu},R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}$). We can regroup them as the Gauss-Bonnet combination $R^{2}+4R_{\mu\nu}R^{\mu\nu}+R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}$, Weyl-squared term $C^{\mu\nu\rho\sigma}C_{\mu\nu\rho\sigma}$ and $R^{2}$. In four dimensions, the Gauss-Bonnet term is a total derivative and hence doesn’t contribute to the equations of motion. The remaining Weyl-squared term and $R^{2}$ will introduce massive spin-2 and massive scalar modes respectively. The absence of the massive scalar mode can lead to critical gravityLu:2011zk . One important property of the quadratic curvature extension in four dimensions is that Ricci-flat metrics continue to be the solutions of the extended theory. This raises a question: does the theory admit new black holes beyond the Schwarzschild black hole? It turns out that the Ricci scalar must be zero owing to the application of the no scalar theorem in higher-derivative gravity Lu:2015cqa . On the other hand, black holes beyond the Schwarzschild black hole can carry massive spin-2 hair and they were constructed in Lu:2015cqa .

Taub-NUT spacetime is another one of the simplest Ricci-flat solutions in Einstein gravity which was constructed in 1960s Taub:1950ez ; Newman:1963yy . Compared to the famous Schwarszchild black hole, Taub-NUT spacetime has an additional integration constant $n$, which is called NUT parameter. Due to this NUT parameter, it has a string-like singularity along the polar axis, which is now known as Misner string singularityMisner:1963fr . Surprisingly, the Taub-NUT metric has no curvature singularity though it possesses an event horizon and two other Killing horizons. Those peculiar properties of Taub-NUT spacetime have been attracting many attentions for years Clement:2015cxa ; Clement:2015aka ; Hennigar:2019ive ; Bordo:2019slw ; Durka:2019ajz ; Awad:2022jgn ; Wu:2019pzr ; Rodriguez:2021hks ; Liu:2022wku ; Liu:2023uqf ; Jiang:2019yzs ; Chen:2023eio ; Siahaan:2022jrl ; Yang:2023hll .

It is clear the Ricci-flat Taub-NUT black hole will continue to be the solution in the extended gravity theory. In this paper, we would like to investigate whether new black holes carrying the NUT parameter beyond the Ricci-flat solution can arise. As in the case of Lu:2015cqa , the no-hair theory will rule out the possibility of having the hair associated with the massive scalar mode. We focus on the construction involving the massive spin-2 mode. The general solution should contain two parameters $(m,n)$. In the Ricci-flat case, any $(m,n\neq 0)$ can lead to an event horizon, but this is no longer true when higher-derivative terms are involved. The parameters $(m,n)$ has to be fine tuned in order to form an event horizon. This fine-tuning leads to difficulties in constructing the numerical solutions, since we have to start at some fiducial point of the parameter space.

Since the massive spin-2 mode of the quadratic extension is relevant for the new solution, we perform a spin-2 tensor mode analysis upon the Ricci-flat Taub-NUT background, which is called Lichnerowicz mode analysis. We find that there exists a Licnerowicz mode for some specific values of $(m,n)$ with a negative eigenvalue, which implies that the Taub-NUT background of this $(m,n)$ is not stable and a more stable and new black hole solution should emerge at this point of the parameter space. With this guidance we construct new Taub-NUT like black hole solutions carrying spin-2 hair in the theory of Einstein gravity added with quadratic curvature terms.

Theroy.- We consider Einstein gravity extended with the most general quadratic terms in four dimensional spacetime

where $C_{\mu\nu\rho\sigma}$ is Weyl tensor and $\kappa,~{}\alpha$ and $\beta$ are constants. For simplicity we shall set $\kappa=1$. The equations of motion can be obtained through the variation of the metric and written in a simple form

where $B_{\mu\nu}=(\nabla^{\rho}\nabla^{\sigma}+\frac{1}{2}R^{\rho\sigma})C_{\mu\nu\rho\sigma}$ is the trace free Bach tensor.

It is worth pointing out that we don’t use the Gauss-Bonnet terms to get rid off the Riemann square but to group the quadratic curvature terms into Weyl square. The advantage is that the Einstein gravity theory with quadratic curvature terms contains three modes in the linear spectrum, while the three terms in the theory (1) correspond to the three modes respectively, a usual massless spin-2 graviton, a massive spin-2 mode with mass $m_{2}=1/\sqrt{2\alpha}$ and a massive scalar mode with mass $m_{0}=1/\sqrt{6\beta}$. It was proved that the scalar hair is forbidden for static black hole solutionLu:2015cqa , since the trace part of the equation of motion leads to the vanishing of the Ricci scalar, $R=0$. However, the massive spin-2 mode can survive, and the theory can admit black hole solutions carrying spin-2 hair . In this paper, we want to explore the situation in Taub-NUT spacetimes.

The Taub-NUT solution was constructed in pure Einstein gravity,

with

The solution is a generalization of the famous Schwarzschild black hole solution. In add to the the mass parameter $m$, the Taub-NUT has one more integration constant $n$ which is usually called NUT parameter. Setting $n=0$, the Taub-NUT solution reduces to Schwarzschild solution.

As is mensioned, the Taub-NUT solution is naturally a solution of theory (1), since the solution is Ricci flat, $R_{\mu\nu}=0$. It is easily to verify that the Ricci scalar also vanishes under the Taub-NUT like spacetime (3) with arbitrary $h(r)$ and $f(r)$ satisfying condition of $h$ and $f$ vanish on the event horizon. The trace of the equations of motion (2) is $\beta(\Box-m_{0}^{2})R=0$. Multiplying $R$ and integrating over the horizon to infinity gives $\int dx^{4}\sqrt{g}R(\Box-m_{0}^{2})R=\int dx^{4}\Big{[}\partial_{r}(\sqrt{g}g^{rr}R\partial_{r}R)-\sqrt{g}(g^{rr}\partial_{r}R\partial_{r}R+m_{0}^{2}R^{2})\Big{]}=0$. Since $g^{rr}$ vanishes on the horizon and $\partial_{r}R$ goes to 0 at spatial infinity, the total derivative term gives no contribution, and the remaining two terms are non positive, thus the Ricci scalar has to be 0. Any Taub-NUT like solution of metric ansatz (3), with undetermined $h$ and $f$, in the theory (1) will have a vanishing Ricci scalar. This property of vanishing Ricci scalar makes a great simplification to the equations of motion.

Lichnerowicz Mode.- Since Taub-NUT black hole is a solution of the theory (1), we want to analyze the linear exited modes under the Taub-NUT background. The Taub-NUT spacetime is locally flat, $R_{\mu\nu}=0$, making a perturbation of Ricci tensor, we find

The trace of above equation is$(6\beta\Box-1)\delta R=0$. With the same analysis in the last section, we can deduce that the $\delta R$ must vanish. It means that any Taub-NUT like black hole which is perturbatively around Taub-NUT spacetime must have a vanishing Ricci scalar at the linearised level. Substituting $\delta R=0$ back to equation (5) gives

where

is called Lichnerowicz operator. For Ricci flat spacetime($R_{\mu\nu}=0$) such as Taub-NUT black hole, firstly, we can deduce that $g^{\mu\nu}\delta R_{\mu\nu}=0$ from $\delta R=0$. Secondly, taking variation of the Bianchi identity $\nabla^{\mu}R_{\mu\nu}-(1/2)\nabla_{\nu}R=0$, we can have $\nabla^{\mu}\delta R_{\mu\nu}=0$. Thus, $\delta R_{\mu\nu}$ is a transverse traceless tensor on the Ricci flat background.

This implies that once the Lichnerowicz operator has a transverse traceless eigenfunction $\psi_{\mu\nu}$ with eigenvalue $\lambda=-1/(2\alpha)$, $\Delta_{L}\psi_{\mu\nu}=\lambda\psi_{\mu\nu}$, there exists a linearized perturbation away from the Taub-NUT solution.

In the large $r$ limit, the Lichnerowicz equations can be solved and two Yukawa like modes, $e^{\pm\sqrt{-\lambda r}}$, emerge, which is under expectation. In general, the Lichnerowicz equation can not be solved analytically, and to solve the equation numerically, one need to use shooting method to seek proper values of the parameters to elude the diverging mode $e^{+\sqrt{\lambda r}}$, left only the normalisable mode $e^{-\sqrt{\lambda r}}$. Without loss of generality, we set $\alpha=1/2$, and find that there exists one mode corresponding to each value of $n$, and relation between the starting point $r_{0}$ of the linear perburbation and $n$ are shown in Fig 1. The technical details can be found in the Appendix. We can see that when $n=0$, $r_{0}=0.876$, which covers the result for Schwarszchild case Lu:2015cqa .

New Taub-NUT Black Hole Solutions.- The Lichnerowicz mode analysis shows that a spin-2 perturbation exists under the analytic Taub-NUT background, which implies that there must exist black hole solutions carrying spin-2 hair. At this stage, we carry out a numerical search of Taub-NUT like black hole solutions which carry spin-$2$ hair in the theory (1). We start with Taub-NUT like metric ansatz (3) and suppose there exists a black hole with event horizon located at radius $r_{0}$, where metric functions $h$ and $f$ both vanish. Then, we take a Taylor expansion around the horizon for $h$ and $f$,

Substituting these back into the equations of motion, the coefficients $h_{i}$ and $f_{i}$ for $i>1$ can be solved and expressed in terms of $f_{1}\,,h_{1}$ and $r_{0}$, the first two coefficients are

We use these near-horizon expansion to set initial data and then integrate the solution from horizon to large $r$. Since Taub-NUT solution corresponds to $f_{1}=1/r_{0}=h_{1}$, we set $f_{1}=\frac{1}{r_{0}}+\delta,~{}h_{1}=\frac{1}{r_{0}}$. And $\delta$ characterizes the deviation of our newly constructed numerical solution from the analytical Taub-NUT solution.

We set the theory parameter $\alpha=1/2$ as we did in the last section, and $n=1/2$, then start the investigation of new solution around $r_{0}=0.852$ as suggested by the result of the Lichnerowicz mode analysis. We find a newly Taub-NUT like black holes carrying spin-2 hair under expectation. We show the metric functions $h$ and $f$ in the Fig.2, it is obvious that the two metric functions $h$ and $f$ are different, unlike that of the analytical Taub-NUT solution in which $h$ and $f$ are the same. The profile of metric functions are qualitatively different for different horizon radius. The Fig.2 shows the plot of metric function with $r_{0}=1$ and $r_{0}=1.5$, corresponding to black hole with a positive mass parameter and black hole with a negative mass parameter, respectively. We read of the mass parameter $m$ from the fall off of metric function $f\sim 1-2m/r+\dots$ in the large $r$. As is for analytic Taub-NUT black hole solution, there is no objection for the mass parameter taking negative value, and the mass of the black hole can have contributions from both mass parameter $m$ and the NUT parameter $n$ Liu:2022wku ; Liu:2023uqf . We also show the pattern of the profile of metric function with negative or positive mass parameter for various NUT parameter $n$, whilst keep $r_{0}$ fixed in Fig.3.

In fact, this pattern also emerges in the analytical Taub-NUT solution, where one can easily derive that the condition for critical value of mass parameter $m=0$ is $r_{0}=n$. Here, for our numerical solution, we can not give an analytical relation of ($r_{0}\,,n$) for $m=0$, but we can plot the curve of the relation ($r_{0}\,,n$) in the $r_{0}-n$ space, which is displayed in Fig.4. We also plot the cure $r_{0}=n$ of analytic Taub-NUT solution as a comparison, the critical curve of the newly found Taub-NUT like solution crosses with the critical curve of analytical Taub-NUT at $n=0.83$. It is worth pointing out that when $n=0$, the critical value of $m=0$ emerges at $r_{0}=1.143$, which is consistent with the result in Lu:2015cqa .

We also plot the mass parameter as a function of horizon radius in Fig.5. The red line is the new found Taub-NUT like black holes, while the blue line represents the analytic Taub-NUT black hole. One can see that the mass parameter decreases as the event horizon $r_{0}$ increases. For $n=0.5$, the mass parameter can be negative when the black hole is larger enough. However, for larger n (like n=1.5), the mass parameter is always negative for the newly found Taub-NUT like black holes. It is consistent with Fig.2, mass parameter decreases as NUT parameter grows. Compared with the Schwarzschild case in Lu:2015cqa , our new Taub-NUT solutions have two different properties. The first is that the new solutions have two branches and the smooth connected. The other is that, for fixed NUT parameter n, there is a minimal value of event horizon radius,$r_{min}=0.721$ for $n=0.5$ and $r_{min}=0.67$ for $n=1.5$.

We now turn to the thermodynamical properties of the newly found Taub-NUT black hole solution. The Hawking temperature can be calculated through the standard method

The temperature of the solution as a function of event horizon radius is shown in Fig.6 for $n=0.5$ and $n=1.5$. As a comparison, we also plot the temperature cure of Taub-NUT solution, the two curves cross at $r_{0}=0.852$ for $n=0.5$ and at $r_{0}=0.716$ for $n=1.5$. It is consistent with the Lichnerowicz mode analysis in the previous section. There is a maximum value of temperature, which is more obvious for $n=1.5$. There can exist two black holes in the same temperature. Furthermore, with $r_{0}$ decreasing, a notable inflection point will occur at $r_{0}=0.721$ in the $n=0.5$ case. As is seen in Fig.6, the two branches of the new solution are smoothly connected in the minimal radius. And the temperature of one branch is always higher than the other one.

The entropy of the newly found black hole is not equal to one quarter of the event horizon area due to the higher derivative terms. According to Wald’s entropy formula Wald:1993nt ; Iyer:1994ys , the expression of entropy with higher curvature correction is

As is mentioned in the beginning of this paper, the Gauss-Bonnet terms is a total derivative term in four dimensional spacetime and do not contribute to the equations of motion. However, the Gauss-Bonnet terms can have a constant contribution to the entropy. We use this freedom to make the entropy is equal to one quarter of the event horizon area for analytic Taub-NUT black hole solution. We plot the $S-r_{0}$ curve for $n=0.5$ and $n=1.5$ in Fig.6, the solid line is the cure for analytic Taub-NUT solution. Here, the two branches of the new solution meet at the minimal radius, and the entropy of one branch is always larger than the other branch. However, it is worth pointing out that the branch with lower temperature has larger entropy, and vice versa.

Conculsion.- We constructed new Taub-NUT-like black hole solutions in Einstein gravity plus quadratic curvature terms. A Lichnerowicz mode analysis was carried out on the Ricci-flat Taub-NUT background, which gives us values of $(m,n)$ where there could exist new black holes carrying the massive spin-2 hair. Then new Taub-NUT like black hole solutions which carry spin-2 hair were constructed numerically, enlarging the black hole family in higher derivative gravity theories. The newly found solutions bifurcate from the Ricci-flat Taub-NUT solutions. Compared with the analytic Taub-NUT solutions, the temperature of one branch new solution has a maximum for a large $n$, and there exists two black holes with the same temperature. Since the Ricci-flat Taub-NUT is a solution, there can be a total of three black holes carrying NUT parameter for an appropriate given temperature. This implies that there could be phase transitions between these black holes and the phase diagram becomes much richer.

As pointed out in Liu:2022wku the mass of the analytic Taub-NUT black hole has extra contribution from the Minser strings which makes the total mass of the spacetime non-negative, though the mass parameter of the Taub-NUT solutions can be negative. The situation here is similar, the mass parameter can be negative, but negativity of the total mass could be avoided due to the extra contribution from the polar strings. However, the spin-2 modes are ghost like, the solution may still suffer from negative mass problem. The definition of mass is a difficult problem, but is very crucial at the same time, and so is the NUT charge. Both the definitions of mass and NUT charge are worth further exploring.

Acknowledgement.- We thanks Ze Li and Hong Lu for useful discussion. We are also grateful for Ling-Yu Zhang’s help in figures. This work is supported in part by NSFC (National Natural Science Foundation of China) Grants No. 12075166.