Black holes in vector-tensor theories and their thermodynamics

In this paper, we study Einstein gravity either minimally or non-minimally coupled to a vector field which breaks the gauge symmetry explicitly in general dimensions. We first consider a minimal theory which is simply the Einstein-Proca theory extended with a quartic self-interaction term for the vector field. We obtain its general static maximally symmetric black hole solution and study the thermodynamics using Wald formalism. The aspects of the solution are much like a Reissner-Nordstr{\o}m black hole in spite of that a global charge cannot be defined for the vector. For non-minimal theories, we obtain a lot of exact black hole solutions, depending on the parameters of the theories. In particular, many of the solutions are general static and have maximal symmetry. However, there are some subtleties and ambiguities in the derivation of the first laws because the existence of an algebraic degree of freedom of the vector in general invalids the Wald entropy formula. The thermodynamics of these solutions deserves further studies.

On the other hand, there exists a no-go theorem which excludes the existence of Einstein-Proca black holes in asymptotically flat space-times [20,21]. However, the theorem is easily evaded. By numerical analysis, it was established in [22] that when Einstein gravity minimally coupled to an even number of real Proca fields, there exist asymptotically flat, stationary, axisymmetric black holes with Proca hair. It was analytically shown in [11,12,13,14,15,16,17,18,19] that the no-go theorem can be avoided in the presence of non-minimal couplings between the curvature and the vector fields.
Yet, there are still some holes left in the literatures which motive our current work. The first is in vector-tensor theories the vector field is as physical as the field strength since the gauge symmetry is explicitly breaking owing to either a nonzero bare mass or nonminimal couplings. A direct consequence of this is one can introduce a non-trivial radial component for the vector field A r when solving black hole solutions in the dual theories [12,14,16]. However, the power of this has not been considered very well. In this paper, we will show that in many cases how one can obtain the most general static maximally symmetric solutions with a nonzero A r in general dimensions. The most simple example we study is a minimal theory which generalizes the free massive Proca theory with a quartic self-interaction term. We also obtain such general static solutions for non-minimal theories with coupling terms of the form R µν A µ A ν and G µν A µ A ν , where R µν , G µν are the Ricci tensor and Einstein tensor of the metric respectively. In the four dimension, the later case has been well studied in [16].
Our second motivation is while people have obtained a series of exact black holes with vector hairs, their thermodynamics was even not studied except for a few papers [11,13].
Here we will adopt the Wald formalism to derive the first law of thermodynamics systematically for all the solutions we obtain. Moreover, we find some subtleties and ambiguities when deriving the first laws for the solutions with a nonzero A r . The underlying reason is A r is a purely algebraic degree of freedom which does not introduce corresponding vector charges in the solutions. However, to govern the validity of Wald entropy formula, one should impose proper boundary conditions on the horizon for both A t and A r , which in general results to a degenerate solution characterized by only one parameter, in contrast to the general two-parameter family solutions. Of course, this does not make sense in the derivation of the first law. Thus, one has to relax the horizon condition for A r but conversely this in general invalids the Wald entropy formula. For more discussions, we refer the readers to section 2 and section 4.
The paper is organized as follows. In section 2, we study a certain type vector-tensor theories. We analyze the structure of the general static maximally symmetric solutions. We also briefly review the Wald formalism, derive explicit formulas for our gravity model and discuss the subtleties in the derivation of the first law. In section 3, we study the minimally coupled theory by introducing a quartic self-interaction term for the vector. We obtain the general static solution with A r = 0 and study various properties of the solutions. From section 4 to section 6, we study non-minimally coupled vector-tensor theories and obtain a lot of exact black hole solutions depending on the parameters of the theories. We also derive the first law using Wald formalism. We conclude in section 7.

.1 Structure of general static solutions
In this paper, we consider Einstein gravity either minimally or non-minimally coupled to a vector field together with a potential V (it should not be confused with the vector field A).
The Lagrangian density is given by where F = dA and ψ ≡ A µ A µ . Note that the effective gravitational coupling constant is inversely proportional to To avoid ghost-like graviton modes, we require κ eff being positive definite throughout this paper. In addition, the γ coupling term can be written more explicitly as This is a special case discussed in [5], where a general construction of vector-tensor theories preserving parity has been well studied.
The covariant equations of motions are where G µν = R µν − 1 2 Rg µν is the Einstein tensor and For later convenience, we denote the Einstein and the vector equations of motions in (4) by E µν = 0 and P µ = 0 respectively.
In particular, we are interested in a vector potential of the type where γ 4 is a coupling constant characterizing the self-interaction of the vector field. Hence, the general theories are characterized by five independent parameters (β , γ , Λ 0 , m 2 , γ 4 ).
For γ = 0, the theories with such a potential were first studied in [13] whilst the γ 4 = 0 case has been studied in [11,12,14,16] for certain coupling constants but most of them are limited to the four dimension. Instead, in this paper we will investigate the theories for general coupling constants and solve the static maximally symmetric solutions in general dimensions.
The most simple solutions of the theories (1) are given by It follows that depending on the sign of the bare cosmological constant, the maximally respectively. Linearizing the equations of motions around the vacuum, we find that the linear fluctuations of the equations are described by a massless graviton and a Proca which has an effective mass where n denotes the space-time dimensions. Notice that owing to the existence of the nonminimal couplings, an effective Proca mass can be generated in the vacuum even if the bare mass vanishes. Likewise, even if the bare mass is nonzero, the U(1) gauge symmetry of the vector can be restored at the linear level when the parameters are such that m 2 eff = 0. This is true for any Ricci-flat metric, including Schwarzschild and Kerr black holes.
The most general ansatz for static maximally symmetric solutions is where h , f , A t , A r are all functions of r and dΩ 2 n−2 ,k is the metric of the codimension-2 space with spherical/hyperbolic/toric symmetries, corresponding to k = 1 , −1 , 0, respectively. It is easy to see that the vector equation P r is purely algebraic for A r . We find where the dotted term is composed of the functions h , f , A t and their derivatives with respect to r ( this term exactly vanishes for Einstein-Proca theory and hence the solutions with a nonzero A r do not exist in this case. ). It is clear that the above equation has isolated roots A r = 0 and A r = 0, corresponding to different branch solutions. Consequently, in general the solutions with A r = 0 do not have a smooth limit to send A r → 0 and reduce to the solutions with A r = 0. This is also true even if γ 4 = 0, in which case A r in general can not be solved algebraically 1 . Nonetheless, the ansatz (9) is most general for both A r = 0 and A r = 0 solutions. We will study either of the two cases or both of them, depending on whether we can solve exact black hole solutions.
In the near horizon region, the metric functions and the vector fields can be expanded as Taylor series of the form where r 0 denotes the horizon radii and we have set h 1 = 1 owing to the scaling symmetry of the time coordinate. It should be emphasized that unlike the A r = 0 case, for the solutions with a non-vanishing A r the finite norm condition of the vector is insufficient to govern A t vanishes on the horizon. We find that σ = 1 when a 0 = 0 and σ = 1/2 when a 0 = 0.
Both cases are allowed by the equations of motions. Substituting the expansions into the equations of motions, we find that for the minimal theory β = 0 = γ and a certain nonminimal theory with β = γ/2, there are either three independent parameters (r 0 , a 0 , a 1 ) when a 0 = 0 or two parameters (r 0 , a 1 ) when a 0 = 0 on the horizon. In these two cases the coefficient f i , b i are completely fixed because the metric function f and the vector field However, in spite of that a nonzero a 0 is compatible with the equations of motions, it leads to a divergent local diffeomorphism invariant of the vector Aā = E μ a A µ on the horizon, where E μ a is the inverse vielbein. This is something that we do not appreciate 2 and we will not discuss this case further in the remaining of this paper.
The general structure of the asymptotic solutions at infinity heavily depends on the five parameters of the theories as well as the asymptotical structure of the space-times. Here we shall not analyze them in a case-by-case basis since most of the solutions that we obtain contain all the independent integration constants. Nevertheless, it is deserved to show some universal aspects of the general asymptotic solutions. We find where the effective cosmological constant is parameterized by Λ eff = − 1 2 (n−1)(n−2)g 2 , k eff is a function of (k , µ , q 1 , q 2 ) and in general k eff = k (we call it the effective curvature of the codimension-2 space). It is clear that the asymptotic solutions are characterized by three independent integration constants µ , q 1 , q 2 which are associated with the black hole mass and the vector charges respectively. However, only two of the three parameters are truly independent since the boundary conditions on the horizon provide an algebraic constraint for the three parameters. For example, we may take the parametric relation by saying q 1 = q 1 (µ , q 2 ). Then the full solutions are characterized by two independent parameters µ , q 2 , which are analogous to the case of a Reissner-Nordstrøm (RN) black hole.

Wald formalism and thermodynamics
In this paper, we will adopt the Wald formalism to derive the first law of thermodynamics for all the solutions we obtain. The Wald formalism provides a systematic procedure for the derivation of first law of thermodynamics for the solutions of a generic gravity theory.
It was first developed by Wald in [23,24]. Variation of the action with respect to the metric and the matter fields, one finds where φ collectively denotes the dynamical fields and E φ = 0 are the equations of motions.
For our gravity model, the current J µ receives contributions from both the gravity and the vector. We find where G µνρσ is the Wheeler-Dewitt metric, defined by and the current associated with the γ coupling term is Note that we have put the current associated with the non-minimally coupled terms into the vector sector. For a given current J µ , one can define a current 1-form and its Hodge dual as When the variation is generated by an infinitesimal diffeomorphism ξ µ = δx µ , one can define an associated Noether current (n − 1)-form as where i ξ · denotes the contraction of ξ with the first index of the n-form * L it acted upon.
It was shown in [23,24] that the Noether current J (n−1) is closed once the equations of motions are satisfied, namely where e.o.m denotes the terms proportional to the equations of motions. Thus one can further define a charge (n − 2)-form as It was shown in [23,24] that when ξ is a Killing vector, the variation of the Hamiltonian with respect to the integration constants of a specific solution is given by where C is a Cauchy surface, Σ n−2 is its two boundaries, one on the horizon and the other at infinity. For our vector-tensor theories, it is straightforward to derive the various quantities in the Wald formalism though the calculations are a little lengthy. For pure gravity, we For the vector sector, we obtain Notice that the Wald formalism does not explicitly depend on the non-derivative terms of the Lagrangian density. The various quantities have been given in [13] for γ = 0 and in [25,26,27,28,29,30] for β = 0 = γ.
Now we evaluate δH for the general static solutions with maximal symmetry (9). Let and where ω n−2 is the volume factor of the (n − 2) dimensional space and It was shown in [23,24] that evaluating δH on the horizon yields where the temperature and Wald entropy are given by Here κ is the surface gravity on the horizon. Throughout this paper, the Wald entropy is always denoted by S, without any subscript. For our metric ansatz, we have where A = ω n−2 r n−2 0 is the area of the horizon. Evaluating δH at both infinity and on the event horizon yields Thus the first law of thermodynamics is simply This is the standard derivation of the first law when the Wald entropy formula holds.
However, the situation in our case is even more subtle because counterintuitively, the finite norm condition of the vector is not sufficient to govern the validity of the Wald entropy formula (29). The reason is δH + may be non-integrable for general near horizon solutions.
As was discussed in [31], to govern the validity of Wald entropy formula one should require the local diffeomorphism invariant Aā of the vector vanishes on the horizon 3 . However, for our vector-tensor theories such a condition in general turns out to be too strong to be imposed because A r is an algebraic degree of freedom which does not have corresponding vector charges. So we have to relax the condition for Ar(r 0 ) and simply demand a vanishing At(r 0 ). This has fixed the parametric relation between the parameters (µ , q 1 , q 2 ) of the asymptotic solutions but it does not necessarily lead to a vanishing Ar(r 0 ). Consequently, δH + in general becomes non-integrable. We find where Φ ≡ A 2 r (r 0 ) is a dimensionless quantity. The existence of the non-integrable one-form on the r.h.s of the equation invalids the Wald entropy formula as well as a refining entropy defined as δH + ≡ T δS re . Nonetheless, formally one can still write down a "first law" using the Wald equation despite that its physical meaning is unclear. Notice that when γ = 0, one will not encounter the trouble any longer because of δH + = T δS. Furthermore, when β = 1 2 γ, δH + is integrable as well 4 , given by where the improved temperature and entropy are defined by Here comes an intriguing question that how the improved temperature T is interpreted in the thermodynamical content. We leave this as a future direction for research. For generic case, as will be shown later, the above non-integrable one-form may vanish for a certain coupling constant γ.

Minimal theory
In this section, we study a minimally coupled theory described by which generalizes the Einstein-Proca theory with a quartic self-interaction term for the vector field. Despite the simple form of the theory, there are some new interesting and important features in the theory. For instance, although for Λ 0 = 0 the maximally symmetric vacuum of the theory is (A)dS space-times, it also allows a simple solution which is Minkowski space-times supported by a constant vector provided the parametric relation Note that this relation leads to a perfect squared vector potential V = −γ 4 A 2 + m 2 4γ 4 2 and the parameters in the solution (36) are such that V = 0. Reality of the solution naturally requires q 2 1 ≥ m 2 4γ 4 , where the "=" case corresponds to a vanishing A r , which was first studied in [13]. It is worth emphasizing that the above solution (36) is not a vacuum solution because the vector breaks the gauge symmetries explicitly. When the bare cosmological constant deviates from the critical value (37), an effective cosmological constant emerges in the corresponding solutions because the potential now becomes A second new and probably more important feature of the theory (35) is that we can exactly solve its general static maximally symmetric black hole solution with A r = 0.
This is quite surprising since up to now any exact black hole solution has not been found in Einstein-Proca theory. To keep generality, let's discuss how to analytically solve the equations of motions for general parameters.
First, the equations P r and E rr are purely algebraic for A r and f so we can solve the two functions in terms of h , A t and their derivatives The remaining independent equations are E tt and P t , which are second order non-linear ordinary differential equations (ODE) of A t and h. They are in general very difficult to integrate. Fortunately, we find that the two equations become integrable if we parameterize the two functions as where µ is an integration constant. Strikingly, the equation E tt simplifies to a single linear which can be immediately solved by where C = 0 is a new integration constant. Substituting (43) into P r , we find it also reduces to a linear first order ODE which is easy to integrate. We get where q 2 is an integration constant associated with the Coulomb-like charge of the vector, as will be shown later. Plugging Eq.(43) and Eq.(45) into (41), we obtain where the effective cosmological constant (38) is parameterized by Λ eff = − 1 2 (n−1)(n−2)g 2 and q 1,2 are the two vector charges which are analog of the chemical potential/charge density of Reissner-Nordstrøm (RN) black hole. Now it is clear that C is a non-physical parameter which is associated with the scaling symmetry of the time coordinate. Without loss of generality, we set C = 1. Finally, substituting (46) into (40), we find This completes our derivation. To conclude, we obtain Now we are ready to give some comments on the solution. First, as is clear from the derivation, the above solution is the most general static solution with maximal symmetries in the presence of a non-vanishing A r . It contains all the three integration constants µ , q 1 , q 2 which are associated with the black hole mass and the vector charges respectively. Second, formally without A r the solution is simply the RN black hole! In fact, this is easily understood because under the special ansatz h = f , A r = 4γ 4 such that V = 2Λ eff and the stress tensor of the vector reduces to that of a Maxwell: T A µν = T Max µν . In this sense, we may call the solution a stealth Reissner-Nordstrøm black hole (but one should remember the gauge symmetry is breaking).
Moreover, evaluating δH at infinity yields where the black hole mass M is defined by 5 The Wald equation (31) implies the first law is where the temperature and entropy are given by In addition, the Smarr relation is where the thermodynamic pressure P and volume V are defined by as the counterparts of the thermodynamic conjugate (Φ e , Q e ) of a RN black hole. In the latter case, the non-integrable term Φ e dQ e associated to the Maxwell field in the first law is well understood as the working term of the electrostatic force. Thus, in this case the black hole mass is well defined via the first law by using the Wald entropy. Here comes an 5 In this paper, without specification, we always define the mass using the usual fall-off mode 1/r n−3 associated with the condensate of the massless gravitons. It is the standard ADM/AMD mass of the asymptotically flat/AdS solutions for the minimal theory. However, for non-minimal theory the situation is not so simple. Nonetheless, we continue using the concept for all these theories.
intriguing question: for our vector field whether the term Φ p dQ p appearing in the first law Eq.(51) can be interpreted as the work of the force associated to the vector field as well. We find surprisingly, the answer is yes! The reason is for our solution the vector field enjoys a detailed balance condition A 2 = −m 2 /(4γ 4 ) such that its self-interaction terms on the r.h.s of the equation of motion Eq.(4) are exactly cancelled. Thus, in this sense the vector charge Q p is globally conserved and the mass in addition to the ADM definition is well defined via the first law by making use of the Wald entropy. In addition, we also find that the absence of a naked curvature singularity at the origin leads to an upper bound for the vector charge where the bound is saturated for an extremal solution which is asymptotically flat.
Third, in general the limit A r → 0 is not allowed except for the asymptotically flat solution Note that the limit A r → 0 gives rise to an extremal solution [13] because of f ∼ A 2 t . Finally, the reality of A r at any position of the space-times strongly constraints the parameters of the solution (here we do not clearly distinguish the integration constants of the solution and the coupling constants of the theory). We find Note that the limit of a free vector field γ 4 → 0 is not well defined for our solutions.
4 Non-minimal theory: case I From now on, we turn to study the generalized Einstein-Proca theories with non-minimally coupled terms. In this section, we study a simple theory which has a single γ term, namely whilst the theory with a single β term was studied in [13]. It turns out that for this simple theory, we can find different kinds of stealth black hole solutions which satisfy G µν = 0 = T µν , depending on the non-minimal coupling constant.

Stealth black hole: γ = 1
The first case we consider is when γ = 1, we always have T µν = 0 for the special ansatz h = f , A t ∝ f , A r = 0. Therefore, we easily find a stealth Schwarzschild black hole solution where q 1 , µ are two independent integration constants associated with the vector charge and the black hole mass. The solution can be trivially generalized to including a cosmological constant, additional matter fields (such as a Maxwell field) or higher curvature terms (such as the Love-Lock terms ) in the Lagrangian density.
However, the solution is not most general since it contains one less integration constant.
To derive its first law of thermodynamics, we shall first analyze the structure of general asymptotic solutions and derive the corresponding first law. We find It is easy to see that in general h = f and the mass of the black hole does not receive contributions from the back-reaction of the vector. Substituting the asymptotic solutions into the Wald formula, we obtain where the mass is defined by (50). It follows that the first law reads It is interesting to note that there are two new pairs of thermodynamic conjugates: (q 1 , q 2 ) and (µ , q 2 1 ). This is very different from the Einstein-Proca black hole [25] which only has the first pair of conjugates. For the special solution (60), we have q 2 = 2µq 1 and It is straightforward to verify the first law (63) is indeed satisfied. In addition, define a new energy function as the first law can be cast into the form of So we may take E as a function of E = E(S , q 1 ). Note that the coupling constant γ is dimensionless so the above equation contains all the dimensionful quantities in the theory and the solution. It follows that using the scaling dimensional arguments we can derive a Smarr relation Written back in terms of the original mass function, we find This can be easily verified although the solution (60) is degenerate.

The solution and thermodynamics
Interestingly, for generic γ = 0, we can also obtain a stealth Schwarzschild black hole solution which has a nonzero A r . First, the vector equation P r reads while the Einstein equation E rr is no longer algebraic for A r or the metric functions. To proceed, we choose a special ansatz h = f . Then the vector equations dramatically simplify to y ′′ + n−2 r y = 0 , which can be immediately solved as The remaining equations are E tt , E rr , both of which are ODE of A r . Here, one may worry about the two equations are inconsistent with each other. Fortunately, we are able to find an unique solution for A r which satisfies both equations This completes our derivation. In spite of that we do not expect to find the most general solutions of the theory at the very start, the solution we get contains all the three independent integration constants µ , q 1 ,2 . This gives us strong confidence that the above solution is the general static spherically symmetric black hole solution with a nonzero A r . To govern the reality of A r at any position of the space-times, the coupling constant is bounded For γ = 1, we always have A 2 r ≥ 0 so the parameters µ , q 1 , q 2 are free in this case. For γ = (n−3) 2(n−2) , the condition leads to µq 1 (µq 1 − q 2 ) ≥ 0. For generic case, the constant terms in the square bracket of (73) should be nonnegative.
For later convenience, we list the above solution as follows Evaluating δH at infinity yields where the refining mass M is defined by Then the Wald equation suggests the first law is For our solution (75), the temperature and Wald entropy are given by However, by plugging these results into the first law (78), we find that it picks out a special coupling γ = 1 when we impose the boundary condition that Aā vanishes on the horizon. In fact, relaxing the horizon condition for Ar, we find in general δH + becomes non-integrable This is a correct mathematic equation though its physical meaning is not so clear. Define a Legendre transformed energy function the above first law simplified to In addition, we also find a Smarr-like relation which is a natural result of the scaling dimensional arguments.

Euclidean action
Since the non-integrability of δH + invalids the Wald entropy formula, we shall explore whether there exists an alternative approach to define the black hole entropy. It was first proposed in [32] that thermodynamic quantities for black holes can be calculated by means of quantum statistical relation: where F is the free energy, I reg is the regularized Euclidean action of black hole solutions and M , S are black hole mass and entropy, respectively (they should not be confused with the mass and entropy defined from Wald formalism.). The regularized Euclidean action can be defined by subtracting the action of a background solution with µ = 0 from the action of the black hole, However, for our solution (75) the limit µ → 0 is singular for a generic coupling γ. Instead, we derive a proper background solution by taking double scaling limit: µ → 0 , q 2 → 0 with q 2 2 /µ → const. The resulting expression for free energy is very simple In usual cases (such as a Schwarzschild black hole), one can derive both the mass M and entropy S independently as  Non-minimal theory: case II Now let's consider the general non-minimally coupled theory described by The maximally symmetric vacuum is Minkowski space-times. However, the Lorentz symmetry of the vacuum can break down because a constant vector is admitted as well, namely Depending on the coupling constants, we find that there exist significantly different classes of asymptotically flat black hole solutions.

Unconventional black hole
The first class solution we find has an unconventional fall-off at asymptotic infinity. It reads provided the parametric relation Here µ is an integration constant which should not be confused with the usual fall-off mode 1/r n−3 . This type solution was first found in [11] for β = γ/2 and in [13] for γ = 0. Note that the reality of the vector requires γ < 1 and the limit γ → 1 or equivalently β → 0 is singular. In fact, the solution does not exist in the theory (59) which has a single γ term.
As was shown in [13], the unusual fall-off 1/r n− 3 2 in the metric functions corresponds to the longitudinal graviton mode, which is excited by the back-reaction effect of a background vector.
Since the solution has only one integration constant, we shall first analyze the general asymptotic solutions of the theory for generic coupling constants before deriving the first law. Linearing the equations of motions around the background (90), we find where σ > 0 is an under-determined constant. The conventional solution has σ = n − 3 and µ = 4βq 1 q 2 + (1 − βq 2 1 )µ 1 + βq 2 1 . (94) Thus, the general solution is characterized by three independent parameters (µ , q 1 , q 2 ), as expected. However, the unconventional solution with σ = n − 3 also exists provided and the vector charge q 1 has been fixed as a function of (n , β , γ) (the details is irrelevant in our discussion). Hence, one may worry about the existence of this type solutions since it needs a delicate fine tuning of the boundary conditions on the horizon 6 . It is interesting that we do find such a solution when σ = (n − 3)/2 (one can check that in this case µ = µ and the linearized solution (93) is just the exact solution (91) ).
To understand the solution (91) better, we develop its full large-r expansions. We find 7 where m is a new parameter associated with the condensate of the transverse gravitons. The higher order coefficients a i , h i , f i can be solved in terms of functions of the two parameters (µ , m). Substituting (96) into the Wald formula, we obtain Surprisingly, it does not receive any contributions from the transverse graviton mode! This is very different from the conventional black holes. We define the mass as δM uc ≡ δH ∞ , giving rise to For our solution (91), the temperature and entropy are given by It follows that the first law and Smarr relation are given by To end this subsection, we point out that in the weak field limit the unconventional solution (91) predicts a stronger gravitational force than the Schwarzschild black hole. For example, in the four dimension it has 1/r 3/2 -law rather than the well-known 1/r 2 -law. More interestingly, since the general asymptotic solution (96) has the usual fall-off mode 1/r n−3 as well, one can turn on or turn off the unconventional mode freely. This gives rise to new possibilities and candidates how the Newtonian inverse-squared law can be modified in galaxies and may be tested by observational data in astrophysics in the future.

Stealth black hole and beyond
Following the derivation in sec.4.2, we find that there exists an exact stealth Schwarzschild black hole at the critical coupling γ = n−3 2(n−2) whilst β remains free. The solution reads We demand q 1 (µq 1 − q 2 ) ≥ 0 to govern the reality of A r . Formally, the solution is simply the one (75) when γ = n−3 2(n−2) . Likewise, its first law can be studied along the discussions in sec.4.2.
Surprisingly, if we fix the vector charge q 1 = 1/ √ γ − β, a new fall-off mode emerges in the metric function without altering anything else in the above solution. We find where λ is a new independent integration constant. The coupling constant γ is still equal to n−3 2(n−2) and β < γ , β = γ/2. For β = 0 , ξ = 0 and for β = − 1 2 (n − 3)γ , ξ = n − 3. In both cases, the new mode is trivial and can be dropped in the metric function. When γ/2 < β < γ or β < − 1 2 (n − 3)γ, we find ξ > (n − 3), the λ mode falls off faster than 1/r n−3 whilst for − 1 2 (n−3)γ < β < 0, we have 0 < ξ < n−3, the new mode falls off slower than the conventional one. Moreover, when 0 < β < γ/2, we have ξ < 0, implying that the solution is no longer asymptotically flat, although the maximally symmetric vacuum of the theory is Minkowski space-times. In particular, when β = γ/n, we have ξ = −2, the solution becomes asymptotically (A)dS and the cosmological constant emerges as an integration constant, which is totally independent of the parameters in the Lagrangian density. It may be the first time to observe this phenomena in Einstein gravity except for conformal gravity.
Finally, it should be emphasized that the existence of such a new mode is peculiar for the theory we consider because the presence of both non-minimal couplings is essential to govern the existence of this type solution.
6 Non-minimal theory: case III Now we study a certain non-minimally coupled theory which has β = γ/2 and also includes a bare cosmological constant and a bare mass term The theory has been extensively studied in [11,12,14,16] but most of them are limited to the four dimension.

Without Proca mass
First, let's consider a simpler case m 2 = 0. From the discussions in above section, it is immediately to see that when the bare cosmological constant also vanishes, the asymptotically flat solutions (91) and (101) are still valid with the coupling constant γ = (n−3) 2(n−2) in both cases.
For the same coupling constant, the solution (101) can be generalized to non-asymptotically flat space-times when Λ 0 = 0. We obtain 8 where various parameters are specified by Here Λ p is related to the bare cosmological constant .

With Proca mass
With a nonzero Proca mass, it is of great difficult to solve exact black hole solutions in the theory (103). Interestingly, in [14] the author found some exact solutions for certain coupling constants. Furthermore, in [16] the authors developed a nice procedure to derive the general solution in the four dimension. Here we follow the discussions in [16] and generalize the method to general dimensions.

Derivation of the solutions
A neat observation in [16] is that the equations P r and E rr are purely algebraic for the metric function f and the vector field A r . Hence, they can be solved in terms of other functions and their derivatives immediately The two equations (110) encode some universal information about the general solution.
For example, requiring the metric functions behaves standard at asymptotically AdS spacetimes, namely at leading order h = f = g 2 r 2 + · · · at infinity, we find the effective cosmological constant should be proportional to the Proca mass squared where the effective cosmological constant is parameterized by Λ eff = − 1 2 (n−1)(n−2)g 2 . On the other hand, the non-negativity of A 2 r strongly constrains the parameters in the theory as well as those in the solution. For instance, at asymptotic infinity, A 2 r behaves as (at leading order) which implies that Λ 0 should not be bigger than Λ eff , namely These results are universal for the general solution.
To proceed our derivation, we parameterize the metric function h and the vector field It turns out that the remaining two independent equations P t , E tt are integrable for H , F at the critical coupling constant γ = (n−3) 2(n−2) . First, the vector equation P t simplifies to which can be solved immediately as where C 1 is an integration constant. Its physically meaning will be explained later. Then the Einstein equation E tt reduces to a linear first order ODE of H which is easy to integrate (we do not list it in the following due to its lengthy expressions). We obtain where C 2 is a new integration constant. Substituting (116) and (117)  The results depend on the topological parameter k as well as the space-time dimension n.
In the following, we will discuss the solutions in a case-by-case basis. Before this, we point out the above two integration constants C 1 , C 2 are related to the vector charge q 1 and the effective cosmological constant, respectively. We find Something else is we will not list the explicit expressions for f and A r since they are lengthy and we have too many solutions in the following. The readers in need can easily find them through the equation (110).

Planar black holes
For simplicity, let us consider planar black hole solutions at first. For n = 5 dimensions, we The event horizon is determined by the largest real root of h(r 0 ) = 0 and its effective curvature is nonzero, given by The black hole mass M is given by Eq.(50). Notice that the solution (119) is singular in the n → 5 limit. As a matter of fact, a logarithmic term emerges in the original fall-off mode 1/r 2 . We find For above solutions, there is an another singular limit Λ 0 → −Λ eff = (n−2)m 2 n−3 . However, in this case we do not find any physically interesting solutions. Notice that the above solutions do not have a regular limit for a vanishing Proca mass m 2 → 0. Thus, one cannot find planar black holes when m 2 = 0.
Evaluating the Wald formula for above solutions in general n ≥ 4 dimensions, we obtain where M is defined by As discussed earlier, when β = 1 2 γ, δH + is given by Eq. (33-34). Then the Wald equation implies the first law is For our solutions, the temperature T and the quantity Φ are given by where h 0 is the effective curvature of the horizon which in general is not equal to unity.
For even dimensions, the black hole mass is given by Eq.(50) whist for odd dimensions which receives contributions from the vector field as well.
However, in the n = 5 and n = 7 dimensions, the solution (126) becomes singular. As a matter of fact, for the n = 5 dimension we find The black hole mass is given by (129) with In the n = 7 dimension, we find A t = q 1 + Notice that for all these solutions the limit m 2 → 0 is singular. Thus, one cannot recover the solution (104) by simply sending m 2 → 0 from these solutions. Instead, one should follow the derivation in sec.6.2.1 (there the integration constants C 1 ,2 in (118) become singular in the limit m 2 → 0 and should be chosen properly again) and then re-derive the solution (104).
Finally, for all the solutions above, the first law is still given by Eq. (124)  .
This is expected since the Wald formula given in Eq. (24)(25)(26) do not explicitly depend on the topology of the space-times.

Conclusions
In this paper, we study generalized Einstein-Proca theories in general dimensions by introducing either a quartic self-interaction term for the vector or non-minimally coupled terms between the curvature and the vector. In general, the gauge symmetry of the vector is explicitly breaking but can be restored at the linear level around any Ricci-flat metric, depending on the parameters of the theories.
We find that there are two distinct class solutions, both of which are general static and have maximal symmetry, depending on whether A r vanishes. In particular, the solutions with a nonzero A r have some attractive features that we do not find for the solutions with a vanishing A r . The first is in many cases (for example the minimal theory which is simply the Einstein-Proca theory extended with a quartic self-interaction term for the vector field) we can analytically solve all the equations of motions and exactly obtain the general static maximally symmetric black hole solutions. This is quite surprising since it is known that the Einstein equations are highly non-linear and one has not found any analytical solutions in the standard Einstein-Proca theory. The underlying reason is the equations P t , E tt are purely algebraic for the metric function f and the vector component A r , which can be immediately solved in terms of functions of h , A t and their derivatives. It turns out that under certain parametrizations of h , A t , the remaining independent equations P t , E tt are greatly simplified to first order ODEs, which are easy to integrate.
Second, the reality of A r provides strong constraints on the parameters of the solutions as well as those in the Lagrangian density. Third, we adopt the Wald formalism to derive the first law of thermodynamics for all of the solutions. However, the situation is subtle for the solutions with a nonzero A r because to govern the validity of Wald entropy formula, we need impose proper boundary conditions that the local diffeomorphism invariant of the vector Aā vanishes on the horizon, which unfortunately turns out to be too strong for this type solutions. The reason is A r does not have corresponding vector charges since it is a purely algebraic degree of freedom. Thus, we have to relax the horizon condition for Ar but this conversely results to a non-integrable δH + which invalids the Wald entropy formula.
The thermodynamics of such solutions deserves further studies.
Finally, we also obtain some exact black hole solutions with vanishing A r . In particular, one of the solutions has an unconventional fall-off mode, which is interpreted as the longitudinal gravitons excited by the vector field. In the weak field limit, the solution has a stronger gravitational force than the usual Newton's 1/r 2 -law. This is particularly interesting in astrophysics since it provides new candidates to modify the Newton's inverse-squared law.