Testing Symmergent gravity through the shadow image and weak ﬁeld photon deﬂection by a rotating black hole using the M87 ∗ and Sgr. A ∗ results

,


I. INTRODUCTION
In the Wilsonian sense, quantum field theories (QFTs) are characterized by a classical action and an ultraviolet (UV) cutoff Λ.Quantum loops lead to effective QFTs with loop momenta cut at Λ.The effective QFTs suffer from UV oversensitivity problems: The scalar and gauge boson masses receive O(Λ 2 ) corrections.The vacuum energy, on the other hand, gets corrected by O(Λ 4 ) and O(Λ 2 ) terms.The gauge symmetries get explicitly broken.The question is simple: Can gravity emerge in a way restoring the explicitly broken gauge symmetries?Asking differently, can gravity emerge in a way alleviating the UV oversensitivities of the effective QFT?This question has been answered affirmatively by forming a gauge symmetry-restoring emergent gravity model [1][2][3].This model, briefly called as symmergent gravity, has been built by the observation that, in parallel with the introduction of Higgs field to restore gauge symmetry for a massive vector boson (with Casimir invariant mass) [4][5][6], spacetime affine curvature can be introduced to restore gauge symmetries for gauge bosons with loop-induced (Casimir non-invariant) masses proportional to the UV cutoff Λ [1][2][3].Symmergent gravity is essentially emergent general relativity (GR) with a quadratic curvature term.It exhibits distinctive signatures, as revealed in recent works on static black hole spacetimes [7][8][9][10].In the present work, we use observational data on M87* and Srg.A* black holes to study symmergent gravity.We focus on rotating black hole spacetimes as these black holes can act as a richer test and analysis laboratory.What we are doing can be viewed as "doing particle physics via black holes" as it will reveal salient properties of the QFT sector through various black hole features.
In 1919 Arthur Eddington led an expedition to prove Einstein's theory of relativity using gravitational lensing and since then lensing has become an important tool in astrophysics [11][12][13][14][15][16][17][18][19][20].In astrophysics, distances are very important when determining the properties of astrophysical objects.But Virbhadra showed that just the observation of relativistic images without any information about the masses and distances can also accurately give value of the upper bound on the compactness of massive dark objects [21].Furthermore, Virbhadra proved that there exists a distortion parameter such that the signed sum of all images of singular gravitational lensing of a source identically vanishes (tested with Schwarzschild lensing in weak and strong gravitational fields [22]).
In 2008, Gibbons and Werner used the Gauss-Bonnet theorem (GBT) on the optical geometries in asymptotically flat spacetimes, and calculated weak deflection angle for first time in literature [23].Afterwards, this method has been applied to various phenomena [24][25][26][27][28][29][30][31][32][33][34][35][36][37].One of the aims of the present paper is to probe symmergent gravity through the black hole's weak deflection angle using the GBT and also the shadow silhouette as perceived by a static remote observer.In essence, by the very nature of the symmergent gravity, we will be studying the numbers of fermions and bosons and related effects in the vicinity of a spinning black hole, which has the potential to affect the motions of null and time-like geodesics as well as their spin parameter a. (By definition, a = cJ/GM 2 where J and M are the angular momentum and mass of the rotation black hole.)Recently, shadow coming from different black hole models has been analyzed extensively in the literature as a probe of the imprints of the astrophysical environment affecting it  Calculation of the shadow cast of a non-rotating black hole was pioneered by [87] and [88], and later on extended by [89] [90] to an axisymmetric spacetime.Recently, after the detection of gravitational waves in 2015 [91], the first image of the black hole in M87* was formed by the Event Horizon Telescope (EHT) using its electromagnetic spectrum [92].At the time of writing of this paper, new milestone has been achieved since the EHT revealed the shadow image of the black hole in our galaxy, Sgr.A* [93].This new result indicates that black hole theory and phenomenology is a hot research topic in vivew of the rapidly developing observational techniques.
The paper is organized as follows.In Sec.II, we give a detailed discussion of the symmergent gravity.We give in Sec.III the rotating metric, indicating dependence on the symmergent parameters.In Sec.IV, we study the effects of symmergent gravity on null geodesics, which in turn influence the shadow radius and observables associated with it.We also examine the said effect on time-like orbits in the same section.We devote Sec.V to investigation of the weak deflection angle at finite-distance, and in Sec.V A we calculate the center-of-mass energy (CM) of two particles and study particle acceleration near rotating symmergent black hole background.In Sec.VI, we conclude and give future prospects.

II. SYMMERGENT GRAVITY
In this section, we give a brief description of the symmergent gravity in terms of its fundamental parameters.The starting point is quantum field theories (QFTs).Quantum fields are endowed with mass and spin as the Casimir invariant s of the Poincaré group.Fundamentally, QFTs are intrinsic to the flat spacetime simply because they rest on a Poincaré-invariant (translation-invariant) vacuum state [94,95].Flat spacetime means the total absence of gravity.Incorporation of gravity necessitates the QFTs to be taken to curved spacetime, but this is hampered by Poincaré breaking in curved spacetime [95,96].This hamper and absence of a quantum theory of gravity [97] together lead one to emergent gravity framework [98][99][100] as a viable approach.
In general, loss of Poincaré invariance could be interpreted as the emergence of gravity into the QFT [101].In a QFT, curvature can emerge at the Poincaré breaking sources.One natural Poincaré breakings source is the hard momentum cutoff on the QFT.Indeed, an ultraviolet (UV) cutoff Λ [102] limits momenta p µ within −Λ 2 ≤ η µν p µ p ν ≤ Λ 2 interval as the intrinsic validity edge of the QFT [102].Under the loop corrections up to the cutoff Λ, the action S[η, φ, V ] of a QFT of scalars S and gauge bosons V µ receives the correction (with (+, −, −, −) metric signature appropriate for QFTs) in which η µν is the flat metric, V O is the vacuum energy which is not power-law in Λ, m i stands for the mass of a QFT field ψ i (summing over all the fermions and bosons), and c O , c m , c S and c V are respectively the loop factors describing the quartic vacuum energy correction, quadratic vacuum energy correction, quadratic scalar mass correction, and the loopinduced gauge boson mass [103,104].As revealed by the gauge boson mass term c V Λ 2 V µ V µ , the UV cutoff Λ breaks gauge symmetries explicitly since Λ is not a particle mass, that is, Λ is not a Casimir invariant of the Poincaré group.The loop factor c V (as well as c S ) depends on the details of the QFT.(It has been calculated for the standard model gauge group in [1,2].)In Sakharov's induced gravity [98,99], the UV cutoff Λ is associated with the Planck scale, albeit with explicitly-broken gauge symmetries and Planckian-size cosmological constant and scalar masses.In recent years, Sakharov's setup has been approached from a new perspective in which priority is given to the prevention of the explicit gauge symmetry breaking [1].For this aim, one first takes the effective QFT in (1) to curved spacetime of a metric g µν such that the gauge boson mass term is mapped as agreement with the fact that the Ricci curvature R µν (g) of the metric g µν can arise only in the gauge sector via the covariant derivatives [1][2][3].One next inspires from the Higgs mechanism to promote the UV cutoff Λ to an appropriate spurion field.Indeed, in parallel with the introduction of Higgs field to restore gauge symmetry for a massive vector boson (Poincare-conserving mass) [4][5][6], one can introduce spacetime affine curvature to restore gauge symmetries for gauge bosons with loop-induced (Poincare-breaking mass) masses proportional to Λ [1,2,105,106].Then, one is led to the map in which R µν (Γ) is the Ricci curvature of the affine connection Γ λ µν , which is completely independent of the curved metric g µν and its Levi-Civita connection [105][106][107].This map is analogue of the map M 2 V → φ † φ of the vector boson mass M V (Poincare conserving) into the Higgs field φ.Under the map (2) the effective QFT in (1) takes the form in which R(g) ≡ g µν R µν (Γ) is the scalar affine curvature [1].This metric-Palatini theory contains both the metrical curvature R(g) and the affine curvature R(Γ).From the third term, Newton's gravitational constant G is read out to be where M 2 is the mass-squared matrix of all the fields in the QFT spectrum.In the one-loop expression, str[. . .] stands for super-trace namely str[M 2 ] = i (−1) 2si (2s i + 1)tr[M 2 ] si in which tr[. . .] is the usual trace (including the color degrees of freedom), s i is the spin of the QFT field ψ i (s i = 0, 1/2, . . .), and [M 2 ] si is the mass-squared matrix of the fields having that spin (like mass-squared matrices of scalars (s i = 0), fermions (s i = 1/2) and so on).One keeps in mind that tr[. . .] encodes degrees of freedom g i (like color and other degrees of freedom) of the particles.It is clear that known particles (the standard model spectrum) cannot generate Newton's constant in (4) correctly (in both sign and size).It is necessary to introduce therefore new particles.Interesting enough, these new particles do not have to couple to the known particles since the only constraint on them is the super-trace in (4) [1,2].
The action (3) remains stationary against variations in the affine connection provided that such that Γ ∇ λ is the covariant derivative of the affine connection Γ λ µν , and is the field-dependent metric.The motion equation (5) implies that D µν is covariantly-constant with respect to Γ λ µν , and this constancy leads to the exact solution in which g Γ λ µν is the Levi-Civita connection of the curved metric g µν .The Planck scale in (4) is the largest scale and therefore it is legitimate to make the expansions and so that both Γ λ µν and R µν (Γ) contain pure derivative terms at the next-to-leading O (G) order [2,3].The expansion in (8) ensures that the affine connection Γ λ µν is solved algebraically order by order in G despite the fact that its motion equation (5) involves its own curvature R µν (Γ) through D µν [105,106].The expansion (9), on the other hand, ensures that the affine curvature R µν (Γ) is equal to the metrical curvature R µν (g) up to a doubly-Planck suppressed remainder.In essence, what happened is that the affine dynamics took the affine curvature R from its UV value Λ 2 ℘ in (2) to its IR value R in (9).This way, the GR emerges holographically [108,109] via the affine dynamics such that loop-induced gauge boson masses get erased, and scalar masses get stabilized by the curvature terms.This mechanism renders effective field theories natural regarding their destabilizing UV sensitivities [1,2].It gives rise to a new framework in which (i) the gravity sector is composed of the Einstein-Hilbert term plus a curvature-squared term, and (ii) the matter sector is described by an M S-renormalized QFT [1,2].We call this framework gauge symmetry-restoring emergent gravity or simply symmergent gravity to distinguish is from other emergent or induced gravity theories in the literature.
It is worth noting that symmergent gravity is not a loop-induced curvature sector in curved spacetime [99,110].In contrast, symmergent gravity arises when the flat spacetime effective QFT is taken to curved spacetime [1,2] in a way reviving the gauge symmetries broken explicitly by the UV cutoff.All of its couplings are loop-induced parameters deriving from the particle spectrum of the QFT (numbers and masses of particles).It is with these loop features that the GR emerges.In fact, the metric-Palatini action (3) reduces to the metrical gravity theory after replacing the affine curvature R µν (Γ) with its solution in (9).Of the parameters of this emergent GR action, Newton's constant G was already defined in (4).The loop factor c S depends on the underlying QFT.(It reads c S 0.29 in the standard model.)The loop factor c O , which was associated with the quartic (Λ 4 ) corrections in the flat spacetime effective QFT in (1), turned to the coefficient of quadratic-curvature (R 2 ) term in the symmergent GR action in (10).At one loop, it takes the value in which n B (n F ) stands for the total number of bosonic (fermionic) degrees of freedom in the underlying QFT (including the color degrees of freedom).Both the n B bosons and n F fermions contain not only the known standard model particles but also the completely new particles.As was commented just above (5), it is a virtue of symmergence that these new particles do not have to couple to the known ones, non-gravitationally.
The last parameter of the symmergent GR action (10) is the vacuum energy density V O .It belongs to the non-power-law sector of the flat spacetime effective QFT in (1).At one loop, it takes value (tr[. . .] involves all degrees of freedom g i of particles, like color) after discarding a possible tree-level contribution.Being a loop-induced quantity, Newton's constant in (4) involves supertrace of (masses) 2 of the QFT fields.In this regard, the potential energy V O , involving the super-trace of (masses) 4 of the QFT fields, is expected to expressible in terms of G. To see this, it proves useful to start with mass degeneracy limit in which each and every boson and fermion possess equal masses, m b = m f = M 0 , for all b and f .Needless to say, M 0 is essentially the characteristic scale of the QFT.(Essentially, M 0 is the mean value of all the field masses.)Under this degenerate mass spectrum the potential V O can be expressed as follows: where use have been made of the G formula in (4) and c O formula in (11).Now, the problem is to take into account realistic cases in which the QFT fields are not all mass-degenerate.For a QFT with characteristic scale M 0 but with no detailed knowledge of the mass spectrum, realistic cases might be represented by parametrizing the potential energy as in which the new parameter α is introduced as a measure of the deviations of the boson and fermion masses from the QFT characteristic scale M 0 .Clearly, α = 1/2 corresponds to the degenerate case in (13).Alternatively, α = 1 represents the case in which B m 4 B = F m 4 F in (13).In general, α > 1 (α < 1) corresponds to the boson (fermion) dominance in terms of the trace (masses) 4 .
Symmergence makes gravity emerge from within the flat spacetime effective QFT.Fundamentally, as follows from the action (1), Newton's constant G in (4), the quadratic curvature coefficient c O in (11), and the vacuum energy V O in (12) transpired in the flat spacetime effective QFT from the matter loops.(The loop factor c S and similar parameters couple the matter fields.)Once the gravity emerges as in (10), however, G, c O and V O gain a completely new physical meaning as the curvature sector parameters (not involving the matter fields).In this sense, they are symmergent gravity parameters in the sense that they became curvature sector parameters via the symmergence.
A glance at the second line of (10) reveals that symmergent gravity is an R+R 2 gravity theory with non-zero cosmological constant.In fact, it can be put in the form after leaving aside the scalars S and the other matter fields, after switching to (−, +, +, +) metric signature (appropriate for the black hole analysis in the sequel), and after introducing the f (R) gravity function comprising the vacuum energy V O and the quadratic-curvature term proportional to c O .The Einstein field equations arising from the action (15) , and possess the solution [111,112] after using the vacuum energy formula in (14) in the second equality.It is clear that this R 0 value would vanish if the vacuum energy were not nonzero.Indeed, for a quadratic gravity with f (R) = bR 2 one gets the solution R 0 = 0.But for a quadratic gravity like The Einstein field equations ( 17) possess static, spherically-symmetric, constant-curvature solutions of the form [111,112] in which is the lapse function following from (18).In the sequel, we will analyze this constant-curvature configuration (R = R 0 = 0) in both the dS (V O > 0 namely α < 1) and AdS (V O < 0 namely α > 1) spacetimes in the case of rotating black holes.We will use rotating black hole properties to determine or constrain the model parameters c O and α.

III. ROTATING BLACK HOLES IN SYMMERGENT GRAVITY
Appropriate for a rotating black hole geometry, we use Boyer-Lindquist coordinates with a metric free of the coordinate singularities in the spacetime both exterior to the black hole and interior to the cosmological horizon.The metric is given by [111,[113][114][115][116][117]] with the various quantities Here, a is the black hole spin parameter, and α and c O are the parameters of the symmergent gravity.The physical energy E and angular momentum J of the black hole are related to the parameters M and a via the relations [118,119] in which M is the root of the equation ∆ r (r + ) = 0 namely Variations of the functions ∆r (left panel) and gtt (right panel) with the radial coordinate r for the black hole spin parameter a = 0.9 M .The plots show the null boundaries, with solid (dashed) lines mimicking the behavior of AdS (dS) black holes.We consider only negative values of cO in agreement with [120].
where r + is the radius of the event horizon.In this way, one can express energy and angular momentum in terms of the parameters r + , a, and G. (One here notes that 1 l 2 = (â−1) 24πG cO , where the curvature radius l is related to the negative cosmological constant as Λ = −3l −2 ).As a result, one finds as the explicit expressions for the definitions in (26).
Numerically, horizons can be determined by analyzing the lapse function, that is, by solving ∆ r = 0 using its definition in (22).Besides, the location of the ergoregions can be determined numerically by plotting the metric component g tt and looking for points satisfying g tt = 0.In fact, plotted in Fig. 1 are ∆ r and g tt as functions of r for the integration constant α of 0.90 and 1.10, black hole spin parameter a = 0.9 M , and the various values of c O .In general, symmergent black hole can mimic the dS (α > 1) or AdS-Kerr (α < 1) black holes depending on the sign of α.But, as already shown by Pogosian and Silvestri [120], consistent f(R) gravity theories require dF (R) dR > 0 in order to remain stable the (non-tachyonic scalaron).This constraint imposes the condition c O < 0 on the loop-induced quadratic curvature coefficient, and restricts viable solutions to the AdS-Kerr type.In view of the underlying QFT, c O < 0 implies that n F > n B -a mostly fermionic QFT [1][2][3].In the extreme case of no new bosons (in agreement with n F > n B ), one finds out that the new particle sector remains inherently stable thanks to the fact that fermion masses do not have any power-law sensitivity to the UV cutoff.
Usually, value of the cosmological constant is scaled to a higher value to see its overall effect.In this sense, assumed values of α and c O are also rescaled similarly.As is seen from Fig. 1, null boundaries are greatly affected by the values of the symmergent parameters α and c O .Indeed, α > 1 leads to three horizons where, in the context of the dS geometry, the third horizon corresponds to the cosmic outer boundary.It shifts farther as c O gets more negative.The two remaining null boundaries are the Cauchy and event horizons.In constrast to the dS solution, the symmergent AdS solution gives only two horizons such that symmergent effects are nearly vanishingly small on the inner horizon compared to the outer horizon.As for the ergoregions, α > 1 leads to three ergoregions, and it can be seen that the farthest one turns out to be even beyond the cosmic horizon.
The Hawking temperature is given by such that the surface gravity involves the null Killing vectors χ ν .The two Killing vectors ξ ν = ∂ t and ζ ν = ∂ ϕ in the metric are associated with the time translation and rotational invariance, respectively.Consequently, we take and determine Ω under the condition that χ ν is a null vector.This leads to the constraint from which one gets which reduces to at the event horizon ∆ (r + ) = 0.For this Ω value one then gets for the surface gravity, for the Hawking temperature, and for the Bekenstein-Hawking entropy.Needless to say, all thermodynamic quantities are necessarily non-negative.
The rotating symmergent black hole satisfies therefore the first law of thermodynamics with IV. GEODESICS AROUND ROTATING SYMMERGENT BLACK HOLES In this section we give a detailed analysis of the geodesics and orbits around the symmergent black holes.

A. Null geodesic and shadow cast
We start the analysis with null geodesics.The Hamilton-Jacobi equation gives where S is the Jacobi action in terms of the affine parameter λ (proper time) and coordinates x µ .The Hamiltonian is given by in the GR so that as follows from (40) above.Using the separability ansatz for the Jacobi function with the particle mass µ, one is led to the following first-order motion equations [121] after introducing the functions From the third equation in (44) above, constants of motion can be correlated via the relation K = Ξ 2 (aE − L) 2 , which is a consequence of a hidden symmetry in the θ-coordinate [121,122].
First, we study null geodesic and shadow cast for massless particles (µ = 0).Photon circular orbits, which are always unstable, must then satisfy the following condition The null geodesic is important is an important quantity for symmergent black holes.We need, in particular, the photon region in constant r 0 , the so-called photonsphere, which is an unstable orbit.The shadow cast ultimately depends on the photon region.With µ = 0, it proves convenient to define these two impact parameters These parameters are found to possess the following explicit expressions after a lengthy algebra.The photonsphere radius r ph can then be determined by solving η(r) = 0 for r.The analytical solutions are well-known for both Schwarzschild and Kerr black holes.In our case, we plot (49) in Fig. 2 to get the qualitative impression about the location of r ph .The left panel of Fig. 2 concerns location of the photonsphere along the equatorial plane.As we would expect, the Schwarzschild case gives r ph = 3M .The extreme Kerr case gives two possible locations, that is, the prograde orbit at r ph = M and the retrograde orbit at r ph = 4M .The plot shows the locations for a = 0.90M , which fall near the said extreme values.For the symmergent effect, the dS type (dashed lines) gives slightly lower values in the retrograde case while the AdS type gives a larger one.It is clear that the more negative the c O the closer the photon ring to the Kerr case.The symmergent effect is also present in the prograde case though the deviation is smaller than in the retrograde case.In the right panel of Fig. 2, we consider photons with zero angular momentum, that is, those that traverse the equatorial plane in a perpendicular manner (the so-called nodes).The retrograde case gives higher values for such an orbit than the prograde case.The deviation caused is barely evident, especially in the prograde case (the inset highlights the deviation).As a final remark, we note that the symmergent effect mimics dS or AdS cases, where such a behavior does not occur in the Schwarzschild case.The coupling between the spin and symmergent parameters made it possible to affecy the behavior of the photonsphere.Escaping photons in the unstable orbit gives the possibility for remote observers to backward-trace and obtain a shadow cast by using the celestial coordinates (r obs , θ obs ).Such an observer is also known as the Zero Angular Momentum Observer (ZAMO).Also, in this type of co-moving frame, without loss of generality, one makes the approximation r obs → ∞ and takes θ obs = π/2.The celestial coordinates are defined as [123] and the condition r obs → ∞ leads to the simplified relations These expressions further simplify to X = 0 and Y = ± √ η when θ obs = π/2.If, furthermore, a = 0 then shadow cast of a Schwarzschild black hole (a circle) is obtained.The plot of Y vs. X is shown in Fig. 3 for the black hole spin parameter value of a = 0.90M .Overall, we observe the D-shaped nature of the shadow cast due to the spin parameter.Considering first the effect of the inclination angle, we see that as we decrease the value of θ (going up from the equatorial plane), we see that the shape changes.Remarkably, in the upper right figure, the shape becomes more compressed in the α axis, and as one continues to go up (lower left figure), it goes back to an almost spherical shape.Finally, near the pole, a drastic change in the shadow shape occurs due to stretching of the α-axis.For the symmergent effects, we note that the dS type tends to decrease the photonsphere radius.Nevertheless, as photons travel in the intervening space under the effect of the symmergent parameter the shadow size is seen to increase.We can see the contrast in the AdS type symmergent effect in that while the photonsphere increases the shadow size tends to decrease.
In the azimuthal plane, we see how the shadow becomes "D-shaped" when the balck hole spin parameter a is near extremal.The numerical value of the shadow radius associated with this shape [59] can be calculated [40] by using which we plot in Fig. 4 (upper left), along with the shadow's angular radius θ sh (upper right), where the latter is defined by after taking the black hole mass M in units of M and D in parsecs.We emphasize that these behaviors are consistent with Fig. 3. Other observables that can be derived from the shadow are the distortion parameter δ s and the energy emission rate d 2 E dωdt , which are defined as such that the energy absorption cross-section can be approximated as Π ilm ∼ πR 2 s for an observer at r obs → ∞.These observables are plotted in Fig. 4 in the lower two panels.From these panels we are able to see how the spin parameter increases the distortion of the shadow.We can also see that even if the dS symmergent effect tends to increase the shadow, the distortion it gives is smaller than the AdS symmergent effect.In the lower right panel of Fig. 4, we comparatively show a = 0.25M and a = 0.90M for the Kerr, dS and AdS symmergent black holes.For the Kerr case, we see that low values of a leads to higher values for energy emission rate and peak frequency.Adding the symmergent effect at low values of a and α = 0.90, we see that a less negative c O provides the highest peak frequency and energy emission rate.For α = 1.1, this very same c O gives the lowest emission rate and peak frequency.For a = 0.90M , we observe in the inset plot that the roles are reversed near the peak of the curve.We see that there is a point in the plot where the behavior flips as the frequency σ increases whilst the emission rate decreases.
We now investigate observational constraints on the symmergent parameter c O using the black hole shadow data collected from M87* [124] and Sgr.A* [93].The data is tabulated in Table I.With these data, one can determine diameter of the  shadow size in units of the black hole mass with In this sense, the diameter of the shadow image of M87* and Sgr.A* are d M87* sh = (11 ± 1.5)m and d Sgr.A* sh = (9.5 ± 1.4)m, respectively.Meanwhile, the theoretical shadow diameter can be obtained via d theo sh = 2R sh with the use of equation (52).Based on Fig. 5, we can conclude that the parameters we used in describing the shadow cast fall within the upper and lower bounds at the 1σ level [38,93,124,125].Where this happens, the range in c O is more extensive in M87* than in Sgr.A*.We also observe that less negative c O introduces more sensitivity to deviation than more negative c O .

B. Symmergent gravity effects on time-like orbits: Effective potential and ISCO radii
We now turn our attention to time-like orbits.Here, we will still use the equations of motion from the Hamilton-Jacobi approach but this time set µ = 1.We are interested in determining the radius of the innermost stable circular orbit (ISCO),   II.Bounds on c0 as read off from the curves in Fig. 5.
which is very important in the study of matter accretion disks and in the study also of the effective potential for getting a qualitative description of bound and unbound orbits under the symmergent effects.
To determine the qualitative behavior of massive particle's motion around the symmergent black hole (i.e.bound, stable and unstable circular orbits), we use the effective potential formula [126] written in terms of the inverse metric g µν and angular momentum L. Shown in Fig. 6 is the effective potential for L = ±3.50M .In this figure, maxima and minima correspond to the unstable and stable circular orbits, respectively.One notes that for a = 0 (Schwarzschild) case, the maxima of the effective potential curve are way lower compared to when a = 0.90M (Kerr case).We can see that the energy for an unstable circular orbit is higher in the AdS symmergent effect, where less negative c 0 giving the highest energy.Its minima, which represents the stable circular orbit, has a greater radius for more negative c 0 .Peculiar as it may seem, particles with energies higher that than the peak energy tend to get deflected back when it reaches at some radial point near the cosmic horizon (one recalls that we scaled the Symmergent dS and AdS parameters).As for the dS symmergent effect, we can see that the least negative c O gives the lowest peak energy in the unstable circular orbit.Particles with low energy are deflected back to r → ∞ when they reach low values of r.While the AdS type can have elliptical bound orbits the same way as in the Kerr case, we see that dS type does not admit such orbit as far as the value of the angular momentum used herein is concerned.In Fig. 6 (right-panel), we plotted also the effective potential in which the particle attains a negative angular momentum.It is known that it does not happen in the Schwarzschild case.As we can see in the right-panel, particles with negative energies can be created near the rotating black hole's event horizon.Clearly, symmergent gravity affects the particle's energy, which in turn contributes to the Penrose process and black hole evaporation.
The ISCO radius occurs when the maximum and minimum of the effective potential curve merge in an inflection point for a given value of L. It is clear that we have two separate cases for ISCO determination since we need to consider both the prograde and retrograde motion of the time-like particles.Besides, such a circular orbit is only marginally stable, which means that any inward perturbation will lead the particle to spiral toward the event horizon.In locating the ISCO radius, we use equation (46) and solve for the energy E cir .One starts with [121,127], Differentiating equation ( 58) with respect to r and after some algebra we get (E cir = E isco ) where B = −8a 2 + r (r∆ r − 5∆ r ) + 8∆ r .The ISCO radius can then be found by solving the equation [128] η(r which redues to the Kerr case when ∆(r) = r 2 − 2mr + a 2 .The upper and lower signs in equation ( 60) correspond to the prograde and retrograde orbit of the massive particle at the ISCO radius, respectively.We determine ISCO radius numerically by plotting equation (60) in Fig. 7. First, we note that in the Schwarzschild case, the ISCO radius is located at r = 6M .For the extreme Kerr case, on the other hand, the prograde ISCO is at r = M (which coincides with the prograde photonsphere) and retrograde ISCO is at r = 9M .Thus, in Fig. 7, the left-panel represents the retrograde ISCO where we can see clearly how the symmergent parameters affect the radius.The innermost radius, which determines the accretion disk's innermost region, is smaller in the AdS type whilst the dS type is slightly larger in the Kerr case.We note that the first root of η isco is unphysical since those radii are smaller than the photonsphere in the retrograde case.Finally, we see the same symmergent effects in the prograde ISCO (right-panel of Fig. 7).Before we close this section, we remark that we used lower values of c 0 in time-like geodesics than when analyzing the null boundaries and geodesics.Thus, there must be some different constraints to be considered in studying time-like geodesics.If we used the same parameters for the null case, the results would be unphysical.Nonetheless, even if c 0 are more negative, we see that time-like particles are more sensitive to symmergent effects.

V. WEAK DEFLECTION ANGLE FROM ROTATING SYMMERGENT BLACK HOLES USING GAUSS-BONNET THEOREM
In this section, we investigate the deflection angle of light rays from the rotating symmergent black hole in the weak field limit using the methodology given in [29,31].Below, we briefly discuss the finite distance method.First of all, the deflection angle is defined as in which φ RS is the longitude separation angle, Ψ R is the angle between light rays and radial direction, and Ψ S is the angle between the observer and the source.Unit tangential vector e i is used to write the above angle as follows (e r , e θ , e φ ) = dr dφ , 0, 1 where is a radial quantity.For a stationary spacetime described by the metric one straightforwardly finds where b is the impact parameter.We can write conserved quantities energy and the angular momentum per unit mass in the form where the overdot denotes the derivative of the coordinates t and φ with respect to the affine parameter λ (proper time).
In the equatorial plane (θ = π 2 ), the line element reduces to with the spatial metric γ ij .It is then possible to define as the angle between the direction e i and the radial vector R j = ( 1 √ γrr , 0, 0).Having cos Ψ at hand, we use equations ( 63) and (68) to get Now, using this methodology we can calculate the deflection angle of the symmergent black hole in the weak field limit.Indeed, with the impact parameter b = L E and radial coordinate 1/u in place of r we get where we considered only the retrograde solution.For the prograde solution, we just flip the sign of the spin parameter a.It is convenient to start by computing the separation angle integral.This can be done by writing du dφ 2 = F (u) and extracting the angle where u S (u R ) is inverse distance to the source (observer), and u 0 is the inverse of the closest approach r 0 .Considering the weak field and slow rotation approximations, the impact parameter can be related to u 0 as Then, performing the necessary calculations we find where the Kerr term is given by In this equation, the Schwarzschild term is To get the remaining terms appearing in the light deflection angle expression, one should identify the Ψ terms.We find This relation produces where one has and where one has found Combining the above equations, we get an expression of the light deflection angle given by This form can be reduced to a simplified one using certain convenable approximations.Taking u S b << 1 and u R b << 1, we can get an expression involving divergent terms coupled to the cosmological contributions.These terms should be existed to show the cosmological background dependence.The desired deflection angle of light rays is found to be An examination shows that this expression recovers many previous findings.Without the cosmological contributions, we get the results of the ordinary Symmergent black hole solutions investigated in [9].
Plotted in Fig. 8 are the equations (80) (left-panel) and ( 81) (right-panel).The solid curves represent photons co-rotating with the black hole, while the dashed curves represent the counter-rotating photons.For the left-panel, which holds for finite distance, the results show that the co-rotating photons produce weak deflection angle at lower impact parameters compared to the counter-rotating photons.If we compare θ for a specific orbital direction of photons, we see that the symmergent AdS type gives a higher value for θ at large impact parameters.However, AdS type in the counter-rotation case gives a higher θ than the co-rotation case.Finally, the right-panel stands for the case where we used the large distance approximation.Since the symmergent parameters are scaled, the location of the receiver in this case is comparable to the cosmic horizon.In this case, the symmergent effect is strong, giving the effects shown in the plot.We could see that there are cases where θ is repulsive.In general, weak deflection angle is seen to have a strong sensitivity to the symmergent parameters.

A. Particle acceleration near rotating Symmergent black hole backgrounds
As Banados, Silk and West (BSW) show that collisions in equatorial plane near a Kerr black hole can occur with an arbitrarily high center of mass-energy and Kerr black holes can serve as particle accelerators [129].Then the BSW mechanism has been studied different black hole spacetimes [130][131][132][133][134][135][136].In the present work, we find that the symmergent parameter c 0 has an important effect on the result.To do that first we consider the particle motion on the equatorial plane (θ = π 2 , ρ 2 = r 2 ).In fact, generalized momenta P µ can be written as (µ, ν = t, r, φ, θ) where the dot stands for derivative with respect to the affine parameter λ.Then, one can write the generalized momenta P t (test particle's energy per unit mass E) and P φ (the angular momentum per unit mass L parallel to the symmetry axis) where one keeps in mind that P t and P φ are constants of motion.
Our goal is to study the CM energy of the two-particle collision in the background spacetime of the rotating symmergent black holes.To this end, we start with the CM energy E CM [129] in which u µ 1 , u ν 2 are the 4-velocity vectors of the two particles (u = ( ṫ, ṙ, 0, φ)).Correspondingly, particle i = 1, 2 has angular momentum per unit mass L i and energy per unit mass E i .We also consider particles with the same rest mass m 0 .Initially the two particles are at rest at infinity (E 1 /m 0 = 1 and E 2 /m 0 = 1) and then they approach the black hole and collide at a distance r.Then, their CM energy takes the form with the functions Collisions occur at the horizon of the black hole so that ∆ r = 0 in equation ( 86) so CM energy could diverge when the particles approach the horizon.It is not difficult to see that the denominator of E 2 CM is zero there.Evidently, the maximal energy of the collision occurs if L 1 and L 2 are opposite (such as head-on collision).Plotted in Fig. 9 is E CM as a function of the radial distance r.As is seen from the figure, CM energy blows up at the horizons.It demonstrates the radial dependence of the CM energy of the particles moving along circular orbits with the different symmergent parameters.The symmergent parameter decreases the CM energy.

VI. CONCLUSION
In this paper, we have studied rotating black holes in symmergent gravity, and used deviations from the Kerr black hole to constrain the parameters of the symmergent gravity.Symmergent gravity generates parameters of the gravitational sector from flat spacetime loops.It induces the gravitational constant G and quadratic-curvature coefficient c O .In the limit in which all fields are degenerate in mass, the vacuum energy V O can be expressed in terms of G and c O .We parametrize deviation from this degenerate limit by the parameter α.The black hole spacetime is dS for α < 1 and AdS for α > 1.In constraining the symmergent parameters c O and α, we utilize the EHT observations on the M87* and Sgr.A* black holes.
We have analyzed symmergent gravity parameters for different spinning black hole properties.We first investigated the effects of symmergent gravity on the photonsphere, black hole shadow, and certain observables related to them.Our findings show that location of the photonsphere depends on the symmergent effects due to its coupling to the black hole spin parameter.We have seen how the instability of the photon radii depends on and gets affected by c O and α.Also, symmergent gravity effect is evident, even for observers deviating from the equatorial plane.Interesting effects are also seen in the energy emission rate, which may directly affect the black hole's lifetime.Constraints on symmergent parameters from the shadow radius data coming from recent observations are also discussed.For a low and high spin parameter a values (the parameter α dictates whether the symmergent effect mimics the dS and AdS type), we found that the parameter c 0 fits within ±2σ at 95% confidence level for M87* than in Sgr.A*.It means that the true value is somewhere within such confidence level.The plot in Fig. 5 reveals that the detection of symmergent effects can be achieved more easily in M87* than in Sgr.A*, for the obvious reason that the farther the galaxy is the stronger the symmergent effects are.Finally, for comparison, we also briefly analyzed such an effect on the time-like geodesic.Our findings indicate that massive particles are more sensitive to deviations caused by the symmergent parameter than null particles.
In addition, we investigated the symmergent effects on photons' weak deflection angle as they traverse near the black hole.We also considered the finite distance effects.We have found that the deflection angle depends on α but the deviations are also caused by how far the receiver is from the black hole.This effect is clearly shown in Fig. 8.It implies that the weak deflection angle can detect symmergent effects more easily for black holes that are significantly remote compared to our location.
Lastly, we study the particle collision near the rotating symmergent black hole background and we analyze the possibility that rotating symmergent black holes could act as particle accelerators.We show that rotating symmergent black holes can serve as particle accelerators because the center-of-mass energies blow up at the horizons.As a result, the BSW mechanism depends on the value of the Symmergent parameter c 0 of the black hole.
Future research direction may include investigation of the shadow, and showing its dependence on the observer's state.One may also study the spherical photon orbits or the stability of time-like orbits.
FIG.1.Variations of the functions ∆r (left panel) and gtt (right panel) with the radial coordinate r for the black hole spin parameter a = 0.9 M .The plots show the null boundaries, with solid (dashed) lines mimicking the behavior of AdS (dS) black holes.We consider only negative values of cO in agreement with[120].

FIG. 2 .
FIG.2.Location of the photonsphere (left panel), and location of the photonsphere crossing at φ = 0 or along the equatorial plane (right panel).We take a = 0.90M in both plots (a is the black hole spin parameter).

1 FIG. 3 .
FIG.3.Shadow plot for different inclination angles for the black hole spin parameter a = 0.90M (aspect ratios of figures are set to obtain the actual view).

2 FIG. 4 .
FIG. 4. Shadow observables associated with the upper left panel of Fig. 3. Shown here are shadow radius (upper-left panel), angular diameter in µas for M87* (upper-right panel), distortion parameter (lower-left panel), and energy emission rate (lower-right panel).

FIG. 6 .
FIG.6.The effective potential as a function of r/M for L > 0 (left-panel) and L < 0 (right-panel).In both plots, the black hole spin parameter is set to a = 0.90M .

FIG. 7 .
FIG.7.Location of ISCO radius for counter-rotation of a massive particle (left-panel) and co-rotation of a massive particle (rightpanel).We set a = 0.90M in both panels.

1 0 1 FIG. 8 .
FIG.8.Weak deflection angle of photons as a function of b/m for various values of the spin and symmergent parameters.The solid lines correspond to a > 0 (prograde) and dashed lines correspond to a < 0 (retrograde).In making these plots we take u S = u R = 0.5b −1 (left), u S = u R = 0.01b −1 (right), and α = 0.90.The horizontal dotted lines represent critical impact parameters in each case.

TABLE I .
Black hole observational constraints.