Observational appearances of magnetically charged black holes in Born–Infeld electrodynamics

In this paper, we investigate the observational appearances of magnetically charged black holes in Born–Infeld (BI) electrodynamics. We examine the effects of the magnetic charge and the BI parameter on the geodesics with different impact parameters. Using the backward ray tracing method, we investigate how spherically symmetric accretions interact with black hole shadows and photon spheres. The shadows of infalling accretion are darker than that of static ones. Moreover, the radius of the photon sphere is an intrinsic property of the spacetime independent of accretions. We then study how the thin disk models affect the black hole shadows. After obtaining the transfer functions, we divide photons emitted from the thin disk into three categories: direct emission, lens ring, and photon ring. Applying three emission models, we find that the width of the shadow is dominated by the direct emission, the photon ring and the lens ring can hardly be identified by changing the emission models.


I. INTRODUCTION
Recently, the Event Horizon Telescope (EHT) Collaboration processed images of the supermassive black hole (BH) Sgr A* at the center of the galaxy [1][2][3][4][5][6].Similar to the images of M87* [7][8][9][10][11][12], these images are relatively dark in the middle and have a bright ring.The central dark region is known as the shadow of the BH, and the bright ring corresponds to the photon sphere [13].The black hole's gravity bends the trajectories of photons emitted by a distant light source behind the black hole, forming a shadow and photon sphere.Modern astrophysical observations suggest that the emissions from BHs are mainly disk-like accretions [14,15].The times of photons intersecting with the accretion disk play an important role in the optical appearances of shadows.The photon ring consists of light rays intersecting the accretion disk three or more times [16].
Observing BH shadows is an important way to detect BHs directly, which can enhance our understanding of their properties [17][18][19].Furthermore, the structures of shadows can be used to examine general relativity in strong gravitational fields [20].As a result, analyzing the images of BH shadows and photon spheres has become one of the important topics.Schwarzschild BH's shadow was first studied in [21].Sooner, Bardeen et al. studied Kerr BH's shadow [22].Through analytical calculation, the first figures of BH shadow with a rotating accretion disk are proposed in [23].In some modified gravity theories, the shadows of rotating BHs were also investigated [24,25].Compared to the EHT observations for M87* and SgrA*, there have been studies scrutinizing visible shapes of these black holes for the brightest point in accretion disk [26,27].Physical origin of the dark spot at the SgrA* image was studied [28].For invisible black holes, their silhouettes was also investigated [29].The research on BH shadows has been expanded to many aspects in recent years [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47].
There are some pieces of evidence suggesting that magnetic fields may exist near BHs.After two years of in-depth research and data processing, the EHT Collaboration unveiled the shadow of M87* with polarized light [48,49].The image clearly showed that the shadow of the BH is affected by magnetic fields, implying there probably exists a strong magnetic field surrounding the BH.Just recently, the motion of the hot spot orbiting Sgr A* is predicted to be relevant to the magnetic fields surrounding the BH [50].In recent decades, there have been many theoretical studies on magnetic black hole properties, like formation, evolution and classical instability [51][52][53].
To eliminate the singularity of point-like charge in Maxwell electrodynamics, Born and Infeld proposed the BI electrodynamics [54], inspired by special relativity and modified the Lagrangian of electromagnetic fields.In addition to BI electrodynamics, many other theories have been proposed to modify Maxwell electrodynamics (MED) to avoid singularity [55][56][57][58].These Maxwell's theory extensions are usually known as nonlinear electrodynamics (NLED).Among NLEDs, BI is distinguished for it is the only one that ensures no birefringence [59].Surprisingly, BI Lagrangian can be exactly derived from low energy string theory [60].In addition, BI has been studied as coupling with general relativity [61].In recent years, Einstein-Born-Infeld gravity has been studied in various papers [62][63][64][65][66][67][68][69][70][71][72][73][74][75][76].Magnetically charged BH has been studied with the presence of BI NLED [77], as well as other NLEDs [58,78].Considering two NLED models and comparing the shadows of these models to the observation of M87*, an upper bound on the magnetic charge of BHs was given [79].The potential astrophysical signatures of magnetically charged BHs are also studied [80].
According to [81], the magnetic charge strongly affects the outer horizons and the ergospheres of rotating BHs, but the circular photon orbits are not explicitly dependent on magnetic charge.Furthermore, the effects of strong magnetic fields on shadows are examined [82].Theoretically, photons move along the null geodesics of the spacetime governed by the BH.Being affected by NLED, photons move along the null geodesics of the effective metric rather than the background metric [83,84].Besides, it was shown that magnetically charged BH can have regular solution [85] when NLED becomes Maxwell electrodynamics at weak fields.And a good example is the rational NLED model first introduced by Kruglov [86], which has not only limited electric field but also the possibility to regard the mass of the electron as pure electromagnetic energy.After that, BH shadows' radii with the presence of rational NLED are investigated, which results are compatible with M87* data [84].We can compare the basic properties and observational appearance of BI BHs with magnetic regular BHs in rational NLED.
This paper investigates the observational appearances of magnetically charged BI BH.The shadows and photon spheres of BHs are derived using the backward ray tracing method [23].Additionally, we investigate the impact of the accretion's dynamics and shapes on the shadows.This paper is organized as follows.In section II, we investigate the metric and the trajectories of photons deflected by a magnetically charged BI BH.Additionally, we discuss the range of the magnetic charge and the BI parameter.Section III studies the shadows and photon spheres with spherical accretion.In section IV, the shadows produced by thin accretion disks are investigated.In section V, we discuss and conclude our works.Moreover, we compares the basic properties and shadows of BI BHs with magnetic regular BHs in rational NLED in Appendix A.

A. The metric
In this section, we derive the metric of BI BHs with magnetic charge in static spherically symmetric spacetime.The Einstein-Born-Infeld action is defined as follows [87] where R is the Ricci scalar, F = F µν F µν .The BI Lagrangian L(F ) takes the form as where β is called the BI parameter with the dimension of mass.When β → ∞, L(F ) degenerates into the standard electromagnetic form.Varying the action with the gauge field A µ , the BI equation yields Combining with the Einstein equation, the energy momentum tensor takes the form as Restricting to static and spherically symmetric spacetime, one can have the ansatz and the electromagnetic tensor generated by magnetically charged BH takes the form as where P is a parameter representing the magnetic charge.
The metric function f (r) for magnetically charged BI BHs is given in [88] where M is the mass of the BH, and 2 F 1 is the hypergeometric function.When β → +∞, f (r) degenerates into the form of magnetically charged Reissner-Nordström (RN) BH.
In addition, the equation f (r) = 0 has two roots.The larger one r h represents the radius of the event horizon.Since it is not easy to obtain analytical solutions of r h , its numerical solutions of different β are listed in Table 2. Setting the minimum of f (r) to zero will yield the critical value of P with respect to different β, as Fig. 1 and Table 1 show.

B. The innermost stationary circular orbit
Next, we derive the innermost stationary circular orbit (ISCO) for massive particles.The Lagrangian of a particle with unit mass, L m , has the form as where the dot on x µ represents the derivative with respect to the proper time τ .For massive particles, there are two conserved quantities: energy E m and angular momentum L m , which are obtained by Substituting Eq. ( 10) into Eq.( 9), noting that L m = −1, we derive the equation of radial motion We define the effective potential of massive particles to be The radius of the ISCO satisfies Solving Eq. ( 13) numerically, we obtain the radii of the innermost stationary circular orbit r ISCO for different BI parameters β and magnetic charges P in Tables 2 and 3.
We calculate the values of r ISCO − 3r h for different parameters β and P in Tables 2 and  3 to compare the spacetime of magnetically charged BI BH to that of Schwarzschild BH, for r ISCO = 3r h in Schwarzschild spacetime.The non-monotonicity of r ISCO − 3r h shown in Tables 2 and 3 implies that there is a critical value of β that determines at which r ISCO − 3r h switches from negative to positive.Fig. 2 is the contour plot of r ISCO − 3r h in β-P space.We start from β = 0.05 rather than β = 0 since numerical instability occurs at β close to 0. Changes in β and P have contrary impacts on r ISCO − 3r h on each side of the critical values.In the red region, r ISCO − 3r h increases with P and β.In the blue region, the impact of P on r ISCO − 3r h is opposite to that of β.In other words, the impact of P is contrary to that of β if it is less than the critical value and vice versa.When P = 0, a magnetically charged BI BH degenerates into a Schwarzschild BH, regardless of the BI effect, as the leftmost bolded line in Fig. 2 shows.

C. The trajectories of photons
In observations, for black hole shadows in nonlinear electrodynamics, photons dispersion relations need to be revised due to NLED effects.From a different perspective, the dispersion relations can be interpreted as that of the null vector with the geometry described by the effective metric.It means that the motion of lights can be viewed as electromagnetic waves propagating through a classical dispersive medium in a nontrivial vacuum [89].The effective metric for magnetic BI black holes takes the form as [88,90,91]  where The BI Parameter β characterizes the strength of NLED effects.In the limit β → ∞, both the metric function in Eq. ( 7) and the effective metric function in Eq. ( 14) reduces to the ordinary one for the RN black hole.Therefore, for magnetic black holes if one uses ordinary metric without NLED effects, the observation appearance reduces to the one for magnetic RN black holes, which has been extensively studied [92][93][94][95].Thus, we can derive the trajectories of photons.The Lagrangian of a photon is where G µν are the covariant components of the effective metric, and the dot on x µ represents the derivative with respect to the affine parameter λ.We simply focus on the equatorial plane (θ = π/2 and θ = 0) because the metric is spherically symmetric.Moreover, since the Lagrangian does not explicitly depend on coordinates t and φ, there are two conserved quantities along the geodesics For photons L = 0, one can derive the following equations representing the motion of photons The equations of photons' motion are derived by redefining the affine parameter λ into λ/L and setting the impact parameter b = L/E.
where ± in Eq. ( 20) represents the counterclockwise and clockwise motion of the photon respectively.We define the effective potential V (r) = f (r)/h(r).Specially, when ṙ| r=r ph = 0 and r| r=r ph = 0, it suggests that the photon can circle the BH for infinite times in an unstable orbit called photon sphere, which indicates where r ph is the radius of the photon sphere, b ph is the corresponding impact parameter.Since it is not easy to derive an analytical expression of r ph and b ph , we list their numerical solutions in Tables 2 and 3.As Table 2 shows, r h , r ISCO , r ph , and b ph decrease monotonically when β increases.Despite that r ph drops steeper than r h as β increases, r h < r ph is always satisfied because a BI BH degenerates into an RN BH when β → ∞.The effective potentials with respect to different parameters β and P are shown respectively in Figs.3(a) and 4(a).The trajectories of photons can be derived by combining Eqs. ( 20) and ( 21), yielding where ± represents the photons moving radially inward and outward, respectively.We can obtain the null geodesics by integrating Eq. ( 23).For b > b ph , photons approach the BH before reaching the closest approach r = r 0 , then escape towards infinity.The closest approach r 0 can be obtained by solving F (r) = 0.In the backward ray tracing method [23], r 0 is significant for determining the photons' trajectories as it can be the lower integral limit.Meanwhile, for b < b ph , photons approach the BH continuously before falling into it.We set the photons incident from infinity on the x-axis.Firstly, we study the geodesics for fixed P = 0.6 and M = 1 as β varies.The numerical solutions of r h , r ph , and b ph are shown in Table 2. Interestingly, r ph and b ph decrease dramatically with the increase of β when it is less than 0.2.When β increases, r ph and b ph slowly decline and gradually converge to the case of RN BH.Therefore, we select β = 0.05, 0.1, 0.2 to investigate the photons' trajectories in Fig. 3.
In Figs.3(b), 3(c), and 3(d), the black disk represents the BH.And the orange and grey trajectories represent photons originating with b > b ph and b < b ph , respectively.In the case of photons coming from b = b ph , they continually approach the BH and keep rotating around it in an unstable orbit located at r = r ph , i.e., the photon sphere, which is shown by red dashed lines.The BI parameter significantly affects the radius of the photon sphere and the curvature of the geodesics.However, as β increases, the decrease in r h becomes so insignificant that it is scarcely visible in those figures.
Then, we study the influence of magnetic charge on geodesics with fixed β = 0.1 and M = 1.Using the same method, we list the numerical solutions of r h , r ph , and b ph in Table 3.The geodesics of different magnetic charges are plotted in Fig. 4(b), 4(c), and 4(d).As the BH gets more magnetically charged or the BI effect becomes stronger, the range between the photon spheres and the event horizon extends.

III. SHADOWS AND PHOTON SPHERES WITH SPHERICAL ACCRETIONS A. Stationary spherical accretions
In this subsection, we study the BH shadows of static and spherically symmetric accretions.We employ the backward ray tracing method to study the specific intensity received by a distant observer [23].The emissivity of photons can be expressed as [38] j(ν e ) ∝ ρ(r)P (ν e ), where ρ(r) is the density of photons at a given radius r, and P (ν e ) is the probability of photons emitting at a given frequency ν e .For spherical accretions, we assume that the density of photons follows a logarithmic normal distribution [96] where γ is a constant that affects how fast the density of photons decays with respect to the distance from the BH, r m is the median of the logarithmic normal distribution.Since the photon density reaches the maximum at the photon sphere, we select r = r ph as the highest value of the probability density of the lognormal distribution.Solving dρ(r) dr r=r ph = 0, one can obtain We introduce a non-monochromatic light source satisfying a normal distribution with central frequency ν c , which takes the form as In this paper, we only take the frequency within ν c − √ 2σ < ν < ν c + √ 2σ into account and ignore the rest.The probability of a photon whose frequency falls within the aforementioned range is Combining Eqs. ( 24), ( 25), (26), and (28) together, one can derive Integrating Eq. ( 29) along the geodesics and setting the coefficient of proportion to be unit, we obtain the specific intensity where g = ν o /ν e is the redshift factor, ν o is the frequency received by the observer, dl prop is the infinitesimal proper length observed in the instantaneous rest frame of the photon.The redshift factor g is given by where p α and p β are the 4-momenta of the photons received and those emitted, respectively.And u α o is the 4-velocity of the observer, u β e is the 4-velocity of the static accretions.Setting the observer at r = ∞, from G µν u µ e u ν e = 1 we derive Substituting p t = −E into the Hamiltonian of photons H = G µν p µ p ν = 0, and redefining the affine parameter λ to λ/L, we derive where the first positive and negative sign, respectively, represent photons moving radially inward and outward.And ±b represents photons moving counterclockwise and clockwise.Substituting Eqs. ( 32) and (33) into Eq.( 31), we get And the infinitesimal proper length can be obtained as follows where dφ/dr = 1/F (r), can be obtained from Eq. ( 23).Substituting Eqs. ( 29) and (35) into Eq.( 30), we obtain the specific intensity I of the impact parameter b.For different decay factors γ, we plot the specific intensity in Fig. 5.For convenience, we set γ = 1 in the following discussions.We set the inclination angle of the observer to be zero due to the spherically symmetric of the spacetime.The impact parameter b determines the light ray's trajectory and thus determines the observed intensity [97].Thereon, we can obtain the specific intensity I(b) from Eq. ( 30) and then plot the simulated shadows in Figs. 6 and 7.The specific intensity is represented by several colors, which are brighter if the specific intensity is greater, and vice versa.The specific intensity and the shadows with static accretions when P = 0.6 and M = 1 with respect to different β are shown in Fig. 6.The photon spheres are the brightest rings in Figs.6(b), 6(c), and 6(d).And the radii of the photon spheres decrease continuously with the increase of β, but the decreasing speed gradually slows down, which is consistent with Table 2. Also, we discuss the specific intensities and shadows in the case of β = 0.1, M = 1 as P varies in Fig. 7 by applying the same method above.The radii of the photon spheres obviously increase as the magnetic charge increases, as described in Table .3. In summary, the average radial distance between the photons received by the observer and the BH would increase, causing a rise in the specific intensity, if the BH gets more magnetically charged or the effect of BI becomes more significant.
Additionally, We summarized and compared our results for magnetic BI BH with magnetic regular BH in rational NLED [86] (Appendix A).Observational appearances of regular BH in rational NLED were also investigated with different magnetic charge and rational NLED parameter β r .We infer that BI and the rational NLED model have similar effects on the observational appearances of magnetic black holes with static spherical accretions.

B. Infalling spherical accretions
In this subsection, we study the shadows of free-falling spherically symmetric accretions in a static and spherically symmetric spacetime.The accretions have radial velocity pointing to the BH in this model.Similarly, we also employ Eq. (30) to describe the intensities of the shadows.In infalling spherical accretions, the redshift factor g can be calculated through Eq. ( 31).The 4-velocity of a spherically symmetric infalling accretion influenced by the magnetic charged BI BH is [98] u For simplicity, we set A(r) = k(r)f (r) and B(r) = 1/[k(r)f (r)].The observer's position is fixed as the same as that in section III A, through the same method in obtaining Eqs.(32) and ( 33), we derive where the ± sign represents the photons moving radially inward and outward respectively.Therefore, we obtain the specific intensity I with the impact parameter b in Fig. 8 by combining Eqs. ( 30) and (37).The specific intensity in infalling accretion decreases slower compared with the static one.We extend the plot range to b ≤ 25 to display the properties of shadows more comprehensively.
As shown in Fig. 8, the specific intensities of infalling accretion are obviously darker than those of the static one.Compared to static accretion, the infalling one has a smaller redshift factor g. As a result, the impact of redshift on infalling accretion is more obvious than that on static accretion.If the photons' impact parameter b is small, their trajectories will get close to the BH.As a consequence, the integration along the geodesics contributes less to the specific intensity, resulting in a significantly lower value.Interestingly, the specific intensity of infalling accretion decreases inapparently as b → 0 when b < b ph , compared to that of static accretion.That is because the average radial distance between photons and the BH increases as the impact parameter increases, and the impact of the redshift factor on photons attenuates with their distance from the BH.
Additionally, the specific intensity of either static or infalling accretion, converges to zero as b approaches infinity.Firstly, the average distance between photons and the BH positively correlates to the impact parameter b.And the photon density decreases with the distance from the BH.As a result, the density of photons is less when the impact parameter is greater.Secondly, as illustrated in the end of the previous paragraph, the effect of the redshift factor on photons becomes weaker with increasing distance from the BH.
It is also noteworthy that, in both static and infalling accretion, the impact parameter corresponding to the peak of the specific intensity is always the same.This is because, in this paper, we ignore the other factors that may affect the photons' equivalent metric, such as the electromagnetic interaction between the accretion matters and photons [96].
We also investigate the impact of different magnetic charges on the specific intensity with infalling spherical accretion with fixed β = 0.1 in Fig. 9. Due to the slow increase of b ph as the BH gets more magnetically charged (see Table 3), the differences between Figs. 9(a) and 9(b) are not as obvious as that between Figs. 8(a the observational appearance shrinks as the BH becomes more magnetically charged or the influence of BI becomes more significant.The reasons are the same as our discussions in the end of the penultimate paragraph, section III A.

A. Direct emission, lens ring, and photon ring
In this section, we study the shadows cast by a thin accretion disk.The accretion disk can be placed on the equatorial plane by defining that the z-axis is perpendicular to the disk plane, for spherically symmetric spacetime.Inspired by [16], the trajectories of photons can be classified into three categories according to the total laps n that orbit the BH.
The trajectories with n < 3/4 are defined as direct emission, where photons intersect the accretion disk only once.The lens ring, where photons pass through the accretion disk twice, corresponds to the orbit's laps within 3/4 < n < 5/4.Photons on the photon ring travel across the accretion disk at least three times, satisfying n > 5/4.The ranges of direct emission, lens ring, and photon ring with respect to the impact parameter b are numerically solved in Table 4 and Fig. 11 with different parameters β and P .Fig. 10 is a visualization for Table 4.   4. The grey, red, and yellow ranges are direct emission, lens ring, and photon ring, respectively.
increase of P , but the ranges of the lens ring and photon ring both decrease.Interestingly, all these three ranges are shifted outwards with the increase of P or the decrease of β.Hence, in the presence of magnetic charge, the increase of BI effect obviously strengthens the effect of magnetic charge, which is similar to the result of the blue region of Fig. 2 in section II B. Changes in the BI parameter and the magnetic charge thus have contrary impacts on these three regions.We plot the photons' trajectories close to the BH in Fig. 12 to illustrate how the magnetic charge and BI effect affect these three ranges.From Fig. 12 we note that when β is fixed, the increase of P significantly extends the photon ring, and slightly reduces the radius of the event horizon.Meanwhile, the range of the lens ring becomes narrower.If the magnetic charge is fixed, the increase of the BI parameter would expand the range of the lens ring.Interestingly, the impact of weakening BI effect is similar to that of reducing the magnetic charge of BH, which is consistent with our previous analysis.

B. Observed specific intensities and transfer function
We assume that the thin accretion disk emits isotropically in the rest frame of static worldlines, and the accretion matter is set to be fixed.The accretion disk is placed at the equatorial plane (the y-axis in Fig. 12), and the observer is located at infinity above the north pole (x = +∞ in Fig. 12).The specific intensity received by the observer with emission frequency ν e is [16] I obs (r, ν o ) = g 3 I emi (r, ν e ), g = k(r)f (r), (38) where I emi (r, ν e ) is the emitted specific intensity by the accretion disk, ν o is the observed frequency, and g is the redshift factor derived in Eq. (34).The total specific intensity I o (r) can be obtained by integrating all frequencies of I obs (r, ν o ), which is denoted as where we have defined I e (r) = +∞ 0 I emi (r, ν e )dν e as the total emitting specific intensity.We assume the disk is optically thin, i.e., no photons are absorbed by the disk.When the photons are emitted from the accretion disk, the light rays are bent by the BH and may intersect with the accretion disk for several times.For the light rays in the regions of direct emission, they intersect with the accretion disk only once.The photons with 3/4 < n < 5/4 cross the equatorial plane twice.The light rays intersecting the accretion disk three or more times form the photon ring.
The sum of the observed specific intensities from each intersection determines the total received intensity, for photons contribute the specific intensity received by the observer once every time they pass through the accretion disk [99], yielding The r n (b) above is the transfer function, where n is a positive integer.It is the radial position where the photons cross the accretion disk for the n-th time.The slope of a transfer function is called the demagnification factor [16].At given impact parameter b it represents how much the corresponding part of the accretion disk is demagnified in the observer's vision at infinity.For different β and P , the first three transfer functions r n (b) are plotted in Figure 13.Here, we ignore the transfer functions of n ≥ 4 due to limited numerical calculation precision.It is noteworthy that the demagnification factors dr/db of these transfer functions are quite different.The first transfer function corresponds to the direct emission, which looks like a straight line with a demagnification factor of 1 approximately.So r 1 (b) is related to the redshifted source profile.The second transfer function relates to the lens ring (including the photon ring).The slope of r 2 (b) is much greater than that of r 1 (b), which indicates the lens ring is a significantly demagnified image of the back side of the accretion disk.The third transfer function is a line almost perpendicular to the horizontal axis and only contains the photon ring.Hence, the image of the photon ring is an extremely demagnified image of the front side of the disk.Thereon, the demagnification factors of further transfer functions approach infinity, and contribute negligibly to the observed intensity.Ignoring the higher-order transfer functions, therefore, has little effect on the total observed intensity.

C. Observational appearances
In this subsection, we study the specific intensity with some numerical toy models based on Eq. (40).The first emissive function I e1 is assumed to be an exponential distribution that starts from r = r ISCO .
where I 0 is the maximum intensity, similarly hereinafter.It is plotted in the top-left figures in Figs.14, 15, and 16 for different parameters β and P .It reaches its peak at r = r ISCO , and plummets to zero as r < r ISCO .
The second emissive function I e2 is an inversely proportional function multiplied by an exponential distribution, which starts from r = r ph .We call it power law exponential model for convenience in the following text.The third emissive function is assumed to be a bell-shaped function starting from r = r h .
The third emissive function is plotted in the top-right figures in Figs.14, 15, and 16 for different parameters β and P .After reaching its maximum at r = r h , it decreases slowly at first and then rapidly.By employing Eq. ( 40) to these emissive functions, we numerically calculate the observed intensities, which are then plotted in the middle row in those figures corresponding to these three emissive functions.The panels in the second row of Figs.14, 15, and 16 correspond to the observed intensities of different parameters and emissive functions.For the exponential model, the corresponding impact parameter for the first peak of the observed intensity is about b ≈ 6.1 for β = 0.1 and P = 0.6 (b ≈ 5.3 for β = 0.2 and P = 0.6; b ≈ 5.5 for β = 0.2 and P = 0.9), which is located at the photon ring.The lens ring causes the second, considerably higher peak that occurs immediately after the first one.In the middle-left panels of Figs. 14 and 15, the observed intensities reach their maximum after the second peak.For β = 0.1 and P = 0.6, the third peak of the observed intensity's corresponding impact parameter is approximately b ≈ 6.6 (b ≈ 6.1 for β = 0.2 and P = 0.6), representing the peak observed intensity produced by direct emission.In the case of β = 0.2 and P = 0.9, there are only two peaks in the observed intensity, but the second peak is higher than those in the other cases.The reason is that the lens ring is merged with the peak corresponding to the direct emission.In the bottom-left panels of Figs. 14 and 15 we can find a narrow ring just within the brightest ring, which is the lens ring.One may zoom in to find the extremely thin photon ring just within the lens ring.The brightness and width of the whole shadow are dominated by direct emission in the exponential model.For the power law exponential model, all the observed intensities shoot up to their first peak at b ≈ 4.9, b ≈ 4.1, and b ≈ 4.3 in the central panel of Figs.14, 15, and 16 respectively.These local maximums are related to direct emission.The observed intensities then steadily decline until they reach the maximum produced by the combination of the photon ring and the lens ring.The bottom-center panels of those three figures show that the lens ring and photon ring are combined to generate a halo with a highly luminous but extremely narrow profile.And the direct emission has a low observed intensity but a wide range.Among these figures, the lens ring and the photon ring are indistinguishable.For this model, the brightness of the whole shadow is dominated by the combination of the photon ring and the lens ring, but the direct emission contributes the most to the width.
In the middle-right panels of Figs.14, 15, and 16, three peaks appear from left to right.The luminosity climbs to the first maximum produced by the direct emission.After the first peak, the observed intensities drop gradually except for a dramatic rise to their maximum and fall immediately.The maximum is approximately b ≈ 6.1 for β = 0.1 and P = 0.6 (b ≈ 5.3 for β = 0.2 and P = 0.6; b ≈ 5.5 for β = 0.2 and P = 0.9), which corresponds to the lens ring.The photon ring is a lower peak very close to the lens ring.Although we can identify the lens ring and the photon ring through I o3 (b) curves, it is hard to distinguish them in the simulated shadow images.As a result, the lens ring and photon ring produce a thin and bright halo, and the direct emission contributes a relatively dark but wide ring.For the bell-shaped model, the lens ring, combined with the photon ring, dominates the maximum brightness of the whole shadow, while the direct emission contributes most to the width.
Combining Figs.14, 15, and 16, we arrive at some interesting conclusions.Firstly, it is easy to notice the differences in the shadows of magnetically charged BI BHs cast by the accretion disks starting from various radii.The central dark area's radius strongly correlates with the model.Secondly, the BI effect and the magnetic charge both have a considerable effect on the shadows.It is noticeable that the increase of the magnetic charge intensifies the luminosity produced by direct emission apparently, and the increase of the BI parameter weakens the effect of the magnetic charge.Lastly, the third is the most consistent with astronomical observations.Different emission models have a large impact on the observed intensity of direct emission but have little impact on the lens ring and the photon ring.This is due to the small range of the impact parameter related to the lens ring and photon ring.The range of the radial distance from the BH of photons emitted within the aforementioned range of the impact parameter is also narrow.The shapes and locations of the lens ring and the photon ring are almost unchanged with the emission models, for the emissivity of the accretion disk is almost constant in such a short range of the radial distance.The direct emission, however, relates to a wide range of the impact parameter.Photons with an impact parameter in this range have a wide range of radial distance from the BH when emitted from the accretion disk.The shapes and locations of observed intensity corresponding to the direct emission are also changing because the emission intensity of the accretion disk varies with the change in radial distance.Actually, the lens ring and the photon ring of the three models are difficult to identify, for the ranges of the impact parameter corresponding to the transfer functions of n = 2 and n = 3 are very narrow.In this case, regardless of the accretion model, the emission intensity is nearly constant in such a narrow range of radial distance.Therefore, it is difficult to distinguish the photon ring and the lens ring only by changing the emission model.

V. CONCLUSION AND DISCUSSION
In this paper, we investigated the observational appearances of magnetically charged Born-Infeld black hole.Firstly we presented the metric for BI BH with magnetic charge, and studied the existence of black holes with unit mass in the parameter space of the BI parameter β and magnetic charge P .The critical value of the magnetic charge increases with the increasing effect of BI nonlinear electrodynamics.The radii of the event horizon, the innermost stationary circular orbit and the photon sphere, i.e., r h , r ISCO , and r ph , and their corresponding impact parameter b ph decrease as β increases.r ph and r ISCO monotonically decrease as P increases, but r ph and b ph increase consistently.
To compare the spacetime of a magnetically charged BI BH with that of the Schwarzschild BH, we further investigated the values of r ISCO − 3r h and found that β has critical values for each magnetic charge.The value r ISCO − 3r h increases monotonically as the BI effect weakens, i.e., the increase of β.But the increase of the magnetic charge decreases the aforementioned value if β is below the critical value, and vice versa.The geodesics are then derived from the equivalent metric of photons.The photon trajectories are null geodesics in the effective metric rather than the background metric because of the effects of BI.We found that the photon geodesics bent more as the BI parameter decreases or the magnetic charge increases.
BH shadows are produced by photons that are deflected in curved spacetime.We first consider static spherically symmetric accretion, assuming that the density of photons follows a lognormal distribution with a maximum at the radius of the photon sphere.Also, we supposed that the frequency of the photons emitted is not monochromatic but has a normal distribution with central frequency ν c .The specific intensity received by a distant observer was obtained by integrating along the geodesics using the backward ray tracing method.And we plotted the specific intensity on a two-dimensional plane with the observer located at infinity above the north pole to display the BH shadows.The intensity grows slightly as b < b ph , then rises sharply to the maximum as b approaches b ph , and asymptotically drops to zero when b > b ph .The specific intensity reaches the maximum when the impact parameter of photons is equal to the critical value, as a result of the photons orbiting the BH infinite times.
We also took infalling spherically symmetric accretion models into account.The shadows of infalling accretion are obviously darker than the static ones with the same BI parameter β and magnetic charge P because the impact of the redshift factor g on specific intensity becomes significant when the photons get close to the BH.The brightest ring in each shadow represents its photon sphere in both infalling and static accretion.The radii and intensities of the shadows vary depending on the BI parameter and the magnetic charge.It is also noteworthy that the impact parameter b corresponding to the peak of the specific intensity is always the same in both static and infalling accretions, indicating that the critical value b ph is only related to the effective metric.Moreover, the radius of the photon sphere is independent of the accretions.
We then studied the shadow produced by accretion disks.The photons emitted by the disk are divided into direct emission, lens ring, and photon ring based on their numbers of intersections with the disk.The regions of direct emission, lens ring, and photon ring with parameter changes are similar to the previous analysis of r ISCO .We introduced transfer functions to describe the relationship between the radial positions where the photons intersect the accretion disk and their impact parameters.From the demagnification factors of the transfer functions, we illustrated that the slope of direct emission is almost fixed.However, the lens ring and the photon ring are noticeably demagnified.Then the observational appearances are simulated using three simple emission models of the accretion disk.Among all the simulations, the direct emission dominates the width of the shadow.The photon ring and the lens ring can hardly be distinguished only by changing the emission model.
Moreover, we compare our findings to regular BH in rational NLED in Appendix A. The comparison of regular magnetic BH in rational NLED with magnetically charged BI BH is summarized in Table 5. Observational appearances of regular BH in rational NLED with different magnetic charge and rational NLED parameter β r were also investigated.We found that the observational appearances of a magnetic BH with spherical accretions are inferred to be similar to the one for BI BH.
With increasing resolutions of astrophysical observations, the nuances between the inner radii of shadows and the theoretical solutions for r ISCO may be found, implying the existence of magnetic charge combined with the BI NLED is probable.We expect this paper to provide a reference for the investigation of shadows of BI BH and the effect of magnetic fields on BI BH.where Φ(r) is given in Eq. (A3), and u(x) = x 8 (x 4 + 1) 2 , v(x) = Through the same process as in deriving Eq. ( 34), the redshift factor of regular magnetic BH in rational NLED is given by g = u(r)Φ(r).

(A8)
We apply lognormal distribution to photons' distribution again.Combining Eqs. ( 29), (A7), and (A8), integrate the specific intensity along the geodesics like Eq. ( 30), the specific intensity I(b) received by an observer at infinity with respect to the impact parameter b for different magnetic charges P and rational NLED's effects β r are plotted in Fig. 17.
Comparing the tendency of the specific intensity of regular magnetic BH in rational NLED (Fig. 17) to that of magnetically charged BI BH (Figs. 6(a) and 7(a)), two conclusions can be verified.Firstly, the peak of the specific intensity always appears at b = b ph , for the photons rotate the BH infinite times if they incident from b = b ph .Secondly, b ph is an intrinsic property of the spacetime that is only dependent on the magnetic charge and the effect of NLED.As the tendency of the specific intensity of the two models is similar, the observational appearances of the two models are also similar.It is conceivable that changing the NLED model will not have a great impact on the observational appearances of magnetic  black holes with static spherical accretions.This point is also clarified by comparing Figs. 6, 7, and 18.

FIG. 2 :
FIG.2: Contour plot in β −P space with fixed M = 1.The colors represent the value of r ISCO −3r h , and the contours are the isosurfaces of r ISCO − 3r h .On the bolded curve, r ISCO = 3r h .The impact of the magnetic charge P is opposite to the BI parameter β below the bolded curve, and vice versa.

2 FIG. 3 :
FIG. 3: The effective potential V (r) for different BI parameter β are shown in the top-left panel.The other figures are trajectories of photons for different values of β with fixed P = 0.6 and M = 1.

9 FIG. 4 :
FIG. 4: The effective potential V (r) for different magnetic charges P are shown in the top-left panel.The other figures are trajectories of photons for different values of P with fixed β = 0.1 and M = 1.

2 FIG. 6 :
FIG. 6: The specific intensity I with respect to different impact parameters b in the top-left panel.And the other figures are the BH shadows with static accretions for different BI parameters β with P = 0.6 and M = 1.

9 FIG. 7 :
FIG. 7: The specific intensity I with respect to different impact parameters b in the top-left panel.And the other figures are the BH shadows with static accretions for different magnetic charge with β = 0.1 and M = 1.

FIG. 8 :
FIG. 8: The specific intensity I with respect to different impact parameters b in the left panels with static and infalling accretions for different β with fixed P = 0.6 and M = 1.And the shadows with these parameters for infalling accretions are shown in the right panels.

FIG. 9 :
FIG. 9: The specific intensity I with respect to different impact parameters b in the left panels with static and infalling accretions for different P with fixed β = 0.1 and M = 1.And the shadows with these parameters for infalling accretions are shown in the right panels.

βFIG. 11 :FIG. 12 :
Photon ring FIG.13:The first three transfer functions of a magnetically charged BI BH with a thin accretion disk.The black, red and yellow lines represent the first, second and third transfer functions, respectively.

( 42 )FIG. 14 :
FIG. 14: Observational appearances of a thin disk for different models with fixed β = 0.1 and P = 0.6.The first row shows the emissive functions I e (r), and the figures in the second row are the observed intensities I o (b).In the third row, we plot these observed intensities into two-dimensional disks.

FIG. 15 :
FIG. 15: Observational appearances of a thin disk for different models with fixed β = 0.2 and P = 0.6.The first row shows the emissive functions I e (r), and the figures in the second row are the observed intensities I o (b).In the third row, we plot these observed intensities into two-dimensional disks.

FIG. 16 :
FIG. 16: Observational appearances of a thin disk for different models with fixed β = 0.2 and P = 0.9.The first row shows the emissive functions I e (r), and the figures in the second row are the observed intensities I o (b).In the third row, we plot these observed intensities into two-dimensional disks.

5 :
Summary of the comparisons between a regular magnetic BH in rational NLED and a magnetically charged BI BH.
FIG.18: Observational appearances of regular BH in rational NLED with different magnetic charge and rational NLED's effects.

Black holes can exist No black holes
FIG. 1: Region plot in lg β − lg P space for magnetically charged BI BH with fixed M = 1.Black holes can only exist in the light yellow area.The critical value of the magnetic charge P increases with the decrease of β.The black dashed line (P = 1) represents extreme RN BHs.To display a wider plot range, we use logarithmic coordinates.

TABLE 1 :
Data of critical value of magnetic charge P crit for different β with fixed M = 1.The critical value of magnetic charge decreases as the BI effect weakens until β = ∞, when P crit becomes the same as magnetically charged RN BH.

TABLE 2 :
Data of r h , r ISCO , r ph , and b ph for different β with fixed M = 1, P = 0.5.Numerical solutions of different P are listed in Table3.With the increase of P , r h , and r ISCO decrease, while r ph and b ph increase monotonically.

TABLE 3 :
Data of r h , r ISCO , r ph , and b ph for different P with fixed M = 1, β = 0.1.

Table 4 ,
Figs. 10 and 11 present that the range of direct emission extends with the