arxiv: v2 [gr-qc] 6 Jan 2018

Transcription

1 Smarr Formula for Einstein-Æther-Maxwell Theroy Fei-hung Ho a,c, Shao-Jun Zhang a, Hai-Shan Liu a,, and Anzhong Wang a,b a Institute for Advanced Physics & Mathematics, Zhejiang University of Technology, Hangzhou, 3003, China b GCAP-CASPER, Physics Department, Baylor University, Waco, Texas , USA c Department of Physics, Jimei University, Xiamen, 360, China Dated: January 9, 08) arxiv: v [gr-qc] 6 Jan 08 Using the Wald formalism, we investigate the thermodynamics of charged black holes at both universal and Killing horizons in the Einstein-æther-Maxwell theory. In particular, using the existence of a scaling symmetry of the action, we obtain the Smarr integral formula. Assuming its canonical form among the mass, temperature and electrical potential, we find that the entropy in general takes a very complicated form, and is a function of various parameters. Only in very particular cases, it is given by one fourth of the area of the horizon. I. INTRODUCTION JacobsonandMattinglyintroducedaunittime-like4-vectorfieldu a, theso-called ætherfield coupledtoeinstein s general relativity GR), which is often referred to as the Einstein-æther theory []. The existence of the æther field defines a preferred frame, whereby the Lorentz invariance is locally broken and leads to some novel effects: superluminal gravitational modes with different speeds exist []. Due to the existence of the superluminal gravitational modes, the causality of the theory will be quite different from that of GR. In particular, since the speeds of these particles in principle can be arbitrarily large, so the corresponding light-cones can be completely flat, and the causality of the theory is more like that of Newtonian theory [3, 4]. Then, some natural questions rise: does black holes still exist? how to define a black hole now? To answer these questions, we need first to understand the corresponding causality of the theory. In Newtonian theory, an absolute time t is first introduced, and the causality is guaranteed by assuming that all particles move in the increasing direction of t. With this in mind, Blas and Sibiryakov first realized that the time-like æther field can play the role of the absolute time of Newton [5]. In particular, the æther field naturally defines a globally timelike foliation of a given spacetime, and all the particles are assumed to move alone the increasing direction of the foliation. Once the causality of the theory is defined, a universal horizon can be defined as a one-way membrane of particles with arbitrarily large speeds, that is, a surface that traps all particles: once they get in, they shall be trapped inside of the horizons all the time, even they have infinitely large speeds. In the decoupling limit, Blas and Sibiryakov showed that such a surface indeed exists [5]. Since then, universal horizons have been intensively studied [6]. In particular, to study universal horizons in gravitational theories without the presence of the æther field, one way is to promote the æther field to play the role as that of a Killing vector [7]. With the establishment of the existence of universal horizons in static and stationary spacetimes [6], another natural question is: do universal horizons have a thermal interpretation, similar to Killing horizons in GR? In the framework of the Einstein-æther theory, it was argued that universal horizons should possess such properties, and showed explicitly that the first law of black hole mechanics indeed exists for static and neutral universal horizons [8], provided that the definition of the surface gravity at the universal horizon is replaced by, κ peeling = d dr dr dv = d v) s r = κ UH UH dr s UH uu χ),.) UH which was first obtained by considering the peeling behavior of a null ray near universal horizons [9], and later was shown that the above relation holds for all static universal horizons [0], where u µ is the four-velocity of the æther field, and s µ is a spacelike vector orthogonal to u µ and tangent to the universal horizons. The quantities r and v The corresponding author Electronic address: fei.h.ho@gmail.com Electronic address: sjzhang84@hotmail.com Electronic address: hsliu.zju@gmail.com Electronic address: Anzhong Wang@baylor.edu

2 denote the Eddington-Finkelstein coordinates, in terms of which the static spacetimes are described by the metric [0], ds = Fr)dv +fr)dvdr +r dσ k,.) where k = 0,±, and dθ +sin θd, k =, dσ k = dθ +d, k = 0, dθ +sinh θd, k = -..3) It is interesting to note that universal horizons always exist inside the Killing horizons in the well-known Schwarzschild, Schwarzschild-de Sitter and Reissner-Nordström spacetimes [0], which are also solutions of Hořava gravity [4]. It must be noted that it is still an open question how to generalize such a first law to charged and/or rotating universal horizons [ 5]. In addition, in the case of a neutral universal horizon, it was showed, via the tunneling approach, that the universal horizon radiates as a blackbody at a fixed temperature, even if the scalar field equations alsoviolatelocal Lorentzinvariance[6], wherethe temperaturet UH satisfiesthe canonicalrelationt UH = κ UH /π). However, if the Lorentz invariance is also allowed to be broken in the matter sector, such as in the case of Hořava gravity [4, 7], the dispersion relation of a massive particle will contain generically high-order momentum terms, E = m +c p p + z ) n= ) n p a n,.4) M where E and p are the energy and momentum of the particle considered, and c p and a n s are coefficients, depending on the species of the particle, and z is the dynamical exponent and denotes the power of the leading terms in the ultravioletuv) regime of the dispersion relation. In Hořava gravity, the renormalizability of the theory requires z d, where d is the spatial dimension of the spacetime, and M is the suppression energy scale of the higher-dimensional operators. Then, it was shown that the universal horizon radiates still as a blackbody at a fixed temperature, but remarkably now it depends on z [], T z UH = z ) z κuh π ),.5) where κ UH is the surface gravity defined by Eq..). When z = we obtain the results of [6], which considered only the case z =. Recently, more careful studies of ray trajectories showed that the surface gravity for particles with a non-relativistic dispersion relation.4) is given by [8], κ z UH = z ) z ) κ UH..6) Similar results with a factor difference) were also obtained in [9]. In the study of black hole thermodynamics, specifically the first law, the energy of the system is usually identified as the ADM mass, and the entropy term is identified as the quasi-)local energy. This was done initially by Brown and York, who applied the quasi-local energy to the black hole thermodynamics [0], whereby the first law was obtained. Soon, Wald formulated it by using the Noether theorem []. Later, a new formula of quasi-local energy was proposed, and then applied to the study of black hole thermodynamics []. From the above investigations, it can be seen that the quasi-local energy plays a key role in the understanding of the black hole thermodynamics [3]. Applying the above conceptions to black holes in Einstein-æther theory, Eling first considered, respectively, the Einstein and Landau-Lifshitz pseudo-tensors [4], while Foster investigated the same problem by using the Wald formulasim [5]. From their expressions, one can deduce the ADM mass in asymptotically flat spacetimes with a modification due to the presence of the æther field. Recently, along the same line, Pacilio and Liberati obtained an integral form of the Smarr formula by assuming a scaling symmetry in Einstein-æther theory [4]. They also showed that the cosmological constant term can be also included, although it breaks the scaling symmetry. An important aspect of the form is that it is also valid for spacetimes with rotations. In this paper, our goals are two-fold: ) we would like first to generalize the integral Smarr formula obtained by Pacilio and Liberati in [4] to include an electromagnetic field. ) Using such obtained formula, we study charged black holes obtained in the Einstein-æther-Maxwell theory in [, 3] by paying particular attention to the first law, as mentioned above, the first law of charged black holes has not been properly formulated, yet. The rest of the paper is organized as follows. In Sec. II we give a brief review over the Einstein-æther-Maxwell theory, while in Sec. III, applying Wald s formalism, we derive the integral form of the Smarr formulas for the

3 Einstein-æther-Maxwell theory. In Sec. IV we apply such obtained Smarr formulas to study thermodynamics of both Killing and universal horizons of the two classes of exact solutions of the Einstein-æther-Maxwell theory obtained in [] for spherically symmetric spacetimes. In Sec. V we summarize our main results and present some general remarks and discussions. 3 II. EINSTEIN-ÆTHER MAXWELL THEORY The Einstein-æther theory minimally coupled with an electromagnetic field is described by the action [], S = d 4 x [ ] ) g R Λ+L æ +L M,.) 6πG æ where R is the Ricci scalar, and Λ is the cosmological constant. The æther Lagrangian L æ is given by where the tensor Z ab cd is defined as [] L æ = Z ab cd a u c ) b u d )+λg ab u a u b +),.) Z ab cd = c g ab g cd +c δ a cδ b d +c 3 δ a dδ b c c 4 u a u b g cd,.3) where c i i =,,3,4) are coupling constants of the theory. The Newtonian constant G N is related to G æ via the relation G æ = c 4 /)G N [4 6]. There are a number of theoretical and observational bounds on the coupling constants c i [, 7, 8]. Here, we impose the following constraints, 0 c 4 <, +c 3 +3c > 0, 0 c 3 <,.4) where c 4 c +c 4, and so on. Another useful expression of L æ is related to the irreducible decompositions of the æther field [9], where h ab g ab +u a u b, and a u b = 3 θh ab +σ ab +ω ab u a a b,.5) θ u, a a u u a σ ab a u b) +u a a b) 3 θh ab, ω ab [a u b] +u [a a b],.6) where X X b b. Note that, if the æther satisfies the hypersurface-orthogonal condition, u [a b u c] = 0, then the twist ω ab vanishes identically. Then, we find that where The source-free Maxwell Lagrangian L M is given by L u Z ab cd a u c ) b u d ) = 3 c θθ +c σ σ +c ω ω c a a,.7) c θ c +c 3 +3c, c σ c +c 3, c ω c c 3, c a c +c 4..8) L M = F ab F ab, F ab = a A b b A a,.9) where A a is the four-vector potential of the electromagnetic field. To proceed further, let us introduce the tetrad, u a,s a,m a,n a ), where the space-like unit m a and n a lie on the tangent plane of the two-spheres S, and s a is another space-like vector, which is orthogonal to u a, m a and n a. In terms of these four vectors, the metric takes the form, g ab = u a u b +s a s b +ĝ ab,.0)

4 where ĝ ab m a m b +n a n b is the metric on the unit two-sphere S. Due to the spherical symmetry, any vector V a can be decomposed as, V a = V u)u a +V s)s a..) In particular, the acceleration has only one component along s a and is given by a a = a s)s a. Additionally, any rank-two tensor F ab can have components along the directions of the bi-vectors u a u b, u a s b), u [a s b], s a s b and ĝ ab. Especially, in [8, 30] the following useful relations were given, where ˆK u and ˆK s are given by, a u b = a s)s b u a +K 0 s b s a + ˆK u ĝ ab,.) a s b = a s)u b u a +K 0 u b s a + ˆK s ĝ ab,.3) ˆK u = [m c c +n c c ] u, ˆK s = [m c c +n c c ] s,.4) where superscript u and s represent that the extrinsic curvature operator act on the vector u a and s a respectively. The Maxwell field F ab is usually decomposed into the electric and magnetic fields E a and B a, respectively, E a = F ab u b, B a = eabmn g F mnu b,.5) where e abmn is the Levi-Civita tensor. However, because of the spherical symmetry, we have B a = 0, and then we find F ab = E a u b +E b u a..6) Because the electric field is spacelike, we have E u) = 0 and E a = E s)s a. Thus, F ab = E s)ˆǫ ab, where ˆǫ ab = u a s b + u b s a ). Then, from the vacuum Maxwell equations, we find E s) = Q/r, where r denotes the geometric radius of the two-sphere, and Q is an integration constant, which represents the total charge of the whole spacetime. Therefore, we have F ab = Q r ˆǫ ab..7) The equations of motion, obtained by varying the action.), respectively, with respect to g ab, u a, A a and λ are given by G ab +Λg ab = Tab æ ætab M.8) Æ a = 0,.9) a F ab = 0,.0) g ab u a u b =,.) where G ab is the Einstein tensor, the æther and Maxwell energy-momentum tensors Tab æ and T ab M are given by 4 with T æ ab = c [ a u c ) b u c ) c u a ) c u b )]+c 4 a a a b T M ab = + c X c ab +λu au b g aby c d cu d,.) [ ] 6πG æ 4 g abf mn F mn +F am Fb m,.3) Æ a = b Y b a +λu a +c 4 a u b )a b,.4) X c ab = Y c a u b) u a Yb) c +u c Y ab),.5) Y a b = Zac bd cu d..6)

5 5 III. SMARR FORMULAS FOR EINSTEIN-ÆTHER MAXWELL THEORY For any diffeomorphism-invariant D-form Lagrangian, one can always associate a local symmetry with its corresponding conservation quantities, i.e. a Noether current D )-form J as well as a Noether charge D )-form Q [, 5]. If such a symmetry is represented by a Killing vector field, ξ, then the corresponding Noether current D )-form J is given by where ξ is the Lie derivative along the ξ direction, and J is closed on shell, J[ξ] = Θϕ, ξ ϕ) ξ L, 3.) which implies that corresponding to ξ there is a D )-form Noether charge, Q, In the Einstein-æther Maxwell theory, the Noether charge is given by dj[ξ] = E ξ ϕ 0, 3.) J[ξ] dq[ξ]. 3.3) Q[ξ] = [ a ξ b +u a Y b cξ c +u a Y b c ξ c +Y ab u ξ)+f ab A c ξ c] ǫ ab. 3.4) 6πG For a solution {ϕ} of the theory, when it is invariant along ξ, that is, ξ ϕ 0, we have Θϕ, ξ ϕ) 0. From 3.) and 3.3), we integrate the expression, J[ξ]+ξ L 0, over an hypersurface Σ and its boundary, Σ, respectively, we find that Q[ξ]+ ξ L = ) Σ Σ In [3], it was observed that Eq.3.5) leads to the Smarr formula for sufficiently general diffeoinvariant theories, while in [4] it was shown that scaling symmetry leads the Lagrangian of Einstein-Aether theory to be a total divergence on-shell, i.e. L db, with a D )-form B. Then we can derive the Smarr formula by taking a surface integral under a flow of the Killing vector field ξ. For the action.), the scaling symmetry is g ab γg ab u a γ ua λ γ λ A a γ A a, 3.6) where γ is a constant. Following [4, 3], we find ) R+L u +L M = m u a Y m a u m Y a a F ma A a ). 3.7) D Then, total Lagrangian, L = R+L u +λu )+L M, becomes a total divergence on shell, i.e. L db, where ) B = D 6πG ua Y m a u m Y a a F ma A a )ǫ m. 3.8) Thus, on shell, Eq.3.5) reduces to a surface integral, 0 Σ Q[ξ] ξ B. 3.9) Here, represent the equality when the field equation, E = 0, hold on shell. In [4] the Smarr formulas for the infrared Horava gravity were also studied, and found that they reduce to a surface integral only when the integrals are evaluated in the preferred frame.

6 6 A. The Expression for Q[ξ] ξ B Inserting the expressions of Z ab cd given by.3) and Y a b, X a bc given by.4) into 3.4) and 3.8), we find [ Q[ξ] ξ B = a ξ b c 4 u a a b u ξ)+c 3 u a ξ c b u c) +c c 3 ) a u b u ξ) 6πG D c 4ξ a a b +c 3 u a ξ b u))+ F ab A c ξ c ) ] D Fac A c ξ b ǫ ab. 3.0) Alternatively, in terms of the fluid coefficients of Eq..8), we find [ Q[ξ] ξ B = a ξ b +c σ u a ξ c b u c) [ c3 u a ξ b u)+c a ξ a a b] 6πG D c a u ξ)u a a b +c ω u ξ)ω ab +c σ u ξ)u a a b + F ab A c ξ c ) ] D Fac A c ξ b ǫ ab. 3.) It should be noted that the above expressions hold for any stationary spacetimes, including the ones with rotations. B. Spherically Symmetric Charged Black Holes Restricting the above formulas to static spherically symmetric spacetimes, we first note that the hypersurfaceorthogonal condition is automatically satisfied, so the twist ω ab vanishes identically, while ǫ ab can be written as ǫ ab = ˆǫ ab ǫ. Inserting Eq..7) into Eq. 3.), we find Q[ξ] ξ B = [ a ξ bˆǫ ab +c σ s ξ)s a s b a u b 6πG D [c 3s ξ) u) c a u ξ)a s)] c a u ξ)a s)+c ω u ξ)ω abˆǫ ab + Q r ξ A)ˆǫ ab + ] D A ac cξ bˆǫ )ˆǫ ab ǫ, 3.) whereξ a = u ξ)u a +s ξ)s a andawillbedefinedbelow. Finally, withnulltwistconditionand b u a +u b a a )ˆǫ ab = 0, we obtain [ a ξ bˆǫ ab +c σ s ξ)s a s b a u b c 3 Q[ξ] ξ B = 6πG D 3 D c au ξ)a s)+ Q r When D = 4, the above expression reduces to, Q[ξ] ξ B = 8πG where we had used.) to get the first term, D s ξ) u) ξ A)+ D [u A)u ξ) s A)s ξ)] [ c a )u ξ)a s)+ c σ)s ξ)s a s b a u b + c 3 s ξ) u) Q r[ξ A)+u A)u ξ) s A)s ξ)] )] ǫ. 3.3) ] ǫ, 3.4) a ξ bˆǫ ab = [u ξ)a s) s ξ)k 0 ], 3.5) in r.h.s. of Eq.3.4). With the above expression for Q[ξ] ξ B and by taking the integral of Eq.3.9) over the boundary Σ = Σ Σ BH, Q[ξ] ξ B) Q[ξ] ξ B) ) Σ Σ BH now we are in the position to study thermodynamics of charged black hole horizons either Killing or universal) for asymptotic flat spacetimes found in [].

7 7 IV. SMARR FORMULAS FOR CHARGED BLACK HOLES The metric for the charged Einstein-æther black holes solutions takes the form [], with the æther and time-translation Killing vectors are given, respectively, by ds = er)dv +dvdr+r dσ +, 4.) u a = αr),βr),0,0), u a = eα+β,α,0,0), χ a =,0,0,0), 4.) In asymptotically flat spacetime, the æther field is required to be parallel to the time-like Killing field, u χ) =. In addition, we also have the relations, αr) = s χ) u χ), βr) = s χ), er) = u χ) s χ), 4.3) and the Maxwell potential A a has the following non-vanishing components, u A) = Q α β ), r e s A) = Q β r e, χ A) = Q r [ α β e )u χ) βe s χ) ], 4.4) where A a = u A)u a + s A)s a. The Killing horizon is located at e ) = 0, while the universal horizon is located at [4, 5] u χ) = 0, from which we find that eα + ) UH = ) At the universal horizon, the surface gravity is defined by Eq..), i.e., κ peeling uu χ) = UH UH a s)s χ). 4.6) The surface gravity at the Killing horizon is still defined as in GR, κ inaffinity aχ b ) a χ b ) = u χ)a s) s χ)k 0 KH. 4.7) In [], two classes of exact solutions of the Einstein-æther Maxwell theory were found, depending on the choices of the free parameters c i. For c a = 0, the solutions are given by [], A. Charged Black Hole Solutions with c a = 0 s χ) = r æ r, 4.8) u χ) = r 0 r + Q r + c σ)ræ 4 r 4, 4.9) er) = r 0 r + Q r c σr 4 æ r 4, 4.0) where r 0, Q and r æ are the integration constants. From Eqs.4.5) and 4.9), we find ) r UH = r Q, 4.) r 4 æ = r3 UH c σ r UH r 0 r 0 ), 4.)

8 while from Eq.4.0), we obtain where = r 0 r UH ) +L+ N P + 4Q /r0, 4.3) 4L N L = +P, N = 4Q 3r0, P = /3 I +Q 4 /r0 4) + H 3H I = c ) σ r 4 UH c σ r0 4 r3 UH r0 3, H = J + On the other hand, Eqs.4.6) and 4.7) yield, κ peeling UH = 3 c σ ) κ inaffinity KH = r 0 rkh Q rkh 3 + c σræ 4 rkh 5 3 /3, 4I +Q 4 /r 4 0 )3 +J ) /3, J = 7I 7IQ /r 0 +Q6 /r ) ) ) Q ruh Q ruh Q=0, r UH=3r 0/4 r 0 r KH Q=0, r UH=3r 0/4 8 3r 0 3 c σ ), 4.5) + c σ 7r0 4 c σ 8rKH ) At the universal horizon, we have u χ) = 0, s χ) = ræ/r UH, and then Eq.3.4) becomes, [ c σ )s χ)s a s b a u b + c 3 Σ UH Q[χ] χ B = = G = G Σ UH Q s χ) u)+ 8πG r 3 c σ + c ] 3 )s χ)s χ) r c 3 s χ) r+ Q r [ [ r UH r 0 )+ Q r UH ] ǫ r=r UH ], 4.7) where s a s b a u b = K 0 and u = K 0 + ˆK u with K 0 = s χ),r and ˆK u = s χ)/r. Note that the last equality is obtained from Eq.4.). When Q = 0, we have r UH = 3r 0 /4, and then the above expression reduces to Q[χ] χ B = r ) Σ UH 4G Q=0 At the Killing horizon, we have u χ) = s χ), and the integration of Eq.3.4) over the Killing horizon yields, Σ KH Q[χ] χ B = G = G [ s χ)u χ) r ] )s χ)s χ) r c 3 s χ) r + Q r r= c σ + c 3 [ ] r0 )+ Q Q. 4.9) Note that in the above calculations, we had used the expression, c σ r 4 æ r 4 KH = r 0 + Q r KH at the Killing horizon. It is interesting to note that when Q = 0 the above expression reduces exactly to the expression given by Eq.4.8).

9 When r, we have u χ) =, s χ) = 0, α = and a s) = u χ). Then, we find ) α β e Q Σ [Q[χ] χ B] = G = G c a )u χ)u χ) r c a r u χ) r ) M ADM, 4.0) where M ADM = r 0 / denotes the ADM mass. Generalizing the total energy of the Einstein-æther theory first introduced in [4, 5] to the Einstein-æther-Maxwell theory considered here, we obtain M c ) a M ADM, 4.) which we shall consider as the total mass in the Smarr formulas 3.6). This mass was also used in [8, ]. Without charges, Eling showed that the total mass can be obtained by integrating the Einstein pseudotensor over a sphere at spatial infinity [4], while Foster showed that the same result can be also achieved from the Wald formalism [5]. Here in the Einstein-æther-Maxwell theory, we have shown that the total energy takes the same form as in the uncharged case. At the universal horizon, substituting Eqs.4.7) and 4.0) into Eq.3.6) and assuming that the standard form, holds, where Φ = Q/r UH, we find that where T UH := κ peeling π M = T UH S +ΦQ, 4.), S = 4πr r=ruh UHζ UH, 4.3) 9 ζ UH = 3 c σ) r0 r UH Q r UH ) Q r UH ) Q=0 cσ ). 4.4) r UH=3r 0/4 6 Note that ζ UH is a very complicated function of r 0,r u,q and c i, and in general is quite different from /4, a factor that would be expected at a Killing horizon in GR. In fact, even when Q = 0, it is different from /4. This is similar to the results obtained in [, 4, 5]. It is also very interesting to compare the current results with the ones of uncharged black holes obtained in [8] 3. At the Killing horizon, from Eqs.4.9), 4.0), and the definition given in Eqs.4.6), we find that Eq.3.6) can take the form, where Φ = Q/, and where M = T KH S +ΦQ 4.5) πt KH := κ inaffinity, S = 4πr KH)ζ KH, 4.6) ζ KH = r 0 Q 4 r0 Q + cσr4 æ r 5 KH ) cσ= ) Thus, at the Killing horizon, we can get the desirable factor when c σ = 0, but in general it is not the case. These results again are quite similar to the ones obtained in [, 4, 5]. 3 In [8], the Smarr formula takes M æ = q UH A UH /4πG æ, in which M æ = c 4 /)M ADM, q UH is the linear combination of the corresponding surface gravities and A UH is the area of the universal horizon.

10 0 B. Charged Black Hole Solutions with c 3 = 0 When c 3 = 0, the corresponding solutions are given by [] where r u is given by s χ) = r 0 +r u, u χ) = + r 0 r r, 4.8) er) = r 0 r r ur 0 +r u ) r, 4.9) r u = r 0 c a / 4Q /r 0 c σ. 4.30) In order to have u a well-defined in the whole spacetime, we must keep r u positive, that is Q r 0 c σ c a, c σ c a. 4.3) Then, we find = r 0 +r u, r UH = r 0, 4.3) and the corresponding surface gravity at the universal and Killing horizons are given, respectively, by κ peeling UH = c a r UH cσ Q Q=0 ruh c a /, 4.33) r 0 c σ κ inaffinity KH = + r u r KH Q=0 + r 0 4rKH c a /. 4.34) c σ At the universal horizon, we have u χ) = 0, s χ) = r 0 +r u )/r, r UH = r 0 / and α = /s χ). Then, we find, Q[χ] χ B = ] [ c σ )s χ)s χ) r + Q Σ UH G r r=r UH = [r UH c ] a )+ G Q Q ruh. 4.35) r UH When Q = 0, we have r UH = r0, and the above expression reduces to Q[χ] χ B = c ) a r ) Σ UH 4G Q=0 On the other hand, at the Killing horizon, we have Q[χ] χ B = [ c ] a )s χ)u χ) r c σ )s χ)s χ) r + Q Σ KH G r r= = [ [ c ] ] a r0 )+ G Q Q. 4.37) When Q = 0 we have = r 0 and r UH = r 0 /, and then the above expression reduces precisely to the one given by Eq.4.36) at the universal horizon! At the universal horizon, substituting Eqs.4.35) and 4.0) into Eq.3.6) and comparing the entropy term with Eq.4.33), we find that M = T UH S +ΦQ, 4.38)

11 where Φ = Q/r UH and but now with ζ UH = T UH := κ peeling π c σ ) c ) a Q ruh, S = 4πrUH ζ UH, 4.39) r=ruh Q=0 r UH=r 0/ c σ ) c a ). 4.40) At the Killing horizon, from Eqs.4.37), 4.0), and the definition in Eqs.4.34), Eq.3.6) takes the form, M = T KH S +ΦQ 4.4) where Φ = Q/ and πt KH := κ inaffinity, S = 4πrKHζ KH, 4.4) where ζ KH = c a ) r0 Q + ru ) ) ca=0,cσ=0 Thus, at the Killing horizon, we can get the desirable factor only when c a = 0 and c σ = 0 [, 4, 5]. V. CONCLUSION In this work, we have derived the Wald formalism for the Einstein-æther Maxwell theory, with the fact that a scaling symmetry exists. Applying it to both Killing and universal horizons, we have obtained the Wald formulas at both Killing and universal horizons in the Einstein-æther-Maxwell theory, whereby the Smarr integral formulas were read off for the charged black hole solutions found in []. We have found that the entropy takes a complicated form at both Killing and universal horizons, if we assume that the Smarr integral formulas take their canonical forms []. This is quite similar to neutral black hole solutions in asymptotically anti-de Sitter spacetimes studied recently in [4, 5]. Thus, a desirable Smarr integral formulas are still absent at both Killing and universal horizons in the Einstein-æther-Maxwell theory. Acknowledgments This work was supported in part by National Natural Science Foundation of China NNSFC) with the grant Nos.: A.W.), A.W.), H.S.L.), H.S.L.), and S.J.Z). [] T. Jacobson and D. Mattingly, Phys. Rev. D 63, ); T. Jacobson, Proc. Sci. QG-PH, ). [] T. Jacobson and D. Mattingly, Phys. Rev. D 70, ). [3] J. Greenwald, J. Lenells, J. X. Lu, V. H. Satheeshkumar and A. Wang, Phys. Rev. D84 0) [4] A. Wang, Inter. J. Mod. Phys. D6 07) [5] D. Blas and S. Sibiryakov, Phys. Rev. D 84, ). [6] E. Barausse, T. Jacobson and T. P. Sotiriou, Phys. Rev. D83, ); B. Cropp, S. Liberati, and M. Visser, Class. Quantum Grav. 30, ); A. Mohd, arxiv: ; F.-W. Shu, K. Lin, A. Wang, and Q. Wu, JHEP 04, ); S. Janiszewski, A. Karch, B. Robinson, and D. Sommer, JHEP 04, 63 04); M. Saravani, N. Afshordi, and R.B. Mann, Phys. Rev. D89, ); C. Eling and Y. Oz, JHEP, ); J. Bhattacharyya and D. Mattingly, Int. J. Mod. Phys. D3 04) ; K. Lin, F.-W. Shu, A. Wang, and Q. Wu, Phys. Rev. D 9, ); M. Tian, X.-W. Wang, M.F. da Silva, and A. Wang, Gravitational collapse and formation of universal horizons, arxiv: ; F. Michel and R. Parentani, Phys. Rev. D9 05) 4049; D. Blas and E. Lim, Int. J. Mod. Phys. D3 05) ; Y. Misonoh and K.-i. Maeda, Phys. Rev. D9 05) ; K. Lin, V. H. Satheeshkumar, and A. Wang, Phys. Rev. D93 06) 405; A. Maciel, Phys. Rev. D93 06) 0403; K. Lin, S. Mukohyama, A. Wang, and T. Zhu, Phys. Rev. D95, ).

12 [7] K. Lin, E. Abdalla, R.-G. Cai, and A. Wang, Inter. J. Mod. Phys. D 3, ). [8] P. Berglund, J. Bhattacharyya and D. Mattingly, Phys. Rev. D 85, 409 0). [9] B. Cropp, S. Liberati, A. Mohd, and M. Visser, Phys. Rev. D 89, ). [0] K. Lin, O. Goldoni, M. F. da Silva and A. Wang, Phys. Rev. D9, ). [] C. Ding, A. Wang, and X. Wang, Phys. Rev. D 9, ). [] C. Ding, A. Wang, and X. Wang, and T. Zhu, Nuc. Phys. B 93, ). [3] C. Ding, C. Liu, A. Wang, and J. Jing, Phys. Rev. D 94, ). [4] C. Pacillio and S. Liberati, Phys. Rev. D95, ). [5] S. Liberati and C. Pacilio, Phys. Rev. D 96, ). [6] P. Berglund, J. Bhattacharyya and D. Mattingly, Phys. Rev. Lett. 0, ). [7] P. Hořava, Phys. Rev. D 79, ). [8] C.-K. Ding and C.-Q. Liu, Sci China-PMA 60, ) [arxiv:6.0353]. [9] B. Cropp, Strange Horizons: Understanding Causal Barriers Beyond General Relativity, arxiv: [0] J. D. Brown and J. W. York, Phys. Rev. D 47, 407, ). [] V. Iyer and R. M. Wald, Phys. Rev. D 50, ). [] C.-M. Chen, J. M. Nester and R.-S. Tung, Phys. Lett. A 03, 5 995). [3] L. Szabados, Living Reviews in Relativity, 4 009). [4] C. Eling, Phys. Rev. D 73, ). [5] B. Foster, Phys. Rev. D 73, ). [6] S. M. Carroll and E. A. Lim, Phys. Rev. D 70, ). [7] J.W. Elliott, G.D. Moore, and H. Stoica, JHEP 08, ). [8] K. Yagi, D. Blas, E. Barausse and N. Yunes, Phys. Rev. D89, ). [9] T. Jacobson, Phys. Rev. D 89, ). [30] B. Cropp, S. Liberati and M. Visser, Class. Quantum Grav. 30, ). [3] S. Liberati and C. Pacilio, Phys. Rev. D 93, ).