NAT Black Holes

We study some physical properties of black holes in Null Aether Theory (NAT)--a vector-tensor theory of gravity. We first review the black hole solutions in NAT and then derive the first law of black hole thermodynamics. The temperature of the black holes depends on both the mass and the NAT \textquotedblleft charge"of the black holes. The extreme cases where the temperature vanishes resemble the extreme Reissner-Nordstr\"{o}m black holes. We also discuss the contribution of the NAT charge to the geodesics of massive and massless particles around the NAT black holes.


Introduction
Black holes are of fundamental importance today. This is because of the fact that studies of their properties from both theoretical and observational points of view are being expected to shed much light on the nature of the gravity at strong gravity regimes and at very high energy scales where the gravitational force becomes dominant over the other interactions. For this reason, they have always been at the heart of the theoretical investigations involving gravitational phenomena, especially since the discovery of the four laws of black hole mechanics [1] and Hawking radiation [2] in the context of general theory of relativity (GR). More importantly, the thermodynamic interpretation of the four laws [3] and the attributions of temperature and entropy to black hole horizon have provided useful information about the nature of quantum gravity through holography [4] and its specific realization AdS/CFT correspondence [5]. Observationally, recent GW events [6] and the image taken by Event Horizon Telescope Collaboration [7] have proven the existence of black holes by direct observations, which has also justified the theoretical studies conducted so far.
The event horizon of a black hole is a globally-defined causal boundary which separates the inside of the black hole from the outside. More formally, it is a null surface separating those light rays reaching infinity from those falling to the singularity inside. Since it is defined globally, the determination of the location of the event horizon requires in general the knowledge of the global structure of the spacetime. However, in the case of static spherically symmetric spacetimes, one can introduce convenient coordinate systems in which the determination is made by looking for places where the local light cones tilt over. This implies that the existence of event horizons (and of black holes) has to do with the local Lorentz invariance of the spacetime. Therefore, it is of great importance to explore the properties of black holes in gravity theories that exhibit violations of local Lorentz invariance.
Lorentz symmetry is built in GR which describes gravitation well at low energy scales by assuming the spacetime structure as continuous and smooth, excluding singularities. But this symmetry might be broken at very high energy scales, especially at the Planck or quantum gravity scales, where quantum gravitational effects must be taken into account. In fact, there are theories, such as string theory and loop quantum gravity, contemplating that the quantum fluctuations at or beyond the Planck scale might be so violent that the spacetime ceases to be continuous and has a discrete structure, and thereby the Lorentz symmetry is not valid [8]. This way of reasoning immediately leads to the contemplation of gravity theories in which Lorentz symmetry is broken explicitly.
One way to construct a Lorentz-violating gravity theory is to assume the existence of a vector field of constant norm which dynamically couples to the metric tensor at each point of spacetime. In other words, the spacetime curvature is determined together by the metric tensor and the coupled vector field in spacetime. Such a vector field is referred to as the "aether" because that generally defines a preferred direction in spacetime and breaks the local Lorentz invariance. Eintein-aether theory [9] is such a theory in which the vector field is timelike everywhere and explicitly breaks the boost sector of the Lorentz symmetry. The internal structure and dynamics of this theory have been studied extensively in the literature .
Recently, a new vector-tensor theory of gravity called the Null Aether Theory (NAT) [35] has been introduced into the realm of modified gravities. This theory assumes the dynamical vector field (the aether) inherent in the theory to be null at evey point in spacetime. In the paper [35], we first studied the Newtonian approximation of the theory and showed that it reproduces the Poisson equation at the perturbation order by, in some cases, rescaling the Newton's constant G N . Then we obtained exact spherically symmetric solutions in this theory by properly choosing the null vector field and we showed that there is a large class of solutions depending on the parameters of the theory. Among these, there are Vaidya-type nonstationary solutions because of the null aether behaving as a null matter source, and for some special values of the parameters, stationary Schwarzschild-(A)dS and Reissner-Nordström-(A)dS type solutions with some effective cosmological constant and some "charge" sourced by the aether, respectively. We also discussed the existence of stationary black holes among these exact solutions for arbitrary values of the parameters of the theory. (See [35] for details and explicit structures of these solutions.) To see the effect of the null aether in cosmology, we studied the flat FLRW metric and, taking the spatial component of null aether lying along the x axis, we found all possible perfect fluid solutions of NAT. We also discussed the existence of the Big-bang singularity and the accelerated expansion of the universe in NAT. In addition to these, to better understand the internal dynamics of the theory, we constructed exact wave solutions by specifically considering the Kerr-Schild-Kundt (KSK) class of metrics [36,37] with maximally symmetric backgrounds. After giving the exact AdS-plane wave solutions of NAT in D ≥ 3 dimensions, we also obtained all possible pp-wave solutions of the theory propagating in the flat background spacetime. These exact wave solutions are consistent with the linearized waves of the theory [38].
In this paper, we will continue our explorations in the implications of the exact spherically symmetric solutions and black hole spacetimes found in [35]. After giving a brief review of the Newtonian limit and static spherically symmetric solutions of NAT, we will first discuss the possible effect of the null aether field on the solar system dynamics by extracting the so-called Eddington-Robertson-Schiff parameters β and γ for our solutions, which explicitly appear in the perihelion precession and the light deflection expressions. We will see that, at the post-Newtonian order, there is no contribution from the aether field to the deflection of light rays passing near a massive body; that is the same as in GR! However, there is an explicit contribution, at the post-Newtonian order, from the aether field to the perihelion precession of planetary orbits. This fact can be used to constrain the parameters of the theory from solar system observations. Then we shall present the details of the black hole spacetimes by discussing the singularity structure, the ADM mass of the asymptotically flat solutions, and the thermodynamics in order. In the thermodynamics of NAT black holes especially, it is interesting to note that an appropriate definition of the NAT "charge" reduces the horizon thermodynamics to that of the Reissner-Nordström-(A)dS black hole in GR and the first law takes the standard form if the theory's parameters c 2 and c 3 satisfy a strict condition. Lastly, we will also discuss the circular geodesics of massive and massless particles around the NAT black holes to see the effect of the null aether on the particle trajectories in the spacetime. We will show that the null aether substantially changes the behavior of the circular orbits of massive and massless particles. We will also calculate the perihelion precession of planets and the deflection of light rays explicitly in the case of a nonzero cosmological constant.
The organization of the paper is as follows. In Sec. 2, we give the Null Aether Theory in detail. In Sec. 3 we review the Newtonian approximation of the theory and observe that the results we obtained in this section are consistent with the exact solutions in the next section. In Sections 4 and 5, we discuss exact spherically symmetric solutions and black hole spacetimes in NAT, respectively. In Sec. 6, we obtain the ADM mass of the asymptotically flat NAT black holes. In Sec. 7, we study the first law of black hole thermodynamics. In Sec. 8, we obtain the circular orbits of massive and massless particles around the NAT black holes.

Null Aether Theory
Aether theory is a generally covariant theory of gravity in which the metric tensor (g µν ) of the spacetime dynamically couples, through covariant derivatives, to a vector field (v µ )referred to as the "aether." In the absence of matter fields, the action of the theory can be written as [35] where R is the Ricci scalar, Λ is the bare cosmological constant, and with the dimensionless constant parameters c 1 , c 2 , c 3 , c 4 . From now on, throughout the text, we shall use the shorthand notation c ij = c i + c j for combinations of these constants. When ε = −1, the aether field is timelike and this case corresponds to the Einstein-Aether theory of [9]. In our case, however, ε = 0 and the aether becomes a null vector field. The Lagrange multiplier λ in (1) is introduced into the theory to explicitly enforce the nullity of the vector field; that is, to have v µ v µ = 0 at each point of the spacetime. Therefore the independent variables in the theory are g µν , v µ , and λ. The field equations are then obtained by varying the action (1) with respect to these fields: Varying with respect to λ immediately leads to the null constraint (3) and, making use of it, varying with respect to g µν and v µ respectively yields where we used the identificationsv Obviously, the Minkowski metric (η µν ) together with a constant null vector (v µ = const.) and λ = 0 constitute a solution to NAT. Since being null, the zero ather field (i.e. v µ = 0) with an arbitrary λ reduces the theory to the usual general relativity; however, this trivial case can be distinguished from the nontrivial aether case by imposing certain initial and boundary conditions on the solutions of the Einstein-Aether equations (4) and (5). (See the discussion in [35].) Since the aether field in NAT is null by construction, one can always introduce a scalar degree of freedom into the theory. The reasoning is as follows: First set up, at each point in spacetime, a null tetrad e a µ = (l µ , n µ , m µ ,m µ ), where l µ and n µ are real null vectors with l µ n µ = −1, and m µ is a complex null vector orthogonal to l µ and n µ , and then assume the null aether v µ is proportional to the one null leg of this tetrad, say l µ ; i.e. v µ = φ(x)l µ . Thus this geometric construction enables us to naturally introduce a scalar function φ(x)-the spin-0 part of the aether field-which generally contains the physical meaning of the aether by carrying a nonzero "aether charge."

Newtonian Limit of Null Aether Theory
The Newtonian limit of the theory can be achieved by assuming the gravitational field is weak and static and produced by a nonrelativistic matter field. Also the cosmological constant plays no role in this context so that it can be set equal to zero. Therefore in taking the Newtonian limit, we can write the metric in x µ = (t, x, y, z) as where Φ( x) is the gravitational potential on the order of G, and take the matter energymomentum tensor as where u µ = √ 1 + 2Φ δ 0 µ is the four-velocity of the matter field, ρ m and p m are the mass density and pressure, and t µν is the stress tensor with u µ t µν = 0. Then perturbing also the aether field appropriately, we consider only the zeroth and first order (linear) terms in v µ and g µν in the Eistein-Aether equations (4) and (5). At this point, however, there appear three distinct cases in perturbing the aether field, with the associated Newtonian limits: Case 1: Let us decompose the null aether field as where a µ = (a 0 , a 1 , a 2 , a 3 ) is a constant null vector representing the background aether field and k µ = (k 0 , k 1 , k 2 , k 3 ) is the perturbation which need not necessarily be a null vector. The null constraint (3) then implies that at the perturbation order. Since the metric is symmetric under rotations, we can take, without loosing any generality, a 1 = a 2 = 0 and for simplicity we will assume that k 1 = k 2 = 0. Then one can show that The last equation (15) is in the form of the Poisson equation and implies that Newton's gravitation constant G N is an effective one defined by the scaling Similar scaling also appears in the context of Einstein-Aether theory [22,23]. The constraint c 3 + c 1 = 0 can be removed by taking the stress part t µν into account in the energy momentum tensor, then there remains only the constraint c 2 = c 1 .
Case 2: Take the null aether field as v µ = φ( x)l µ where l µ is a null vector defined by the geometry (9) as with r = x 2 + y 2 + z 2 and i = 1, 2, 3. Note that any multiplicative function of x can be absorbed into the scalar function φ( x). Now assuming the perturbation φ( where φ 0 = const. = 0 and φ 1 is at the same order as G, we obtain Again, the effective value of Newton's constant can be seen from (19) This is, however, a very restricted aether theory because there is only one independent parameter c 1 left in the theory.
Case 3: Take the zeroth order scalar aether field in Case 2 as zero; i.e., φ 0 = 0. This means that φ( x) = φ 1 ( x) and is at the same order as G. Therefore, there is no contribution to the equation (4) from the aether field at the linear order in G, and from the 00 component of (4), we get which is the Poisson equation unaffected by the null aether field at the perturbation order.
On the other hand, from the ith component of the aether equation (5) we obtain, at the linear order in G, after eliminating the Lagrange multiplier λ by using the zeroth order equation.
In the case of spherical symmetry, outside the mass distribution of mass M, the Poisson equation (21) gives and the condition (22) gives where a 1 and a 2 are arbitrary constants on the order of G and we have defined the parameter which is always positive by definition. Therefore, we can immediately see that the three of the parameters of NAT must satisfy the constraint Specifically, when q = 0 (c 1 = −9c 23 /8), we have or when q = 1 (c 1 = −c 23 ), we have

Spherically Symmetric Static Solutions in Null Aether Theory
In this section, we shall review the spherically symmetric static solutions in NAT found previously in the original work [35]. The metric written in the Eddington-Finkelstein coordinates x µ = (u, r, θ, ϕ) is where u is the null coordinate, then taking the null aether field-assumed to be present at each spacetime point in the theory-is aligned with this coordinate, we obtain the solution where a 1 , a 2 ,m, and m are just integration constants and As we will show later, the constantsm and m are the mass parameters of the solutions. At this point, it is also important to note that the exact solution (32) is the same as the linearized one (24) obtained in the previous section. This means that the null aether contribution to the metric [see Eq. (33)] comes in at the order of G 2 . Now performing the coordinate transformation one can bring the metric (30) into the Schwarzschild coordinates where This metric describes the spherically symmetric static solutions in NAT, and interestingly we have lots of them due to the free parameters q, b 1 , and b 2 in the theory. The solution for q = 0 is the usual Schwarzschild-(A)dS spacetime but there are also solutions corresponding to some other specific values of the parameter q which are of special importance; for instance, This solution with A = 0 has been obtained by Mannheim and Kazanas [39] in conformal gravity who also argue that the linear term Br can explain the flatness of the galaxy rotation curves.
Before concluding this section, one last remark must be made on the possible effects of the null aether field on the solar system observations. For this purpose, we will consider the post-Newtonian parameters in the case of a static, spherically symmetric mass distribution like the Sun. Since the cosmological constant is totally negligible in this setting, the metric produced by such a body can be expanded to post-Newtonian order as [41] where M is the mass of the body and β and γ are the so-called Eddington-Robertson-Schiff parameters. These two parameters explicitly appear in the expressions for the perihelion precession of a planetary orbit and the deflection of light rays passing near the body which are respectively given by where a is the semi-major axis and e is the eccentricity of the orbit and b is the impact parameter. In general relativity, from the Schwarzschild metric, it can immediately be seen that β = γ = 1.
In NAT, we have the solutions given by (36) and (37). So taking Λ = 0, for the case q = 0, since we recover the usual Schwarzschild solution, we can immediately have β = γ = 1 just as in GR, but when q is a positive integer, the expanded metric is where we have assumed a 2 = 0 just for simplicity. It should be noted that the terms with q > 1 do not contribute to the post-Newtonian order. In other words, only the term with q = 1 has contribution to the post-Newtonian order. Now, knowing thatm ∼ G and a 1 ∼ G and comparing (41) with (38), we can read off the post-Newtonian parameters as Therefore, we can see from (40) that the null aether does not affect the light deflection at the post-Newtonian order; it is the same as in GR. However, it is obvious from (108) that it does affect the perihelion precessions of planets as This result tells us that, if b 1 > 0, the perihelion advance is greater than that of GR, and if b 1 < 0, it is less than that of GR.

Black Hole Solutions in Null Aether Theory
The metric (36)  However, we can generally make the following points. If q is an integer, h(r) = 0 becomes a polynomial equation which may have at least one positive real root representing the event horizon of the corresponding black hole. And, if q is not an integer, the limits lim r→0 + h(r) and lim r→∞ h(r) may be used to just determine the existence of the real roots; more explicitly, since h(r) is a continuous function of r, when the signs of the limits are opposite, it is certain that there is at least one real root of h(r). For example, in Table ( Black hole solutions may have one or multiple horizons. We call r = r 0 the largest root of h(r) and hence the one corresponding to the event horizon. When there is only one event horizon, the metric function h(r) can be written as where g(r) is a continuous function for r ≥ r 0 and g(r) > 0 because h(r) must be positive for r > r 0 . This means that where g(r) > 0 for r greater than the largest root, say r 0 . Again, due to the continuity of where we assume that all the roots are distinct and the event horizon is at r 0 , the largest root of (46). When some or all of the roots are coincident, we have the extreme case. For example, for two coincident roots, where g(r) > 0 for r > r 0 . Then h ′ (r 0 ) = 0.
From now on, we shall admit this condition as the indicator of an extreme black hole.
To understand the singularity structure of our solutions given in (36) and (37), we shall calculate the two of the curvature scalars; namely, the Ricci and Kretschmann scalars. For q = 0, they are where we made the definitions A 1 ≡ a 2 1 b 1 and A 2 ≡ a 2 2 b 2 . It can be seen that the only singularity is at r = 0. From these, we can also recover the standard Schwarzschild-(A)dS expressions by setting A 1 = 0 and A 2 = 0 simultaneously.

ADM Mass of Asymptotically Flat Solutions
To obtain asymptotically flat solutions, we should immediately take Λ = 0, and the metric (36) becomes where As is obvious, in the q = 0 case, the metric is just the usual Schwarzschild spacetime which is explicitly asymptotically flat. However, in the q = 0 case, to achieve asymptotically flat boundary condiaitions, one should consider the following cases separately: Since q > 0 by definition (see Eq. (34) For stationary spacetimes with the time translation Killing vector χ µ , the ADM and Komar masses are identical. So, the ADM mass can be calculated from where dΣ µν = −u [a s b ] dA, with dA = r 2 sin θdθdϕ, is the differential surface element on a two-sphere B living in a spacelike hypersurface Σ of the spacetime. Here, u µ = − √ hδ t µ and s µ = δ r µ / √ h are the unit timelike and spacelike normals to B, respectively, and B ∞ is a two-sphere at spatial infinity. Regarding the stationary nature of our spacetime (36), the corresponding Killing vector field is χ µ = δ µ t and where h(r) is given by (37) with Λ = 0 and the prime denotes differentiation with respect to r. Then, the ADM mass in (55) reduces to For the case q = 0, the ADM mass reads as but for the case q = 0, we obtain Then one realizes that, for having an asymptotically well defined ADM mass for NAT black holes, In all these cases, the ADM mass is rescaled through the definition of G in the theory; for example, in the Newtonian limit Case 1 of the previous section, G = G N (1 − c 1 a 2 3 ) and .
Although both cases a 2 = 0 and b 2 = 0 give the same ADM mass (60) for q > 0 for an observer at infinity, they differ if one considers the aether field φ by putting different constraints on the parameter q. That is, For both of these cases, the constraints on q parameter guarantees that the aether field is also well behaved at asymptotic region.

Thermodynamics of NAT Black holes
Now we shall study the thermodynamics of NAT black holes that we reviewed in Sec. 5. Here we first consider the case a 2 = 0. Then the metric function h(r) and the scalar aether field φ(r) take the forms The location of the event horizon r 0 is given by h(r 0 ) = 0 and the area of the event horizon is A = 4πr 2 0 . Now let a 1 = GQr , where Q is the NAT "charge." With this identification, (63) and (64) become At the event horizon location r 0 , we then have It is interesting that the horizon condition (67) is independent of the parameter q and, when b 1 ≡ 1 8 [c 3 − 3c 2 + c 23 q] = −1/2, it becomes that of the Reissner-Nordstrom-(A)dS black hole in GR. In addition, the scalar aether field φ(r) resembles the electric potential at r = r 0 . Now assuming the entropy as S = kA, where k is a positive constant which takes the value 1/4 [2], and varying that, we obtain where r 0m = ∂r 0 ∂m , r 0Q = ∂r 0 ∂Q , and r 0Λ = ∂r 0 ∂Λ . This relation can be translated into the form of the first law of thermodynamics as where the temperature T , the NAT charge potential V φ , the event horizon volume V , and the pressure P are given by where b 1 takes −1/2 to get the standard expression for the fist law. By using the discussions in Sec. 5, we can now explicitly see from (71) that T > 0 for the non-extreme cases and T = 0 for all the extreme cases. As a remark, it is worth mentioning the following point. The extremal event horizon r 0 is a radius where h(r 0 ) = 0 and h ′ (r 0 ) = 0, and so, when Λ = 0, the extremal event horizon for (67) can be obtained as which can equivalently be written in terms of mass and aether charge as This relation tells us that b 1 must always be less than zero and particularly for b 1 = − 1 2 , one can obtain the relationm 2 = G 2 Q 2 similar to the one in the case of the Reissner-Nordstrom black hole in Einstein gravity, which is also obvious from (67).
The thermodynamics of the other case a 1 = 0 is similar to the case above in which a 2 = 0. In this case, the metric function h(r) and the scalar aether field φ(r) become This time, defining a 2 = GQr −(q+1)/2 0 , where Q is the NAT "charge" again, we can write (77) and (78) as At the event horizon location r 0 , however, we obtain the same equations (67) and (68) h The rest goes on like in the case of a 2 = 0; the only difference is that b 1 must be replaced by b 2 in all the equations (69)-(76).

Circular Orbits
Here, we study the circular orbits at the equatorial plane, i.e θ = π 2 , for the metric (63) with a 2 = 0. Accordingly, we have two Killing vector fields K µ = (∂ t ) µ = (1, 0, 0, 0) and R µ = (∂ ϕ ) µ = (0, 0, 0, 1) corresponding to the conserved energy E = −K µ dx µ dσ and conserved angular momentum L = R µ dx µ dσ , respectively, where σ is an affine parameter along the geodesics. Then, regarding the metric, the energy and angular momentum magnitude of the orbiting body are given by On the other hand, using the geodesics equation g µν dx µ dσ dx ν dσ = ǫ, where ǫ = 0 and −1 denote the null and timelike geodesics, respectively, we obtain Using the energy and angular momentum (83), we arrive at where E = E 2 2 and the potential V reads as Substituting the metric function h in (63), we find the potential as where the first four terms are the standard terms as in GR [42], and the last four terms are the new correction terms by the cosmological constant and aether field, respectively. In Figure 1, we have plotted the potential function V versus r for some sets of q, L and a 2 1 b 1 parameters for the massive and massless particles, respectively. For each set of parameters, one can see that in general the deviation of the potential V from GR potential for the massive particles is more than for the massless particles. For both the massive and massless cases, by increasing q, the potential tends to GR. However, by increasing L, the potential increases and deviates more from GR. For b 1 > 0, the potential decreases by increasing a 2 1 b 1 values and vice versa. The circular orbits can be obtained as the radii where the potential is flat, i.e dV dr | r=rc = 0. Here r c denotes the circular orbits. Then, the equation governing the circular orbits can be obtained as For the GR limit by turning off the cosmological constant and aether field (Λ = 0 and a 1 = 0), we arrive at −L 2 r c + 3mL 2 − ǫmr 2 c = 0, which admits the following solutions for the massless and massive particles respectively In the presence of the cosmological constant and aether field, the equation (88) for the null geodesics reduces to where one can see that cosmological constant does not contribute for the null geodesics but the aether field does as the last term. Here, one may consider the particular case q = 1. This case has two solutions as Considering 9m 2 ≫ 16a 2 1 b 1 , we have Then, the second solution is a physical orbit only for b 1 < 0. Thus, in contrast to GR which has only one null circular orbit as in (90), in the presence of aether field for b 1 < 0, there are two null circular orbits in which the radius of the outer one is smaller than GR. For b 1 > 0, there is only one null circular orbit greater than the one in GR. For the case of timelike circular orbits, solving Eq. (88) for a generic q is impossible. Thus, one may consider the particular case of q = 1 where the resulted equation will be a 6th order equation for r c as (88) reduces to Finding the general real and positive solutions to this equation is not easy task. However, for realizing the effect of aether field, one may consider Λ = 0 and r c > 2m ≫ 2a 2 1 b 1 in the equation (88) which leads to This equation has two solutions as Following Carroll [42] for large L values, we obtain where we have considered 1 − 2a 2 1 b 1 > 0. Then, one can see that the aether field changes the inner and outer circular orbits of massive particles in GR given by r c − = 3m and r c + = L 2 m , respectively. Accordingly, for b 1 < 0, the outer and inner circular orbits will be larger and smaller, respectively, relative to GR and vice versa. For these orbits coincide at r c = 6m One can see from (99) that the aether field changes the smallest possible circular orbit for the massive particles as r c = 6m in GR.

Perihelion Precession
The perihelion precession represents that non-circular orbits are not perfect closed ellipses.
To derive it, one should obtain the evolution of the radial coordinate r as a function of angular coordinate ϕ, i.e. r = r(ϕ). To do this, using (83) we write the equation (85) in the following form 1 2 For more convenience, we introduce a new variable as x = 1 r . Then, the above equation takes the following form 1 2 whereṼ Then, Eq.(101) for the timelike geodesics becomes This equation is the master equation for the perihelion precession in the context of NAT for generic q and b 1 parameters. Analytically solving this equation for generic q is not an easy task and one may consider specific cases. For the case of q = 1, this equation reduces to Then, in comparison to the Newtonian gravity possessing the equation d 2 x dϕ 2 + x =m L 2 , one can realize the GR, cosmological constant and aether field corrections, respectively. One can show that this equation admits the following solution [43] x(ϕ) =m where the first term is the solution for the Newtonian gravity with the eccentricity parameter e, and the other terms are the corrections by GR, cosmological constant and aether field. Neglecting the higher order terms of the small eccentricity parameter e and using the conditionsm 2 L 2 ≪ 1, where Then, during each orbit of the planet, there is a perihelion advance given by One can rewrite this relation by converting L to the geometric quantities of each orbit. For this end, using the relation governing ordinary ellipses as where a is the semi-major axis, one can obtain the angular momentum as Then, by substituting (110) in (108), we obtain Here, the correction term by the aether field is exactly same as the one we previously obtained in (43) by using the post-Newtonian approximation. It is seen that for Λ > 0 and b 1 < 0, we always have less perihelion precession relative to GR. However, for Λ > 0 and b 1 > 0, depending on the value of contributions by aether field and cosmological constant, we may have more or less precession. Also, there is an interesting case for b 1 > 0 in which the cosmological constant and aether fields cancel out the effect of each other, i.e for Λ(1 − e 2 ) 4 a 4 = 2a 2 1 b 1 , leading to the same precession as in GR.

Light Deflection
To obtain the deflection angle of null geodesics, we set ǫ = 0 in the potential V in (102). Then, the equation governing null geodesics takes the form of which shows that similar to the closed null geodesics in section 8.1, the cosmological constant does not contribute to the light deflection angle. However, the aether field contributes. Considering the case of q = 1, this equation reduces to which has the following solution [43] x(ϕ) = 1 r 0 sin(ϕ) +m where the first term represents a straight line in polar coordinates (x, ϕ), and r 0 denotes the distance of closest approach of the light from the gravitational center. Then, the second and third terms denote the GR and aether field contributions to the light deflection angle, respectively. The light deflection angle, say ξ, can be obtained using the condition x(π + ξ) = 0 in (114) as where we have used the approximation relations sin(π + ξ) ≃ −ξ and cos(π + ξ) ≃ −1 and dropped higher order terms inm and a 2 1 b 1 . Here, one realizes that depending on the sign of the aether field parameter b 1 , the light deflection can be more or less than the GR value given by the above first term. For b 1 < 0, the aether field decreases the light deflection angle relative to the Schwarzschild case in GR. This is similar to the effect of charge in the Reissner-Nordström solution [43,44].

Conclusion
In this work, we investigated the properties of the black hole solutions found in NAT [35] which is a vector-tensor theory of gravity with the vector field being null and defining the aether field at each point of the spacetime. We first reviewed the Newtonian limit of the theory and showed that the Poisson equation is recovered at the linear order in the gravitation constant G of the theory which, depending on the form of the null vector, is related to the Newton's constant G N by a scaling factor. We also reviewed the exact spherically symmetric static solutions in NAT and extracted the post-Newtonian parameters β, γ from these solutions when Λ = 0. In GR, these parameters are β = γ = 1 and in NAT, for q = 0, we have the same values because the solution is the usual Schwarzschild metric in this case. However, for solutions with q > 0, taking a 2 = 0 for simplicity, we found that β = 1 − a 2 1 b 1 m 2 and γ = 1, meaning that, at the post-Newtonian order, the aether does not contribute to the light deflection expression, which is determined only by γ, while it contributes to the perihelion advance expression, which is determined by both β and γ. Since the perihelion advance differs from the GR value by the term ( (43)], the effect of the null aether is such that the GR value for the perihelion advance of planets is increased (for b 1 > 0) or decreased (for b 1 < 0). That is to say, solar system observations can be used to put some constraints on the parameters of the theory.
We also studied the exact static black hole solutions in NAT. We observed that, depending on the parameter q, there is a large class of black hole solutions in the theory and showed, by calculating the curvature scalars Ricci and Kretschmann, that all the these solutions are singular only at r = 0. These black holes possess in general multiple event horizons and the locations of these horizons are dependent on the parameters (q, Λ, a 1 , a 2 , b 1 , b 2 ,m, m) and the relations between them. There are also extreme cases in which some or all of the event horizons coincide. To determine the mass parameters of these solutions, we calculated the ADM mass of the asymptotically flat black holes and showed that, just like the mass parameter m in the case q = 0, the mass parameter in the case q > 0 readsm = GM ADM , where G is the gravitational constant appearing in the theory.
In the thermodynamics discussion of the NAT black holes, we carried out a generic analysis in which the cosmological constant is nonzero. First, defining the NAT "charge" appropriately, we showed that the horizon condition h(r 0 ) = 0 and the scalar aether field φ(r 0 ) at the horizon become similar to the ones of the Reissner-Nordström-(A)dS black hole in GR, independently of the values of the parameter q. Then we obtained the first law of thermodynamics in which the contribution of the aether field appears as V φ δQ, where V φ = −2b 1 φ(r 0 ) with φ(r 0 ) = GQ r 0 and Q is the NAT charge. Therefore, for consistency, it turns out that b 1 = −1/2 to recover the standard form of the first law.
Lastly, we studied both the null and timelike geodesics in the NAT black hole geometries. We explicitly derived the general expression for the effective potential governing the motion of the particles in the gravitational field including the correction terms due to the cosmological constant and the aether field. As is shown in Fig. (1), depending on the values of (q, L, a 2 1 b 1 ), it turns out that the deviation of the potential from the GR value is more in the case of massive particles than in the case of massless particles. In addition, by increasing q, the potential tends to the GR one, while, by increasing L, it deviates more from the GR value for both the massive and massless particles. We also obtained the general equation governing the location of the circular geodesics for both massive and massless particles to which there is no contribution from the cosmological constant for the null geodesics as in GR. For specifically q = 1, we showed that, in contrast to GR possessing only one null circular orbit, there are two circular orbits in the presence of the aether field for b 1 < 0, and of them, the outer one has a smaller radius than that of the one in GR. For b 1 > 0, on the other hand, there is always only one circular orbit the radius of which is greater than the one in GR. In the case of timelike geodesics, again for b 1 < 0, there are two different circular orbits: the outer and the inner ones are, respectively, larger and smaller that the ones in GR. As a particular case, when these circular orbits coincide, the aether field makes the location smaller than in GR if b 1 < 0. We further studied the perihelion advance of massive particles in this context. We explicitly calculated the contributions of the cosmological constant and the null aether field and showed that, when Λ = 0, the aether field contribution is the exactly the same as the one obtained in the post-Newtonian order. Finally, we investigated the issue of light deflection angle. We showed that the cosmological constant does not contribute to the light deflection angle. However, the aether field contributes in which, depending on the sign of the b 1 parameter, the light deflection can be more or less than in GR. Indeed, for b 1 < 0, the aether field decreases the light deflection angle relative to the Schwarzschild solution in GR.
NAT is a new modified theory of gravity recently introduced [35]. So far, we have investigated this theory from various respects: Newtonian limit, spherically symmetric solutions, black holes, thermodynamics, circular geodesics, flat cosmological solutions, exact plane waves, etc. But there are some open problems regarding, for example, the stability of the theory, linearized waves, rotating black holes, generic cosmological solutions, inflationary cosmologies and etc. Therefore, to gain more understanding on the internal structure and dynamics of NAT, one needs to further investigate and pose analytical solutions to the theory.