Exploring tidal force effects and shadow constraints for Schwarzschild-like black hole in Starobinsky–Bel-Robinson gravity

The current manuscript deals with the tidal force effects, geodesic deviation, and shadow constraints of the Schwarzschild-like black hole theorised in Starobinsky–Bel-Robinson gravity exhibiting M-theory compactification. In the current analysis, we explore the radial and angular tidal force effects on a radially in-falling particle by the central black hole, which is located in this spacetime. We also numerically solve the geodesic deviation equation and study the variation of the geodesic separation vector with the radial coordinate for two nearby geodesics using suitable initial conditions. All the obtained results are tested for Sag A* and M87* by constraining the value of the stringy gravity parameter β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} using the shadow data from the event horizon telescope observations. All the results are compared with Schwarzschild black hole spacetime. In our study, we found that both the radial and angular tidal forces experienced by a particle switch their initial behaviour and turn compressive and stretching, respectively, before reaching the event horizon. The geodesic deviation shows an oscillating trend as well for the chosen initial condition. For the constrained value of β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}, we see that the spacetime geometry generated by Sag A* and M87* is effectively same for both Schwarzschild and Starobinsky–Bel-Robinson black hole. Furthermore, we also calculated the angular diameter of the shadow in Starobinsky–Bel-Robinson black hole and compared with the Schwarzschild black hole. It is observed that the angular diameter of shadow for M87* and Sgr A* in Starobinsky–Bel-Robinson black hole is smaller than the Schwarzschild black hole. The calculated results satisfy the event horizon telescope observational constraints.


I. INTRODUCTION
Black holes are one of the most fascinating components in our universe, representing a fundamental prediction within Einstein's theory of general relativity.The scientific community and researchers have been intrigued by black holes for numerous years, driven by various reasons.Their captivating nature makes them curious and exotic entities within the universe.A widely accepted theory is that black holes are formed through the gravitational collapse of massive stars [1].Another possible alternative to the astrophysical origin of black holes are primordial black holes, formed due to over-densities seeded during inflation or by topological defects during radiation dominated era of the universe.[2][3][4][5].The existence of black holes has been strongly sup-ported by several observational characteristics.Notably, measurements of black hole spin in X-ray binaries, the detection of gravitational wave signals from binary black hole mergers by LIGO [6][7][8], and the remarkable image of a black hole at the core of the galaxy M87 captured by the Event Horizon Telescope (EHT) collaboration [9] have provided compelling evidence.Additionally, the discovery of a wide star-black hole binary system through radial velocity measurements [10] has further contributed to our comprehensive understanding of black holes.These achievements collectively reinforce the evidence for the existence and properties of black holes.These enigmatic objects have captured considerable attention across the fields of astrophysics, astronomy and high-energy physics.
Through studying black holes, researchers have made remarkable discoveries in areas such as thermodynamics, quantum effects, and gravitational interactions within curved spacetime.The physical properties of black holes remain a fascinating subject for the scientific community, even after the first discovery of a black hole as a solution to Einstein's equations more than a hundred years ago in General Relativity (GR).This initial discovery gave rise to the well-known Schwarzschild solution [11], characterized by a single parameter -the mass of the black hole.Since then, multiple solutions to Einstein's equations have been unveiled, going beyond the Schwarzschild solution.These encompass the Kerr black hole, exhibiting rotation [12], the Reissner-Nordström black hole, demonstrating electrical charge [13,14], and the Kerr-Newman black hole, possessing both charge and rotation [15].Thus in the context of GR, black holes can be described using three fundamental parameters: mass, charge, and angular momentum.However, when it comes to astrophysical black holes, it is generally anticipated that they are characterized by only two of these parameters -mass and angular momentum.Despite this simplification, these black holes may exist within highly dynamic environments, featuring elements like accretion disks and electromagnetic fields [16].Furthermore, in modified theories of gravity, black holes are linked to slight deviations from their standard counterparts, leading to unique features in their behavior [17].Investigating how these deviations can be observed in electromagnetic and gravitational wave observations is currently a highly active area of research in the scientific literature [18][19][20].In the last decade, modified gravity models have been extensively explored to address various cosmological issues and investigate the properties of black holes.Numerous theories have been proposed, taking diverse approaches and methods.Among these, Einstein Gauss-Bonnet (EGB) gravity, which finds support from string theory and related dual models, including M-theory [21][22][23], has emerged as one of the most extensively studied theory.In the context of EGB gravity, black holes have been investigated with respect to the Gauss-Bonnet gravity parameter.This parameter plays a pivotal role in constructing and studying black holes within these gravity models [24,25].The thermodynamics and optical properties of these black holes have been carefully examined concerning the Gauss-Bonnet parameter [26,27].Notably, it has been observed that variations in the Gauss-Bonnet parameter significantly impact the physical properties of black holes, with the shadow radius decreasing as the parameter increases [28].These findings provide valuable insights into the interplay between gravity models and the intriguing behavior of black holes, furthering our understanding of the cosmos at both macroscopic and microscopic scales.Recently, a novel gravitational model called Starobinsky-Bel-Robinson (SBR) gravity has been introduced, incorporating a new parameter denoted by β [29].This four-dimensional gravity is inspired by M-theory, which, at lower energies, is described by eleven-dimensional supergravity.Specifically, there are suggestions that this modified gravity is inspired from M-theory compactified on a two-sphere factor represented by S 3 × S 4 [29].Various models have been proposed within this novel gravitational framework.Inflation scenarios involving the β parameter have been studied, and associated cosmological observables have been computed [30].These investigations have explored how observational data can impose constraints on the value of β.
Furthermore, Schwarzschild-type black holes in SBR gravity have been constructed, and their thermodynamic properties have been analyzed in terms of the parameter β.Relevant quantities, including entropy and pressure, have been computed, revealing that the presence of β introduces corrections to these quantities [31].These explorations shed light on the influence of the β parameter on the thermodynamic properties of black holes in the context of SBR gravity.In a recent study, Belhaj et al. [32] investigated the deflection angle and shadows produced by black holes within the framework of Starobinsky-Bel-Robinson Gravity inspired from M-theory.
In the current paper, we focus on the analysis of tidal forces in the black holes in Starobinsky-Bel-Robinson modified Gravity from M-theory, recently obtained by Delgado and Ketov [31].It is widely known that tidal forces within a Schwarzschild spacetime cause a body falling towards the event horizon to experience stretching in the radial direction while being compressed in the angular directions [33,34].The stretching and compression are induced by gravity's tidal effect, which is created by a difference in gravity's strength between two neighbouring places.This phenomenon is highly widespread in our cosmos and have sparked public scientific interest for most of the twentieth century.In the last decades, tidal forces were analyzed for Kiselev black hole [35], Hayward black hole spacetime [36], null naked singularity spacetime [37] and some regular black holes [38].The authors of [39] compared the tidal force effects in a stringy charged black hole to those in Schwarzschild and RN black holes.The authors of [40] studied tidal forces in Reissner-Nordström spacetime and discovered that the radial and angular tidal forces switches sign at the event horizon.Also, the tidal forces in presence of a cosmological constant were studied in [41].
The tidal forces in the Kerr Black Hole were investigated in Refs.[42][43][44][45][46][47].Tidal forces play a significant role in astrophysics, leading to intriguing phenomena such as tidal disruption events (TDEs).In these events, a star may be disrupted by the tidal forces exerted by a black hole [48] ,giving rise to bright flares of X-ray [49], ultraviolet [50], and optical [51] radiation.These TDEs offer valuable insights into the interactions between black holes and surrounding celestial bodies.
In this paper, we thoroughly analyze the tidal forces and geodesic deviation for a radially in-falling particle in the Schwarzschild-like black hole spacetime in the Starobinsky-Bel-Robinson Gravity.Furthermore, we also find the constraint of the black hole parameter β with the black hole shadow by the supermassive black holes M 87 * and SgrA * data.The shadow of a black hole is created by photons emitted from the central black hole and eventually reaching an observer.These photons follow unstable orbits, forming light rings around the black hole.When projected onto an observing screen, these light rings give rise to the shadow.Unstable orbits are critical because small perturbations can cause photons to either be captured by the central compact object or escape to infinity.The light rays composing the shadow pass very close to the event horizon and experience significant bending due to the strong gravitational lensing effect.Consequently, observations of the shadow provide a means to test the strong-field properties of gravity.The seminal work of Synge [52] and Luminet [53] was instrumental in developing a formula for calculating the angular radius of photons.This formula, in turn, facilitated the characterization of the region surrounding a Schwarzschild black hole by elucidating the behavior of the light deflection angle as it diverges.Following this, Bardeen's [54] groundbreaking research in 1973 delved into the intricacies of the shadow cast by a spinning Kerr black hole, demonstrating how the black hole's rotation caused distortions in the shadow's shape.Within this context, the photon ring, encircling the black hole's shadow, emerges as a pivotal element that can be harnessed to determine the black hole's parameters and unveil gravitational features within the information field.Subsequently, there has been a surge of activity in both analytical/numerical investigations and observational studies related to the shadows cast by the different kind of black holes [28,[55][56][57][58][59][60][61][62][63].Numerous studies have been conducted in this area, including those by the works of literature [64][65][66][67][68][69][70].
The organization of this paper is as follows.
In section II, we reconsider the nature of the black holes in the SBR gravity.In Section III we discuss the radial and angular tidal forces in SBR black hole spacetime using tetrad formalism.Section IV presents the numerical solution for the geodesic deviation equation and compares the results with Schwarzschild spacetime.In section V, we investigate the shadow behaviors and provide predictions for β by the help of EHT data.We conclude our results in section VI.Throughout the manuscript, we use as h = 1 = c.

II. BLACK HOLES IN SBR GRAVITY
In this section, we discuss about the Schwarzschild like black hole solution in SBR gravity which is theorized to be embedded in M-theory existing in eleven dimensional spacetime.This solution is discussed by the authors in [71].M-theory has a bosonic sector having a metric g M N and a tensor 3 -form C M N P coupled to M2-branes which are dual to M5-branes [72].Using the compactification process coupled with the existence of stringy fluxes, the corresponding 4-D gravity models may be produced [30,31,71].The action considered is given by: .
(1) where g is the metric determinant and R is the Ricci curvature scalar.m is a free mass parameter which could have various interpretations depending on the underlying theory.β is a positive dimensionless coupling whose value is determined by the compactification of M-theory, and could be fixed by studying the optical behaviors of black hole.P 2  4 and E 2 4 are the Pontryagin and the Euler topological densities which are related to the Bel-Robinson tensor T µνλρ in four dimensions by means of the relation [30,31,71] The Bel-Robinson tensor is defined as: It can be seen from Eq.( 1) that the gravity action depends on two parameters (m, β).This gives rise to plethora of applications such as studying Hawking radiation [73], entropy [74] and inflation [75].The black hole solution in SBR gravity depends on the value of β corrected to first order perturbations.The line element of this non-rotating solution has been found to be [31] where the metric component f (r) is given by Here, r s = 2GM is Schwarzschild radius, where M represents the mass of a black hole.Taking β = 0, we recover the Schwarzschild spacetime metric.

III. TIDAL FORCE EFFECTS IN SBR BLACK HOLES
In this section, we will delve into the complexities of the equations involving tidal forces in our spacetime.To investigate the equation for the distance between two infinitesimally close and free falling particles, we use the following equation for the spacelike components of the geodesic deviation vector η µ where R µ αβγ is the Riemann curvature tensor and v µ is the unit tangent vector to the geodesic.It is known that a test body moving in space solely under the influence of gravity follows a geodesic [76].Now, due to the difference in curvature at every point on the geodesic and the fact that each point on the test body tends to follows a unique geodesic, the body experiences a difference in acceleration, leading to a stretching and squeezing effect known as tidal forces.This is the reason that the Riemann curvature tensor is used to study the tidal interactions due to gravity.
To calculate the tidal forces in the falling body's frame, we make use of the tetrad formalism [77].Tetrads are geometric objects that form a set of local coordinate bases, i.e. a locally defined set of four linearly independent vector fields known as tetrads or vierbien.At each point on a geodesic, there is a tetrad frame that forms a local inertial reference frame, where the laws of special relativity apply.Since, Lorentz transformations can connect an infinite number of orthonormal bases at a specific point, they cannot give enough information about the connection on the metric.To resolve this, we use the transformation from the orthonormal basis to the coordinate basis: where − → e µ represents the coordinate basis, − → e î represents the orthonormal basis and êî µ are the tetrad components.The metric tensor components in the tetrad basis are given as [77]: This is the key equation that enables us to use orthonormal basis in curved spacetime.The tetrad components for freely falling frames can be obtained from Eq. ( 7) and Eq.(8).
where the tetrad components follow the following rule: We notice that the component êµ 0 = v µ is the tangent vector to the geodesic i.e. the four velocity of the particle.Also, the geodesic deviation vector follows the transformation rule from global coordinates to local orthonormal coordinates: It is evident from Fig. (1) that the event horizon begins shifting for values of the order ∼ 10 −4 .For larger values of β, the location of the horizon stays the same.The non vanishing independent components of Riemann tensor for spherically symmetric spacetimes, including SBR black hole spacetime are given by [78]: For computing the Riemann curvature tensor components in tetrad basis, we use the tetrad formalism as: Following the Eq. ( 20), the tidal tensor components are given by: A. TIDAL FORCE EQUATIONS Upon obtaining the expressions from Eq. (21)(22), we can obtain the relative acceleration between two nearby particles as: By substituting the values of f (r) and its higher derivatives in Eq. (23)(24)(25) we get:

B. RADIAL TIDAL FORCE
Following the expression in Eq.( 23) and we notice that the radial tidal force vanishes at a point R = R rtf 0 given by solving the below equation: where the function f (β, r) is taken for simplicity in calculation and representation, and is given by: It can be seen from the graph in fig.
(2) that as the value of β is increased, the peak value of tidal force shifts to the right.This suggests that the maximum radial stretching happens far from the horizon as the departure from Schwarzschild geometry is significant.The peak value of radial tidal force can be calculated by equating the first derivative of the term in RHS of Eq.( 23) to zero.
The peculiar trend that can be noticed from the graph is that the nature of force changes from radial stretching to compression for non zero values of β.This is in strict contrast to the Schwarzschild black hole case where the radial tidal force profile shows a monotonically increasing trend and becomes infinity as r → 0. The transition from radial stretching to compression reverts back very close to the event horizon, after which the tidal force increases drastically till infinity.Fig (3) shows the radial force profile for higher order values of M .We notice that as the value of the mass of the black hole increases, the peak value for radial stretching keeps on increasing and is attained much earlier for a radially in-falling particle.We test the behavior of radial tidal force for both Sag A* and M87* and the results are displayed in figures 4 and 5 respectively.For both of these cases, the force profile for SBR black hole coincides perfectly with the Schwarzschild black hole case in the range −0.00038 < β < 0.00089 and −0.002 < β < 0.0003 for Sag A* and M87* respectively.This suggests that the spacetime geometry generated by Sag A* and M87* is best explained by the Schwarzschild geometry.

C. ANGULAR TIDAL FORCE
Just similar to the radial tidal force profile, the angular force profile also vanishes at a point R = R atf 0 .It can be computed by equating the RHS of Eq.( 24) to zero and is given by: For simplicity, the term f ′ (r) is expressed as 2M r 2 + f ′ (β, r).where f (β, r) is given by the Eq.( 30).Angular tidal forces show compressive behavior as the particle falls radially inwards from infinity.The force reaches a peak value and inverts it nature to stretching.The peak value shifts to the right for increasing value of β similar to the case of radial force profile.This suggests that maximum compression happens far from the horizon as the departure from the Schwarzschild geometry increases.Following the convention in Eq.( 30), the peak value of angular tidal force can be calculated by equating the first derivative of the expression in Eq.( 24) to zero.
Angular force profile also shows a very different trend from the Schwarzschild geometry.In Schwarzschild black hole case, the angular tidal force becomes increasingly compressive as the radial coordinate becomes smaller only to reach infinity as r → 0. But, in the case of SBR black hole, the angular tidal force switches its nature from compressive to stretching as the radial coordinate decreases for in-falling particle.The graph in fig.(6) shows the general force profile.The compression to stretching trend again reverts back very close to the horizon, tending to infinity thereafter.Fig (7) shows the angular force profile for higher order values of M .We notice that as the value of the mass of the black hole increases, the peak value for angular compression keeps on increasing, indicating an increase in the tidal force on the radially in-falling particle by the black hole.We test the behavior of angular tidal force for both Sag A* and M87* and the results are displayed in figures 8 and 9 respectively.We find that for both the cases, the SBR and Schwarzschild force profile coincides.This confirms that the spacetime geometry generated by both compact objects resembles Schwarzschild geometry.We take M = 1,M = 2 and M = 3 respectively.We use β = 0.60

IV. GEODESIC DEVIATION
In this section we discuss the variation of geodesic deviation vector with the radial coordinate in SBR black hole spacetime.The deviation vector measures the deformation of a body falling radially in any spacetime geometry.We can convert Eq. ( 23) and (24) in second derivatives w.r.t.r by substituting dr/dτ = − E 2 − f (r) which results from Eq. ( 9).This gives us the following second order differential equations in r: where i = {θ, ϕ}.The analytic solution for Eq. ( 34) and (35) as pointed out in [79] can be given as: where C 1 , C 2 , C 3 , C 4 are constants of integration.In order to find their value, we numerically solve the differential equations by imposing some initial conditions.For the purpose of this study we take the following initial conditions.
where β = {r, θ, ϕ}.η β (b) represents the separation between two nearby geodesics at r = b in the radial and angular directions.The initial condition in Eq. ( 38) represents a mass released from rest at r = b.In the next sections we will discuss the components of radial and the angular deviation vectors in detail.

A. RADIAL COMPONENT
In Fig. (10) we have shown the radial component of the geodesic deviation vector after solving Eq. ( 36) with initial condition.We notice that for non zero values for β, the radial component of the geodesic deviation vector reaches a finite value outside the event horizon, falls rapidly and again peaks just inside the horizon as r → 2M , unlike in the Schwarzschild black hole scenario, where it crosses the event horizon and goes till infinity because of the infinite stretching radial tidal force at the singularity.This peak value attained outside the horizon keeps on decreasing as the value of β is increased.The radial geodesic deviation trend for higher values of mass can be seen in Fig. (11).For M = 1 case the deviation profile shows an oscillating trend.Subsequently, for higher values of M , the curve becomes more steep signifying a rapid increase in the separation of two nearby geodesics as the radially in-falling particle approaches the spacetime singularity.For all the cases, the geodesic deviation approaches infinity as r → 0. By using the constrained value of β, we show the radial geodesic deviation trend for Sag A* and M87* in Fig (12) and Fig. (13).As seen from figures, the SBR black hole and Schwarzschild black hole lines perfectly coincide, leading to the conclusion that Schwarzschild geometry perfectly defines the geodesic deviation for a test particle falling radially in the spacetime generated by Sag A* and M87* for −0.00038 < β < 0.00089 and −0.002 < β < 0.0003 respectively.

B. ANGULAR COMPONENT
In Fig. (14) we have shown the angular component of geodesic deviation vector by constraining the Eq.(37) with the initial condition.The results are compared with the Schwarzschild black hole.It is observed that in the SBR black hole case, the angular deviation decreases for a radially in-falling particle, shows a little spike near the event horizon and then continues to fall beyond the event horizon.This is in contrast to the Schwarzschild black hole scenario in which the deviation decreases throughout the range of the radial coordinate.The graph in Fig. (14) shows the behavior of the geodesics.The peak value is higher for increasing value of β.The graph in Fig. (15) shows the angular geodesic deviation profile for higher values of black hole mass M .It is noticed that for increasing values of mass, the slope of the deviation vector increases, signifying a rapid decline of the separation vector of two nearby geodesics as the in-falling particle approaches the spacetime singularity.
By using the constrained value of β, we show the angular geodesic deviation trend for Sag A* and M87* in Fig (16) and Fig. (17).As seen from the graphs below, the SBR black hole and Schwarzschild black hole lines perfectly coincide, leading to the conclusion that  Schwarzschild geometry perfectly defines the geodesic deviation for a test particle falling radially in the spacetime generated by Sag A* and M87* for −0.00038 < β < 0.00089 and −0.002 < β < 0.0003 respectively.

V. BLACK HOLE SHADOW WITH M87* AND SGR A*
In this section we explore the shadow of the BH SBR gravity from M-theory.We know that for the angular radius of the BH shadow we study [64,80] with where r ph , r obs are the radius of photon sphere and observer distance, respectively.h(r) was early introduced in several references (for example, see [80]).α sh is the angular radius of the BH shadow.
Here we combine Eqs. ( 39) and ( 40), and for a distant observer the Eq. ( 39) takes the following form One can easily find the observable radius of BH shadow for observer at infinite distance using Eq. ( 41) in the following form as [64] R sh ≃ r obs sin α sh ≃ r ph f (r ph ) .
Now we consider that the supermassive BHs M87* and Sgr A* are spherically symmetric static and SBR parameter β from M-theory.Although the observation got by the EHT collaboration does not support the assumption taken here.However, here we explore theoretically the constrain on the parameter β, from the data provided by the EHT project.To constrain this parameter β in SBR gravity we use the observational data released by the EHT project for the BH shadows of the supermassive BHs M87* and Sgr A*.The angular diameter of the shadow, the distance from sun system and the mass of of the BH at the centre of the galaxy M87, are Ω M87* = 42 ± 3 µas, D = 16.8 ± 0.8 Mpc and M M87* = (6.5 ± 0.7)x10 9 M ⊙ , respectively [9].For the Sgr A* the data recently obtained by the EHT project is Ω Sgr A* = 51.8 ± 2.3 µas, D = 8277 ± 9 ± 33 pc and M Sgr A* = 4.297 ± 0.013x10 6 M ⊙ (VLTI) [81].Using this data, we can estimate the diameter of the shadow cast by the BH, per unit mass from the following expression [28,82], Now we can obtain the diameter of the shadow from the expression d theo sh = 2R sh .Thus, the diameter of the BH shadow image is d M87* sh = (11 ± 1.5)M for M87* and d Sgr A* sh = (9.5 ± 1.4)M for Sgr A*.From the data by the EHT collaboration, we obtain the constrain on the parameter β for the supermassive BHs at the centre of the galaxy M87* and the Sgr A*.We present our results obtained here in the Fig. 18.In this figure,it is observed that the angular diameter of BH shadow decreases with the increasing value of parameter β.It is also observed that the angular diameter of shadow for M87* and Sgr A* BH in the context of the SBR black hole BH is smaller than the other ordinary astrophysical BH such as Schwarzschild BH.Consequently, the SBR black hole parameter β has been constrained as −0.00038 < β < 0.00089 for Sgr A* and −0.002 < β < 0.0003 for M87* respectively.

VI. CONCLUSIONS
In this work, we have discussed tidal force effects and constrained the shadow of Schwarzschild type black holes in Starobinsky-Bel-Robinson modified theory of gravity.This work has been done in an attempt to define some observational signatures that can be used to distinguish SBR gravity from Einstein's general theory of relativity.The metric for spherically symmetric and static spacetime in this gravity differs from the Schwarzschild metric by introducing a stringy gravity parameter β which has a positive value.
We see that both the radial and angular tidal force profiles depend on the value of β and approach infinity as r → 0. Both force profiles show a varying trend where the radial force becomes compressive and angular tidal force becomes of stretching nature as a particle radially falls towards the black hole.This is in contrast with Schwarzschild black hole case in general relativity, where the radial tidal force monotonically increases till infinity showing stretching and the angular tidal force monotonically increases till infinity and shows compression.Both tidal forces vanish at two points outside the event horizon for different values of β giving us an opportunity to observe this phenomenon in some cases.In case of SBR gravity, the event horizon of black hole forms prior to reaching a radial distance of r = 2M .Another important observation that can be made from the radial and angular force profiles is that the the peak value is attained farther away when compared with the Schwarzschild black hole case.
In addition to this, we also examined the geodesic deviation for a radially in-falling particle in the Schwarzschild like spacetime in SBR gravity.
We consider the initial condition where a body is released from rest at a far away point where the spacetime is Minkowskian.It is observed that the radial deviation vector acts very similar to that of the Schwarzschild black hole case i.e. increasing as the particle approaches the event horizon, but shows an oscillating trend outside the event horizon.Opposite behaviour can be seen for angular deviation.It mimics the Schwarzschild black hole as the particle approaches the event horizon, but shows a sudden peak just around it and continues to rapidly fall to zero.Also, we have studied the tidal forces and geodesic deviation for Sag A* and M87* in SBR gravity by constraining the value of β using shadow size from EHT observations.We can see that in the particular β range, the tidal force and geodesic deviation profiles and both Sag A* and M87* in SBR gravity perfectly co-incide with the Schwarzschild black hole case.Furthermore, we also calculated the angular diameter of the shadow in SBR black hole and compared it to Schwarzschild BH.It is observed that the angular diameter of shadow for M87* and Sgr A* BH in SBR black hole is smaller than the Schwarzschild BH.Consequently, the SBR black hole parameter β has been constrained as −0.00038 < β < 0.00089 for Sgr A* and −0.002 < β < 0.0003 for M87* respectively.It suggests that such black holes satisfy the EHT observational constraints, and it may be possible to detect the SBR black hole and distinguish it from the other astrophysical BHs in the future.

FIG. 2 .
FIG.2.Radial tidal force as a function of the radial coordinate.We notice that the force vanishes at a single point in contrast to what happens with Schwarzschild spacetime.We take M = 1.Dotted line shows the location of event horizon for Schwarzschild spacetime

FIG. 5 .
FIG. 5. Radial tidal force profile for M87* as a function of the radial coordinate.Dotted red line shows the location of event horizon for Schwarzschild spacetime

FIG. 6 .FIG. 7 .
FIG.6.Angular tidal force as a function of the radial coordinate.We notice that the force vanishes at a single point in contrast to what happens with Schwarzschild spacetime.We take M = 1.Dotted line shows the location of event horizon for Schwarzschild spacetime

FIG. 10 .
FIG.10.Radial geodesic deviation profile for a radially infalling particle in SBR black hole spacetime.We use M=1 to show the trend for different values of β.Dotted line shows the location of event horizon in Schwarzschild spacetime.

FIG. 11 .FIG. 12 .
FIG.11.Radial geodesic deviation profile for a radially infalling particle for higher values of black hole mass M. We use β = 0.0033.

FIG. 13 .
FIG.13.Radial geodesic deviation profile for M87* using the constrained values of β.The dotted red line shows the location of the event horizon of Schwarzschild spacetime.

FIG. 14 .
FIG. 14. Angular geodesic deviation profile for a radially infalling particle in SBR black hole spacetime.We use M=1 to show the trend for different values of β.Dotted line shows the location of event horizon in Schwarzschild spacetime.

FIG. 15 .FIG. 16 .
FIG.15.Angular geodesic deviation profile for a radially infalling particle for higher values of black hole mass M. We use β = 0.0033.

FIG. 17 .
FIG. 17. Angular geodesic deviation profile for M87* using the constrained values of β.The dotted red line shows the location of the event horizon of Schwarzschild spacetime.
) Event horizon location in SBR black holes for different values of stringy parameter β.We take M=1 for computation.