Strong gravitational lensing by loop quantum gravity motivated rotating black holes and EHT observations

We investigate gravitational lensing in the strong deflection regime by loop quantum gravity (LQG)-motivated rotating black hole (LMRBH) metrics with an additional parameter l besides mass M and rotation a. The LMRBH spacetimes are regular everywhere, asymptotically encompassing the Kerr black hole as a particular case and, depending on the parameters, describe black holes with one horizon only (BH-I), black holes with an event horizon and a Cauchy horizon (BH-II), black holes with three horizons (BH-III), or black holes with no horizons (NH) spacetime. It turns out that as the LQG parameter l increases, the unstable photon orbit radius xps\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_{ps}$$\end{document}, the critical impact parameter ups\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$u_{ps}$$\end{document}, the deflection angle αD(θ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _D(\theta )$$\end{document} and angular position θ∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{\infty }$$\end{document} also increases. Meanwhile, the angular separation s decreases, and relative magnitude rmag\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{mag}$$\end{document} increases with increasing l for prograde motion but they show opposite behaviour for the retrograde motion. Using supermassive black holes (SMBH) Sgr A* and M87* as lenses, we compare the observable signatures of LMRBH with those of Kerr black holes. For Sgr A*, the angular position θ∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{\infty }$$\end{document}∈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\in $$\end{document} (16.4, 39.8) μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}as, while for M87* ∈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\in $$\end{document} (12.33, 29.9) μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}as. The angular separation s, for SMBHs Sgr A* and M87*, differs significantly, with values ranging ∈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\in $$\end{document} (0.008–0.376) μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}as for Sgr A* and ∈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\in $$\end{document} (0.006–0.282) μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu $$\end{document}as for M87*. The deviations of the lensing observables Δθ∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \theta _{\infty }$$\end{document} and Δs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta s$$\end{document} for LMRBH (a=0.80,l=2.0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a=0.80,l=2.0$$\end{document}) from Kerr black holes can reach up to 10.22μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10.22\,\upmu $$\end{document}as and 0.241μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.241~\upmu $$\end{document}as for Sgr A*, and 7.683μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$7.683~\upmu $$\end{document}as and 0.181μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.181~\upmu $$\end{document}as for M87*. The relative magnitude rmag\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{mag}$$\end{document}∈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\in $$\end{document} (0.047, 1.54). We estimate the time delay between the first and second relativistic images using twenty supermassive galactic centre black holes as lenses to find, for example, the time delay for Sgr A* and M87* can reach approximately 23.26 min and 33,261.8 min, respectively. Our analysis concludes that, within the 1σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1 \sigma $$\end{document} region, a significant portion of the BH-I and BH-II parameter space agrees with the EHT results of M87* and Sgr A*. The possibility of LMRBH being a BH-III with three horizons has been almost ruled out, except for a small portion of parameter space, by θsh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _{sh}$$\end{document} bounds of Sgr A* and M87* measured by EHT. In contrast, NH without a horizon is completely ruled out. We discover that the EHT results of Sgr A* place more stringent limits on the parameter space of LMRBH black holes than those established by the EHT results of M87*.


I. INTRODUCTION
Einstein's general relativity (GR) has been validated by the Event Horizon Telescope (EHT) observation of shadows cast by supermassive black holes M87* [1] and SgrA* [2].These also serve as further evidence of the remarkable accuracy of GR.The black hole shadow results from strong gravitational lensing of light by the intense gravitational field of the black hole.The extreme gravitational field of the black hole bends the light from the surrounding matter and stars, creating a dark area in the centre known as the "black hole shadow", and it is surrounded by a bright ring of light known as the "photon ring" which is caused by the bending and amplification of light around the black hole.Utilizing gravitational lensing can offer a robust means to investigate gravitation on a large scale.By using stronggravitational lensing by black holes and compact objects, it is possible to conduct gravity tests on a smaller scale, surpassing the limitations imposed by the Solar System [1,2].In its strong field limit, gravitational lensing exhibits intriguing characteristics, making it one of the most remarkable phenomena in astronomy.Gravitational lensing has proven to be a valuable tool for gaining insights into the structure of spacetime [3][4][5][6][7][8].In the case of strong gravitational fields, where light passes close to the source, and the bending angle is much greater, phenomena such as shadow, photon rings, and relativistic images can occur [9][10][11][12][13][14][15][16].Recent discussions have focused on strong-field gravitational lensing by black holes as another avenue to test general relativity [13][14][15][16][17][18][19][20][21][22][23].Virbhadra and Ellis [17] provided a numerical method for studying the large deflection of light rays resulting in strong gravitational lensing.Later, Bozza [13][14][15][16] and Tsukamoto [21] analyzed the strong-field gravitational lensing analytically for general spherically symmetric and static spacetimes.Gravitational lensing by black holes has been critical in quantitative studies of the lensing by Kerr black holes [13,24,25].With current observational facilities, the gravitational deflection of light by rotating black holes has received significant attention [26][27][28][29][30][31][32][33][34][35].It turns out that we can understand the properties of black holes from Parameter space (a/M, l/M ) for LMRBH spacetime [52].The red line corresponds to the extremal black holes with degenerate horizons, where the outer two horizons merge.The blue line also corresponds to the black hole with degenerate horizons, but unlike the red line, the inner two horizons merge.For the green line, ∆(r) = 0 admits two positive and the third root at r = 0, while for the dashed green line, we have r = 0 as the only root.The black dot corresponds to the black hole with three degenerate horizons located at r = 0.5M .
the gravitational lensing effect.The observables in strong gravitational lensing can diagnose the properties of black holes in modified theories of gravity and compare them with their counterparts in GR.Moreover, the photon region of gravitational lensing also provides vital properties in the black hole shadow.In addition, gravitational lensing by regular or non-singular exhibit several exciting features compared to Kerr black holes [33,34,36].The prevailing belief is that singularities are an artifact of classical GR, and they can be resolved through a quantum theory of gravity [37].Although a complete quantum gravity theory is not yet available, we must focus on regular models motivated by quantum arguments.Bardeen suggested the first regular solution for a black hole [38].In Bardeen's model, there are horizons, but there is no curvature singularity.Instead, the center of the black hole develops a de Sitter-like region, resulting in a black hole with a regular center.Since then, several regular black hole models have been proposed based on Bardeen's idea, which mimics the behaviour of the Schwarzschild black hole at large distances.It is a wide belief that loop quantum gravity (LQG) could potentially address singularities in classical general relativity [39][40][41].Because of the complexity of the complete LQG system, research has mainly focused on spherically symmetric black holes [42][43][44][45][46]. Semiclassical polymerization, which preserves the discreteness of spacetime suggested by LQG, has proved to be an effective technique for resolving the singularity issue [44][45][46].As different polymerizations can cause various types of regularized spacetimes, exploring a broader range of models and methods is of great interest.Building upon previous research [42][43][44][45][46][47][48][49], Peltola and Kunstatter [50] used effective field theory and partially polymerized theory arguments to construct a static, spherically symmetric black hole that is asymptotically flat and encompasses the Schwarzschild black hole [50,51].Notably, unlike most regular black holes with two horizons, this LQG-corrected black hole has just one horizon.However, astrophysical observations cannot test non-rotating black holes, as the black hole spin plays a critical role in any astrophysical process.Using the modified Newman-Janis algorithm, prompted us to generalise these regular solutions to the axially symmetric case or the Kerr-like solution -LQG-motivated rotating black holes (LMRBH) [52,53].Testing LMRBH metrics with astrophysical observations, like the EHT observations, is crucial.The spin is essential as it signifies the current-dipole moment of the gravitational field produced by a compact object.It serves as the primary correction to the term mass-monopole.
The prime aim of this investigation is to explore the gravitational lensing properties of recently obtained LMRBH metric [52,53] and compare them to the gravitational lensing by Kerr black holes.In addition, we investigate the observable characteristics of LMRBH versus Kerr black holes when supermassive black holes such as Sgr A* and M87* act as lenses.Notably, although strong deflection lensing effects by LMRBH black holes could be detected using the Event Horizon Telescope (EHT), distinguishing between two black holes is challenging because of deviations being on the order of O(µas).
The paper is organized as follows: We briefly review the horizon structure and calculate the deflection angle by LMRBH spacetime in Sect.II.The strong-lensing observables by the LMRBH, including the image positions θ ∞ , separation s, and magnifications µ n are also part of Sect.II.The time delay between the first and second images on the same side of the source has been calculated for supermassive black holes SgrA*, M87* and those at the centers of 19 other galaxies in Sect.II.A numerical analysis of the observables by taking the supermassive black holes Sgr A* and M87* as the lens is part of Sect.III.The constraints on the LMRBH parameters inferred using black hole shadow observational data of Sgr A* and M87* are discussed in Sect.IV.Finally, we summarize our results to end the paper in Sect.V.
Throughout this paper, unless otherwise stated, we adopt natural units (8πG = c = 1)

II. GRAVITATIONAL LENSING BY LQG-MOTIVATED ROTATING BLACK HOLES
We derived the metric for the LQG-motivated rotating black hole (LMRBH) by building upon a partially polymerized static and spherically symmetric black hole solution [50,54].To construct the rotating spacetime LMRBH, we employed the revised Newman-Janis algorithm (NJA) described in [55,56].This procedure has been successful in generating imperfect fluid rotating solutions in Boyer-Lindquist coordinates from spherically symmetric static solutions, and it can also produce generic rotating regular black hole solutions [55][56][57][58][59].We express the metric of the LMRBH in the Boyer-Lindquist form as [52,53]. where The rotating metric, derived from the NJA, captures important elements of LQG, such as a transition surface at the black hole centre and the global regularity of spacetime.It is noteworthy that LMRBH (Eq.( 1)) encompasses the Kerr spacetime [60] in the limit l → 0 and spherical LQG black hole [50,54] in the limit a → 0. When a = M = l = 0, then Eq. ( 1) gives flat spacetime.The horizons of the LMRBH are determined by the roots of the null surface ∆(r) = 0, which is a coordinate singularity of Eq. (1).Depending on the values of a and l, up to three real roots may exist, with one to three positive roots.Only positive roots correspond to horizons, and we label them as r 1 , r 2 , and r 3 with r 3 ≤ r 2 ≤ r 1 , where r 1 is the event horizon and r 2 , if present, is the Cauchy horizon.The additional root r 3 lies inside the Cauchy horizon.The parameter space (a, l) for the LMRBH is illustrated in Fig. 1.There exists a critical value of l (a) on the red line, denoted by l c (a c ), such that ∆(r) = 0 has a double root, for given a (l), corresponding to an extremal LMRBH with degenerate horizons.For a < a c (l > l c ), ∆(r) = 0 has two simple positive roots, corresponding to LMRBHs with Cauchy and event horizons (BH-II).In contrast, for a > a c (l < l c ), ∆(r) = 0 has no positive roots, corresponding to no-horizon (NH) spacetimes.Similarly, for a given a (l), we can find the critical value of l (a) on other transition lines and dots in Fig. 1.We consider four regions for our study, namely black holes with only one horizon (BH-I), black holes with an event horizon and Cauchy horizon (BH-II), black holes with three horizons (BH-III), and black holes with no horizon (NH).The coloured lines in Fig. 1 denote the boundaries that divide these regions (see Kumar et al. ([52]) for more information).The presence of a black hole has an influential impact on the motion of nearby photons.The null geodesics that describe photon orbits around black holes are crucial for observing gravitational effects caused by the black hole.Two linearly independent killing vectors, η µ (t) = δ µ t and η µ (φ) = δ µ φ , associated with the time translation and rotational invariance, are admitted by the black hole metric (1) [61].The total energy E and the angular momentum L, two conserved quantities corresponding to the Killing vectors, govern the photon's path.We identify u = L/E as the impact parameter and employ the Hamilton-Jacobi approach to establish the relationship between the effective potential V eff and impact parameter u.The relationship between u and V eff determines the orbit's qualitative characteristics.Photons with different impact parameters approach the black hole and are deflected by its strong gravitational field, reaching a minimum distance of r 0 .It is critical to understand that light rays with impact parameters smaller than the critical value u ps fall into FIG.2. Effective potential for photons, in LMRBH spacetime, having different impact parameters u.Black curves correspond to the photons whose impact parameter u is exactly equal to the critical impact parameter ups.These photons revolve around the black hole in unstable photon circular orbits at radial distance x = xps.Photons with u ≈ ups but u > ups make several loops around the black hole and are scattered to infinity after reaching some closest approach distance x0.These photons form the (strong) gravitationally lensed image of the source.Photons with u < ups fall into the black hole.the event horizon, leading to a dark spot on the observer's sky.On the other hand, light rays with impact parameters slightly larger than the critical value make several loops around the black hole and are scattered to reach the observer at infinity (cf.Fig. 2).Finally, photons with impact parameters equal to the critical value follow unstable circular photon orbits with a constant radius r ps (cf.Fig. 2).By considering light rays near the critical impact parameter, we investigate the effect of the parameter l on strong gravitational lensing using the LMRBH model to represent the black holes.The LMRBH metric, similar to the Kerr metric, exhibits reflection symmetry θ → π − θ, allowing for studying light ray motion in the equatorial plane where they initially reside.To explore the strong gravitational lensing effects, we follow the approach in Refs.[15,33,35].We measure the quantities r, a, l, and t in units of M [15] such that and use x instead of radius r to rewrite the LMRBH metric (1) in the equatorial plane (θ = π/2) as where and A photon from a source travels towards a black hole until it is at least at a distance of x 0 away from it before being redirected by the black hole's gravitational field to reach the observer at infinity.The impact parameter, which is the perpendicular distance from the black hole's centre of mass to the initial direction of the photon at infinity, stays constant throughout the trajectory because the approach phase resembles the departure phase with time-reversed due to spacetime symmetry.Furthermore, the minimum approach distance x 0 marks the turning point and hence vanishing effective potential, which gives an expression for the impact parameter u in terms of the closest approach distance x 0 as By restricting the light rays to the equatorial plane, the unstable circular photon orbits radius x ps in terms of metric components, given by Eq. ( 4), is the solution of the equation [16] A which implies The unstable photon orbit radius is the largest root of the Eq. 7. It turns out that the photon sphere depends on the parameter l and the rotation parameter.By fixing the winding of light rays to be counterclockwise, we assign a > 0 (prograde orbits), if the black hole also rotates in the counterclockwise direction and a < 0 (retrograde) if the black rotates in the opposite direction of photon winding i.e., clockwise.Fig. 3 shows the decrease in radius x ps with the rotation parameter a while the opposite behaviour concerning the parameter l and suggests that the x ps of LMRBH is greater than the Kerr black holes.Fig. 3 also suggests that the photons forming prograde orbits can get closer to the black hole than the photons forming retrograde orbits.The critical impact parameter u ps is the parameter for which the closest approach distance x 0 equals the photon orbit radius x ps .We have plotted the critical impact parameter (cf.Fig. 4) with varying a and l and found that it varies similarly to unstable photon orbit radius.
Moreover, the deflection angle of a photon moving in the equatorial plane of the LMRBH spacetime can be obtained by using the null geodesic equations, which are first-order ordinary differential equations [61].The light bending angle in a general rotating stationary spacetime described by the line element (3), for a closest distance approach x 0 is given by [15,33,35,62] This integral is non-trivial to solve, and hence, we expand it near the unstable photon sphere radius [15,17,63] by defining a new variable z = 1 − x 0 /x in strong deflection limit (SDL) [21,64].This technique not only shows the behaviour of photons near the photon sphere but also provides an analytical representation of the deflection angle as [15,65,66] α where ā, b are the lensing coefficients.The details of this calculation can be found in [15,65,66].The deflection angle increases as x 0 approaches x ps , eventually exceeding 2π radians and diverging logarithmically at x 0 = x ps .We can investigate the deflection angle of strong gravitational lensing by LMRBH and compare it with the analogous results of the Kerr black hole.Like Kerr black hole, the deflection angle for LMRBH spacetime in SDL (cf.Fig. 5) is more than 2π.The effect of the LMRBH parameter l and spin a can be seen in Fig. 5.The deflection angle diverges at larger u ps for larger l while at smaller u ps for larger a.As can be seen in Eq. ( 9) if a light beam's impact parameter is sufficiently close to its critical value, it may also approach the unstable photon orbit radius and make one, two, or even more turns around the lens before it reaches the observer.As a result, a strong gravitational field can generate many images, creating what may appear to be an endless series of images.The two infinite sets of relativistic pictures represent the clockwise and anticlockwise winds of the black hole.In Fig. 5, we show that for higher impact parameter values, the deflection angle for some parameters can become negative, which shows that the path of the light or particles is bent opposite the black hole's centre, causing photons to be deflected away from the black hole's location.But for our configuration in which the black hole lies between a light source and observer, we calculated the deflection angle using the expression in Eq. ( 9), which is valid for SDL.However, this expression of the deflection angle in SDL technique is valid only in close vicinity to the critical impact parameter.In reality, a negative deflection angle value is not achievable.Still, as is usually done in the literature, one can utilise the weak deflection limit technique to get the deflection angle at larger impact parameters [13].

II.1. Lens equation and observables
The lens geometry is essential in understanding black hole lensing because it can be used to identify the exact positions and magnifications of relativistic images in particular.Since the light source and observer are placed far  enough from the black hole for the gravitational forces around them to be insufficient, the lens configuration we are interested in positions the black hole between a light source and observer [15,67] (cf.Fig. 6).Assuming that the source and observer are almost aligned, the lens equation reads Instead of using a full deflection angle, we used an offset of deflection angle ∆α n = α − 2nπ with n being an integer with n ∈ N , and 0 < ∆α n ≪ 1.Here, β and θ are the angular positions of the source and image from the optical axis, respectively.The distances of the source and lens from the observer are given by D OS and D OL , respectively (cf.Fig. 6).
Next, we shall estimate the observables for the strong gravitational lensing by LMRBH spacetime as in [15,16,34,35].Using the lens Eq. ( 10) and Eq. ( 9), and following the condition where the source, lens and the observer are aligned, the angular separation between the lens and the n th image is given by [16] where Here θ n 0 is the angular position of the image when a photon encircles complete 2nπ and the second term in Eq. ( 11) is the extra term exceeding 2nπ such that θ n 0 ≫ ∆θ n [15].When β = 0 is entered into Eq.( 11), one can obtain the angular radius of the Einstein rings for spherically symmetric cases [3][4][5][6][7][8].Einstein rings for Sgr A* and M87* are plotted in Fig. 7.
In gravitational lensing, the light is deflected while maintaining the surface brightness, but the appearance of the solid angle changes, enhancing the brightness of the images.The magnification, for the n−loop images, is evaluated as the quotient of the solid angles subtended by the n th image and the source as [15,16]  The magnification µ n is inversely proportional to D 2 OL and so the images are faint.But they can be bright in the limit β → 0, i.e., when the source, lens and observer are perfectly aligned.When β → 0, the Eq. ( 15) diverges, suggesting that the perfect alignment maximises the possibility of the detection of the images.The brightness of the first image is dominant over the other images, as a result, we will focus on the most straightforward scenario, in which just the outermost image, θ 1 , is resolved as a single image, while the subsequent images are all crammed together at θ ∞ .In practice, if the 1-loop image can be distinguished from the rest packed images, we can have three characteristic  Estimates for the lensing observables by considering supermassive black holes at the center of nearby galaxies as LMRBHs.We measure the quantities a and l in units of the mass of the black hole, M .
In the above expression, θ 1 is the angular position of the outermost single image, θ ∞ is the angular position of the rest packed images, s is the angular separation between the θ 1 and θ ∞ , r mag is the ratio of the flux of the first image

and the all other images.
There is a time discrepancy between the two images because the time taken by the light routes corresponding to the various images is different.A further significant observable known as the time delay is the interval between the creation of relativistic images that take the source's varying brightness into account.Finally, the time delay ∆T 2,1 when the images are on the same side of the lens can be tentatively calculated as follows [68] ∆T The source must be variable to measure the time delay, which is not a stringent condition as variable stars are common in all galaxies.It could calculate the time difference between the relativistic images by implicitly assuming these fluctuations, which will manifest themselves with a temporal phase in all images.From an observational standpoint, the accurate measurement of the time delay has a significant advantage of dimensional measurement of the system's scale, which can be used to accurately estimate the black hole's distance.

III. STRONG GRAVITATIONAL LENSING BY SUPERMASSIVE BLACK HOLES, SGR A* AND M87*
Assuming that Sgr A* and M87* as characterised by LMRBH, we study the consequences of strong gravitational lensing by these supermassive black holes.By computing observables viz., θ ∞ , separation s, and relative magnification r mag for different values of parameter l, we contrast the lensing outcomes between Kerr and LMRBH black holes.We depict our results in Fig. 8, Fig. 9 and Fig. 10, while Table I and Table II show the lensing observables and their deviation for various values of a, and l in comparison with Schwarzschild (a = l = 0) and Kerr black hole (l = 0).Our analysis of LMRBH as the lens reveals that the angular position of images for Sgr A* and M87* is consistent with the EHT-measured angular shadow diameters of Sgr A* and M87*.The results in Table I and II show that in the case of LMRBH, the angular positions of images are larger than their corresponding values in GR and vary slowly concerning the position of the source β.In fact, the deviation from the Kerr black hole can go up to 10.23 µas and 7.68 µas, respectively, for Sgr A* and M87*, at a = 0.8 and l = 2.00, an effect too tiny to be observed with current telescopes.Further, the separation s in the case of LMRBH for Sgr A* and M87* range between 0.008-0.376µas and 0.006-0.282µas,respectively.The angular separation s between the first and other packed images due to the LMRBH for Sgr A* and M87* is beyond the threshold of the current EHT observation, and we may have to wait for the next generation event horizon telescope (ngEHT) for this purpose.For higher spin values, the angular separation s decreases with l while it increases with l at lower spin levels.The relative magnification of the first-order images of presented in Table I using Eq. ( 15) for black holes in GR and LMRBH.The first-order images by LMRBH are more highly magnified than the corresponding images of black holes in GR, and the magnification increases slowly with the parameter l.The LMRBH images are brighter than their spherically symmetric equivalents, as indicated by the quick decline in the flux ratio of the first image to all other images with a (cf.Fig. 10).We have also calculated the time delay for different black holes in nearby galaxies in Table III.The time delay of the first image from that of the second image, ∆T 2,1 , for the LMRBH as Sgr A* and M87* can reach up to 23.25 min and 554.36 hrs, respectively, while the deviation from the Kerr black hole for Sgr A* and M87* is 8.93 min and 212.88 hrs, respectively.Observing the time delay in Sgr A* is much shorter and more difficult for measurement.In the case of M87*, the time delay can reach up to a few hundred hours, sufficient for astronomical measurements, provided we have enough angular resolution separating two relativistic images.

IV. CONSTRAINTS FROM EHT OBSERVATIONS SHADOWS OF M87* AND SGR A*
Black holes within the optically thin accreting region are anticipated to display a dark "shadow" surrounded by a bright ring [69,70].The boundary of the shadow is determined by the photon ring, which principally relies on the black hole's parameters [71].The application of the shadow in understanding near-horizon geometry has sparked a flurry of activity in analyzing, both analytically and numerically, shadows for black holes in GR [12,34, and in LQG [53,54,[94][95][96].The EHT observation unveiled the shadows of supermassive black holes Sgr A* and M87* [1,97] whose sizes are within 10 percent of the Kerr predictions, furnishing another tool to explore the nature of strong-field gravity.One can put constraints on potential deviations from the Kerr, i.e., such as LMRBH stemming from LQG, as the quantum effects cannot be overlooked in the strong-field regime.Thus, the EHT observation results of M87* and Sgr A* shadows can restrict the LMRBH-black hole parameters.We confine the parameters (a, l) within the 1-σ level by using the apparent radius of the photon sphere (θ ∞ ) as the angular size of the black hole shadow.We model the M87* and Sgr A* as the LMRBH and use their shadow results to test the viability of our model to explain the astrophysical black holes spacetimes.
IV.0.0.1.Constraints from M87* shadow In 2019, the EHT collaboration produced a ring of diameter θ sh = 42 ± 3 µas as the first image of the supermassive black hole M87* [1].We find that, for all values of a and l, the Kerr black hole, with a mass of M = (6.5 ± 0.7) × 10 9 M ⊙ and distance of D OL = 16.8Mpc, casts the smallest shadow that is contained inside the 1-σ region [1,98,99].Fig. 11 depicts the angular diameter θ sh as a function of (a, l) for the LMRBH as M87*, with the black corresponding to θ sh = 39 µas .The LMRBH metric, when investigated with the 2θ ∞ ( μas) .The EHT observation used three independent algorithms, EHT-imaging, SIMLI, and DIFMAP, to find out that the averaged measured value of the angular shadow diameter lies within the range θ sh ∈ (46.9, 50) µas.The 1-σ interval is ∈ (41.7, 55.6) µas.The average measured value of the angular shadow diameter is within the range of θ sh ∈ (46.9, 50)µas, according to the EHT observation, which employed three independent techniques, EHT-imaging, SIMLI, and DIFMAP.The average strongly constrains the parameters a and l for the LMRBH observed angular diameter of the EHT observation of Sgr A* black hole i.e., θ sh ∈ (46.9, 50)µas, which falls within the 1 − σ confidence region with the observed angular diameter such that 0.128662 ≤ a ≤ 0.27171 and all values of l are allowed for the LMRBH.Thus, within the finite parameter space, LMRBH agrees with the EHT results of Sgr A* black hole shadow (cf.Fig. 11).

V. CONCLUSIONS
LQG faces challenges in testing its predictions because its framework lacks rotating black hole models.We recently addressed this issue, used a nonrotating LQG black hole as a starting point, and applied the modified Newman-Janis algorithm to construct LMRBH with an additional parameter l [52,53].This solution encompasses the Kerr black hole in the limit l → 0. It can describe different regions of the parameter space, including multi-horizon rotating regular black holes [52,53].The spacetime singularities are artefacts of classical GR, as the prominent singularity theorems also predicted them.However, it is a widespread belief that singularities do not exist in nature and are the limitations of general relativity.Therefore, it is instructive to explore gravitational lensing by LMRBH to assess the dependence of observables on the parameter l and compare the results with those for the Kerr black holes.
With this motivation, we have analyzed the strong gravitational lensing of light because of LMRBH, which, besides the mass M and angular momentum a, has an additional deviation parameter l.We have examined the effects of the LQG parameter l on the light deflection angle α D (θ) and lensing observables θ ∞ , s, r mag , u m , in the strong field observation, because of LMRBH and compared them to the Kerr black holes.By studying the observational implications of the gravitational lensing by this LMRBH, we can constrain the fundamental parameter l of LQG.Indeed, we can achieve the analysis of strong fields and the testing of theories of gravity using gravitational lensing, which is a powerful tool.Comparing lensing by different gravities within the strong-field limit is valuable.In this scenario, the LMRBH model provides an exciting opportunity to discuss the observational signatures of LQG over the Kerr black hole of general relativity.Our investigation focuses on the gravitational lensing of light around the LMRBH in the strong deflection limits, which reveals that photons get deflected from their straight path and result in multiple images of a source depending on the impact parameter u of the photon.At u = u ps , photons follow circular orbits around the black hole, and the deflection angle diverges.Intrestingly, the unstable photon orbit radius x ps , the critical impact parameter u ps increase with l.We also found an increase in deflection angle with l such that the LMRBH leads to a larger deflection angle than the Kerr black hole, and the deflection angle diverges at larger u ps for larger l while at smaller u ps for larger a.
By considering the spacetime to be defined by the LMRBH, we estimated the lensing observables θ ∞ , s and r mag of the relativistic images for supermassive black holes, specifically Sgr A* and M87*.LMRBH exhibit potentially apparent distinctions from the black holes in general relativity in its predictions for gravitational lensing caused by supermassive black holes.The presence of parameter l rapidly increases θ ∞ when compared to the Kerr black hole.We observe that θ ∞ for Sgr A* ∈ (16.40, 39.80) µas and its deviation from its GR counterpart can reach as much as 10.22 µas.While for M87* ranges between 12.32 -29.91 µas, and deviation is as high as 7.68 µas.The separation s decreases, with increasing l, for the prograde motion while it shows opposite behaviour for the retrograde motion.Moreover, the deviation in separation s with l for retrograde motion is much smaller than the prograde motion.The separation s due to LMRBH black holes for Sgr A* and M87* range between 0.008-0.376µas and 0.006-0.282µas, respectively with the deviation from Kerr black hole can go as much as 0.241 µas for SgrA* and 0.181 µas for M87*.Compared to the analogous images of black holes in GR, the first-order images of LMRBH are significantly magnified.With increasing l, the magnification increases for the prograde motion but shows opposite behaviour for the retrograde motion with the relative magnitude r mag ∈ (0.047, 1.54).Also, the deviation in magnification with l for retrograde motion is smaller than the prograde motion.Finally, the time delay ∆T 2,1 for Sgr A* and M87* deviate significantly from the Kerr black hole by 8.93 min and 212.8 hrs, respectively, and for the LMRBH as Sgr A* and M87*, ∆T 2,1 can reach up to 23.25 min and 554.3 hrs, respectively.The time delay in Sgr A* is smaller and more challenging to measure.If we have enough angular resolution to distinguish between two relativistic images, the time delay with M87* can be as little as a few hundred hours, which is sufficient for astronomical measurements.We show how the shadow-size measurements significantly constrain the deviation parameters for the LMRBH.From M87* results of EHT, we found that 0 < a ≤ 0.03684 and all values of l are allowed.The EHT results of SgrA* strictly allow 0.128662 ≤ a ≤ 0.27171 and all values of l.Thus, the LMRBH parameters a and l are constrained by the EHT bounds on θ sh of SgrA* and M87*, within the 1σ region, such that a significant portion of the BH-I and BH-II parameter space and small portion of BH-III agrees with the EHT results of SgrA* and M87*, whereas the NH spacetime is entirely ruled out.We have found that the results gained by the Event Horizon Telescope (EHT) for Sgr A* impose stricter constraints on the parameter space of LMRBH than the limits established by the EHT results for M87*.
The primary constraint of our approach is that the LMRBH metric is not a result of the direct loop quantization of the Kerr spacetime.The LMRBH has thrilling properties, e.g., a transition surface substitutes the classical ring singularity in the Kerr black hole and provides a singularity solution of the Kerr black hole.Thus, we can expect that LMRBH captures some description of LQG.Our analysis reveal that LMRBH being BH-I with one horizon and BH-II with Cauchy and the event horizons are also observationally favoured with no restriction LQG parameter l.
The results presented here generalize previous discussions on black hole lensing in GR and our results go over to Kerr and Schwarzschild black holes in the limits, l → 0, and a, l → 0 respectively.Although resolving the order estimated in SDL is challenging, the outlook for future observations looks bright.

VI. ACKNOWLEDGMENTS
FIG.1.Parameter space (a/M, l/M ) for LMRBH spacetime[52].The red line corresponds to the extremal black holes with degenerate horizons, where the outer two horizons merge.The blue line also corresponds to the black hole with degenerate horizons, but unlike the red line, the inner two horizons merge.For the green line, ∆(r) = 0 admits two positive and the third root at r = 0, while for the dashed green line, we have r = 0 as the only root.The black dot corresponds to the black hole with three degenerate horizons located at r = 0.5M .

FIG. 3 .
FIG. 3. Behavior of the photon sphere radius xps with respect to the parameter a (left) for different l and with respect to the parameter l (right) for different a for LMRBH spacetime.Negative values of a in the left diagram correspond to the retrograde motion of the photon.

FIG. 4 .FIG. 5 .
FIG. 4. Behavior of the critical impact parameter ups with respect to the parameter a (left) for different l and with respect to the parameter l (right) for different a for LMRBH spacetime.Negative values of a in the left diagram correspond to the retrograde motion of the photon.

FIG. 6 .
FIG.6.Formation of primary images of source space (S) in case of gravitational lensing.Light rays are deviated by the black hole to be observed at an angular position θ by the observer space (O).

FIG. 8 .
FIG.8.Behavior of strong lensing observables θ∞ in strong field limit, as a function of the parameters a and l by considering that the spacetime around the compact objects at the centers of Sgr A*(left panel) and M87*(right panel) is LMRBH spacetime.

FIG. 9 .FIG. 10 .
FIG. 9. Behavior of strong lensing observables s in strong field limit, as a function of the parameters a and l by considering that the spacetime around the compact objects at the centers of Sgr A*(left panel) and M87*(right panel) is LMRBH spacetime.

40 FIG. 11 .
FIG.11.Shadow angular diameter θ sh = 2θ∞ of LMRBH as a function of (a, l).The black and black dashed lines correspond to the Sgr A* black hole shadow at θ sh = 46.9µas and θ sh = 50 µas, respectively, such that the region between these lines satisfies the Sgr A* shadow 1-σ bound (left).M87* shadow angular diameter when considered as a LMRBH.The black line is θ sh = 39 µas, and the region contained within it meets the M87* shadow 1-σ bound (right).

TABLE III .
Estimation of time delay for supermassive black holes at the center of nearby galaxies in the case Kerr (a = 0.8) and LMRBH (a = 0.8 and l = 2.0).Mass (M ) and distance (DOL) are given in the units of solar mass and Mpc, respectively.Time Delays are expressed in minutes.
TABLE II.Deviation of the lensing observables of LMRBH black holes from Kerr black hole for supermassive black holes at the center of nearby galaxies for a = −0.40, a = 0.40 and a = 0.80.Here ∆(X) = XLMRBH − XKerr