Probing geometric proca in metric-palatini gravity with black hole shadow and photon motion

Extended metric-Palatini gravity, quadratic in the antisymmetric part of the affine curvature, is known to lead to the general relativity plus a geometric Proca field. The geometric Proca, equivalent of the non-metricity vector in the torsion-free affine connection, qualifies to be a distinctive signature of the affine curvature. In the present work, we explore how shadow and photon motion near black holes can be used to probe the geometric Proca field. To this end, we derive static spherically symmetric field equations of this Einstein-geometric Proca theory, and show that it admits black hole solutions in asymptotically AdS background. We perform a detailed study of the optical properties and shadow of this black hole and contrast them with the observational data by considering black hole environments with and without plasma. As a useful astrophysical application, we discuss constraints on the Proca field parameters using the observed angular size of the shadow of supermassive black holes M87∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^*$$\end{document} and Sgr A∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^*$$\end{document} in both vacuum and plasma cases. Overall, we find that the geometric Proca can be probed via the black hole observations.


Introduction
In the present work, our goal is to probe a minimal form of Weyl gravity by the use of existing observational data on black holes.The work can be divided in two parts: The model and its black hole candidates (Sec. 2 and Sec. 3) and observational implications of the candidate black holes (Sec.4, Sec. 5 and Sec. 6).
The question of if the general relativity (GR) is the sole theory of the gravitation is at the heart of the ongoing research in cosmology and astrophysics.In this regard, studying testable extensions of the GR proves particularly useful.One class of extensions concerns higher-order curvature invariants as in, for example, the f (R) gravity [1].One other class involves extension of the metrical geometry of the GR to non-Riemannian geometries based on metric-incompatible connections [2][3][4].Our framework is a special case of such extensions.
The simplest non-Riemannian extension is the Palatini formulation [5], which is characterized by metric g µν and the Ricci curvature R µν (Γ ) of a general symmetric affine connection Γ λ µν (a torsion-free connection independent of the metric g µν and its Levi-Civita connection g Γ λ µν ).This formulation gives the Einstein field equations with no need to extrinsic curvature [6,7].With general curvature invariants, it leads to the GR along with geometrical scalars, vectors and tensors [8].The relevance of the Palatini gravity for emergent gravity theories [9,10] is it's another application area.It was shown that the extension of the Palatini gravity with fundamental scalars like the Higgs field leads to natural inflation [11,12].Higher-curvature terms were also studied in the Palatini formalism [4,13,22] and their certain effects in astrophysics and cosmology were analysed in [14].
One step further from the Palatini formulation is the inclusion of a term like R [µν] (Γ )R [µν] (Γ ), where R [µν] (Γ ) is the anti-symmetric part of the affine Ricci tensor R µν (Γ ) .What is important about this inclusion is that it leads dynamically to the GR plus a purely geometric massive vector field Q µ [13,15].This vector field, a geometric Proca field as we will call it henceforth, is defined as Q µ ≡ 1  4 Q ν µν where Q λµν ≡ − Γ ∇ λ g µν is the non-metricity tensor [13,[15][16][17][18]. (We call non-metricity vector as geometric Proca to distinguish it from the generic Einstein-Proca system as well as the generic Z ′ bosons in the literature.)With a symmetric affine connection (torsion-free), one is left with a special case of non-Riemannian geometries in which non-metricity Q λµν is the only source of the deviations from the GR.This special case is the Weyl gravity [19][20][21][22] (see also gauge invariance analysis in [23]).The geometric Proca is a direct signature of the Weyl gravity.More specifically, it is signature of metric-incompatible symmetric connections (torsion-free).It is not something put by hand.It is not something that comes from gauge theories either.It is a geometrical massive vector field that characterizes the Weyl nature of the geometry [20,21].It has been studied as vector dark matter in [13].Its couplings to fermions (quarks and leptons) were explored in [24] in regard to the black hole horizon in the presence of the Proca field [18].
The gravity theory we explore in the present work is the metric-Palatini gravity extended with the R [µν] (Γ )R [µν] (Γ ) invariant and a negative cosmological constant (CC).Indeed, as it was made clear in [24], in the presence of the geometric Proca Q µ , the CC is a necessity for having black hole solutions.Our framework, which we call extended metric-Palatini gravity (EMPG), possesses these four basic properties: (iv ) It contains a negative CC.
Given these properties, the EMPG action takes the schematic form (whose exact form will be discussed in Sec. 2) in which R µν ( g Γ ) and R µν (Γ ) are, respectively, the metrical and affine Ricci curvatures.The EMPG is an Einstein-geometric Proca-Anti de Sitter (AdS) system.We will constrain its parameters and dynamics by using the observational data on black holes.In fact, in the literature, Einstein-Proca system (not the geometric Proca in the present work) has been analyzed for finding Reissner-Nordström type spherically-symmetric vacuum solutions [31][32][33][34], for determining the role of the Proca field [35][36][37], for obtaining spherically-symmetric static solutions [38][39][40], and for revealing the structure of the horizon radius [41][42][43].In view of this rich literature, construction and analysis of the Einsteingeometric Proca-AdS black holes in the present work can probe affine curvature, with various applications in other astrophysical and cosmological phenomena.In essence, in the present work, we are probing, for the first time in the literature, non-metricity tensor via the Einstein-geometric Proca-AdS black holes.One of the main features of the metric theories of gravity is that the electromagnetic wave propagation will be affected by the spacetime curvature.The light deflection near the compact object due to the strong gravity can be used to test the corresponding metric theory.The first observational test of the general relativity proposed by Einstein in 1915 was performed using the gravitational lensing effect observed during the solar eclipse in 1919 by Arthur Eddington [44].Later the gravitational lensing and electromagnetic wave propagation have been intensively studied by numerous authors in Refs.[45][46][47][48].Particularly, one may distinguish the gravitational lensing in weak field [49,50] and strong field [51][52][53][54][55][56][57][58] regimes, where the deflection of the light from the initial line is not large and has large values, respectively.
One of the consequences of the light propagation in curved spacetime is the appearance of the shadow of the black hole (BH).If the BH is located between the light source and the observer, then due to the capturing the part of the light by the central object observer detect the black spot in the bright background on the celestial plane.This black spot is referred to as BH shadow and first theoretically predicted by Synge [59] and further developed by Luminet [60] and Bardeen [61].Even before the first ever observation of the BH shadow in 2019 by the Event Horizon Telescope collaboration [62,63] various authors have studied the shadow of the black hole/compact objects within general relativity and modified/alternative theories of gravity (see, e.g.[64-68, 68-72, 72-84, 84-99]) and also the authors studied the effect of plasma on the BH shadow in [100][101][102][103][104][105][106].Observation data on the shadow of supermassive BHs at the center of M87 [62] and Sgr A* [107] are used to get estimations and/or constraints of BH parameters within different gravity models [108].
In fact, despite large uncertainties of black hole shadow size including the mass and distance measurements, astrophysical observation related to black holes shadows in strong and weak gravity regimes is helpful in testing gravity theories.The first black hole shadow has been observed by Event Horizon Telescope observing at a wavelength of 1.3 mm, in 2019, the image of M87* (which has 6.5±0.7 billion solar masses, located at 16.8 kpc) with the angular size 42±3 µas with the observational resolution 20 µas [62].Three years later, in 2022, the image of the central part of the Milky Way galaxy, known as Sgr A* (with 4 +1.1 −0.6 million solar masses, located at about 8 kpc), has also been observed with shadow size 48.7±7µas together with the radiation ring with the size 51.8±2.3 µas [107].The appearance of the radiation ring may be due to the presence of a plasma medium around Srg A*.Therefore, studies of the gravitational lensing effects of black holes including plasma is important.
The present work is organized as follows: In Sec. 2, we give a detailed analysis of the EMPG, starting with a ghost-free Lagrangian in the AdS background.We show that the EMGP gravity dynamically reduces to the GR plus a geometric Proca field (the non-metricity vector from the rank-3 non-metricity tensor).In Sec. 3, we give static, spherically-symmetric solutions of the EMGP and apply our results to determine the Einsteingeometric Proca-AdS black hole solutions.In Sec. 4, we study photon motion around a compact object in the EMGP model to determine the optical properties of black holes and the black hole shadow.In Sec. 5, we study plasma effects on the black hole shadow.In Sec. 6, we determine constraints on the EMPG parameters using the observational image size of the supermassive black holes M87* and Sgr A*.Finally, we discuss our results in Sec. 7. (Throughout the work we use a system of units in which G = c = 1 and signature (−, +, +, +)) for the spacetime metric.)

Einstein-Geometric Proca Model in AdS background
The EMPG model, defined schematically in (2), takes its exact form [13,24] with an action based on the metric tensor g µν like the GR and on a symmetric affine connection Γ λ µν = Γ λ νµ unlike the GR.The affine connection is independent of the Levi-Civita connection generated by the metric g µν .This connection sets the covariant derivative ∇ µ such that ∇ α g µν = 0.It also sets the Ricci curvature R µν ( g Γ ), which contracts to give the scalar curvature R(g) ≡ g µν R µν ( g Γ ).The affine connection Γ λ µν , on the other hand, is independent of the metrical connection g Γ λ µν , defines the covariant derivative Γ ∇ µ , and sets the Riemann curvature with R µ ανβ (Γ ) = −R µ αβν (Γ ).Its contractions lead to two distinct Ricci curvatures R µν (Γ ) = R λ µλν (Γ ) and R (Γ ) = R λ λµν (Γ ).The latter is actually the antisymmetric part of the former and does necessarily vanish when the affine connection Γ λ µν is replaced with the metrical one g Γ λ µν .The total affine Ricci curvature leads to the affine scalar curvature R(g, Γ ) ≡ g µν R µν (Γ ) [8,11].
The action (3) is composed of physically distinct parts [24].The first part proportional to M 2 would be the usual Einstein-Hilbert action if M were equal to the Planck mass M P l .The second term proportional to M 2 is the standard Palatini action [5], which leads to the Einstein field equations with no need for extrinsic curvature [6,7].The third term proportional to ζ was considered in both [15] and [13].The fourth term V 0 is the vacuum energy not considered in [24] (It is proportional to the CC in (2) in Introduction).The matter Lagrangian L m ( g Γ, ψ) governs the dynamics of the matter fields ψ (matter sectors involving Γ λ µν (not g Γ λ µν ) have been analyzed in [13]).Our setup differs from the socalled metric-Palatini setup [25][26][27] by the third term proportional to ζ (and dropping of the higher powers of R µν ( g Γ ) and R µν (Γ ) in view of gravitational ghosts).
In non-Riemannian geometries with symmetric connections (Γ λ µν = Γ λ νµ ) the torsion vanishes identically.Then, the difference between such geometries and the metrical geometry of the GR [15,17,18] is sourced by the non-metricity tensor as a measure of the metric-incompatibility of the affine connection Γ λ µν .In fact, non-metricity is the main feature of the Weyl geometry [20,21,23].The use of the affine connection in (7) in the action (3) and the use also of the equation of motion for the non-metricity tensor in (8) leads one to the reduced action [13,15,24] in which 2 ) is Newton's gravitational constant, and is the squared-mass of the Y µ (the geometric Proca field).We now bring the EMPG action (9) into a more compact form in which κ = 8πG N , Λ = κV 0 is the CC in (2), and Ŷµ ≡ √ κY µ is the canonical dimensionless Proca field.
The actions ( 9) and ( 11) are readily recognized to belong to the Weyl geometry with Weyl vector Y µ [20,21,23].It should nevertheless be kept in mind that neither the schematic action (2) nor the exact action (3) nor the reduced Weyl action (9) are conformalinvariant.The reason is that the EMPG, beginning from (2), is based on the Ricci curvatures and the CC not the Weyl tensor W λ µνρ whose quadratic is conformalinvariant.(Conformal transformations in Weyl gravity have been discussed in [8,109].)Variation of the action (11) with respect to the metric g µν leads to the Einstein field equations and its variation with respect to Ŷµ gives the Proca equation describing Ŷµ as a free dimensionless massive vector field.

Static Black Hole Solutions of the Einstein-Geometric Proca Model
In an attempt to find spherically-symmetric static solutions of the field equations in space coordinates (r, θ, φ), we put forth the ansatz Having done with the metric, the geometric Proca field Y µ obeying the equation of motion ( 13) can be taken as a purely time-like field in agreement with a spherically-symmetric background.With this time-like vector, the gravitational and geometric-Proca parts of the EMPG model are described by three real functions h(r), f (r) and φ(r).
In order to have a dimensionless equation's system, we define the following dimensionless quantities: So that along with Ŷµ the entire system gets expressed in terms of the dimensionless quantities.
The subtraction of the (µν = 00) and (µν = 11) components of the Einstein equations accordingly leads to the following equation The (µν = 22) component of the Einstein field equations reads as The Proca equation of motion (13) becomes In order to have an idea of what solution to expect we first solve the Proca equation ( 19) in the flat spacetime limit (f (r) = h(r) = 1).In this particular limit we find that which is what is expected of a Yukawa potential (massive vector field) in the flat spacetime.Now, we turn to curved geometry in the weak field limit.In this limit, we take the metric functions in the form after letting Λ = − 3 l 2 in (11) in which l stands for the AdS radius.In this case, the Proca equation ( 19) takes the form with the exact solution where The Breitenlohner-Freedman window of non-tachyonic negative mass-squared values lies in the following range where is the Breitenlohner-Freedman mass bound which is required to avoid tachyonic run-away instabilities in the AdS background.The range of parameter σ accordingly is Expanding φ(r) around infinity the Proca field takes the following form It is clear that in the Maxwell limit (σ → 1) this Proca field behaves as φ(r) → q 1 + q2 r and it means that q 2 has the meaning of an electromagnetic-like charge while q 1 represents a uniform potential.
Corresponding to the Proca field in (28), the metric potentials are expected to get modified as follows in which the higher inverse powers of r have been ignored.En passant, one notices that the AdS geometry makes the Proca field to have the Coulombic form in (28) in place of the Yukawa form in (20) of the flat spacetime.The parameters n 1 , n 2 and m 1 are found in terms of q 1 , q 2 and m 2 by substituting ( 28) and ( 29) into the equations of motion where q 1 , q 2 and m 2 are free parameters.
The ADM mass for this solution can be expressed as in which γ is the coefficient of the surface term for the geometric Proca.Now, for simplicity, we take after normalizing m 2 and q 1,2 in terms of the ADM mass M ADM (M ADM = 1).For small r namely when r ≪ l, the metric functions in (29) take the form after using the relations in (30) and (32).Here, one recalls that r = κ −1/2 r is the dimensionless radius.
In order to study the horizon structure, first we need to find the allowed range of parameters.(To simplify the notation, here on we drop hats on parameters namely we use r, M Y and φ to mean r, MY and φ.We also set γ = 1.)We find it by setting the condition 0 < h(r) < 1. Fig 1 shows this region for σ = 0.8.The blue region shows where h(r) > 0 and on the blue borders h(r) = 0.In the orange region h(r) < 1 and on its borders h(r) = 1.The intersection of these two regions is where 0 < h(r) < 1, and we chose the value of q 1 and q 2 from this region.Fig 2 shows how the horizon radius varies according to the parameters q 1 and q 2 .When q 1 and q 2 both increase or decrease, the horizon radius decreases, but when one of q 1 or q 2 increases and the other one decreases, the horizon radius increases.We see in both plots when q 1 = 0, the horizon radius equals 2, which corresponds to the Schwarzschild case as we expected.It is also visible by setting q 1 = 0 in equation (33).In Fig 3, we can see how the horizon radius varies with respect to σ.In all cases, horizon radius increases as σ increases except the case of q 1 = 0, which represents the horizon radius for Schwarzschild black holes.Fig 4 shows that horizon radius decreases as q 1 or q 2 increases for selected values of σ.The lines intersect where q 1 = 0 in the right plot, but they descend and spread as q 1 grows.
Fig. 1: The allowed region for choosing the values of the parameters q 1 and q 2 .The blue region is where h(r) > 0 and the orange region is where h(r) < 1.The values of q 1 and q 2 should be chosen from the intersection of these two regions.The values of γ and σ have been set to 1 and 0.8, respectively.
The Kretschmann scalar (KS) R µνσρ R µνσρ and the Ricci-squared (RS) scalar R µν R µν possess lengthy expressions as invariant functions of r, σ, q 1 and q 2 .It is clear from (33) that, for q 1 = 0, the geometry reduces to Schwarzschild geometry namely f (r) = h(r) = 1 − 2/r with a horizon at r = 2.In this sense, to reveal the ge-Fig.2: The dependence of horizon radius on q 1 and q 2 for two values of σ = 0.5 and σ = 0.9.ometry of the EMPG black hole near the Schwarzschild limit we expand the KS and RS scalars around q 1 = 0 up to quadratic order to find and  Fig. 3: Variation of horizon radius with increase of the parameter σ for selected values of q 1 and q 2 .In the left plot q 1 is fixed and q 2 takes three different values while in the right plot, q 2 is fixed and q 1 takes three different values.If q 1 → 0, the radius of horizon is the same as the horizon of Schwarzschild BH.

Horizon Radius
Fig. 4: The dependence of horizon radius on q 1 and q 2 for three values of σ.The solid black lines, dashed red lines and dotted green lines correspond to σ = 0.5, σ = 0.7 and σ = 0.9, respectively.In the left plot q 1 is fixed and in the right plot q 2 is fixed.
These invariants exemplify the KS and RS scalars near the physically-interesting domain of the Schwarzschild black hole.The KS scalar, as revealed by (34), fully agrees with the Schwarzschild limit at the zeroth order in q 1 .The RS scalar, on the other hand, vanishes up to the quadratic order q 1 , in accordance with the Schwarzschild limit.In general, KS and RS behave differently in different geometries.In Schwarzschild solution, for instance, the spacetime is Ricci-flat and KS remains as the indicator of the black hole singularity.
In the EMPG this picture changes as a funtion of q 1 (and other parameters).
From hereon, we want to dwell on the singularity structure of the black hole by analyzing the KS in detail.The EMPG black hole develops two true singularities, one at the origin r = 0 and the other at the horizon r = r H .We call them primary and secondary singularities, respectively.(In the literature singulari-ties at different r values have already been discussed for the KS for the 5-dimensional Schwarzschild-AdS spacetimes [110], for the brane solutions of supergravity theory [111] and for spherically-symmetric solutions of general relativity with scalar fields [112].)Near the Schwarzschild limit q 1 = 0, the KS expression in (34) shows explicitly that the secondary singularity occurs at r = 2, which is the event horizon of the Schwarzschild black hole.The primary singularity remains just as in the Schwarzschild case.The secondary singularity, however, changes with changing parameters.Indeed, as we illustrate in Fig. 5, the secondary singularity varies with different sets of q 1 and q 2 .(Setting σ = 0 and γ = 1).
It is worth nothing that the secondary singularity is a signature of the Proca field.Indeed, as σ → 1 (for which M Y → 0) it is found that the KS reduces to (KS) EMP G = (KS) Sch + (12q 1 q 2 (q 1 q 2 − 4))/r 6 having a sole singularity at the origin just like Schwarzschild case (with a q 1 and q 2 dependent residues, though).In this limit one recovers the electromagnetism-like geometric vector field.It is in this sense that the secondary singularity at the horizon r = r H is a direct signature of the Proca field.In fact, the plots in Fig. 5 illustrate σ = 0, a point far from the electromagnetism-like limit of σ = 1.Basically, sensitivity of the EMPG black hole on the parameters σ, q 1 and q 2 introduce a certain form of "hair".It turns out the EMPG black hole possesses structures beyond the Schwarzschild case.
The stability of the EMPG black hole is another point to dwell on.In the literature, stability conditions for the spherically symmetric solutions of the extended gravity theories have been discussed in [113,114].The EMPG black hole admits the solution as follows from the metric potentials in ( 29) and (30) after ignoring the matter contribution by taking σ = 1 and after setting m 2 = −2M BH .As was shown in [113,114], near the extremal limit one is led to the Nariai black hole metric in which where τ = arccos(cosh t) −1 .Clearly, −∞ < t < ∞ corresponds to −π/2 < τ < π/2 in the Nariai solution.
It is necessary to study especially the perturbations δφ about the Nariai solution to determine stability of its spherically-symmetric structure: For δφ, the perturbed Einstein field equations give the motion equation in terms of the cosmological time t and the mass parameter where the parameter α is set by the curvature sector.
It is clear that perturbations remain stable for 0 < α < 1/2 for which 2α − 1 < 0. For f (R) gravity is given by [113,114] at constant curvature R = R 0 = 4Λ.In this expression f R denotes derivative of f with respect to R and f ′ denotes derivative with respect to the radial coordinate r.
At large times t → ∞ tanh t → 1 and the perturbation in (40) behave as with the initial value φ 0 and the Lyapunov exponent It is linear and it has therefore α = 0 (44) as follows from (42).This α value is at the edge of the stability interval.It requires m 2 → −∞ and this corresponds to a fully stable perturbation as follows from (40).The corresponding Lyapunov exponent implies periodic time-dependence for δφ and this dependence ensures stability.This stable solution of the perturbation shows that there arise no metric instabilities in the EMPG black hole (in the absence of matter).The perturbation does not diverge.It is in this regard that the EMPG model is devoid of any metric instabilities.
It is known that the critical impact parameter b C and the critical photon-trajectory are determined by the time-varying horizon radius which depends on the Lyapunov exponent as [113] b (45) so that for the EMPG model one expects no instability in the photon trajectory.Indeed, as expected from the imaginary value of λ g for α = 0, the photon radius shows small variations about the EMPG horizon radius r H (0). and different values of q 1 with σ = 0 and γ = 1.The KS is seen to develop two singularities, one at the origin r = 0 (independent of parameters) and the other at the horizon r = r H (dependent of the parameters).
The black holes can be destabilized also by nonminimal couplings of the Proca field [115].This has been shown explicitly for the Schwarzschild solution in the presence of the non-minimally coupled Proca fields.But the EMPG model is a minimal Proca model (with no higher-order couplings to the self and to the curvature), and no instability is expected in this sense.

Photon motion
In this part we explore photon motion around BH in the EMPG model.We will use the Hamiltonian approach to investigate the photon motion.The Hamiltonian of the photons reads as where p α is the 4-momentum of the photons.The components of the four velocity for the photons in the equatorial plane (θ = π/2, p θ = 0) are given by ṫ where we used the relationship ẋα = ∂H/∂p α .From Eqs. ( 48) and ( 49), we obtain a governing equation for the phase trajectory of light.
Using the constraint H = 0, one can rewrite the above equation as [100] where we defined The radius of a circular orbit of light, particularly the one which forms the photon sphere of radius r ph , is determined by solving the following equation [98] d(γ 2 (r)) dr The solution of this equation is depicted in Figs. 6 and  7.The plots for photon sphere radius exhibit exactly the same pattern as the horizon radius, as shown in Figs. 3  and 4. One may see that photon sphere enlarges as σ increases except the case of q 1 = 0, for which the photon sphere radius is r ph = 3 and represents Schwarzschild black holes.

Black hole shadow
In this subsection we study the shadow of the BH described by the EMPG spacetime.For the angular radius of the BH shadow we explore [98,116] sin  Fig. 6: Variation of photon sphere radius with increase of the parameter σ, for selected values of q 1 and q 2 .In the left plot q 1 is fixed and q 2 takes three different values while in the right plot, q 2 is fixed and q 1 takes three different values.If q 1 → 0, the radius of photon orbit is the same as with Schwarzschild BH's orbit.Photon Sphere Radius Fig. 7: The dependence of photon sphere radius on q 1 and q 2 for three values of σ.The solid black lines, dashed red lines and dotted green lines correspond to σ = 0.5, σ = 0.7 and σ = 0.9, respectively.In the left plot q 1 is fixed and in the right plot q 2 is fixed. with α sh is the angular radius of the BH shadow, r obs is observer distance and r ph is already introduced in the previous subsection.The Observer distance is very large but finite, and it has the value r obs = D = 8.3 kpc for the Sgr A* [107] or r obs = D = 16.8Mpc for the M87* [62].Now we combine Eqs. ( 54) and ( 55), and for an observer the Eq. ( 54) takes the following form One can find the radius of BH shadow for an observer at large distance using Eq. ( 56) as [116] R sh ≃ r obs sin α sh ≃ r ph h(r ph ) h(r obs ).( 57) Fig. 8 shows how the shadow radius changes with respect to σ for different values of q 1 and q 2 .Increasing the value of q 1 or q 2 , when one of them held fixed, decreases the shadow size.It is also worth mentioning that if we reverse the signs of q 1 and q 2 at the same time in Figs. 3, 6 and 8, we get exactly the same plots.This is because in the equation (33), in the lapse function, we have q 2  1 in the second term and the product q 1 q 2 in the third term.Then, when q 1 and q 2 are both positive or negative, we get the same result.When they have opposite signs, negative q 1 and positive q 2 give the same result as a positive q 1 and negative q 2 provided that |q 1 | and |q 2 | remain unchanged.Fig 9 shows the change of shadow size according to q 1 and q 2 .Again, in the right plot, we see the shadow size of Schwarzschild black holes where the lines intersect at q 1 = 0.

Photon motion in the presence of plasma
In this subsection we investigate the effect of a plasma environment on the photon motion around the BH.The Hamiltonian in the presence of a plasma environment is written as [59,100,103,105] where ω p = 4πe 2 me N e (x) is the plasma frequency [100].N e is the number density of electrons, e and m e are the charge and mass of electron, respectively.
The components of the four velocity for the photons and the governing equation for the phase trajectory of light are the same as equations ( 47) to (51) with the new γ 2 (r) modified as [100,103] where ω 0 = −p t .Then photon sphere is defined as a solution of the following equation Here we consider two cases: uniform plasma and nonuniform plasma.Uniform plasma is defined as ω 2 p /ω 2 0 = constant.For the non-uniform plasma, for the simplicity, we choose the following radial dependence of plasma frequency [103,117]: z 0 is a free constant parameter [103,117].Fig 10 shows the photon sphere size for uniform and non-uniform plasma.In both cases (left panels) photon sphere grows as ω 2 p /ω 2 0 or z 0 /ω 2 0 grow but the effect of ω 2 p /ω 2 0 on increasing the size of photon sphere in uniform case is stronger than the effect of z 0 /ω 2 0 in non-uniform case.The right panels show the effect of q 1 on the photon sphere size in the plasma environment.In contrast to the vacuum case, photon sphere enlarges as q 1 grows and becomes bigger than photon sphere of the Schwarzschild black holes.

Black hole shadow in plasma
Now we investigate the radius of the shadow of the BH with the EMPG model in a plasma medium.The radius of the BH shadow in plasma can be written as [100,103] In the vacuum case ω p (r) ≡ 0, we recover the radius of the shadow of the BH in the EMPG model without the plasma medium.
Fig 11 shows the shadow size for uniform and nonuniform plasma.In the left panels, we see the shadow size decreases as the plasma frequency increases, but we have seen that, for the same values of parameters, the photon sphere grows as plasma frequency grows in Figs 10.In other words, plasma brings photon sphere and shadow close to each other.Again one can see the effect of q 1 on the shadow radius in the right panels and compare our result with the Schwarzschild case, which is shown in the figure coresponding to the case of q 1 = 0.It is also worth mentioning that for all the cases we studied till now, photon sphere and shadow had the same behaviour with respect to σ, for the same parameter choices.
Fig. 12 compares photon sphere and shadow for three cases: vacuum, uniform plasma, and non-uniform plasma.Here we can clearly see that plasma increases the photon sphere radius but decreases the shadow radius.Fig. 8: Variation of shadow radius with increase of the parameter σ, for selected values of q 1 and q 2 .In the left plot q 1 is fixed and q 2 takes three different values while in the right plot, q 2 is fixed and q 1 takes three different values.If q 1 → 0, the radius of BH shadow is the same as with Schwarzschild BH.Shadow Radius Fig. 9: The dependence of shadow radius on q 1 and q 2 for three values of σ.The solid black lines, dashed red lines and dotted green lines correspond to σ = 0.5, σ = 0.7 and σ = 0.9, respectively.In the left plot q 1 is fixed and in the right plot q 2 is fixed.

Astrophysical Constraints from the EHT observations
In general, a BH shadow has two observable parameters: (i) shadow radius, and (ii) its distortion from the circle shape.Observations of the supermassive black holes (SMBHs) M87* and Sgr A* have revealed that their shadow shape are almost circular with a very small distortion parameter (with about 1 − 5 % variation).Also revealed is that shadows of both SMBHs have some common features like for example shadow's center has some brightness depression and both are nearlycircular rings.As a result, the EHT collaboration has determined sizes of these BHs by assuming ring-shaped structures.Here, in this section, we will explore potential constraints on the Proca field parameters using the Hioki& Maeda method and the EHT observations of the shadows cast by SMBHs M87* and Sgr A*.Through-out our analysis will be based on the Einstein-geometric Proca BH solution.
In Figs. 13 and 14, depicted are the constraints on possible values of the parameters q 1 and q 2 (between red-dashed and blue solid curves) using the observational angular size of SMBHs M87* and Sgr A* for σ = 0.5 and σ → 1.It is seen from the figure that the constrained range increases slightly with the increase in the angular size.The ranges of q 1 and q 2 considerably increase in the presence of a uniformly distributed plasma with the frequency ratio ω p /ω 0 = √ 0.3 (corresponding to a plasma frequency of about 127 GHz).In the limit σ → 0, constrained values of q 1 and q 2 obey power-law behavior: q 1 q 2 = const, where the constant depends on the plasma frequency and the angular size.It is seen from both figures that the value of the constant decreases as the plasma frequency and the angular size increase.However, in σ = 0.5 case, q 1 and q 2 parameters have a complicated relationship depend- Fig. 10: The dependence of photon sphere radius on σ for selected values of q 1 and q 2 in uniform (upper panels) and non-uniform (lower panels) plasma.In the left panels q 1 and q 2 are fixed while plasma frequency is changing.
In the right panels, plasma frequency and q 2 are fixed but q 1 is changing.If q 1 → 0, the radius of photon orbit is the same as with Schwarzschild BH case.
ing on the size and plasma frequency.Moreover, in this case, the constrained ranges of q 1 and q 2 parameters for M87* are bigger than that which for Sgr A*, due to its large mass.

Discussions and Conclusion
In this paper, we have performed a systematic study of the EMPG model [24] in the AdS background in the linear ghost-free limit in which quadratic and higherorder curvature terms are all dropped.The model is in the class of linear, torsion-free, metric-Palatini gravity theories [25][26][27], with the extensions that a term quadratic in the antisymmetric part of the affine curvature [4,13,22,24] exists.We performed a detailed investigation of shadow and photon motion around this black hole for the main purposes of probing the properties of the geometric Proca field.(As mentioned before, we use "geometric Proca" to distinguish the massive non-metricity vector of the metric-Palatini gravity from generic Einstein-Proca systems as well as the Z ′ gauge boson literature.) The results obtained in this report can be summarized and discussed as follows: -We have shown that the metric-Palatini gravity theory in the AdS background reduces to the GR plus Fig. 11: The dependence of shadow radius on σ for selected values of q 1 and q 2 in uniform (upper panels) and non-uniform (lower panels) plasma.In the left panels q 1 and q 2 are fixed while plasma frequency is changing.In the right panels, plasma frequency and q 2 are fixed but q 1 is changing.If q 1 → 0, the radius of BH shadow is the same as with Schwarzschild BH case.a geometrical massive vector theory, which we call the EMPG model [13,24].-The EMPG system provides a novel geometro-dynamical framework.It can set the stage for diversely different physical phenomena.It would certainly have effects in various astrophysical media such as neutron stars, BHs, magnetars, and other compact objects.
In the present work, we study BHs with asymptoti-cally AdS solutions which are determined by solving the EMPG field equations ( 12) and ( 13).-Using the horizon structure of a gravitational compact object, the range of EMPG gravity parameters are also obtained in detail.It is shown that orbits enlarge or shrink depending on the value of the EMPG parameter σ.However, the effects of parameter q 1 and q 2 are opposite to each other.Fig. 13: Constraints on possible values of the parameters q 1 and q 2 for σ = 0.5.Here, we used the observational angular size of SMBHs M87* and SgrA* as θ sh = 42 ± 3µas and θ sh = 48.7 ± 7µas, respectively.
-We study the singularity structure of the EMPG black hole by analyzing the Kretschamann scalar (34) and the Ricci-squared scalar (35).Since q 1 → 0 leads to the Schwarzschild solution, we simplify KS and RS expressions by expanding q 1 around zero up to the second order.The physical properties of the KS is studied and its behaviour is shown in Fig. 5 for various parameter values q 1 at fixed values of q 2 and σ.We show that there are basically two singularities of the KS at r = 0 and at the horizon r = r H .One can see from Fig. 5 that the place of the singularity at the horizon depends on the model parameters and changes accordingly.Near q 1 = 0, as expected, the singularities in (34) occur at r = 0 and r = r H = 2, corresponding to the event horizon for the Schwarzschild black hole.It is exactly seen from the KS (34) that the Proca field causes the KS to deviate from the Schwarzschild limit by its charge (q 1 = 0 and q 2 = 0 ).The KS in the EMPG for q 1 = q 2 = 0 with σ namely 0 ≤ σ < 1 (ensuring M 2 Y = 0) implies a Proca "hair" of the black hole.Given the singularity at the horizon the black hole is expected to have a "hair", and that hair is provided by the Proca field.
-We show that the EMPG model is devoid of any metric instabilities since for our model the stability related parameter is zero which makes the Lyapunov exponent imaginary.The imaginary Lyapunov ex- ponent accordingly prevents the perturbation of the metric to diverge.The critical impact parameter b c and the critical photon-trajectory depend on the time varying horizon radius which is dependent on the imaginary Lyapunov exponent so that for the EMPG model one expects no instability in the photon trajectory.The photon radius shows small variations about the EMPG horizon radius r H (0).Moreover, since in the EMPG model the Proca field couples minimally, one expects no instabilities arising normally from non-minimal couplings of the Proca field.
-Additionally, we have analyzed photon orbits and the influence of parameter σ.Photon sphere radius increases with σ when q 1 or q 2 is fixed.Moreover, plasma effects on photon orbits have been studied for uniform and non-uniform cases.In both cases, the photon sphere grows as the plasma frequency grows, but this effect is stronger in the uniform case.-We have discussed the observable quantity which is the shadow of a black hole studied by the EMPG spacetime metric in vacuum and plasma.In the vacuum case, increasing the value of q 1 or q 2 , when one of them held fixed, decreases the shadow size.In the presence of plasma, the shadow size decreases as the plasma frequency increases, but we have seen that, for the same values of parameters, the photon sphere grows as plasma frequency grows.In other words, plasma brings photon sphere and shadow close to each other.
-All the results (photon orbits, BH shadow) obtained above are compared to the case of usual Schwarzschild black hole spacetime, which corresponds to case q 1 → 0, parametrically.It has been shown that under the effects of the parameter q 1 the size of BH shadow gets smaller.-Finally, we have determined the constraints on possible values of the parameters q 1 and q 2 using the observed shadow sizes of M87* and Sgr A* for σ = 0.5 and σ → 1.It is shown that the constraining ranges of the parameters increase with the increase of shadow size and plasma frequency.

q 1 = 0. 1 q 1 = 0. 2 q 1 =Fig. 5 :
Fig.5: Variation of the Kretschmann scalar (KS) of the EMPG black hole with the radial distance r for q 2 = 1 and different values of q 1 with σ = 0 and γ = 1.The KS is seen to develop two singularities, one at the origin r = 0 (independent of parameters) and the other at the horizon r = r H (dependent of the parameters).
Photon Sphere Radius