QED and accretion flow models effect on optical appearance of Euler-Heisenberg black holes

Taking the quantum electrodynamics (QED) effect into account, we investigate the geometrical-optics appearance of the Euler-Heisenberg (EH) black hole (BH) under the different accretion flows context, which depends on the BH space-time structure and different sources of light. The more significant magnetic charge leads to the smaller shadow radius for the EH BH, while the different values of the EH parameter do not ruin it. Different features of the corresponding two-dimensional shadow images are derived for the three optically thin accretion flow models. It is shown that the total observed intensity in the static spherical accretion flow scenario leads than that of the infalling spherical accretion flow under same parameters, but the size and position of the EH BH shadows do not change in both of these accretions flows, implying that the BH shadow size depends on the geometric space-time and the shadows luminosities rely on the accretion flow morphology. Of particular interest is that a thin disk accretion model illuminated the BH, we found that the contribution of the lensing ring to the total observed flux is less than 5\%, and the photon ring is less than 2\%, indicating that the direct emission dominates the optical appearance of the EH BH. It is also believed that the optical appearance of the BH image depends on the accretion disk radiation position in this scenario, which can serve as a probe for the disk structure around the active galactic nucleus (AGN) of M87^{*} like.

Abstract Taking the quantum electrodynamics (QED) effect into account, we investigate the geometrical-optics appearance of the Euler-Heisenberg (EH) black hole (BH) under the different accretion flows context, which depends on the BH space-time structure and different sources of light. The more significant magnetic charge leads to the smaller shadow radius for the EH BH, while the different values of the EH parameter do not ruin it. Different features of the corresponding two-dimensional shadow images are derived for the three optically thin accretion flow models. It is shown that the total observed intensity in the static spherical accretion flow scenario leads than that of the infalling spherical accretion flow under same parameters, but the size and position of the EH BH shadows do not change in both of these accretions flows, implying that the BH shadow size depends on the geometric space-time and the shadows luminosities rely on the accretion flow morphology. Of particular interest is that a thin disk accretion model illuminated the BH, we found that the contribution of the lensing ring to the total observed flux is less than 5%, and the photon ring is less than 2%, indicating that the direct emission dominates the optical appearance of the EH BH. It is also believed that the optical appearance of the BH image depends on the accretion disk radiation position in this scenario, which can serve as a probe for the disk structure around the active galactic nucleus (AGN) of M87 * like.

Introduction
A black hole (BH) is the solution to Einstein's field equations of general relativity (GR) describing regions of space-time that have undergone gravitational collapse. The first convincing evidence to prove the existence of BHs in astronomical observations is the Laser-Interferometer Gravitational Wave-Observatory (LIGO) found that the gravitational wave emission from the coalescence of two BHs [1]. Another big advance with the BH observations is the Event Horizon Telescope (EHT) collaboration debuted the near horizon image around the supermassive BH in the center of the Messier (M) 87 * elliptical Galaxy (M = 6.5 ± 0.7 × 10 9 M ). The BH image reports a compact asymmetric ring-like source morphology, which present a bright ring-shaped lump of radiation surrounding a circular dark region of an estimated 6.5 billion solar masses [2][3][4][5][6][7]. It convincingly confirms the existence of BHs in our Universe. The EHT measured the corresponding linear polarimetric shadow image of the M87 * , crucial to understanding the launching of energetic jets from its core. It carries the information about the structure of the magnetic fields responsible for the synchrotron emission and found that magnetically arrested accretion disks surround the M87 * (the mass accretion rate onto the BH of (3 − 20) × 10 −4 M yr −1 ) [8,9]. More excitingly, EHT obtained the first horizon-scale radio observations of the Sagittarius (Sgr) A * in our own Milky Way, which shows the measured ring size of the Sgr A * is consistent with the shadow critical curve predicted in GR within 10% [10][11][12][13][14][15].
An astrophysical BH provides a constant space-time structure, but it can be illuminated by external sources of the luminous accretion material, leading the BH to have a variety of shapes and emit a variety of colors. The light from the accretion material is bent toward the singularity by the BH's gravitational field if the light passes close to a BH. This makes it possible to investigate the optical ap-pearance of the BH from accretion flow. Possible observational characteristics of the BHs shadows surrounded by different accretion flows were studied for a long time. By establishing a thin-disk accretion flow model, Luminet presented that the emergence of the BH shadow and ring depends on the accretion flow position since the relative intrinsic intensity of radiation emitted by the disk accretion is a function of the BH radius [16]. Falcke et al. created a ray-tracing code and proposed that the BH shadow is equivalent to the gravitational lensing effect by considering the radiation of a hot optically-thin accretion flow surrounding a supermassive BH in the center of our Galaxy [17]. Assuming that the Schwarzschild BH is surrounded by a thin and heavy disk accretion flow, Cunha et al. investigated the BH gravitational lensing effect in this model and displayed an almost equatorial observer who could observe different patches of the sky near the equatorial plane [18]. Narayan et al. studied the Schwarzschild BH shadow under a simple spherical accretion flow model, revealing that the optical appearance of the spherically accreting BH is independent of the inner radius at which the accreting gas stops radiating [19]. By considering an optically/geometrically thin disk accretion flow model, Gralla et al. investigated the shadow and ring of the Schwarzschild BH shadow and found the image brightness diverges logarithmically at the BH photon ring [20].
These experiments for the shadow image would be robust probes GR, and can further understands the properties of BHs and test the other modified theories. Considering the spherical accretion flow, Zeng et al. investigated the shadow of the Gauss-Bonnet BH and discussed the influence of the quintessence dark energy on the spherical BH. They found that the BH shadow depends on the physical properties of the accretion flow [21,22]. In the framework of the Einstein-Gauss-Bonnet-Maxwell gravity, Ma et al. explored the photon sphere and shadow of the static spherically symmetric charged BH, which found that the shadow and photon sphere satisfy the sequence of inequalities to relate a BH's horizon and mass [23]. Peng et al. argued that the BH shadow feature may provide observational evidence for the quantum effect of GR by investigating the quantum corrected Schwarzschild BH [24]. Gao et al. studied the gravitational lensing effect of a hairy BH in Einstein-scalar-Gauss-Bonnet gravity and found that the BH shadow cast in this gravity is consistent with the measurement by the EHT, while other lensing observables are beyond the present capacity [25]. Li et al. investigated the observational appearances of the global monopole BH illuminated by various accretions under the f (R) gravity theory. They found that the BH shadows and the related rings with some different observable features can be used to distinguish BHs from different gravity theories [26]. Guo et al. considered a charged BH surrounded by a perfect fluid radiation field within the framework of the Rastall gravity. They investigated the shadow and pho-ton sphere of this BH with the static/infalling spherical accretion background and obtained that the shadow luminosity of this BH with infalling spherical accretion is dimmer than that of the static spherical accretion, but the photon sphere luminosity is brighter than the static one [27]. Notably, they also showed that the luminosities of both the shadows and rings of the Hayward BH are affected by the accretion flow property and the BH magnetic charge [28]. The BH shadow characteristics in the various gravity are also investigated in [29][30][31][32][33][34][35][36][37].
One of the biggest problems in GR is the singularities that lay at the beginning of the Universe and also at the center of a BH. Maxwell's equations are known to exhibit singularities which cause the divergence problems in Maxwell's theory. Based on Dirac's positron theory, Euler and Heisenberg proposed a new approach to describe the electromagnetic field, taking into account the one-loop corrections to quantum electrodynamics (QED) and explaining the vacuum polarization in QED [38]. Using this new method, the nonlinear electrodynamic models can explain the inflation of the universe in the early times [39]. By considering the one-loop effective lagrangian density coupled with the Einstein field equation, Yajima et al. obtained the Euler-Heisenberg (EH) BH solution [40]. Kruglov explored the vacuum birefringence effect of the nonlinear electromagnetic field and calculated the field effective Lagrangian, which found that the improvement in the PVLAS experiment would correct the relation between two parameters in the effective Lagrangian [41]. Considering classical models of nonlinear electrodynamics with polynomial self-interaction of the electromagnetic field, Costa et al. proposed that the collation with QED results in the total field energy of a point elementary charge about twice the electron mass [42]. By investigating the EH-type model of nonlinear electrodynamics with two parameters, Kruglov further obtained the charged black hole solution in the framework of nonlinear electrodynamics. It also obtained the corrections to Coulomb's law at r → ∞ and studied the energy conditions [43]. It is worth mentioning that Kruglov further derived the effective geometry induced by nonlinear electrodynamics corrections and determined the shadow's size of thee regular non-rotating magnetic BHs [44].
Although the properties of BH in the nonlinear electrodynamic context have been extensively studied, the impact of the non-linear electrodynamics on the optical appearance of the EH BH shadow within the framework of different accretion flow models is an opening question. Meanwhile, it is unclear whether the accretion form affects the EH BH shadow optical appearance since the BH shadows exhibit various interesting observation characteristics under different accretion flow models. This paper focuses on this issue. The BH shadow is a powerful tool to investigate the observed characteristics of the BHs, it is interesting to in-vestigate the optical appearance of the EH BH and know the effect of the QED effect on it. Moreover, to test the effects of the space-time structure on the EH BH shadow, observation effects can also be investigated with a varying magnetic charge. We investigate the optical appearance of the EH BH within three different optically thin accretion flows context and analysis the EH BH observation characteristics under varying magnetic charges by taking the QED effect into account.
This paper is organized as follows. Section 2 briefly discuss the effective geometry of the EH BH and derives the light ray trajectory by using the ray-tracing method. In section 3, we present the shadows and rings as well as the corresponding optical appearance based on the three accretion flow models. We draw the conclusions in section 4.

EH BH the effective geometry and light deflection
The four-dimensional action of the GR coupled with the non-linear electrodynamics can be described as [45] where g is the metric tensor determinant, R is the Ricci scalar. L (F , G) is the non-linear electrodynamics Lagrangian, depending on electromagnetic invariants, where F = 1 4 F µν F µν , G = 1 4 F * µν F µν with F µν denoting the electromagnetic field strength tensor in which * F µν = ε µνσ ρ F σ ρ /(2 √ −g) its dual. The completely antisymmetric tensor is ε µνσ ρ , which satisfies ε µνσ ρ ε µνσ ρ = −4!. The Lagrangian density of the EH non-linear electrodynamics is given by [38] where a = 8α 2 45m 4 is the EH parameter, regulating the intensity of the non-linear electrodynamics contribution. α is the fine structure constant and m is the electron mass, hence the EH parameter is written as α/E 2 c . For a = 0, we recover the Maxwell electrodynamics, i.e. L = −F .
In this paper, we focus on the static and spherically symmetric metric of the 4D compact objects, the line element is where dΩ 2 = dθ 2 + sin 2 θ dφ 2 , f (r) is the EH BH metric potential and it can be written as [40], It was derived originally as an exact solution to the semiclassical Einstein equations with one-loop quantum corrections to QED. M is the BH mass, Q is the BH magnetic charge and the EH parameter a satisfies 0 ≤ a ≤ 32 7 Q 2 [46]. For the a → 0, the standard Maxwell electrodynamics is recovered, the EH BH will degenerate into the Reissner-Nordström (RN) BH. Note that, Bretón proposed a rotating BH solutions in the EH theory [47], we only consider the optical appearance of a static spherically symmetric EH BH solution in this analysis since the behavior of photons around the spin BH are rather differently from the static spherically symmetric BH.
In order to investigate the light deflection near the EH BH, the movement of light ray in this space-time needs to be understood. Note that, the effective geometry induced by the EH non-linear electrodynamics effects, which correct the background geometry along the null geodesics of which photons would usually propagate [48,49]. Taking this approach into account in our analysis, the photons propagate along the null geodesics of this effective geometry. The null geodesics of the photon paths are described by the effective space-time, that is [44,[48][49][50] where L F ≡ ∂ L ∂ F . The effective geometry seen by photons on the background of the EH BH can be rewritten as [44]. Therefore, we have and Note that the function H(r) and h(r) must be positive in order that the effective geometry does not change its signature during the photon motion. Hence, the Lagrangian of the EH BH space-time is whereẋ ≡ dx/dτ, the τ is the affine parameter. Since we only consider the photons that move on the equatorial plane (θ = π/2,θ = 0 andθ = 0), and the EH BH Lagrangian is independent explicitly on time t and azimuthal angle φ , hence, one can obtain the two conserved constants: The four-velocity of the time, the azimuthal angle, and the radial components can be obtained, where the symbol "±" denotes that the counterclockwise (−) and clockwise (+) direction for the motion of photons. The b is the impact parameter, satisfying b = |L|/E. The equations of motion for the null geodesic is Based on the method of Refs. [44,50] and using Eqs. (12) - (15), the effective potential V e of the EH BH can be written as Utilizing the effective potential function Eq. (16) and the condition of unstable circular orbit dV dr = V = 0 [44,50], we have According to above expression and Eqs. (4) (7) (8), the numerical results of the EH BH event horizon radii, shadow radii and critical impact parameters for different parameter values are listed in Tab.1. It is found that the increase of Q value leads to the decrease of r + , r ph and b ph , implying that the BH photon ring is shrunk inward the BH by increasing the magnetic charge. However, these key quantities are not sensitive to the EH parameter a. Based on Eq. (16), we can obtain that By introducing a parameter u ≡ 1/r, the expression of Π (u) ≡ du dφ is obtained. Utilizing the ray-tracing code, the trajectory of the light ray for different parameter values are shown in Fig. 1, showing that the radius of the black disk is smaller and the light rays are more gentle near the BH for a larger Q. The light density received by a distant observer increases with the increase of Q. Note that the deflection of the light ray is insensitive to the EH parameter. The magnetic charge as a free parameter, the BH shadow radius depends on it. Fig. 2 shows the BH shadow diameter Mpc, and its BH mass M = (6.5 ± 0.9) × 10 9 M . The diameter of its shadow is given as d M87 * ≈ 11.0 ± 1.5 [50,51]. Fig. 2 illustrates our result is consistent with that derived from the EHT observations within the observational uncertainty. Using the 1σ and 2σ confidence intervals of d M87 * , g can be constrained as Q 1.28 at 1σ and Q 1.64 at 2σ .
3 Optical appearance of the EH BH under the three accretion flow models

Static spherical accretion flow model
When the BH is surrounded by a static, optically thin, and geometrically thin spherical accretion flow, the observed intensity (ergs −1 cm −2 str −1 Hz −1 ) with a frequency υ s o is given by [52] where g s ≡ υ s o /υ s e is the red-shift factor, υ s e is the intrinsic photon frequency, dl prop is the infinitesimal proper length, and j(υ s e ) is the emissivity per unit volume in the rest frame of the emitter. For the EH BH, the red-shift factor can be regarded as g s ≡ H(r) 1/2 f (r) 1/2 . Note that, we take into account the effective geometry induced by the EH non-linear electrodynamics here. We consider a simple case of the emission is monochromatic with rest-frame frequency υ t . It emissivity has a radial profile as 1/r 2 , one can get   Under the EH BH context, the proper length measured in the rest frame of the emitter can be written as Utilizing Eqs. (20)- (22), the total observed intensity of the EH BH on the static spherical accretion background is obtained, i.e. Fig. 3 shows that the observed intensity I s (b) as a function of the impact parameter b for several representative parameter values. The observed intensity first ascended with the impact parameter, and reached the peak at b ph , which corresponds to the photon ring of the BH. For optically thin sources, the path length of a ray through the source is a proxy for the intensity picked up by that ray. Due to the light ray rotating around BH many times in the BH photon ring orbit, the optical path is infinite. Thus, a distant observer sees the strongest luminosity at the critical impact parameter. The peak value of the observed intensity increase with an increase of the BH magnetic charge when the EH parameter is a constant, and the corresponding b ph get smaller. The EH parameter does not affect the luminosity change, and the intensity curve is always constant whether a increases or decreases. Fig. 4 expounds that the two-dimensional shadows cast by this BH in the celestial coordinates. One can observe that a bright ring surrounds a central dark area. The shadow is not a totally dark region with zero intensity since the perspective of the shadows seen by a distant observer is occupied by the photons radiated from the bright accretion flow. The luminosities of the smaller magnetic charge BH shadows and photon rings are darker than that of the larger magnetic charge BH since the BH magnetic charge makes the BH photon ring can capture more photons from the accretion flow. Thus, the EH shadows and photon rings' luminosities increase gradually as the magnetic charge increases. Consistent with the results in Tab. 1, the change of a does not affect the size and luminosity of the shadow, which indicates that the luminosity of the EH BH is independent of the one-loop corrections to QED effect.

Infalling spherical accretion flow model
We consider the EH BH is surrounded by a radial infalling and optically/geometrically thin spherical accretion flow. In this case, the equation of the observed intensity on the static spherical accretion flow is still valid (Eq. 20). However, the redshift factor is different from the above situ-  ation, it can be written as [53] in which k µ is the four-velocity of the photons, u µi o is the four-velocity of an observer, and u µi e is the four-velocity of the accretion flow. According to Eqs. (12)- (14), one can obtain the k t is a constant, that is, k t ≡ 1/b. The k r comes from k β k β , we have where the symbol "±" corresponds the photon is approaching (+) or away (−) from the BH. The four-velocity of the accretion flow as u ri e = − 1 − H(r) f (r).
The redshift factor g i in this scenario can be written as The infinitesimal proper length is Hence, one can obtain that the total observed intensity of the EH BH under the infalling spherical accretion flow context: Fig. 5 shows that the total observed intensity function, the function curves shows similar feature to Fig. 3, but the observed intensity has an extremely sharp rise before the peak. Fig. 6 illustrates that the two-dimensional shadows cast on the celestial coordinates in this scenario. Our result  indicates that the size and position of the EH BH shadows do not change in both of these accretion flows.

Summary:
The accretion flow property is crucial to the optical appearance of the BH. For the static and infalling spherical emitters, the luminosities of the photon rings with different magnetic charges are listed in Tab. 2. One can see that the total observed intensity in the static spherical accretion flow scenario leads than that of the infalling spherical accretion flow under same parameters.
Compared with the static and infalling spherical emitters, we suggest that: i) the size and position of the BH shadow does not change in case of static and infalling spherical accretion flows, implying that the BH shadow is a signature of space-time geometry; ii) the BH shadow with infalling spherical accretion flow is darker than that of the static spherical accretion flow in the central region since the doppler-effect caused by the radial infalling accretion flow; iii) the influence of the magnetic charge Q on luminosity is  basically the same under the two spherical accretion flows, i.e. the peak value of intensity increase with an increase of the BH magnetic charge when the EH parameter is a constant; iv) QED effect affects the optical appearance of BH. Its existence makes the optical morphology of the EH BHs completely different from the RN BHs.
Additionally, we blur the two-dimensional image and correspond roughly to the EHT resolution. This simple blur does not correspond to the EHT image reconstruction and can only offer a rough illustration of the EHT resolution. From Fig. 7, one can observe that although the observed shadow luminosity is different, the size of the BH shadow does not change for both static and infalling spherical accretion flow. These results suggest that the BH shadow size depends on the geometric space-time and the luminosity of the BH shadow relies on the accretion flow models. Compared with the results of the EHT, the shadow of the EH BH presents a complete circle, and the photon rings are evenly distributed around the shadow. According to Refs. [2][3][4][5][6][7][8][9], the M87 * image shows that a supermassive rotating BH space is illuminated by a magnetically arrested accretion disk. This rotation results in the accumulation of brightness in the southwest of the EHT image. In the follow-up works, we will further describe the optical appearance of the rotating BH solutions within the framework of the EH theory. Inspired by these results, we investigate the optical appearance of the EH BH surrounded by an optically and geometrically thin disk accretion flow in the next section.

Thin disk accretion flow model
Assuming that the BH is illuminated by an accretion disk in the equatorial plane, the disk emits isotropically in the rest frame of static worldlines, and the observer is at the north pole. According to the definition of the total number of light orbits, the light trajectories emitted from the north pole direction can be subdivided into direct emission, lensing ring, and photon ring [20]. When the light trajectories intersect the accretion disk just once and finally fall on the front, the corresponding light ray forms the direct emission. When the light trajectories intersect the accretion disk twice (the light breaks through thin disk) and finally fall on the back of the disk, which is the lensing ring. In this path, the light picks up additional brightness from the second intersection between the light trajectories and the accretion disk. The light arrives at the front side of the thin accretion disk once again for the photon ring, leading to additional brightness from the three intersections between the light and the disk. Hence, the total observed intensity should be the sum of those intensities.
Based on the Liouville's theorem, I em /(υ em ) 3 is conserved in the direction of light propagation, where I em and υ em delegate the radiation intensity and frequency of the accretion disk, respectively. An observer in infinity receive the specific intensity I d o with red-shifted frequency υ d o ≡ √ f υ em . Based on these results, we have For the EH BH under the thin disk accretion flow context, the observed intensity for a specific frequency can be given by The total intensity is an integral over all frequencies: where I e (r) ≡ I em (r)dυ em is the total radiation intensity of the thin disk accretion flow. Hence, the total observed intensity of the EH BH in this situation can be written as where r n (b) is the so-called transfer function which represents the radial coordinate of the n th intersection between the light with impact parameter b and the accretion disk. The slope of the r n (b) -dr/db -is defined as the (de)magnification factor. Fig. 8 shows that r n (b) as a function of b for different parameter values. For the case of n = 1, the r n (b) is a linear function with a slope of 1, indicating that r n is proportional to b. This corresponds to the "direct emission" scenario. The case of n = 2 is for the "lensing ring". The r n (b) function illustrates as a asymptotic curve, and b is limited in a very narrow range around b 4.73M −5.89M (Q = 0.5, a = 0.5). Therefore, the lensing ring shows up as a thin ring in the shadow image. The case of n = 3 is for the "photon ring". In this scenario, r n (b) is almost a vertical line, suggesting that the photon ring should be an extremely thin ring at a given parameter value. It is also found that the increase of the BH magnetic charge leads to a decrease of the impact parameter, the EH parameter do not influence it. Fig. 9 illustrates that the light trajectories of different rings in the polar coordinate by utilizing the ray-tracing procedure. One can see that the radius of the black disk is smaller and the light rays are more curved for the large BH magnetic charge. The radii of the direct emission, lensing ring, and photon ring are shrunk if the Q value increases. Note that the thickness of the direct emission, lensing ring, and photon ring is not sensitive to the QED parameter.
Next, we investigate the optical appearance of the EH BH in this scenario. We main purpose is to analyze the influence of the accretion disk radiation position on the observation characteristics of the BH. As well known that the radiation of accretion disk in the universe satisfies Gaussian distribution [54], hence, we parameterize the radiations intensity of the accretion disk as a gaussian function, that is where r in is the innermost radiation position of the accretion disk. Our analysis is for three different scenarios: (A) r in = r isco , where r isco is the radius of the BH innermost stable circular orbit; (B) r in = r ph at which r ph is the radius of the BH photon ring; (C) r in = r + , where r + is the radius of the BH event horizon.

Case A: r in = r isco
The innermost stable circular orbit (r isco ) is one of the relativistic effects, representing the bounder between test particles circling the BH and test particles falling into the BH. It is given by [28] When the Q = 0.5 and a = 0.5, the innermost stable circular orbit radius of the EH BH is r isco = 5.61r g . According to Eqs. (34) and (35), the total radiation intensity I A e as a function of the radius, the total observed intensity I A o as a function of the impact parameter, and the two-dimensional image in celestial coordinates are displayed in the left panels of the Fig. 10. It is found that the regions of the direct emission, lensing ring, and photon ring are separated. The direct emission starts at b 6.52M and peaks at b 8.18M. Its maximum intensity is 0.485. The lensing ring is limited in a small range of b 5.24M ∼ 5.65M. The photon ring appears at b 4.98M. In the two-dimensional image, the border of the black disk corresponds to r isco . A bright lensing ring is shown within the black disk, and the dim photon ring is in the inner of the lensing ring.  Fig. 10 The total radiation intensity as a function of radius, the total observed intensity as a function of the impact parameter, and the twodimensional images of the shadows for the EH BH. The Q = 0.5, a = 0.5, and the BH mass as M = 1. bright, extremely narrow ring. An completely dark region as shown in the BH photon ring.

Case C: r in = r +
When the innermost of the accretion disk radiation position is in the radius of the BH event horizon (r + ), the radiation peak at the BH event horizon radius r + 1.86r g . Our results are shown in the right panels of Fig. 10. This scenario is similar to Case B. However, the black area that can be observed is reduced to the BH event horizon. The direct emission starts at b 2.76M and the maximum intensity is 0.22. The lensing ring is constrained to b 4.78M ∼ 5.62M, and the photon ring generated at b 5.24M. It is also found that the positions of the lensing ring and photon ring is almost unchanged by compared with Case B. Meanwhile, the region of the photon ring of the BH becomes slightly wider since the attenuation of the accretion disk radiation in the third case is slightly slow.

Summary:
This section shows that the optical appearance of the EH BH shadow depends on the innermost radiation position of the accretion disk, which applies to the model of the BH surrounded by a thin accretion disk. In this scenario, the optical appearance of the BH shadows do not only depend on the structure of space-time, which is different from the static and infalling spherical accretion flow models. The contribution of the different rings to the total observed flux with different parameters are listed in Tab. 3. It is found that the contribution of the lensing ring to the total observed flux is less than 5%, and the photon ring is less than 2%, indicating that the direct emission dominates the optical appearance of the BH image. We also blur the two-dimensional images as shown in Fig. 11. One can observe that the blurring washes out the lensing ring and photon ring features. Their observational appearances reply on instrument resolution. It is difficult to obtain the ring information with the current resolution of EHT. Thus, we suggest that the optical appearance of the BH image depend on the accretion disk radiation position in this scenario.

Conclusions and Discussion
The optical appearance of the EH BH in the different accretion flows framework have been revealed in this analysis. Taking the QED effect into account, we obtained that the effective geometry induced by the EH non-linear electrodynamics effects. The photons propagate along the null geodesics for this effective geometry is discussed. According to derived the effective potential function, we found that the increase of the magnetic charge leads to the decrease of the EH BH event horizon radius, shadow radius and critical impact parameter, implying that the BH photon ring is shrunk inward the BH by increasing the magnetic charge. By investigating trajectory of the light ray, we found the radius of the black disk is smaller and the light rays are more gentle near the BH for a larger magnetic charge and the light density increases with the increase of magnetic charge.
For the EH BH is illuminated by the static and infalling spherical accretion flows, the total observed intensity function is calculated. We found that the observed intensity first ascended with the impact parameter, and reached the peak at the BH photon ring. The peak value of intensity increase with an increase of the BH magnetic charge, and the corresponding b ph get smaller. The change of a does not affect the size and luminosity of the shadow, which indicates that the luminosity of the EH BH is independent of the one-loop corrections to QED effect. We also found the total observed intensity in the static spherical accretion flow scenario leads than that of the infalling spherical accretion flow under same parameters. The size and position of the EH BH shadows do not change in both of these accretion flows, implying that the BH shadow size depends on the geometric space-time and the shadows luminosities relies on the accretion flow models.
For the EH BH is surrounded by a thin disk accretion flow, we found that the radii of the direct emission, lensing ring, and photon ring are dramatically decreased as the magnetic charge increases. In case of that the innermost of the accretion disk radiation position equals to r isco , the regions of the direct emission, lensing ring, and photon ring are separated. A bright lensing ring is shown within the r isco , and the dim photon ring is in the inner of the lensing ring. In case of the innermost of the accretion disk radiation position equals to r ph , the two rings are so close that it cannot Fig. 11 Blurred two-dimensional images utilizing a Gaussian filter with a standard deviation of 1/12 the field of view for the thin disk accretion flow in three cases. be almost separated. In case of the innermost of the accretion disk radiation position shrinks to the event horizon of the BH, the observable black area is shrunk to the BH event horizon, but the radii of the lensing ring and photon ring is almost unchanged. We also found that the contribution of the lensing ring to the total observed flux is less than 5%, and the photon ring is less than 2%, indicating that the direct emission dominates the optical appearance of the BH image. Blurring the two-dimensional image in this scenario, we observed the blurring washes out the lensing ring and photon ring features, which replies on instrument resolution. It is difficult to obtain the ring information with the current resolution of EHT. Thus, we believe that the optical appearance of the BH image depend on the accretion disk radiation position.
Based on our analysis, we argue that the optical appearance of the EH BH depends on the accretion flows morphology and BH space-time structure. When the spherical accretion flows illuminated the BH, the BH shadow is shown as a geometric feature of space-time, the size of the BH shadow does not change with the position of the spherical accretion flows. When the thin disk accretion flow illuminated the BH, the observable characteristics of the BH shadow rely to the position of the radiating accretion disk with respect to the BH.