This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

thanks: Corresponding author, email: dingchikun@163.com; dck@hhtc.edu.cn

High dimensional AdS-like black hole and phase transition in Einstein-bumblebee gravity

Chikun Ding1,2,3    Yu Shi1    Jun Chen1    Yuebing Zhou1    Changqing Liu1 1Department of Physics, Huaihua University, Huaihua, 418008, P. R. China
2Department of Physics, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
3Key Laboratory of Low Dimensional Quantum Structures and Quantum Control of Ministry of Education, and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha, Hunan 410081, P. R. China
Abstract

Abstract

In this paper we obtain an exact high dimensional anti-de Sitter (AdS) black hole solution in Einstein-bumblebee gravity theory. This AdS-like black hole can only exist with a linear functional potential of the bumblebee field. We find that the Smarr formula and the first law of black hole thermodynamics can still be constructed in this Lorentz symmetry breaking black hole spacetime as long as its temperature, entropy and volume are slightly modified. We find also that there exist two kinds of phase transition: small-large black hole phase transition and Hawking-Page phase transition, like those of Schwarzschild AdS black hole. After Lorentz symmetry breaking, the black hole mass at divergent point of heat capacity becomes small, and the Gibbs free energy of the meta-stable large black hole is also smaller, showing that the large stable black hole can be more easily formed.

pacs:
04.50.Kd, 04.20.Jb, 04.70.Dy

I Introduction

Local Lorentz invariance (LI) is the most fundamental spacetime symmetry of the special/general relativity (SR/GR) and the standard model (SM) of particle physics ashtekar . LI means that the equations of motion for particles and fields are invariant under all Lorentz transformations(global for SR, local for GR) and diffeomorphism transformation for a particle. However, LI may fail to be an exact principle at very high energy scales mattingly . GR and SM are the best theories describing the four fundamental forces and there are no known experimental conflicts. Nonetheless, they are fundamentally different in that SM is a quantum theory disregarding the gravitational effects of all particles, while GR is a classical geometrical theory ignoring all the quantum features of particles. For particles with energies on the order of 103010^{30} eV (above the Planck scale), the gravitational interactions estimated by GR are very powerful, and so one shouldn’t ignore it camelia . Thus, on the very high energy scales, one has to reconsider combining the SM with GR in a unified theory, i.e., ”quantum gravity”. An effective field theory combining GR and SM on low energy scales is called standard model extension (SME), which couples the SM to GR, and involves extra items embracing information about the violations of Lorentz invariance(Lorentz violation, LV) happening on the Planck scale kostelecky2004 .

Studying the Lorentz violation (LV) is a useful approach toward investigating the foundations of modern physics. Besides SME, string theory kostelecky198939 , noncommutative field theories carroll , spacetime-varying fields bertolami69 , loop quantum gravity theory gambini , brane world scenarios burgess03 , massive gravity fernando and Einstein-aether theory jacobson are other proposals of Lorentz violation. The SME is experimentally accessible for studying possible observable signals of violation the particle Lorentz symmetry. In this model, the local Lorentz violating is provided by a spontaneous symmetry breaking potential due to self-interacting tensor fields, whose vacuum expectation value (VEV) yields to background tensor fields. A specific theory is the bumblebee field BμB^{\mu}, an self-interacting tensor field has a nonzero VEV bμb^{\mu} which defines a privileged direction in spacetime and spontaneously breaks the local Lorentz symmetry.

There are many forms of the bumblebee potential. One is the smooth functional bumblebee potential V(X)=kX2/2V(X)=kX^{2}/2 which has a minimum when the bumblebee field equals its VEV, i.e., V=0V=0 and V(X)=kX=0V^{\prime}(X)=kX=0 if X=0X=0, where kk is a constant bluhm and X=BμBμ±b2X=B^{\mu}B_{\mu}\pm b^{2}. With this potential, the four dimensional static black hole solutions were found by Bertolami et al. bertolami and by Casana et al. casana ; the higher DD-dimensional static black hole was found by Ding et al. ding2022 ; the rotating black hole solutions were built by Ding et al. ding2020 and by Jha et al. jha . An other form of bumblebee potential is a linear multiplier function V(X)=λX/2V(X)=\lambda X/2, where λ\lambda is a Lagrange-multiplier. With this potential, the Schwarzschild-anti-de Sitter-like and Schwarzschild-de Sitter-like black hole solutions were found by Maluf et al. maluf , wherein the radial and tangential pressures are different, like as an anisotropic fluid.

The quantum field theory in a spacetime with a horizon exhibits thermal behavior which must be thought of as an emergent phenomena padmanabhan . A black hole event horizon has the conceptions of temperature and entropy hawing1975 , and the first law of thermodynamics dU=TdSPdVdU=TdS-PdV can be built, where PP is the radial pressure of the matter, UU is the internal energy and VV is the volume of the black hole. There also has transfer of heat energy and the conception of phase transition hawking . Then the gravity system can be mapped to a thermodynamic system and studying on black hole thermodynamics can help to construct a quantum theory of gravity.

We will study the high dimensional black hole solutions and its thermodynamics and phase transition with a negative cosmological constant in the theory of Einstein-bumblebee gravity. The rest paper is organized as follows. In Sec. II we give the background for the Einstein-bumblebee theory and derive the black hole solution by solving the gravitational field equations. In Sec. III, we study its thermodynamical properties. In Sec. IV, we study the phase transitions and find some effects of the Lorentz breaking constant \ell. Sec. V is for a summary.

II High dimensional AdS-like Black hole in Einstein-bumblebee gravity

In the bumblebee gravity model, one introduces the bumblebee vector field BμB_{\mu} which has a nonzero vacuum expectation value, leading to a spontaneous Lorentz symmetry breaking in the gravitational sector via a given potential. In the higher D4D\geq 4 dimensional spacetime, the action of Einstein-bumblebee gravity is ding2021 ,

𝒮=dDxg[R2Λ2κ+ϱ2κBμBνRμν14BμνBμνV(BμBμb2)+M],\displaystyle\mathcal{S}=\int d^{D}x\sqrt{-g}\Big{[}\frac{R-2\Lambda}{2\kappa}+\frac{\varrho}{2\kappa}B^{\mu}B^{\nu}R_{\mu\nu}-\frac{1}{4}B^{\mu\nu}B_{\mu\nu}-V(B_{\mu}B^{\mu}\mp b^{2})+\mathcal{L}_{M}\Big{]}, (1)

where RR is Ricci scalar, Λ\Lambda is an negative cosmological constant. κ=8πGD/c4\kappa=8\pi G_{D}/c^{4}, where GDG_{D} is the DD-dimensional gravitational constant and has the relation to Newton’s constant GG as GD=GΩD2/4πG_{D}=G\Omega_{D-2}/4\pi boulware , ΩD2=2πD1/Γ[(D1)/2]\Omega_{D-2}=2\sqrt{\pi}^{D-1}/\Gamma[(D-1)/2] is the area of a unit (D2)(D-2)-sphere. Here and hereafter, we take GD=1G_{D}=1 and c=1c=1 for convenience.

The coupling constant ϱ\varrho dominates the non-minimal gravity interaction to bumblebee field BμB_{\mu}. The term M\mathcal{L}_{M} represents possible interactions with matter or external currents. The potential V(BμBμb2)V(B_{\mu}B^{\mu}\mp b^{2}) triggers Lorentz and/or CPTCPT (charge, parity and time) violation, where the field BμB_{\mu} acquires a nonzero VEV, Bμ=bμ\langle B^{\mu}\rangle=b^{\mu}, satisfying the condition BμBμ±b2=0B^{\mu}B_{\mu}\pm b^{2}=0. The constant bb is a real positive constant. Another vector bμb^{\mu} is a function of the spacetime coordinates and has a constant value bμbμ=b2b_{\mu}b^{\mu}=\mp b^{2}, where ±\pm signs mean that bμb^{\mu} is timelike or spacelike, respectively. It gives a nonzero vacuum expectation value (VEV) for bumblebee field BμB_{\mu} indicating that the vacuum of this model obtains a prior direction in the spacetime. The bumblebee field strength is

Bμν=μBννBμ.\displaystyle B_{\mu\nu}=\partial_{\mu}B_{\nu}-\partial_{\nu}B_{\mu}. (2)

This antisymmetry of BμνB_{\mu\nu} implies the constraint bluhm

μνBμν=0.\displaystyle\nabla^{\mu}\nabla^{\nu}B_{\mu\nu}=0. (3)

Varying the action (1) with respect to the metric yields the gravitational field equations

Gμν+Λgμν=κTμνB+κTμνM,\displaystyle G_{\mu\nu}+\Lambda g_{\mu\nu}=\kappa T_{\mu\nu}^{B}+\kappa T_{\mu\nu}^{M}, (4)

where Gμν=RμνgμνR/2G_{\mu\nu}=R_{\mu\nu}-g_{\mu\nu}R/2, and the bumblebee energy momentum tensor TμνBT_{\mu\nu}^{B} is

TμνB=BμαBνα14gμνBαβBαβgμνV+2BμBνV\displaystyle T_{\mu\nu}^{B}=B_{\mu\alpha}B^{\alpha}_{\;\nu}-\frac{1}{4}g_{\mu\nu}B^{\alpha\beta}B_{\alpha\beta}-g_{\mu\nu}V+2B_{\mu}B_{\nu}V^{\prime}
+ϱκ[12gμνBαBβRαβBμBαRανBνBαRαμ\displaystyle+\frac{\varrho}{\kappa}\Big{[}\frac{1}{2}g_{\mu\nu}B^{\alpha}B^{\beta}R_{\alpha\beta}-B_{\mu}B^{\alpha}R_{\alpha\nu}-B_{\nu}B^{\alpha}R_{\alpha\mu}
+12αμ(BαBν)+12αν(BαBμ)122(BμBν)12gμναβ(BαBβ)].\displaystyle+\frac{1}{2}\nabla_{\alpha}\nabla_{\mu}(B^{\alpha}B_{\nu})+\frac{1}{2}\nabla_{\alpha}\nabla_{\nu}(B^{\alpha}B_{\mu})-\frac{1}{2}\nabla^{2}(B^{\mu}B_{\nu})-\frac{1}{2}g_{\mu\nu}\nabla_{\alpha}\nabla_{\beta}(B^{\alpha}B^{\beta})\Big{]}. (5)

The prime denotes differentiation with respect to the argument,

V=V(x)x|x=BμBμ±b2.\displaystyle V^{\prime}=\frac{\partial V(x)}{\partial x}\Big{|}_{x=B^{\mu}B_{\mu}\pm b^{2}}. (6)

Varying instead with respect to the the bumblebee field generates the bumblebee equations of motion (supposing that there is no coupling between the bumblebee field and M\mathcal{L}_{M}),

μBμν=2VBνϱκBμRμν.\displaystyle\nabla^{\mu}B_{\mu\nu}=2V^{\prime}B_{\nu}-\frac{\varrho}{\kappa}B^{\mu}R_{\mu\nu}. (7)

The contracted Bianchi identities (μGμν=0\nabla^{\mu}G_{\mu\nu}=0) lead to conservation of the total energy-momentum tensor

μTμν=μ(TμνB+TμνMΛgμν)=0.\displaystyle\nabla^{\mu}T_{\mu\nu}=\nabla^{\mu}\big{(}T^{B}_{\mu\nu}+T^{M}_{\mu\nu}-\Lambda g_{\mu\nu}\big{)}=0. (8)

We suppose that there is no matter field and the bumblebee field is frozen at its VEV like in Refs casana ; bertolami , i.e., it is

Bμ=bμ.\displaystyle B_{\mu}=b_{\mu}. (9)

For the forms of potential VV, there are both classes: the smooth functionals of xx, i.e., κx2/2\kappa x^{2}/2, that minimized by the condition x=0x=0 and, the Lagrange-multiplier form λx/2\lambda x/2, where λ\lambda is a non-zero constant and deserves as a Lagrange-multiplier field bluhm . The smooth potential does not admit new black hole solution with nonzero cosmological constant Λ\Lambda maluf , one can also see the proof in the appendix for more details.

In order to generate a black hole solution in the case of non-zero cosmological constant Λ\Lambda, we should apply a linear function form of the Lagrange-multiplier potential

V=λ2(BμBμb2).\displaystyle V=\frac{\lambda}{2}(B_{\mu}B^{\mu}-b^{2}). (10)

This potential is V=0V=0 under the condition (9) and its derivative is V=λ/2V^{\prime}=\lambda/2 which can modify the Einstein equations. However, this additional degree of freedom of the λ\lambda field is auxiliary. Then the first two terms in Eq. (II) are like those of the electromagnetic field, the only distinctiveness are the coupling items to Ricci tensor and this λ\lambda term. Under this condition, Eq. (4) leads to gravitational field equations ding2021

Gμν+Λgμν=κ(λbμbν+bμαbνα14gμνbαβbαβ)+ϱ(12gμνbαbβRαβbμbαRανbνbαRαμ)+B¯μν,\displaystyle G_{\mu\nu}+\Lambda g_{\mu\nu}=\kappa(\lambda b_{\mu}b_{\nu}+b_{\mu\alpha}b^{\alpha}_{\;\nu}-\frac{1}{4}g_{\mu\nu}b^{\alpha\beta}b_{\alpha\beta})+\varrho\Big{(}\frac{1}{2}g_{\mu\nu}b^{\alpha}b^{\beta}R_{\alpha\beta}-b_{\mu}b^{\alpha}R_{\alpha\nu}-b_{\nu}b^{\alpha}R_{\alpha\mu}\Big{)}+\bar{B}_{\mu\nu}, (11)

with

B¯μν=ϱ2[αμ(bαbν)+αν(bαbμ)2(bμbν)gμναβ(bαbβ)].\displaystyle\bar{B}_{\mu\nu}=\frac{\varrho}{2}\Big{[}\nabla_{\alpha}\nabla_{\mu}(b^{\alpha}b_{\nu})+\nabla_{\alpha}\nabla_{\nu}(b^{\alpha}b_{\mu})-\nabla^{2}(b_{\mu}b_{\nu})-g_{\mu\nu}\nabla_{\alpha}\nabla_{\beta}(b^{\alpha}b^{\beta})\Big{]}. (12)

The static spherically symmetric black hole metric in the DD dimensional spacetime has the form

ds2=e2ϕ(r)dt2+e2ψ(r)dr2+r2dΩD22,\displaystyle ds^{2}=-e^{2\phi(r)}dt^{2}+e^{2\psi(r)}dr^{2}+r^{2}d\Omega_{D-2}^{2}, (13)

where ΩD2\Omega_{D-2} is a standard D2D-2 sphere. In this static spherically symmetric spacetime, the most general form for the bumblebee field would be bμ=(bt,br,0,0)b_{\mu}=(b_{t},b_{r},0,0), where btb_{t} and brb_{r} are functions of rr subject to the constraint bt2e2ϕ+br2e2ψ=b2-b_{t}^{2}e^{-2\phi}+b_{r}^{2}e^{-2\psi}=b^{2}, here bb is a positive constant. In this general case, the bumblebee field has both radial and a time component for the vacuum expectation value. In the purely radial case bt=0b_{t}=0, the authors in Refs. bertolami ; casana obtained new black hole solutions indeed. But for the general case(temporal and radial), the authors in Ref. bertolami obtain a slightly perturbed metric, where one cannot constrain the physical parameters from the observed limits on the PPN(parameterized post-Newtonian) parameters. Hence here we consider only the purely radial case to get a black hole solution and let the general case for the future work.

We pay attention to that the bumblebee field has a radial vacuum energy expectation because that the spacetime curvature has a strong radial variation, on the contrary that the temporal changes are very slow. Now the bumblebee field is supposed to be spacelike as that

bμ=(0,beψ(r),0,0,,0).\displaystyle b_{\mu}=\big{(}0,be^{\psi(r)},0,0,\cdots,0\big{)}. (14)

Then the bumblebee field strength is

bμν=μbννbμ,\displaystyle b_{\mu\nu}=\partial_{\mu}b_{\nu}-\partial_{\nu}b_{\mu}, (15)

whose components are all zero. And their divergences are all zero, i.e.,

μbμν=0.\displaystyle\nabla^{\mu}b_{\mu\nu}=0. (16)

From the bumblebee field motion equation (7), we have the projection of the Ricci tensor along the bumblebee field is

bμRμν=κλϱbν.\displaystyle b^{\mu}R_{\mu\nu}=\frac{\kappa\lambda}{\varrho}b_{\nu}. (17)

As to gravitational field equation (11), one can obtain the following three component equations

(D2)(1+)[2rψ(D3)]+e2ψ[(D2)(D3)2Λr2]=0,\displaystyle(D-2)(1+\ell)\big{[}2r\psi^{\prime}-(D-3)\big{]}+e^{2\psi}\big{[}(D-2)(D-3)-2\Lambda r^{2}\big{]}=0, (18)
2r2(ϕ′′+ϕ2ϕψ)2(D2)r(ψ+ϕ)2(D2)rϕ\displaystyle 2\ell r^{2}(\phi^{\prime\prime}+\phi^{\prime 2}-\phi^{\prime}\psi^{\prime})-2\ell(D-2)r(\psi^{\prime}+\phi^{\prime})-2(D-2)r\phi^{\prime}
+e2ψ[(D2)(D3)+2κλb2r22Λr2](1+)(D2)(D3)=0,\displaystyle\qquad\qquad+e^{2\psi}\big{[}(D-2)(D-3)+2\kappa\lambda b^{2}r^{2}-2\Lambda r^{2}\big{]}-(1+\ell)(D-2)(D-3)=0, (19)
(1+)[r2(ϕ′′+ϕ2ϕψ)+(D3)(D4)2+(D3)r(ϕψ)]\displaystyle(1+\ell)\Big{[}r^{2}(\phi^{\prime\prime}+\phi^{\prime 2}-\phi^{\prime}\psi^{\prime})+\frac{(D-3)(D-4)}{2}+(D-3)r(\phi^{\prime}-\psi^{\prime})\Big{]}
+e2ψ[Λr2(D3)(D4)2]=0,\displaystyle\qquad\qquad+e^{2\psi}\Big{[}\Lambda r^{2}-\frac{(D-3)(D-4)}{2}\Big{]}=0, (20)

where we have redefined the Lorentz-violating parameter =ϱb2\ell=\varrho b^{2} and, the prime is the derivative with respect to the corresponding argument, respectively. From the Eq. (18), one can obtain a metric function as

e2ψ=1+f(r),\displaystyle e^{2\psi}=\frac{1+\ell}{f(r)}, (21)

where

f(r)=116πM(D2)ΩD2rD32Λ(D1)(D2)r2,\displaystyle f(r)=1-\frac{16\pi M}{(D-2)\Omega_{D-2}r^{D-3}}-\frac{2\Lambda}{(D-1)(D-2)}r^{2}, (22)

and MM is the mass of the black hole. In order to recover the 4-dimensional Schwarzschild-like solution casana in the cases of Λ=0\Lambda=0, then we let another function to

e2ϕ=f(r).\displaystyle e^{2\phi}=f(r). (23)

It is easy to prove that the above results of e2ϕ=f(r),e2ψ=(1+)/f(r)e^{2\phi}=f(r),e^{2\psi}=(1+\ell)/f(r) can meet the gravitational Eq. (20). Then the bumblebee field is

bμ=(0,b(1+)/f(r),0,0,,0).\displaystyle b_{\mu}=\big{(}0,b\sqrt{(1+\ell)/f(r)},0,0,\cdots,0\big{)}. (24)

With the bumblebee field motion equation (17) or the gravitational Eq. (19), a solution of this type will be possible if and only if

Λ=(D2)κλ2ϱ(1+),\displaystyle\Lambda=\frac{(D-2)\kappa\lambda}{2\varrho}(1+\ell), (25)

which is the constraint on the additional field λ\lambda coming from the potential (10). It is easy to see that if the Lagrange-multiplier field λ=0\lambda=0, then Λ=0\Lambda=0 and the metric (22) is the DD-dimensional Schwarzschild-like black hole ding2022 . So λ\lambda isn’t a new freedom and we can define the effective cosmological constant Λe\Lambda_{e} as following to absorb it,

Λe=(D2)κλ2ϱ.\displaystyle\Lambda_{e}=\frac{(D-2)\kappa\lambda}{2\varrho}. (26)

Then Λ=(1+)Λe\Lambda=(1+\ell)\Lambda_{e} and, the present black hole metric is

ds2=f(r)dt2+1+f(r)dr2+r2dΩD22,\displaystyle ds^{2}=-f(r)dt^{2}+\frac{1+\ell}{f(r)}dr^{2}+r^{2}d\Omega_{D-2}^{2}, (27)

with the metric function f(r)f(r) as

f(r)=116πM(D2)ΩD2rD32(1+)Λe(D1)(D2)r2.\displaystyle f(r)=1-\frac{16\pi M}{(D-2)\Omega_{D-2}r^{D-3}}-\frac{2(1+\ell)\Lambda_{e}}{(D-1)(D-2)}r^{2}. (28)

When LV constant 0\ell\rightarrow 0, the effective cosmological constant ΛeΛ\Lambda_{e}\rightarrow\Lambda and it recovers the Schwarzschild-AdS black hole. When rr\rightarrow\infty, the metric function

f(r)=12(1+)Λe(D1)(D2)r2,\displaystyle f(r)=1-\frac{2(1+\ell)\Lambda_{e}}{(D-1)(D-2)}r^{2}, (29)

and the metric (27) becomes a DD-dimensional AdS-like spacetime. The black hole horizon locates at the largest real root of the following equation

116πM(D2)ΩD2rhD32(1+)Λe(D1)(D2)rh2=0.\displaystyle 1-\frac{16\pi M}{(D-2)\Omega_{D-2}r_{h}^{D-3}}-\frac{2(1+\ell)\Lambda_{e}}{(D-1)(D-2)}r_{h}^{2}=0. (30)

It is easy to see that it can give the black hole mass MM from the horizon radius,

M=(D2)ΩD216πrhD3(1+)Λe8π(D1)ΩD2rh(D1).\displaystyle M=\frac{(D-2)\Omega_{D-2}}{16\pi}r_{h}^{D-3}-\frac{(1+\ell)\Lambda_{e}}{8\pi(D-1)}\Omega_{D-2}r_{h}^{(D-1)}. (31)

It shows that the bumblebee field affects the location of the event horizons, contrary to the Schwarzschild-like solution in which the horizon radius is the same as of the Schwarzschild black hole.

III Thermodynamics of Einstein-bumblebee AdS-like black hole

A black hole is not only a gravity system, but also a special thermodynamic system due to that its surface gravity κ\kappa and horizon area AA are the close similarity to the temperature TT and entropy SS, T=κ/2πT=\kappa/2\pi, S=A/4S=A/4. Four laws of black hole thermodynamics were established in Ref. bardeen . LV can modify the geometry of a black hole, so in this kind of black hole, can the thermodynamical laws still hold? In this section we study the Smarr formula and the first law of this LV black hole.

III.1 Smarr formula

Suppose that \mathcal{M} is this DD-dimensional spacetime satisfying the Einstein equations, ξa=(1,0,0,0,,0)\xi^{a}=(1,0,0,0,\cdots,0) is a Killing vector on \mathcal{M}, timelike near infinity. In \mathcal{M}, there is a spacelike hypersurface SS with a co-dimension 2-surface boundary S\partial S, and ξa\xi^{a} is normal to the SS. The boundary S\partial S has two components: an inner boundary at the event horizon Sh\partial S_{h} and an outer boundary at infinity S\partial S_{\infty}. In order to construct a Komar integral relation for non-zero cosmological constant, one should introduce Killing potential ωab\omega^{ab} kastor which can be obtained according to relation ξb=aωab\xi^{b}=\nabla_{a}\omega^{ab}. For the static Killing vector ξa\xi^{a}, we have

ωrt=ωtr=rD1.\displaystyle\omega^{rt}=-\omega^{tr}=\frac{r}{D-1}. (32)

We can integrate the Killing equation b(bξa)=Rbaξb\nabla_{b}(\nabla^{b}\xi^{a})=-R^{a}_{b}\xi^{b} over a hypersurface SS,

SbξadΣab=SRbaξb𝑑Σa,\displaystyle\int_{\partial S}\nabla^{b}\xi^{a}d\Sigma_{ab}=-\int_{S}R^{a}_{b}\xi^{b}d\Sigma_{a}, (33)

where dΣabd\Sigma_{ab} and dΣad\Sigma_{a} are the surface elements of S\partial S and SS, respectively. With the metric (13) in mind, we have the relation

Rbaξb=2D2(Λe,0,0,0,,0)=2ΛeD2ξa,\displaystyle R^{a}_{b}\xi^{b}=\frac{2}{D-2}(\Lambda_{e},0,0,0,\cdots,0)=\frac{2\Lambda_{e}}{D-2}\xi^{a}, (34)

then the Eq. (33) can be rewritten as

D28πS(bξa+2D2Λeωab)𝑑Σab=0,\displaystyle\frac{D-2}{8\pi}\int_{\partial S}\left(\nabla^{b}\xi^{a}+\frac{2}{D-2}\Lambda_{e}\omega^{ab}\right)d\Sigma_{ab}=0, (35)

which is multiplied by the normalization factor (D2)/8π(D-2)/8\pi and called the Komar integral relation. The non-vanishing components of the tensor bξa\nabla^{b}\xi^{a} are given by

rξt=tξr=1(1+)(D2)[(D3)8πMΩD2rD22(1+)Λer(D1)].\displaystyle\nabla^{r}\xi^{t}=-\nabla^{t}\xi^{r}=\frac{1}{(1+\ell)(D-2)}\left[\frac{(D-3)8\pi M}{\Omega_{D-2}r^{D-2}}-\frac{2(1+\ell)\Lambda_{e}r}{(D-1)}\right]. (36)

The closed 2-surface S\partial S has two parts, horizon Sh\partial S_{h} and infinite S\partial S_{\infty}, so Eq. (35) can be rewritten as

D28πS(bξa+2D2Λeωab)𝑑Σab=D28πSh(bξa+2D2Λeωab)𝑑Σab.\displaystyle\frac{D-2}{8\pi}\int_{\partial S_{\infty}}\left(\nabla^{b}\xi^{a}+\frac{2}{D-2}\Lambda_{e}\omega^{ab}\right)d\Sigma_{ab}=\frac{D-2}{8\pi}\int_{\partial S_{h}}\left(\nabla^{b}\xi^{a}+\frac{2}{D-2}\Lambda_{e}\omega^{ab}\right)d\Sigma_{ab}. (37)

If we use the 2-surface element dΣrt=1+rD2dΩD2/2d\Sigma_{rt}=-\sqrt{1+\ell}r^{D-2}d\Omega_{D-2}/2, which is slightly modified by the factor 1+\sqrt{1+\ell}, the left hand side of this integral is

D28πS(bξa+2D2Λeωab)𝑑Σab=(D3)M1+;\displaystyle\frac{D-2}{8\pi}\int_{\partial S_{\infty}}\left(\nabla^{b}\xi^{a}+\frac{2}{D-2}\Lambda_{e}\omega^{ab}\right)d\Sigma_{ab}=-\frac{(D-3)M}{\sqrt{1+\ell}}; (38)

the second term of its right hand side is

Λe8πSh2ωab𝑑Σab=2PV1+,\displaystyle\frac{\Lambda_{e}}{8\pi}\int_{\partial S_{h}}2\omega^{ab}d\Sigma_{ab}=\frac{2PV}{\sqrt{1+\ell}}, (39)

where the pressure PP and the thermodynamic volume VV are

P=Λe8π,V=(1+)ΩD2D1rhD1;\displaystyle P=-\frac{\Lambda_{e}}{8\pi},\;\;V=(1+\ell)\frac{\Omega_{D-2}}{D-1}r_{h}^{D-1}; (40)

and the first term of its right hand side is

D28πShbξadΣab=(D2)TS1+,\displaystyle\frac{D-2}{8\pi}\int_{\partial S_{h}}\nabla^{b}\xi^{a}d\Sigma_{ab}=-\frac{(D-2)TS}{\sqrt{1+\ell}}, (41)

where the temperature TT and entropy SS are

T=4(D2)1+rh[(D3)MΩD2rhD3(1+)Λerh24π(D1)],S=A4=1+4ΩD2rhD2.\displaystyle T=\frac{4}{(D-2)\sqrt{1+\ell}r_{h}}\left[\frac{(D-3)M}{\Omega_{D-2}r_{h}^{D-3}}-\frac{(1+\ell)\Lambda_{e}r_{h}^{2}}{4\pi(D-1)}\right],\;\;S=\frac{A}{4}=\frac{\sqrt{1+\ell}}{4}\Omega_{D-2}r_{h}^{D-2}. (42)

Lastly, the integral (37) can give the Smarr formula

(D3)M=(D2)TS2PV.\displaystyle(D-3)M=(D-2)TS-2PV. (43)

One can see that the Smarr formula can still be constructed in this LV black hole spacetime as long as its temperature, entropy and volume are slightly modified.

Why the thermodynamic volume VV (40) and area AA (42) of the 2-surface S\partial S is different from those in Ref. kubiznak2017 (see the Eqs. 2.25 and 2.28 in it for 4-dimensions) by the factor 1+\sqrt{1+\ell} and (1+)(1+\ell)? Because these two quantities characterize a spacetime and should be entirely derived from thermodynamic considerations. So one should firstly determine its temperature. The temperature TT should be given by the black hole’s thermodynamic process— its Hawking radiation process ding2008

T=f(rh)g(rh)4π=14π1+rh[(D3)2D2(1+)Λerh2],\displaystyle T=\frac{\sqrt{f^{\prime}(r_{h})g^{\prime}(r_{h})}}{4\pi}=\frac{1}{4\pi\sqrt{1+\ell}\;r_{h}}\left[(D-3)-\frac{2}{D-2}(1+\ell)\Lambda_{e}r_{h}^{2}\right], (44)

where g(r)=f(r)/(1+)g(r)=f(r)/(1+\ell).

For the thermodynamic volume VV, it is the quantity thermodynamically conjugate to the pressure PP, then it should be (see the formula 3.1 in Ref. kubiznak2017 )

V(MP)S=(1+)ΩD2D1rhD1,\displaystyle V\equiv\Big{(}\frac{\partial M}{\partial P}\Big{)}_{S}=(1+\ell)\frac{\Omega_{D-2}}{D-1}r_{h}^{D-1}, (45)

with the Eq. (31) for the mass MM. For the thermodynamic area AA, it is the quantity thermodynamically conjugate to the temperature TT, then it should be

A=4S4dMT=1+ΩD2rhD2.\displaystyle A=4S\equiv 4\int\frac{dM}{T}=\sqrt{1+\ell}\Omega_{D-2}r_{h}^{D-2}. (46)

Therefore, one can see that the LV has affected the known spacetime geometry by redefinition of the geometrical area and volume.

III.2 The first law

Nextly, we derive the first law of black hole thermodynamics via Hamiltonian perturbation method. Note that in Ref. kastor , the authors used this method to lead to the first law and they found it was equivalent to the derivative of MM method which will be appeared in the appendix B, where MM is determined by Eq. (31).

In DD-dimensional manifold \mathcal{M} with metric gabg_{ab} of signature (,+,,+)(-,+,\cdots,+), let Σ\Sigma be a family of (D1)(D-1)-dimensional spacelike submanifolds with unite timelike normal field nan_{a} and induced metric sabs_{ab}, i.e.,

gab=nanb+sab,ncnc=1,ncscb=0.\displaystyle g_{ab}=-n_{a}n_{b}+s_{ab},\qquad n_{c}n^{c}=-1,\qquad n_{c}s^{cb}=0. (47)

The Hamiltonian variables are sabs_{ab} and its conjugate momentum πab\pi_{ab}. The energy density ρ=Tabnanb\rho=T_{ab}n^{a}n^{b} and momentum density Ja=TbcnbsacJ_{a}=T_{bc}n^{b}s^{c}_{\;a} must satisfy the Hamiltonian and momentum constrain equations,

H=16πGρ,Ha=16πGJa,\displaystyle H=16\pi G\rho,\qquad H_{a}=-16\pi GJ_{a}, (48)

where

H=2Gabnanb=[(D2)(D3)(1+)r2+2Λe],Ha=2Gbcnbsac=0\displaystyle H=-2G_{ab}n^{a}n^{b}=-\left[\frac{(D-2)(D-3)\ell}{(1+\ell)r^{2}}+2\Lambda_{e}\right],\qquad H_{a}=-2G_{bc}n^{b}s^{c}_{\;a}=0 (49)

for the given case. Let the Killing vector ξa=Fna+βa\xi^{a}=Fn^{a}+\beta^{a} with ncβc=0n_{c}\beta^{c}=0. The Hamiltonian density \mathcal{H} for evolution along ξa\xi^{a} is given by

=s{F[H+2Λe+(D2)(D3)(1+)r2]+βaHa}.\displaystyle\mathcal{H}=\sqrt{s}\left\{F\left[H+2\Lambda_{e}+\frac{(D-2)(D-3)\ell}{(1+\ell)r^{2}}\right]+\beta_{a}H^{a}\right\}. (50)

We assume that the linear approximation metric g~ab=gab+δgab\tilde{g}_{ab}=g_{ab}+\delta g_{ab} is another nearby solution to the Einstein equations with a perturbed cosmological constant Λ+δΛ\Lambda+\delta\Lambda. The induced spatial metric and momentum for this perturbed metric are s~ab=sab+hab\tilde{s}_{ab}=s_{ab}+h_{ab} and π~ab=πab+pab\tilde{\pi}_{ab}=\pi_{ab}+p_{ab}, where hab=δsab,pab=δπabh_{ab}=\delta s_{ab},p_{ab}=\delta\pi_{ab}. From Hamilton’s equations for the zeroth order, the linearized constrain operators δH\delta H and δHa\delta H_{a} combine to form a total derivative,

FδH+βaδHa=DcBc,\displaystyle F\delta H+\beta^{a}\delta H_{a}=-D_{c}B^{c}, (51)

where DaD_{a} is the covariant derivative operator on Σ\Sigma compatible with sabs_{ab}, and 111 In the general case of non-vanishing extrinsic curvature, the boundary vector BaB^{a} includes the term kastor2014 1sβb(πcdhcdsba2πachbc2pba)\frac{1}{\sqrt{s}}\beta^{b}(\pi^{cd}h_{cd}s^{a}_{\;b}-2\pi^{ac}h_{bc}-2p^{a}_{\;b}).

Ba=F(DahDbhab)hDaF+habDbF.\displaystyle B^{a}=F(D^{a}h-D_{b}h^{ab})-hD^{a}F+h^{ab}D_{b}F. (52)

The linear constraints (48) take the form of a Gauss law 222From the Eq. (49), the perturbation δH=16πGδρ2δΛeδ(1/r2)2δΛe\delta H=16\pi G\delta\rho\sim-2\delta\Lambda_{e}-\delta(1/r^{2})\sim-2\delta\Lambda_{e} at large distance. So the contribution of this additional matter field source, i.e., the purely radial bumblebee field is only included in the term δΛe\delta\Lambda_{e}.

DcBc=2FδΛe.\displaystyle D_{c}B^{c}=2F\delta\Lambda_{e}. (53)

Using Killing potential ωab\omega^{ab}, we have F=naξa=Dc(naωca)F=-n_{a}\xi^{a}=-D_{c}(n_{a}\omega^{ca}), and rewrite Eq. (53) in the integral form,

I=Σ𝑑ac(Bc+2ωcdndδΛe)=0,\displaystyle I=\int_{\partial\Sigma}da_{c}(B^{c}+2\omega^{cd}n_{d}\delta\Lambda_{e})=0, (54)

which is taken to extend from a boundary Σh\partial\Sigma_{h} at the bifurcation sphere of the horizon to a boundary Σ\partial\Sigma_{\infty} at infinity, i.e., IIh=0I_{\infty}-I_{h}=0,

Σ𝑑ac(Bc+2ωcdndδΛe)Σh𝑑ac(Bc+2ωcdndδΛe)=0.\displaystyle\int_{\partial\Sigma_{\infty}}da_{c}(B^{c}+2\omega^{cd}n_{d}\delta\Lambda_{e})-\int_{\partial\Sigma_{h}}da_{c}(B^{c}+2\omega^{cd}n_{d}\delta\Lambda_{e})=0. (55)

We choose Ff,na=Fat,dar=1+rD2dΩD2/fF\simeq\sqrt{f},\;n_{a}=-F\nabla_{a}t,\;da_{r}=\sqrt{1+\ell}\;r^{D-2}d\Omega_{D-2}/\sqrt{f}, and

hrrδ(1+f),δf=[16πδM(D2)ΩD2rD3+2(1+)δΛer2(D1)(D2)],\displaystyle h_{rr}\simeq\delta\left(\frac{1+\ell}{f}\right),\qquad\delta f=-\left[\frac{16\pi\delta M}{(D-2)\Omega_{D-2}r^{D-3}}+\frac{2(1+\ell)\delta\Lambda_{e}r^{2}}{(D-1)(D-2)}\right], (56)

then the integral (54) can be simplified by

I=11+ΣrD2𝑑ΩD2[D2rδffδA1+A+2(1+)rδΛe(D1)],\displaystyle I=\frac{1}{\sqrt{1+\ell}}\int_{\partial\Sigma}r^{D-2}d\Omega_{D-2}\left[\frac{D-2}{r}\delta f-\frac{f^{\prime}\delta A}{\sqrt{1+\ell}A}+\frac{2(1+\ell)r\delta\Lambda_{e}}{(D-1)}\right], (57)

where the prime denotes the derive with the argument, and AA is the area of the boundary. So I=16πδM/1+I_{\infty}=-16\pi\delta M/\sqrt{1+\ell} since 1/A01/A\rightarrow 0 at infinite, and

Ih=16π1+[TδSVδP],\displaystyle I_{h}=\frac{16\pi}{\sqrt{1+\ell}}[-T\delta S-V\delta P], (58)

since f(rh)=0,f(rh)=4π1+T,Λe=8πP,A=4Sf(r_{h})=0,\;f^{\prime}(r_{h})=4\pi\sqrt{1+\ell}T,\;\Lambda_{e}=-8\pi P,\;A=4S and V=(1+)ΩD2rhD1/(D1)V=(1+\ell)\Omega_{D-2}r_{h}^{D-1}/(D-1) 333For the general form of the volume VV, see the formula (22) in Ref. kastor . Lastly, the first law is obtained,

δM=TδS+VδP,\displaystyle\delta M=T\delta S+V\delta P, (59)

which still also holds in this LV black hole spacetime.

IV Phase transition of Einstein-bumblebee AdS-like black hole

A phase transition is a discontinuous change in the properties of a substance, as its environment is changed only infinitesimally. An isolated Schwarzschild black hole is thermodynamically unstable due to its negative heat capacity. However for a Schwarzschild AdS black hole, besides the negative heat phase, a phase of positive heat capacity will occur. At a given temperature and pressure, the stable phase is always the one with the lower Gibbs free energy. In an AdS space, the Gibbs energy of the thermal radiation phase is zero wei . So the negative free energy of the positive heat capacity black hole phase is stable phase. In this section, we study this two phase transitions.

From the Eq. (30), we can express the mass MM as

M=(D2)ΩD216πrhD3+PV.\displaystyle M=\frac{(D-2)\Omega_{D-2}}{16\pi}r_{h}^{D-3}+PV. (60)

MM can be interpreted as an enthalpy kastor H=U+PVH=U+PV, so the black hole internal energy UU is

U=(D2)ΩD216πrhD3.\displaystyle U=\frac{(D-2)\Omega_{D-2}}{16\pi}r_{h}^{D-3}. (61)

Then the first law (84) can be rewritten as following

dU=TdSPdV,\displaystyle dU=TdS-PdV, (62)

which is the same as the fundamental thermodynamical equation of a simple gas system. The temperature (44) can be rewritten as

T=14π1+rh[(D3)+16πPD2(1+)rh2],\displaystyle T=\frac{1}{4\pi\sqrt{1+\ell}\;r_{h}}\left[(D-3)+\frac{16\pi P}{D-2}(1+\ell)r_{h}^{2}\right], (63)

which is shown in Fig. 1 when D=4D=4.

Refer to caption
Refer to caption
Figure 1: Black hole temperature and heat capacity CPC_{P} at constant pressure PP vs. black hole mass MM with D=4D=4 and different coupling constant \ell. The temperature has a minimum value and the the heat capacity is divergent at the point of black hole mass M=MdM=M_{d}, where Md=1/18(1+)πPM_{d}=1/\sqrt{18(1+\ell)\pi P}.

There exist a minimum temperature

Tmin=2(D3)P(D2)π,\displaystyle T_{min}=2\sqrt{\frac{(D-3)P}{(D-2)\pi}}, (64)

when rh=(D2)(D3)/[16(1+)πP]r_{h}=\sqrt{(D-2)(D-3)/[16(1+\ell)\pi P]}. Note that TminT_{min} isn’t dependent on \ell. From the left panel of the Fig. 1, when T>TminT>T_{min}, there are two possible black hole masses hawking which can be in equilibrium with thermal radiation, i.e., M>MdM>M_{d} or M<MdM<M_{d}, where MdM_{d} is the black hole mass whose heat capacity is divergent (66). Fig. 1 shows that the temperature of the large black hole (M>MdM>M_{d}) increases when the LV coupling constant \ell becomes bigger. At a fixed pressure PP, the black hole heat capacity CPC_{P} is

CP=(HT)P=(D2)ΩD24rhD216π(1+)Prh2+(D2)(D3)16π(1+)Prh2(D2)(D3),\displaystyle C_{P}=\Big{(}\frac{\partial H}{\partial T}\Big{)}_{P}=\frac{(D-2)\Omega_{D-2}}{4}r_{h}^{D-2}\frac{16\pi(1+\ell)Pr_{h}^{2}+(D-2)(D-3)}{16\pi(1+\ell)Pr^{2}_{h}-(D-2)(D-3)}, (65)

which is also shown in Fig. 1 when D=4D=4. From the right panel of the Fig. 1, one can see that the small black hole mass has negative specific heat, so it is unstable and will decay. While the large black hole mass has positive heat capacity and therefore it is thermally stable. It is easy to see that the heat capacity at constant pressure is divergent at the black hole mass MdM_{d},

Md=(D2)2ΩD28π(D1)[(D2)(D3)16(1+)πP](D3)/2.\displaystyle M_{d}=\frac{(D-2)^{2}\Omega_{D-2}}{8\pi(D-1)}\left[\frac{(D-2)(D-3)}{16(1+\ell)\pi P}\right]^{(D-3)/2}. (66)

It is easy to see that the LV constant \ell decreases the mass MdM_{d}, showing that the phase transition from a small black hole to large one will more easily occur. At this divergent point, the black hole entropy is

Sd=4(1D)ΩD21+(3D)[(D2)(D3)πP](D2)/2.\displaystyle S_{d}=4^{(1-D)}\Omega_{D-2}\sqrt{1+\ell}^{(3-D)}\left[\frac{(D-2)(D-3)}{\pi P}\right]^{(D-2)/2}. (67)

Fig. 2 shows the entropy with the temperature at unit pressure when D=4D=4.

Refer to caption
Figure 2: Black hole entropy SS vs. black hole temperature TT with D=4D=4 and different coupling constant \ell. They all have both branches, the above one is corresponding to the large black hole M>MdM>M_{d}, the below one to the small black hole M<MdM<M_{d}. At the temperature T=TminT=T_{min}, the entropy is Sd=1/81+PS_{d}=1/8\sqrt{1+\ell}P. At the HP temperature THP=8P/3π,SHP=3/81+PT_{HP}=\sqrt{8P/3\pi},S_{HP}=3/8\sqrt{1+\ell}P for the large black hole, Ssmall=1/241+PS_{small}=1/24\sqrt{1+\ell}P for the unstable small black hole.

The phase transition between these two phases is called small-large black hole phase transition.

There is another kind of phase transition called Hawking-Page(HP) phase transition, which is between a metastable large black hole phase and stable large black hole phase in an anti-de Sitter spacetime. Witten witten explained that this HP phase transition is a confinement/deconfinement phase transition in the AdS/CFT correspondence. It can also be understood as a solid/liquid phase transition in black hole chemistry kubiznak2015 . HP phase transition occurs when Gibbs free energy is zero. The black hole Gibbs free energy can be obtained

G=HTS=ΩD28π(D1)[(D2)2π1+Trh]rhD3,\displaystyle G=H-TS=\frac{\Omega_{D-2}}{8\pi(D-1)}\left[(D-2)-2\pi\sqrt{1+\ell}Tr_{h}\right]r_{h}^{D-3}, (68)

which is shown in Fig. 3 when D=4D=4. It can be seen that there are two branches, the above parts corresponding to small black hole phase whose free energy are positive for any temperature; the below ones corresponding to metastable large black hole phase when G>0G>0 and stable large black hole phase when G<0G<0. When the temperature T<TminT<T_{min}, the spacetime has only one phase—thermal radiation, whose free energy is zero.

Refer to caption
Figure 3: Black hole Gibbs free energy vs. temperature with D=4D=4 and different coupling constant \ell. They vanish at the point THP=8P/3πT_{HP}=\sqrt{8P/3\pi}. They all have both branches, the above one is corresponding to the small black hole M<MdM<M_{d}, the below one to the large black hole M>MdM>M_{d}. At the temperature T=TminT=T_{min}, the free energy has a maximum Gmax=1/122(1+)πPG_{max}=1/12\sqrt{2(1+\ell)\pi P}.

At the temperature T=TminT=T_{min}, its free energy has a maximum value GmaxG_{max},

Gmax=ΩD2(D1)π2(32D)[(D2)(D3)(1+)πP](D3)/2.\displaystyle G_{max}=\frac{\Omega_{D-2}}{(D-1)\pi}2^{(3-2D)}\left[\frac{(D-2)(D-3)}{(1+\ell)\pi P}\right]^{(D-3)/2}. (69)

The temperature at Hawking-Page phase transition point can be obtained when G=0G=0,

THP=2(D2)P(D1)π,\displaystyle T_{HP}=2\sqrt{\frac{(D-2)P}{(D-1)\pi}}, (70)

which is also independent on \ell like the minimum temperature TminT_{min}. At this point rh=(D1)(D2)/[16(1+)πP]r_{h}=\sqrt{(D-1)(D-2)/[16(1+\ell)\pi P]} and the black hole mass is,

MHP=(D2)ΩD28π[(D1)(D2)16(1+)πP](D3)/2,\displaystyle M_{HP}=\frac{(D-2)\Omega_{D-2}}{8\pi}\left[\frac{(D-1)(D-2)}{16(1+\ell)\pi P}\right]^{(D-3)/2}, (71)

its entropy is

SHP=4(1D)ΩD21+(3D)[(D1)(D2)πP](D2)/2.\displaystyle S_{HP}=4^{(1-D)}\Omega_{D-2}\sqrt{1+\ell}^{(3-D)}\left[\frac{(D-1)(D-2)}{\pi P}\right]^{(D-2)/2}. (72)

We summarize the three phases: unstable small black hole, metastable large black hole and stable large black hole in Tab. 1; stable large black hole, metastable large black hole and thermal radiation in Fig. 4.

Table 1: Quantities of heat capacity, entropy, temperature and Gibbs free energy for three phases of the AdS black hole in Einstein-bumblebee gravity.
Black hole mass Heat capacity Entropy Temperature Free energy
Unstable small black hole 0<M<Md0<M<M_{d} Cp<0C_{p}<0 S<SdS<S_{d} T>TminT>T_{min} G>0G>0
Metastable large black hole Md<M<MHPM_{d}<M<M_{HP} Cp>0C_{p}>0 Sd<S<SHPS_{d}<S<S_{HP} Tmin<T<THPT_{min}<T<T_{HP} 0<G<Gd0<G<G_{d}
Stable large black hole M>MHPM>M_{HP} Cp>0C_{p}>0 S>SHPS>S_{HP} T>THPT>T_{HP} G<0G<0
Refer to caption
Figure 4: Phase diagram (TPT-P figure at unit mass MM) for D=4D=4 AdS black hole in Einstein-bumblebee gravity. The red solid and green dashed curves respectively correspond to the black hole minimum temperature and the HP phase transition temperature. Both curves extend to infinity.

From the figures 1 to 4 and the table I, one can see the following picture for the formation of a AdS black hole: i), in an AdS space, the thermal radiation with zero Gibbs free energy has a very high temperature, accidentally, a tiny mass black hole begins to form due to quantum fluctuation; ii), the tiny black hole begins to grow, its temperature decreases quickly, its free energy increases slowly; iii), when the black hole mass grows to MdM_{d}, its free energy approaches to a maximum, its temperature approaches to a minimum, and the small-large black hole phase transition occurs; iv), the black hole mass continues to grow, its positive free energy decreases, its temperature increases slowly and it is a metastable black hole; v), when the mass grows to MHPM_{HP}, the HP phase transition occurs; vi), finally, the stable large black hole has formed.

Figure 3 also shows that the LV constant \ell decreases the Gibbs free energy of the meta-stable large black hole with the mass Md<M<MHPM_{d}<M<M_{HP}, showing that the HP phase transition from a meta-stable large black hole to large stable one will be more easily to occur.

Now we consider some relationships between the both kinds of phase transition. Wei and Liu wei found that a novel exact dual relation between the minimum temperature of the D+1D+1-dimensional black hole and the HP phase transition temperature in DD dimensions, reminiscent of the holographic principle. Here we will examine wether it still hold or not to this LV black hole. The first relationship is for the temperature and entropy,

THPTmin=D2(D1)(D3),SHPSd=(D1D3)(D2)/2,\displaystyle\frac{T_{HP}}{T_{min}}=\frac{D-2}{\sqrt{(D-1)(D-3)}},\;\frac{S_{HP}}{S_{d}}=\left(\frac{D-1}{D-3}\right)^{(D-2)/2}, (73)

which are only dependant on the dimension number DD and the same as those in Refs. wei ; belhaj . The second relationship is between the DD dimensional Hawking-Page temperature THPT_{HP} and the D+1D+1 dimensional minimum temperature TminT_{min},

THP(D)=Tmin(D+1),\displaystyle T_{HP}(D)=T_{min}(D+1), (74)

which is the novel dual relation found in Ref. wei . This relation is independent of the pressure and the bumblebee coupling constant, and is like the AdS/CFT correspondence. TminT_{min} is a physical quantity in the bulk, then THPT_{HP} can be treated as the dual physical quantity on the boundary.

V Summary

In this paper, we have derived the high dimensional AdS-like black hole solution in Einstein-bumblebee gravity theory. It can exists only under the condition that the bumblebee potential has a linear functional form with a Lagrange-multiplier field λ\lambda. This additional field is strictly constrained by the bumblebee motion equation and can be absorbed by the definition of an effective cosmological constant Λe\Lambda_{e}. Like Schwarzschild-like black hole casana , it can’t be asymptotically to anti-de Sitter spacetime, so it is a Schwarzschild-AdS-like black hole solution maluf . Contrary to Schwarzschild-like black hole, the bumblebee field will affect the locations of the black hole horizon.

By using black hole chemistry, we have studied the thermodynamics and phase transitions of this high dimensional Schwarzschild-AdS-like black hole. Via Komar integral method, we find that the Smarr formula can exist as long as its temperature, entropy and volume are slightly modified, and the LV has affected the known spacetime geometry, such as area and volume. Via Hamiltonian perturbation method, there still exist the first law of black hole thermodynamics.

Its temperature has a minimum value TminT_{min}, corresponding to the black hole mass MdM_{d}. Small black hole with mass 0<M<Md0<M<M_{d} has a negative heat capacity CpC_{p}, while the large black hole with mass M>MdM>M_{d} has a positive one. At the point MdM_{d}, the CpC_{p} is divergent. The LV constant \ell decreases the mass MdM_{d}, showing that the phase transition from a small black hole to large one will more easily occur.

Its Gibbs energy has a maximum value GmaxG_{max} at the point of mass MdM_{d} and has a zero point G=0G=0 at the point of mass MHPM_{HP}. Small black hole with mass 0<M<Md0<M<M_{d} and large black hole with mass Md<M<MHPM_{d}<M<M_{HP} have a positive Gibbs free energy, while the large black hole with mass M>MdM>M_{d} has a negative one. The LV constant \ell decreases the Gibbs free energy of the meta-stable large black hole with the mass Md<M<MHPM_{d}<M<M_{HP}, showing that the HP phase transition from a meta-stable large black hole to large stable one will be more easily to occur.

The dualities found in Ref. wei still holds when the Lorentz symmetry is spontaneously broken: the ratio of temperature between small/large black hole phase transition and HP phase transition THP/TminT_{HP}/T_{min}, and the equivalence of the temperature of DD-dimensional HP phase transition to the temperature of D+1D+1 dimensional small/large black hole phase transition THP(D)=Tmin(D+1)T_{HP}(D)=T_{min}(D+1).

Acknowledgements.
This work was supported by the Scientific Research Fund of the Hunan Provincial Education Department under No. 19A257 and No. 19A260, the National Natural Science Foundation (NNSFC) of China (grant No. 11247013), Hunan Provincial Natural Science Foundation of China grant No. 2015JJ2085.

Appendix A Proof for the potential V=κx2/2V=\kappa x^{2}/2 admitting no black hole solution with Λ0\Lambda\neq 0

If one uses the smooth potential V=κx2/2V=\kappa x^{2}/2, then under the condition (9) it becomes V=V=0V=V^{\prime}=0, and the motion equation (17) becomes brRrr=0b^{r}R_{rr}=0. In a four dimensional spacetime, the gravitational equations (18,19,20) become,

(1+)(2rψ1)+e2ψ(1Λr2)=0\displaystyle(1+\ell)(2r\psi^{\prime}-1)+e^{2\psi}(1-\Lambda r^{2})=0 (75)
r2(ϕ′′+ϕ2ϕψ)2r(ψ+ϕ)2rϕ+e2ψ(1Λr2)=0\displaystyle\ell r^{2}(\phi^{\prime\prime}+\phi^{\prime 2}-\phi^{\prime}\psi^{\prime})-2\ell r(\psi^{\prime}+\phi^{\prime})-2r\phi^{\prime}+e^{2\psi}(1-\Lambda r^{2})=0 (76)
(1+)[r2(ϕ′′+ϕ2ϕψ)+r(ϕψ)]+e2ψΛr2=0,\displaystyle(1+\ell)[r^{2}(\phi^{\prime\prime}+\phi^{\prime 2}-\phi^{\prime}\psi^{\prime})+r(\phi^{\prime}-\psi^{\prime})]+e^{2\psi}\Lambda r^{2}=0, (77)

and Rrr=0R_{rr}=0 becomes

Rrr=2rψ(ϕ′′+ϕ2ϕψ)=0.\displaystyle R_{rr}=\frac{2}{r}\psi^{\prime}-(\phi^{\prime\prime}+\phi^{\prime 2}-\phi^{\prime}\psi^{\prime})=0. (78)

by using Eq. (78), the Eqs. (76) and (77) become

(1+)(2rϕ+1)+e2ψ(1Λr2)=0\displaystyle-(1+\ell)(2r\phi^{\prime}+1)+e^{2\psi}(1-\Lambda r^{2})=0 (79)
(1+)r(ϕ+ψ)+e2ψΛr2=0.\displaystyle(1+\ell)r(\phi^{\prime}+\psi^{\prime})+e^{2\psi}\Lambda r^{2}=0. (80)

Combining Eq. (79) and (75), one can obtain ϕ+ψ=0\phi^{\prime}+\psi^{\prime}=0, which will lead Eq. (80) to the result e2ψΛr2=0e^{2\psi}\Lambda r^{2}=0, i.e., Λ=0\Lambda=0. Therefore, it is proved that there is no black hole solution with nonzero Λ\Lambda for the above gravitational equations (75) to (77) under the smooth potential V=κx2/2V=\kappa x^{2}/2.

Appendix B An equivalent way to the first law

By using the thermodynamical volume VV (45) and the entropy SS (46), one can rewrite the mass formula (31) as

M=M(S,P)=(D2)S(rh)4π1+rh+PV(rh).\displaystyle M=M(S,P)=\frac{(D-2)S(r_{h})}{4\pi\sqrt{1+\ell}r_{h}}+PV(r_{h}). (81)

Then its total differential is

dM=(MS)PdS+(MP)SdP.\displaystyle dM=\Big{(}\frac{\partial M}{\partial S}\Big{)}_{P}dS+\Big{(}\frac{\partial M}{\partial P}\Big{)}_{S}dP. (82)

In it, the partial derivative (MP)S=V(\frac{\partial M}{\partial P})_{S}=V, another partial derivative is

(MS)P=(MS)P,V,rh+(Mrh)P,S(drhdS),\displaystyle\Big{(}\frac{\partial M}{\partial S}\Big{)}_{P}=\Big{(}\frac{\partial M}{\partial S}\Big{)}_{P,V,r_{h}}+\Big{(}\frac{\partial M}{\partial r_{h}}\Big{)}_{P,S}\Big{(}\frac{dr_{h}}{dS}\Big{)}, (83)

whose result is the same as the Eq. (44). Lastly Eq.(82) gives the first law of black hole thermodynamics,

dM=TdS+VdP.\displaystyle dM=TdS+VdP. (84)

References

  • (1) A. Ashtekar and V. Petkov, Springer handbook of spacetimeP485, Springer Dordrecht, (2014).
  • (2) D. Mattingly, Living Rev. Rel.  8, 5 (2005).
  • (3) G. Amelino-Camelia, Liv. Rev. Rel. 16, 5 (2013).
  • (4) V. A. Kostelecký, Phys. Rev. D 69, 105009 (2004).
  • (5) V. A. Kostelecký and S. Samuel, Phys. Rev. D 39, 683 (1989).
  • (6) S. M. Carroll, J. A. Harvey, V. A. Kostelecky, C. D. Lane, and T. Okamoto, Phys. Rev. Lett. 87, 141601 (2001); I. Mocioiu, M. Pospelov, and R. Roiban, Phys. Lett. B 489, 390 (2000); A. F. Ferrari, M. Gomes, J. R. Nascimento, E. Passos, A. Yu. Petrov, and A. J. da Silva, Phys. Lett. B 652, 174 (2007).
  • (7) O. Bertolami, R. Lehnert, R. Potting, and A. Ribeiro, Phys. Rev. D 69, 083513 (2004); V. A. Kosteleckyý, R. Lehnert, and M. Perry, Phys. Rev. D 68, 123511 (2003); R. Jackiw and S.-Y. Pi, Phys. Rev. D 68, 104012 (2003).
  • (8) R. Gambini and J. Pullin, Phys. Rev. D 59, 124021 (1999); J. R. Ellis, N. E. Mavromatos, and D. V. Nanopoulos, Gen. Relativ. Gravit. 32, 127 (2000).
  • (9) C. P. Burgess, J. Cline, E. Filotas, J. Matias, and G.D. Moore, J. High Energy Phys. 03 (2002) 043; A. R. Frey, J. High Energy Phys. 04 (2003) 012; J. Cline and L. Valcárcel, J. High Energy Phys. 03 (2004) 032.
  • (10) A. Fernando, T. Clark, Gen. Relativ. Grav. 46, 1834 (2014).
  • (11) T. Jacobson and D. Mattingly, Phys. Rev. D 64, 024028 (2001); T. Jacobson, Proc. Sci. QG-PH (2007) 020.
  • (12) R. Bluhm, Shu-Hong Fung and V. A. Kostelecký, Phys. Rev. D 77, 065020 (2008).
  • (13) O. Bertolami and J. Páramos, Phys. Rev. D 72, 044001 (2005).
  • (14) R. Casana and A. Cavalcante, Phys. Rev. D 97, 104001 (2018).
  • (15) C. Ding, X. Chen and X. Fu, Nucl. Phys. B 975, 115688 (2022).
  • (16) C. Ding, C. Liu, R. Casana and A. Cavalcate, Eur. Phys. C 80, 178 (2020).
  • (17) S. K. Jha and A. Rahaman, arXiv: 2011.14916, (2021).
  • (18) R. V. Maluf and J. C. S. Neves, Phys. Rev. D 103, 044002 (2021).
  • (19) T. Padmanabhan, Gravitation: Foundations and Frontiers, Cambridge University Press, (2010).
  • (20) S.W. Hawking, Commun. Math. Phys. 43, 199 (1975); J. Bekenstein, Lett. Nuovo Cimento 4, 737 (1972); J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973).
  • (21) S. W. Hawking and D. N. Page, Cummum. Math. Phys. 87, 577 (1983).
  • (22) C. Ding and X. Chen, Chin. Phys. C 45, 025106 (2021).
  • (23) D. G. Boulware and S. Deser, Phys. Rev. Lett. 55, 2656 (1985).
  • (24) J. M. Bardeen, B. Carter, and S. Hawking, Commun. Math. Phys. 31, 161 (1973).
  • (25) D. Kastor, S. Ray and J. Traschen, Class. Quantum Grav. 26, 195011 (2009).
  • (26) D. Kastor, S. Ray and J. Traschen, J. High Energy Phys. 11 (2014) 120.
  • (27) D. Kubizňák, R. B. Mann and M. Teo, Class. Quantum Grav. 34, 063001 (2017).
  • (28) C. Ding and J. Jing, Class. Quantum Grav. 25, 145015 (2008).
  • (29) Shao-Wen Wei and Yu-Xiao Liu and R. B. Mann, Phys. Rev. D 102, 104011 (2020).
  • (30) E. Witten, Anti-de Sitter space, Adv. Theor. Math. Phys. 2, 505 (1998).
  • (31) D. Kubizňák and R. B. Mann, Can. J. Phys. 93, 999 (2015).
  • (32) A. Belhaj, A. El Balali, W. El Hadri and E. Torrente-Lujan, Phys. Lett. B 811, 135871 (2020).