Gravitational lensing around Kehagias-Sfetsos compact objects surrounded by plasma

We study the optical properties of the Kehagias-Sfetsos (KS) compact objects, characterized by the"Ho\v{r}ava"parameter $\omega_{_{KS}}$, in the presence of plasma, considering its homogeneous or power-law density distribution. The strong effects of both"Ho\v{r}ava"parameter $\omega_{_{KS}}$ and plasma on the shadow cast by the KS compact objects are demonstrated. Using the weak field approximation, we investigate the gravitational lensing effect. Strong dependence of the deflection angle of the light on both the"Ho\v{r}ava"and plasma parameter is explicitly shown. The magnification of image source due to the weak gravitational lensing is given for both the homogeneous and inhomogeneous plasma.


I. INTRODUCTION
Hořava proposed a field-theory approach to quantum gravity inspired by Lifshitz's ideas of solid state physics, based on an anisotropic scaling of space and time [1,2].The Lagrangian of the Hořava theory is Lorentz invariant at low energies, but the invariance is violated at high energies.Later, within a slightly modified theory, the spherically symmetric, asymptotically flat solution has been found by Kehagias and Sfetsos [3].The Kehagias-Sfetsos (KS) solution is compatible with the Minkowski vacuum and includes an extra new parameter ω KS reflecting the quantum effects.It coincides with the Schwarzschild solution in the limit of large values of ω KS M 2 , for the source with total gravitational mass M .Considering ω KS as an universal constant, one may consider the spacetime to be regulated only by the mass of the object.For the case ω KS M 2 ≥ 1/2, the solution describes a black hole with an event horizon, while for ω KS M 2 < 1/2 it describes a naked singularity.The limits/constraints on ω KS obtained by using the observational tests do not exclude the existence of the compact objects described by KS solution [4][5][6].For example, the Solar system test gives the limit of ω KS > 3.2 × 10 20 cm 2 and implies that the total mass of the object cannot exceed 2.6 × 10 4 M [6].In the present paper, we are motivated to consider possibility of testing the KS solution using optical properties of the spacetime.The properties of the KS spacetime have been studied by various authors, see, e.g.[7][8][9][10].The particle motion around KS spacetime have been studied in [11][12][13].
One of the basic features of the metric theories of gravity is the gravitational lensing or light deflection effect due to gravitational interaction.It was first discovered by Einstein within the General Relativity and now it is considered as an useful tool to study either source or lens system.The effect of gravitational lensing is reviewed in [14][15][16][17].Beside the gravitational force, the plasma surrounding the compact object may also significantly affect the photon motion.The effect of plasma on photon motion in various spacetimes and plasma configurations has been studied by number of authors .
Recently, image of supermassive black hole in the center of galaxy M87 has been disclosed [39][40][41][42][43][44].This observation is due to the Event Horizon Telescope (EHT) based on the very large interferometry (VLBI) technique promise to get deep understanding of the strong gravitational field regime around supermassive black hole (SMBH) and test the theories of gravity.The image of the SMBH or so-called shadow of the black hole has been theoretically studied by many authors [8,25,34,.
Here we study the effect of the "Hořava" parameter on the image of the shadow of the KS compact objects and in the weak-field limit its influence on the gravitational lensing, both in the presence of plasma.
The paper is organized as follows.In Sect.II we introduce the notion of the shadow of the KS compact object in vacuum.In Sect.III we investigate the influence of plasma on the shadow of the KS compact object.Then in Sect.IV we review the photon motion around compact object in the presence of plasma.We apply the general formalism to KS spacetime and study the gravitational lensing effect around KS compact object in the presence of plasma.In the next Sect.V we consider the magnification of image source due to lensing in the presence of plasma.Finally, in Sect.VI we summarize our results.Throughout the paper we use space-like signature (−, +, +, +), the geometric system of units in which G = 1 = c and we restore them when we need to compare our results with observational data.Greek indices run from 0 to 3, Latin indices from 1 to 3. The metric of the Kehagias-Sfetsos (KS) spacetime, expressed in the standard Boyer-Lindquist coordinates and geometric units can be written as [3] where the lapse function reads and ω KS is the "Hořava" parameter.

B. Equations of geodesic motion
We treat the equations of motion by following the Hamilton-Jacobi formalism.The Hamilton-Jacobi equation reads where S is the Hamilton-Jacobi action, and λ is the affine parameter that changes along the geodesic.The four momentum of a test particle is related with the action as p α = ∂S/∂x α .Because of the symmetries of the KS spacetime, we can apply separation of the variables and the action can be written as where m is the mass of the test particle (m = 0 in the case we are dealing with, i.e. for photons), E is the energy of the particle (photon), L is the axial angular momentum of the particle (photon), S r (r) is a function of r and S θ (θ) is a function of θ.
Considering the metric of KS spacetime given in Eq.( 1), we put the Hamilton-Jacobi action given in expression (4) into Eq.( 3).Due to the separation of variables, one can find easily the equations of motion of photons in the KS spacetime in the integrated and separated form (see also [8]), where, R(r) = E 2 − Qf (r)/r 2 , Q is the separation "Carter" constant [77] having in the spherically symmetric spacetime direct meaning of square of the total angular momentum, and f (r) is given by Eq.( 2).We introduce two dimensionless impact parameters ξ = L/E and η = Q/E 2 .To obtain the expressions of the impact parameters of the photon circular orbit (being the boundary of unstable circular orbits), we have to solve simultaneously equations R(r) = 0 = dR(r)/dr.We get that the dimensionless impact parameter η corresponding to the photon circular orbit is determined by the relation C. Circular geodesics We calculate the radius, energy and axial angular momentum of circular geodesics at the equatorial plane (θ = π/2 and p θ = 0) by solving the Hamiltonian-Jacobi equation (3).Let us consider that the test particle has unit mass.In this case the Hamilton-Jacobi equation (3) takes the form where, . Using the condition giving the circular orbits, p r = 0 = ṗr , we get from Eq. ( 10) the expressions for the energy (E), and the axial angular momentum (L) in the form Photon circular orbit can be defined as the orbit where energy and angular momentum diverge.We can calculate the radius of photon circular orbit by solving the equation

D. Radius of the black hole shadow
In order to analyse the apparent shape of black hole's shadow for distant observers, it is useful to introduce the celestial coordinates (see [78][79][80] for reference) defined by the relations where r 0 is the distance between the observer and the black hole and θ 0 is the inclination angle between the normal of observer's sky plane and observer-lens axis.We demonstrate definition of the celestial coordinates in Fig. 1.
The celestial coordinates can be expressed in terms of the impact parameters determining the photon equations of motion (6,7,8) by the relations We have to plot 'γ' vs 'ζ' in order to visualize the apparent shape of the image -we can see from expressions (16) and ( 17) that ζ 2 + γ 2 = η, which implies that the apparent shape of the image is a circle of radius √ η.After solving equation (13), plugging the value of radius of photon circular orbit into expression of η in (9) and calculating the square root of η, one can get the radius of apparent shape of black hole's shadow.The plot at the top left corner of Fig. 2 depicts the apparent shape of the black hole's shadow in vacuum, given for typical values of the "Hořava" parameter.

III. SHADOW OF BLACK HOLE IN PLASMA
We consider a static distribution of plasma with refractive index n, which dependce on the photon frequency ω(x i ) is given by the relation Frequency of a photon (ω(x i )) depends on the spatial coordinates (x i ) as a result of the gravitational redshift.In Eq. ( 18), N = N (x i ) is the electron number density in plasma, m is the electron mass and e is the electron charge.

A. Photon motion in plasma
The Hamilton-Jacobi equation in a static distribution of plasma background having refractive index 'n' is given by the relation [18,81,82] Again, considering the general form of metric (1), we put the Jacobi action (m = 0 in the case of photon) as given by Eq. ( 4) into Hamilton-Jacobi equation (19).Due to the possibility of separation of variables, one can easily arrive at the equations of photon motion in presence of plasma in the integrated and separated from where, G(r) = n 2 E 2 − J f (r)/r 2 , J is the separation "Carter" constant [77], keeping the meaning of the total angular momentum in spherically symmetric backgrounds, and the lapse function f (r) is given by Eq.( 1).We follow the same treatment as in the vacuum case, i.e., solving simultaneously equations G(r) = 0 = dG(r)/dr to obtain the impact parameter η of photon circular orbit given by the relation

B. Radius of photon circular geodesic
Following the same procedure as in the vacuum case, we use the Hamiltonian-Jacobi equation ( 19) for photons and Eq. ( 21) of their radial motion that gives the conditions of the circular photon motion G(r) = 0, dG(r)/dr = 0 implying the relations determining the radius of photon circular orbit

C. Radius of the black hole shadow
Using the equations of motion (21,22,23), we can get the relations giving the celestial coordinates (16,17) of the black hole shadow in terms of the photon impact ω KS = 1.2 Shadows KS black holes.The top left corner image is for vacuum/homogeneous case and rest of the figures are plotted for various combinations of the "Hořava" parameter and the plasma parameter for inhomogeneous distribution of plasma.
To visualize the apparent shape of the shadow image, we need to plot 'γ' vs 'ζ'.In this case, we can see from Eqs ( 26) and ( 27) that ζ 2 +γ 2 = η/n 2 , which implies that the apparent shape of the black hole shadow in plasma is a circle of radius √ η/n.After solving Eq. ( 25), plugging the value of radius of the photon circular orbit into expression of η in (23) and calculating √ η/n, one can get the radius of apparent shape of black hole's shadow in the plasma.Fig. 2 depicts the apparent shape of image of black hole's shadow in plasma in various cases discussed below.

D. Homogeneous distribution of plasma
In this case the electron number density is constant throughout the plasma distribution, i.e.
So, from Eq. ( 18) we get where, κ h = K e N h /ω 2 (x i ).We call κ h as plasma parameter for homogeneous plasma distribution.Therefore, the refractive index 'n of the plasma is also constant and Eq. ( 25) is reduced to Eq. ( 13).This means that the radius of photon circular orbit remains the same as in the vacuum case.We also see that the radius of apparent image of black hole is the same as in the case of vacuum, i.e. √ η/n = r 2 f −1 (r) .This case is similar as vacuum case.

E. Inhomogeneous power-law distribution of plasma
Here we consider that the number density of electrons in plasma is given by the relation where N in is the number density of electrons at r = r 0 .According to Eq. ( 18), the refractive index of the medium (n) depends on position (r) as where a new plasma parameter κ in = K e N in r 0 /ω 2 (x i ) is introduced.We call κ in as plasma parameter for inhomogeneous distribution of plasma.Using Eq. ( 31), we can calculate the radius of photon circular orbit by solving Eq. (25).Plugging the value of radius of photon circular orbit into Eq.( 24), we find the radius of the black hole shadow to be given given by √ η/n = r 2 f −1 (r) .
In Fig. 2, we plot the 'γ' vs 'ζ' relations giving the black hole shadow for various combinations of the plasma coefficient κ in , and the dimensionless "Hořava" parameter.We use 'Schw' abbreviation to refer the case of Schwarzschild limit.From top left corner plot and the second row plots, we see that the radius of the shadow is increasing with increasing plasma coefficient and it coincides with the Schwarzschild case for large values of the plasma coefficient.We also see from the figures of third and fourth row that radius of shadow is increasing with increasing dimensionless "Hořava" parameter.

IV. DEFLECTION OF LIGHT NEAR MASSIVE BODY SURROUNDED BY PLASMA IN THE WEAK FIELD LIMIT
In Ref. [29] the authors introduced a special formalism for treating the gravitational lensing in the weak-field limit, for compact objects surrounded by plasma.In the present paper we use this formalism for calculation of the deflection angle of photons in weak gravitational field around KS black holes.
We consider a static space-time with a metric Assuming the weak-filed limit, we are allowed to write the metric as here η αβ = diag(−1, 1, 1, 1) is the Minkowski metric of the flat spacetime, and h αβ is a small perturbation satisfying conditions |h αβ | 1 and h αβ → 0 for x i → ∞.There is η αβ = η αβ and h αβ = h αβ .
We can express the deflection angle as [29] αi = 1 2 we obtain the relation for i = 1, 2 .We can express the deflection angle given by Eq.( 34) in terms of the impact parameter b as where r = √ b 2 + z 2 .For detailed exploration of this section, see [29].Now, we calculate the deflection angle of photon by using the binomial approximation for the KS metric (1).Considering 4M/ω KS r 3  1, and the binomial expansion we can write Considering the binomial approximation and neglecting the higher order terms, the temporal and radial metric components can be written as and 1 .
Considering the weak-field approximation (33), the component of the metric tensor h αβ in the Cartesian coordinates have the following form: where s i correspond to the unit vectors along the coordinate axes:

A. Deflection angle in vacuum
Now, we get the expression for the deflection angle of photon in vacuum by inserting the electron plasma frequency (ω 0 ) and the electron number density (N (x i )) equal to zero in Eq. (37).We then arrive to where r = √ b 2 + z 2 .Introducing the metric approximation and performing the above integral, we get We can see that the result in expressed by Eq. ( 44) is in agreement with the deflection angle in the Schwarzschild geometry, i.e., α b = −4M/b when ω KS → ∞.
In Fig. 3, we plot dependence of the photon deflection angle on the dimensionless impact parameter for typical values of the dimensionless "Hořava" parameter.We can see that the "Hořava" parameter demonstrate notable influence on photon deflection, if the impact parameter is low enough to allow motion close enough to the compact object.We see that the deflection angle coincides with the Schwarzschild case far from the compact object.In the first two rows, we show the dependence of the deflection angle of photon with dimensionless impact parameter (b/M ) for various combinations of dimensionless "Horava" parameter (ω KS ) and plasma parameter.In the third row, the dependence of the deflection angle of photon with plasma parameter for various combinations of dimensionless "Horava parameter" (ω KS ) and dimensionless impact parameter (b/M ).

B. Deflection angle in homogeneous distribution of plasma
Now we consider simple case when plasma distribution follows Eq. ( 28).Using Eq. ( 37), we get Performing the above integral in the approximate form of the metric, we arrive at Using Eq. ( 29), we can write the above expression as In the first two rows, we show the dependence of the deflection angle of photon with dimensionless impact parameter (b/M ) for various combinations of dimensionless "Horava" parameter (ω KS ) and plasma parameter.In the third row, the dependence of the deflection angle of photon with plasma parameter for various combinations of dimensionless "Horava parameter" (ω KS ) and dimensionless impact parameter (b/M ).
here we are neglecting the redshift of photon.
In the first and second rows of Fig. 4, we demonstrate that the deflection angle is increasing as a result of increasing plasma parameter.We can see that this phenomenon is significant near the KS compact object, and the deflection is increasing with increasing "Hořava" parameter.In the third row of Fig. 4, we demonstrate that plasma parameter has significant influence near the compact object; increasing of the plasma parameter implies increasing of the deflection angle.The effect of the plasma parameter becomes weaker far from the KS compact object.So, if the plasma density near KS compact object increases, the deflection angle increases.

C. Deflection angle in inhomogeneous distribution
of plasma(N = N0/r) We assume that the number density of electrons in plasma varies with position due to Eq. (30).Using Eq. ( 37) with the metric coefficient in the approximate form, and performing the integral, we get the deflection angle Taking into account ω 0 = ω e (∞) = 0, and by using (31), we rewrite the above expression as In Fig. 5, we show how the deflection angle is influenced by the plasma distribution near the compact object.In this case, deflection angle is significantly lower than in the case of homogeneous plasma distribution.Deflection angle decreases monotonically with increasing plasma parameter.This behaviour is opposite to the homogeneous case.

V. THE MAGNIFICATION OF THE SOURCE IMAGE
In this section, we consider magnification of the brightness of the source image due to the weak gravitational lensing.We start from well known lens equation given in [83] where β and θ are the source angle and the image angle, respectively, being related to the observer-lens axis, D s is the distance between the source and the observer, and D ls is the distance between the lens and the source.
Manipulating the lens equation ( 50), we obtain where is the distance between the lens and the observer.In Fig. 6, we present a schematic plot of the source-observer-lens system.Positions of images (θ d ) formed due to the lensing can be found by solving Eq. ( 51), where d is the number of the image.Radius of the Einstein ring is defined by R 0 = D l θ 0 , where θ 0 is the solution of Eq. ( 51) when the source is located on the line of sight i.e. β = 0, and θ 0 is defined as the Einstein angle.The magnification of the image brightness is defined by the formula where f is the number of images, I tot is the total brightness of images, I is the brightness of the source and d is the index of a concrete image.
The above equation has three real roots depicting positions of the images.The solutions are given by where B = β 2 + 4θ 2 E .Now we can calculate expressions for the Einstein angle by putting β = 0 into Eqs (57), (58) and (59).We get the relations for the three images in the form The total magnification of the image brightness can be calculated by using Eq.(52).In Fig FIG. 9. Homogeneous case: In the first row, the variation of total magnification (µtot) for homogeneous distribution of plasma with dimensionless "Hořava" parameter (ω KS ) for different combinations of angle of source from observer-lens axis (β) and plasma parameter keeping θE = 0.4472 constant.In the second row, the variation of total magnification (µtot) for homogeneous distribution of plasma with plasma parameter for different combinations of angle of source from observer-lens axis (β) and dimensionless "Hořava" parameter (ω KS ) keeping θE = 0.4472 constant.
the Schwarzschild limit (at high values of ω KS ), the upper two Einstein angles coincide to one, and lower one goes to zero.In Fig 7, on the right panel, we demonstrate variation of the total magnification of the image brightness with the dimensionless "Horava parameter", for various values of the inclination angle of source from the observer-lens axis (β).We see that the magnification is increasing with the inclination angle decreasing.

B. Image magnification in homogeneous plasma
By using Eqs ( 47) and ( 51), neglecting the terms involving ω 2 KS , we obtain the lens equation in the case of homogeneous distribution of plasma in the form where, θ E and θ F are given by Eqs ( 54) and ( 55), respectively.We get three real roots corresponding to three different image positions by solving Eq. ( 63).The solu-tions are given by where, We can obtain the expressions for the Einstein angle in similar way as in the vacuum case, by putting β = 0 in the above expressions.
In the left panel of Fig. 8, we present for homogeneous distribution of plasma the dependence of the Einstein angle(θ 0 ) on dimensionless "Hořava" parameter for different values of θ E .We can see three Einstein rings for each value of θ E , when the dimensionless "Hořava" parameter is low enough.As in the vacuum case, we obtain only one Einstein ring in the Schwarzschild limit( at large values of ωKS ).The total magnification of the image brightness can be calculated by using the Eq. ( 52).In the first row of Fig 9, we present the variation of the total magnification of the image brightness with the dimensionless "Hořava" parameter, and with the plasma parameter, for FIG.10.Inhomogeneous: In the first row, the variation of total magnification (µtot) for inhomogeneous distribution (1/r) of plasma with dimensionless "Hořava" parameter (ω KS ) for different combinations of angle of source from observer-lens axis (β) and plasma parameter keeping θE = 0.4472 constant.In the second row, the variation of total magnification (µtot) for inhomogeneous distribution (1/r) of plasma with plasma parameter for different combinations of angle of source from observer-lens axis (β) and dimensionless "Hořava" parameter (ω KS ) keeping θE = 0.4472 constant.
various values of β.From these plots we observe that the total magnification is decreasing if we increase plasma parameter, and the total magnification is increasing if we decrease β.This means, the total magnification decreases if plasma concentration gets denser.From the second row of Fig 9, we can conclude that the total magnification is increasing if the dimensionless "Hořava" parameter decreases.So, the total magnification increases with strength of the KS spacetime.

C. Image magnification in inhomogeneous plasma
We get the lens equation by using Eqs (51)(49) -neglecting terms involving ω 2 KS ) we arrive at We get three images whose positions are given by where . We can obtain the Einstein angles by setting β = 0 in the above expressions.In the left panel of Fig. 8, we show for inhomogeneous distribution of plasma the dependence of the Einstein angle (θ 0 ) on the dimensionless "Hořava" parameter for different values of θ E .This case is similar to the case of homogeneous distribution of plasma.We can see three Einstein rings for each value of θ E , if the dimensionless "Hořava" parameter is low enough, and we see only one Einstein ring at the Schwarzschild limit (for large values of ωKS ).
The total magnification of the image brightness can be calculated by using Eq. ( 52).In the Fig. 10, we demonstrate the variations of total magnification of the image brightness with the dimensionless "Horava parameter" and with the plasma parameter, for various values of β.
In this case, the profile has similar character as in the homogeneous case, but the magnification increases significantly due to inhomogeneity in the plasma distribution.

VI. CONCLUSION
In this work we explore the optical properties of the KS black holes in the presence of plasma.We have studied the photon motion using the Hamilton-Jacobi equation, modified due to the presence of plasma.
The radius of the shadow of the black hole is increasing when both the "Hořava" parameter and the plasma parameter increase.It was shown that "Hořava" parameter and the plasma parameter have significant influence on photons deflected near the compact object.In case of inhomogeneous distribution of plasma, the deflection angle is lower than in the homogeneous case, and it decreases monotonically with increase of plasma parameter.
We also see that the image magnification is increasing when both the inlination angle of the source from the observer-lens axis and the dimensionless "Hořava" parameter decrease.We can also conclude that if plasma concentration gets denser, the total magnification decreases.In case of injection of an inhomogeneity in the plasma distribution, the magnification increases.Now we have the observational data of the Event Horizon Telescope (EHT) [39][40][41][42][43][44] which we could use to verify results of our theoretical models with the observational data.If this or any future observation will match with our theoretical findings of deflection angle, magnification, number of image, radius of shadow then we can comment that the compact object is a Kehagias-Sfetsos black hole, and we could put limit on the "Hořava" parameter.We could also be able to distinguish the effects of "Hořava" parameter and plasma on the optical phenomena.Next, we plan to study the effect of "Hořava" parameter on the optical phenomena as gravitational lensing in strong gravitational field with anisotropic distribution of plasma.

FIG. 4 .
FIG. 4. Homogeneous case:In the first two rows, we show the dependence of the deflection angle of photon with dimensionless impact parameter (b/M ) for various combinations of dimensionless "Horava" parameter (ω KS ) and plasma parameter.In the third row, the dependence of the deflection angle of photon with plasma parameter for various combinations of dimensionless "Horava parameter" (ω KS ) and dimensionless impact parameter (b/M ).

FIG. 5 .
FIG. 5. Inhomogeneous case:In the first two rows, we show the dependence of the deflection angle of photon with dimensionless impact parameter (b/M ) for various combinations of dimensionless "Horava" parameter (ω KS ) and plasma parameter.In the third row, the dependence of the deflection angle of photon with plasma parameter for various combinations of dimensionless "Horava parameter" (ω KS ) and dimensionless impact parameter (b/M ).
7, on the left panel, we demonstrate for the three Einstein angles their dependence on the dimensionless "Hořava" parameter.At