Unveiling the unconventional optical signatures of regular black holes within accretion disk

The optical characteristics of three types of black holes (BHs) surrounded by a thin accretion disk are discussed, namely the Schwarzschild BH, Bardeen BH, and Hayward BH. We calculate the deflection angle of light as it traverses the vicinity of each BH using numerical integration and semi-analytical methods, revealing that both approaches can effectively elucidate the deflection of light around the BH. We investigate the optical appearance of the accretion disk and its corresponding observational images at various viewing angles, discovering that the luminosity in the region near the BH on the inner side of the accretion disk is higher than that on the outer side owing to higher material density in closer proximity to the BH. We observe a significant accumulation of brightness on the left side of the accretion disk, attributed to the motion of matter and geometric effects. Our findings emphasize the significant influence of the observation inclination angle on the observed outcomes. An increase in the observation inclination angle results in the separation of higher-order images. With the improvement in EHT observation accuracy, we believe that the feature of a minimal distance between the innermost region of the direct image of the Hayward BH and the outermost region of the secondary image can be used as an indicator for identifying Hayward BHs.


Introduction
Black hole (BH), often regarded as the most mysterious object in the universe, are theoretically predicted to arise from the general relativity (GR).Despite the formidable challenge of directly observing BHs, a substantial body of astronomical evidence supports their presence in the cosmos.The Laser Interferometer Gravitational-Wave Observatory (LIGO) has successfully detected gravitational waves generated by the merger of BHs, not only corroborating Einstein's GR but also furnishing a novel tool for BH research [1].Subsequently, numerous other gravitational wave events have been detected, yielding a wealth of data to further substantiate the existence of BHs [2][3][4][5].The most compelling evidence for BH existence stems from direct observations of properties associated with the event horizon.Utilizing the capabilities of the Event Horizon Telescope (EHT), the image of the supermassive BH in Messier (M) 87 * elliptical galaxy is obtained, which provided the first direct observational evidence for the existence of BHs [6].Recently, the EHT collaboration released an image of the BH at the core of our galaxy, Sagittarius (Sgr) A * [7].These achievements not only demonstrate the spacetime structure of the central BH but also encode dynamic information about the distribution of surrounding matter [8].
An observable astrophysical BH is typically encircled by a luminous accretion disk, which forms a ring-shaped region where material accumulates around the BH and is considered one of the primary energy sources for the BH.Through numerical simulation, the magnetic field and jet characteristics on the BH accretion disk were investigated, revealing the influence of angular momentum on the accretion-ejection theory [9].To elucidate the formation of seed BHs in the early universe, a super-exponential accretion disk evolution model has been proposed.It has been demonstrated that a dynamic mechanism can trigger hyper-exponential accretion when a BH seed is trapped in a star cluster that is fed by ubiquitous dense, cold gas streams [10].By establishing analytic force-free solutions of rotating stars and BHs within the magnetic field, the influence of BH accretion disk and spin on the extraction of electromagnetic energy has been investigated [11].It is noteworthy that a substantial research team employs numerical methods to study the accretion flow of BHs.In particular, Porth et al. provided a list of some of the most advanced codes still in use for this purpose as of 2019 [12].Furthermore, Noble et al. obtained numerical results when the accumulated heat energy is fully radiated.Their findings revealed that the total luminosity of the accretion flow could be approximately ≃ 20% greater than the value predicted by the Novikov-Thorne model [13].
In addition to numerical simulations investigating the characteristics of BH accretion disks, there have been notable studies focusing on imaging these disks.For instance, Luminet employed a semi-analytical method to examine the image of a standard accretion disk within the framework of a Schwarzschild BH and obtained direct and secondary images of the disk [14].Bambi utilized ray-tracing code to fit X-ray data from current and future BH candidates, effectively constraining possible deviations from Kerr geometry in the spin-deformation parameter plane [15].Grenzebach et al. derived an analytical formula for the shadow in the Kerr-Newman-NUT-(anti-)de Sitter BH spacetime, allowing them to visualize the photon sphere regions and shadow radii under different parameter values [16].The EHT Collaboration has successfully conducted rigorous tests on BH image general relativistic radiative transfer codes, resulting in a minor 1% flux error, primarily attributed to software camera setup discrepancies.
These findings provide assurance that numerical uncertainties are unlikely to hinder parameter estimation in current EHT observations [17].Prather et al. assessed various general relativistic radiative transfer codes used to predict polarized emissions from BH accretion systems and found them to be highly consistent.Whether employing an analytic accretion model or a magnetohydrodynamic simulation, these codes consistently produced images with minimal deviations, ensuring reliable results for interpreting polarized BH images [18].Other intriguing studies on imaging BH accretion disks can be found in references [19][20][21][22][23][24][25].Gralla et al. proposed a refined description of the BH shadow, lensing ring, and photon ring within the framework of optically thin emission from a geometrically thin or thick disk.They observed that the brightness of the photon rings exhibits a logarithmic divergence at the shadow radius and that the observable characteristics of the BH shadow are influenced by the details of the accretion disk [26].Subsequently, researchers have investigated image formation in BH accretion disks within various gravity spacetime backgrounds [27][28][29][30][31][32][33][34][35][36][37][38][39].
These studies predominantly center on optically thin emission from geometrically thin or thick accretion disk, despite the fact that most observed accretion disks in astrophysics exhibit optically thick.In one particular study, an optically thick and geometrically thin accretion disk surrounding a BH is derived under the influence of Einstein-Gauss-Bonnet gravity.The results reveal that the image of the Gauss-Bonnet BH bears similarity to that of the Schwarz-schild BH, while the naked singularity exhibits distinct characteristics [40].The accretion disk model employed in this context adheres to the Novikov-Thorne model, which is a standard model characterized by geometric thinness and optical thickness.By examining the observational characteristics of strongly naked static Janis-Newman-Winicour singularity in the Novikov-Thorne model context, Gyulchev et al. provided insights into properties of the resulting multiring structure and elucidated the reasons for the appearance of this ring structure by unveiling its underlying physical formation mechanism [41].Furthermore, Liu et al. discussed the physical properties and optical appearance of the electromagnetic radiation emitted from a thin accretion disk around a charged and slowly rotating BH within the context of the Einstein-aether theory [42].Accretion disk images have also been explored in other background spacetimes, including γ-metric BH [43], wormholes [44], Schwarzschild BH penetrated by cosmic string [45], string cloud Schwarzschild BH [46], Schwarzschild-MOG BH in scalar-tensor-vector gravity [47].
In addition to studying the characteristics of BH images in singular spacetimes, exploring BH images in regular spacetimes represents an intriguing area of research.Different from irregular BHs, which feature inherent singularities at the origin of spacetime, Bardeen introduced a BH solution unaffected by spacetime singularities [48].Building upon Bardeen's idea, Hayward proposed a regular Hayward BH solution by combining nonlinear electrodynamics with Einstein's field equation [49].In our prior research, we examined the observable attributes of the shadow and photon rings of the Hayward BH, uncovering the influence of accretion morphology and the BH's magnetic charge on its optical appearance [30].Another study focused on a thin-shell wormhole with a Hayward profile [50].Additionally, we examined the shadow images and observed luminosity of the Bardeen BH under various accretion scenarios [51].Regrettably, our work only consider regular BHs surrounded by an optically thin emission from geometrically thin accretion disk.Our findings indicated that the disappearance of singularities does not significantly impact the optical appearance of the regular BH, and the lensed ring and photon ring do not reflect spacetime singularity.
In this analysis, our objective is to investigate the appearance of regular BH within the Novikov-Thorne accretion disk context.This investigation aims to provide a comprehensive understanding of the optical morphology of regular BHs, enabling us to distinguish them from singular BHs by investigating their observable characteristics.
To achieve this, we employed the semi-analytical approach initially introduced by Luminet [14], complemented by the ray tracing technique as outlined by Falcke and collaborators [52].Subsequently, we utilized these methodologies to develop our own Mathematica code for the Novikov-Thorne accretion disk model.This enabled us to delve into the investigation of both the primary and secondary images originating from three distinct categories of BHs enveloped by an optically thick accretion disk.We analyze the optical morphology of these BHs, taking into account redshift and flux distribution, and provide a comparative analysis.The structure of this paper is as follows: In Section 2, we briefly review three types of BHs and discuss their effective potential.In Section 3, we derive the direct and secondary images of the target BHs using numerical integration and semi-analytical methods.In Section 4, we present the optical appearance of the accretion disk as well as the corresponding observational images under different observation angles.Finally, we draw conclusions in Section 5.

Effective potential and ray-tracing
The static spherically symmetric BH metric can be described as follows: where A(r) is BH metric potential, B(r) is defined as the reciprocal of A(r), and C(r) is defined as r 2 .In this article, we consider Schwarzschild BH, Bardeen BH, and Hayward BH, whose metric descriptions are [48,49] Schwarzschild BH : Hayward BH : Bardeen BH : in which M represents the mass of BH, while g denotes the BH magnetic charge.As the magnetic charge approaches zero, both the Hayward and Bardeen BHs transition into a Schwarzschild BH.Consequently, the magnetic charge in the regular BH does not alter the spacetime causal structure or the Penrose diagram.The discussion of the causal structure of regular BHs has been elaborated upon in several sources [53][54][55][56].However, since the causal structure of BHs does not provide insights into their observable features, we will not delve into Penrose diagrams in this context.Instead, our attention is directed toward examining the optical characteristics of BHs.
To study the motion trajectory of photons around a BH, it is necessary to analyze the evolution of particles in the BH spacetime.Our analysis begins with the Lagrangian of these particles where λ denotes an affine parameter and ẋµ represents the four-velocity of the particle.The Lagrangian L governing the motion of a particle depends on its mass.For massive particles, the Lagrangian is negative, whereas for massless particles, such as photons, it vanishes, i.e., L = 0.The generalized momentum p µ associated with the particle can be obtained from the Lagrangian by applying the Euler-Lagrange equations of motion, which leads to four equations of motion for a particle with energy E and angular momentum L, we have We only consider the particle that move on the equatorial plane (θ 0 = π/2, θ0 = 0 and θ = 0).In addition, we assume that the test particle is a photon with zero angular momentum (L = 0).Using the four equations mentioned above, we can obtain the following results: where V eff is the effective potential.For the static spherically symmetric BH (g tϕ = 0), one can get in which b is the impact parameter, defining as b ≡ L/E = r 2 φ A(r) ṫ .The effective potential can be further expressed as By utilizing Eq. ( 15), we can determine both the radius of the photon sphere and the critical impact parameter by setting and V ′ eff = 0. Fig. 1 illustrates the effective potential of a BH as a function of radius, showcasing three different BH solutions.It is noticeable that the peak of the effective potential for a Schwarzschild BH occurs at approximately ≃ 3r g , where r g signifies the Schwarzschild radius.For the same set of parameters, the effective potential for a Bardeen BH is greater, and the radius corresponding to the peak of the effective potential is smaller in comparison to that of a Hayward BH.To better illustrate the distinctions among the three types of BHs, we derived the deflection angles for various scenarios.Utilizing Eqs. ( 8) and (9), we obtained the following results: By introducing a parameter u ≡ 1/r, the above equation is written as For a distant observer, the deflection angle of light is where the upper and lower limits of the integral represent the photon's emission point u source and the observer's position u obs , respectively.It's essential to note that in some instances, rays with an impact parameter exceeding the critical value may possess a radial inflection point u o .In such scenarios, the motion of photons can be described by splitting the integral into two separate integrals.
in which u o is photon turn point, satisfying Hence, Eq. ( 19) effectively guarantees the tracking of all light near a BH.Apart from employing numerical integration algorithms to investigate light deflection angles, it's also valuable to delve into the semi-analytical algorithm introduced by Luminet [14].In this article, we also employ a semi-analytical approach to derive the deflection angle of light, with the goal of comparing it with numerical methods.Using Eq. ( 17), one can get Schwarzschild BH : Hayward BH : Bardeen BH : where the lower indicators u s , u h , and u b represent the Schwarzschild, Hayward, and Bardeen BHs, respectively.Using the Cardano formula, we can rewrite Eqs. ( 21)- (23) as follows: Hayward BH : Bardeen BH : where the cubic polynomial G(u) contains two positive roots and one negative root, denoted by u 1 ≤ 0 < u 2 < u 3 .The specific forms of these roots can be determined using the method proposed by Luminet [14] and expressed in terms of the periastron distance P as follows in which the Q 2 ≡ (P − 2M )(P + 6M ).The impact parameter in three different BH scenarios can be written as Schwarzschild BH : Hayward BH : Bardeen BH : Based on the equations provided above, it is possible for an observer located at infinity can capture the light that arrives here at a given periastron distance P .As a result, the bending angle of the light ray can be calculated, which is By introducing the elliptic integral, Eq. ( 31) can be transformed into where F (Ψ i , k) represents an incomplete elliptic integral of the first kind, with amplitudes as parameter variables.The parameter variables are Therefore, the total change in bending angle can be simplified as follows in which K(k) is the complete elliptic integrals of the first kind.We computed the deflection angles of light in three types of BHs using two different methods, both of which, despite their distinct forms, all yield good results.Fig. 2 shows the trajectory of light near a BH using the ray-tracing code, where the yellow line represents b < b c , the red line represents b > b c , and the blue line represents b = b c .It's important to note that in this scenario, photons rotate multiple times around the BH.We can observe that, in comparison to the Schwarzschild BH, the black disk radius of the Hayward BH and Bardeen BH are slightly smaller, and the light density received by a distant observer decreases.

Direct and secondary images of BH
We consider a BH enveloped by an optically thick and geometrically thin accretion disk, as described by the Novikov-Thorne thin disk model.This model originates from research conducted in the previous century [57,58].The details of this model are not discussed here (see Refs. [14,57,58] for details).Our primary focus is on contrasting the images of three distinct types of BHs within the framework of this model.The radiation originates from a point with coordinates (r, ϕ) on the emitting plane M and propagates to a point with coordinates (b, α) on the observation plane m.Two images of the accretion disk can be obtained on the observation plane: the direct image at coordinates (b (d) , α) and the secondary image at coordinates (b (s) , α + π) [14].In this section, we derive the direct and secondary images of BHs surrounded by an accretion disk using a numerical integration algorithm and semi-analytical methods, respectively.
Following the references [40][41][42][43], we assume that there is a relationship between the azimuth angle ϕ of the photon emission point, the celestial angle η, and the observer inclination angle θ 0 , which is where the celestial angle η ranges from 0 to 2π.Using these two trigonometric functions and and Eq. ( 19), we can determine the impact parameters of all photon trajectories that are emitted from a specific circular orbit with radial coordinates r = r source .It should be noted that these trajectories can reach the observer whose observation inclination angle is θ 0 and position is r = r obs .Therefore, we obtain the following expression: This formula only applies to a specific ray of light and cannot account for all possible deflections that may occur.It is possible that some of the light emitted from the disk could be deflected at a greater angle and could even orbit the BH multiple times before reaching a distant observer.To account for these more complex scenarios, the above equation must be extended to a more general situation where k is a positive integer that can represent the image order.k = 0 describes the direct image of the accretion disk, while k = 1, 2, 3... representing second order, third order, and higher order situations.The same effect can also be demonstrated using Luminet's semi-analytic method.In this approach, we assume that the deflection angle from the source to the observer is denoted by γ, and the observer's inclination angle is denoted by θ 0 .This derivation also involves the use of elliptic integrals.By applying trigonometric function relationships, we obtain The Eq. ( 21) can be rewritten as where F (ζ r , k) and F (ζ ∞ ) represent the elliptical integrals, the parameter variables , and length k 2 = Q−P +6M 2Q [14].Meanwhile, the radius r can be written as a function of α and P , which is where sn is Jacobi elliptic function.Depends on this equation, one can present the iso-radial curves for a specific angle θ 0 , which can be used to obtain a direct image of the accretion disk.To obtain the (n + 1)th order image of the accretion disk, the Eq. ( 39) satisfies in which K(k) is the complete elliptic integral, n is a positive integer, which can represent the image order.Figure 3 present the direct (k = 0) and secondary (k = 1) images of circular rings in orbit around BHs.The image displays three different types of BHs: Schwarzschild BH, Hayward BH, and Bardeen BH.The ray trajectories are positioned at r = 6M , r = 10M , r = 15M , and r = 20M (from innermost to outermost orbits), with inclination angles of θ 0 = 17 • , 53 • , and 75 • .As the observation angle decreases, the secondary image becomes embedded inside the direct image, forming a structure similar to a "photon ring".This can be considered a viable alternative for a photon ring.As the inclination angle of the observation increases, the degree of separation between the direct and secondary images increases.Upon comparing the results of the three types of BHs, it can be observed that the orbital size of the Hayward BH is smaller than that of the Bardeen BH under the same parameters.Notably, the difference between the Bardeen BH and Schwarzschild BH is not significant.

Optical appearance of three types of BHs
In this section, we examine the observable characteristics of three types of BHs and explore how observation inclination angles affect them.Our goal is to distinguish these characteristics from their optical appearance.Assuming an optically thick and geometrically thin accretion disk surrounding the BH, we can calculate the electromagnetic radiation flux (ergs −1 cm −2 str −1 Hz −1 ) emitted from a specific radial position r on the accretion disk using the following formula [58]: where Ṁ refers to the mass accretion rate, I is the determinant of the induced metric in the equatorial plane, and r in denotes the inner edge of the accretion disk.The parameters E, Ω, and L represent the angular velocity, energy and angular momentum of the particles in a circular orbit, respectively, and they can be expressed as follows: It is worth noting that the Novikov-Thorne model is specifically designed for Kerr metrics.However, if the metric is changed to that of a regular BH, the structure of the thin disk should also be adjusted.Jaroszyński et al. calculated the radiation spectrum observed by observers at various positions relative to the equatorial plane of the disk [59].They investigated the emission line profiles from self-gravitating toroids around BHs [60].Some literature on Kerr metric accretion disk images and ray-tracing can be found in references [61][62][63].The pioneering work on the four-dimensional spherical symmetry metric originated from reference [14].It is important to note that the spherical symmetry metric component g tϕ = 0 leads to changes in angular velocity, energy, and angular momentum forms.A discussion on spherically symmetric ray-tracing and accretion disk images can be found in references [64][65][66][67][68][69][70][71].In the case of spherical and axisymmetric scenarios, the received light intensity is inconsistent for a distant observer.
Bardeen BH : where 4 illustrates the radiant intensity of the accretion disk in three different scenarios involving BHs.As the observation tilt angle increases, the direct and secondary images gradually separate.Notably, the luminosity in the region near the BH on the inner side of the accretion disk is higher than that on the outer side due to higher material density in closer proximity to the BH.This proximity makes it more susceptible to the BH's gravity, resulting in more intense heating and radiation.Furthermore, the appearance of the three types of BHs can be distinguished from one another.In comparison to the Schwarzschild BH image, the Hayward BH exhibits a feature of minimal distance between the innermost region of the direct image and the outermost region of the secondary image.The appearance of the Bardeen BH closely resembles that of the Schwarzschild BH, but there is a significant difference in radiance.Consequently, our results suggest that observing such an effect can be beneficial when considering an external metric consistent with the regular BH solution.Figure 5 shows that the total observed intensity function, the observed intensity has an extremely sharp rise before the peak.Fig. 4 illustrates that the two-dimensional shadows cast on the celestial coordinates.Our result indicates that the size and position of the BH shadows remain unchanged in different accretion flows, implying that shadows are the spatiotemporal geometric characteristics of BHs.The primary distinction between regular BHs and Schwarzschild BHs is the presence of magnetic charges.To better illustrate the correlation between radiant flux and magnetic charge, we have employed three functional models: the Fourier Function (Case 1), the Gaussian Function (Case 2), and the Exponential Function (Case 3).Section 4.4 provides detailed information on the form and accuracy of these functions.Our objective is to represent the accurate value of F as a function of g in a more analytical manner.Fig. 5 displays the red dot, which represents the precise value of F as a function of g, indicating an increase in F .With the Hayward BH, the slope of the fitting curve grows with an increase in g, whereas the fitting curve of the Bardeen BH exhibits a linear relationship with the magnetic charge.Additionally, we have observed that, under specific fixed parameters, the exponential function provides the best fit to the numerical results.

Redshift factor
To obtain the observable radiant flux for a remote observer, it is necessary to consider the correction for the redshift effect.It is important to note that the redshift effect in this context encompasses both gravitational redshift and the Doppler effect.This correction involves taking into account the redshift factor, denoted as d, which is expressed as follows: where E em and E obs respectively represent the projection of photon four-momentum k µ on the four-velocities of source p µ source and observer p µ obs at a distance.The E em is in which p t and p ϕ represent the photon 4-momentum.For observers at a large distance, the ratio p t /p ϕ represents the impact parameter of the photons relative to the z-axis.This can be related to the trigonometric functions as follows sin θ 0 cos α = cos γ sin β, the p t /p ϕ can be rewritten as Therefore, the redshift factor d is obtained, one can get For the three types of spherical symmetry metrics we are considering, their redshift factors are written as Schwarzschild BH : Hayward BH : Bardeen BH : We assume that the inner and outer boundaries of the accretion disk consist of stable circular orbits at r = 6M and r = 20M , respectively.As illustrated in Fig. 6, we have observed a significant asymmetry between the redshift area (z > 0) and blueshift area (z < 0).Furthermore, the redshift area gradually increases as the observer's inclination angle increases.In comparison to Schwarzschild BH, Hayward BHs exhibit a smaller redshift region, whereas Bardeen BHs display a larger redshift region.
Unveiling the unconventional optical signatures of regular black holes within accretion disk16

Observation images
In this section, we will derive observational images of an accretion disk under three different BH spacetime backgrounds.Building on the discussions in the previous subsections, we can determine the observed flux of the accretion disk Figure 7 presents the results obtained for the three different BH scenarios.It is evident that the inclination angle of observation has a significant impact on the final result.As the inclination angle increases, the observable image of the accretion disk takes on a hat-like shape.It is worth noting that according to the Novikov-Thorne model and Eq. ( 42), we obtain the bolometric luminosity.Additionally, it is noteworthy that, regardless of whether the dip angle increases or decreases, there is a prominent accumulation of brightness on the left side of the accretion disk.This phenomenon may have two contributing factors.Firstly, as the light passes through the material surrounding the BH, it gets absorbed and re-emitted, resulting in the observed brightness distribution.The velocity and density distribution of matter as it moves towards the BH along the accretion disk can influence the probability of light being absorbed and re-emitted, leading to the observed brightness accumulation.On the other hand, the light may pass through the top and bottom of the disk, but due to the line-of-sight effect, we are more likely to observe the light emitted from the bottom of the BH, as it is closer to our line of sight.Therefore, the brightness accumulation often appears on the left side, which is a geometric effect dependent on the path and angle of light emission.
Furthermore, we can observe that the Bardeen BH more closely resembles the Schwarzschild BH, whereas the Hayward BH exhibits a feature of minimal distance between the innermost region of the direct image and the outermost region of the secondary image, which can serve as a significant characteristic for distinguishing Hayward BHs from other types of BHs.Notably, the accretion disk surrounding the Bardeen BH is brighter than that around the Schwarzschild BH.This is due to the fact that the accretion disk around the Bardeen BH experiences stronger gravitational forces, leading to a closer accumulation of matter and the release of more energy.The BH mass is taken as M = 3M ⊙ and magnetic charge is g = 1.

Data analysis
In order to better illustrate the observational characteristics of regular BHs, we aim to establish a relationship between magnetic charges and observed flux.To achieve this, we propose three hypothetical functional models and fit them to our data Fourier function : Gaussian function : Exponential function : where a − g with subscripts are the coefficient in the fitting function.
Tables. 1 -2 present the numerical relationship between the radiation flux F R , the observed flux F obs , and the magnetic charge g at various inclination angle θ 0 .In Fig. 8 shows the results obtained from fitting the data using three hypothetical functions: the Fourier function represented by the black solid line, the Gaussian function represented by the blue solid line, and the exponential function represented by the green solid line.It is evident from the figure that Bardeen's observed flux F obs exhibits a more linear relationship with magnetic charge, showing an increasing trend.Our current analysis relies on intuitive methods.However, for more accurate data fitting, we should employ the Mean Squared Error (MSE).This metric quantifies the difference between a model's predicted values and the actual values.A smaller MSE indicates a higher level of similarity between predicted and actual values, leading to a more accurate model.Conversely, a larger MSE indicates greater disparity and a less accurate model.To determine the best-fit function for our data, we computed the MSE for three hypothesis functions and compared them to the standard results.Tables.3 -4 present the MSE values for different fitting functions.Notably, the MSE for the exponential function is significantly smaller than that of the Fourier and Gaussian functions, by approximately four orders of magnitude.This suggests that the exponential function provides the closest fit to the theoretical results.Consequently, we report the best fit achieved with the exponential function in Tables.5 to 6.It is well known that the magnetic field of a BH accretion disk arises from the plasma within it.The size of the magnetic field depends on the magnetic charge density and the motion speed of the plasma.In this study, we establish a functional relationship between a BH's magnetic charge and the observed flux, which may guide further exploration of the effects of a BH's magnetic charge on optical appearance.

Magnetic charge constraint
Note that the regular BH shadow radius depends on the magnetic charge parameter g.However, from an observational perspective, the diameter of the BH shadow, denoted as d sh , can be measured using data from the EHT.It is well-known that the angular size of the shadow of M87 * is δ = (42±3) µas, and its distance is D = 16.8 +0.8  −0.7 Mpc.The BH mass is estimated to be M = (6.5±0.9)× 10 9 M ⊙ .Using these observational data, we can calculate the diameter of the shadow as d M87 * = Dδ M ≃ 11.0 ± 1.5 [72,73].The EHT has not only measured the angular size of the emission ring of the Sgr A * BH, which is δ d = (51.8±2.3)µas, but also estimated the angular size of its shadow to be δ = (48.7±7)µas.Additionally, the distance and mass of Sgr A * are separately given as D = (8.15±0.15)kpc and M = (4.0+1.1 −0.6 ) × 10 6 M ⊙ [7].Using this observational data, we can calculate the diameter of the shadow, denoted as d SgrA * .By applying the constraints obtained from EHT observations, we can further investigate the magnetic charge parameter g.
As depicted in Fig. 9, the result obtained for Hayward and Bardeen BHs are consistent with those derived from EHT observations within the observational uncertainty.Using the 1σ and 2σ confidence intervals of d M87 * , the magnetic charge parameter of the Hayward BH can be constrained as g ≤ 1.72 at the 1σ level and g ≤ 2.11 at the 2σ level.The magnetic charge parameter of the Bardeen BH can be constrained as g ≤ 0.92 within 1σ and g ≤ 1.33 within 2σ.When applying the confidence intervals of d sgrA * , the magnetic charge parameter of the Hayward BH can be constrained as g ≤ 1.89 within 1σ and g ≤ 2.02 within 2σ.For the Bardeen BH, the corresponding constraints are g ≤ 1.09 at 1σ and g ≤ 1.24 at 2σ.It is important to note that the accretion flow surrounding both Sgr A * and M87 * is known to be optically thin and geometrically thick, while the models analyzed in this analysis assume optically thick and geometrically thin disks.However, our constraints on the magnetic charge are based only on the diameter of the central brightness depression, which can be reasonably expected to remain unchanged in both regimes.Hence, our results remain consistent across both the optically thin and geometrically thick model and the optically thick and geometrically thin model.

Anti degeneracy of mass and magnetic charge
Upon reviewing our previous discussion, it becomes evident that distinguishing between regular BHs and Schwarzschild BHs based on direct images presents a challenge.Moreover, a circular structure at the center of the image contributed by the secondary image can resemble the features of BH photon rings but with smaller observation angles.We previously discussed that the distance between the direct image and the secondary image of a BH can help distinguish BH types.Another distinguishing factor is the impact of luminosity.However, we have yet to delve into the topics of BH mass and magnetic charge.Previous research has indicated that the mass of BHs is directly proportional to the shadow radius, while the effect of magnetic charges is the opposite.It is pertinent to explore whether the mass and charge parameters of BHs can be adjusted to strike a balance between these two factors.If such an adjustment is possible, it might render it challenging to differentiate between the three types of BHs solely based on these parameters.In conclusion, differentiating between regular BHs and Schwarzschild BHs based solely on direct images is a daunting task, but a feature of minimal distance between the innermost region of the direct image and the outermost region of the secondary image can serve as a distinguishing feature.Luminosity also plays a significant role.Further research into the impact of BH mass and magnetic charge parameters is essential to determine their utility in distinguishing between the three types of BHs.
We simulated the apparent images of three types of BHs: the Schwarzschild BH with a mass of 1.732 × 10 6 M ⊙, the Hayward BH with a mass of 1.068 × 10 6 M ⊙, and the Bardeen BH with the same mass as the Hayward BH.In all simulations, the magnetic charge was set to g = 1.5.The celestial coordinates, δ(ξ), were obtained from X(Y ) [47] where δ represents the mass ratio of the BH to the Sun, while D represents the fixed distance of 5kpc.The spin of a rotating BH has a significant impact on the narrowing of the gap between the photon ring and the innermost stable circular orbit.This is due to changes in the BH's geometry and gravitational field.Importantly, this effect can potentially introduce degeneracies with magnetic charge.In certain situations, spin and magnetic charge can have similar effects, while our study primarily focuses on non-rotating BHs.In practical astronomical observations, the presence of rotating BHs needs to be taken into account.Understanding how spin affects the constraints on magnetic charge becomes crucial for accurately measuring BH properties.While our study does not delve into a detailed analysis of rotating black holes, further investigation into the influence of spin on this matter may prove beneficial for achieving a more comprehensive understanding of the problem and its astronomical implications.

Conclusions and Discussion
In this analysis, we have investigated the optical appearance of three distinct types of BHs surrounded by an optically thick accretion disk.Our differentiation among these BHs was based on their observable features and calculated their respective effective potential.Our analysis revealed that for a Schwarzschild BH, the peak of the effective potential occurs at approximately ≃ 3r g .The effective potential for a Bardeen BH is larger, and the radius corresponding to the peak of the effective potential is smaller than that of a Hayward BH.Furthermore, we determined the deflection angle for each type of BH using both numerical integration and semi-analytical methods.Our results demonstrated that although there were differences between the two methods, both successfully depicted the deflection of light around the BH.Furthermore, we employed a ray-tracing code to provide photon trajectories around these three types of BHs, providing further insights into their behavior in an optically thick accretion disk environment.
We employed the Novikov-Thorne thin disk model to investigate BHs surrounded by an accretion disk, utilizing both numerical integration algorithms and semi-analytical methods.Our findings indicate that for smaller observation angles, the secondary image becomes enclosed within the direct image, giving rise to a structure resembling a "photon ring".Additionally, we observed that the degree of separation between the direct and secondary images increases with the inclination angle of observation.Upon comparing three types of BHs, we found that the Hayward BH has a smaller orbital size than the Bardeen BH under similar parameters.
Based on our investigation into the optical appearance of the accretion disk and the corresponding observation images at various observation angles, we have discovered that the luminosity in the region near the BH on the inner side of the accretion disk is higher than that on the outer side.This phenomenon can be attributed to the higher material density in closer proximity to the BH, leading to a greater susceptibility to the BH's gravitational influence.Consequently, this results in more intense heating and radiation.Our findings also indicate that the Bardeen BH and Schwarzschild BH share a similar appearance.However, there is a significant difference in radiance between them.This discrepancy suggests that the spatiotemporal structure of a regular BH enables the surrounding accretion disk to emit more photons.This effect is driven by the more pronounced spatiotemporal curvature generated by magnetic charges.Furthermore, our analysis of gravitational redshift reveals that the redshift area gradually expands with increasing observation inclination angle.Specifically, our results demonstrate that Hayward BHs exhibit a smaller redshift region, while Bardeen BHs exhibit a larger redshift region.
Based on the results presented in Fig. 7, we can conclude that the inclination angle of observation has a significant impact on the observed outcome.Notably, regardless of whether the dip angle increases or decreases, a distinct accumulation of brightness occurs on the left side of the accretion disk.This phenomenon can be attributed to the influence of matter movement along the accretion disk towards the BH.The velocity and density distribution of the matter affect the probability of light absorption and reemitted, ultimately resulting in the observed brightness accumulation.Additionally, the brightness accumulation is a geometric effect that depends on the angle and path of light emission.Our research reveals that the relationship between the Bardeen BH and the Schwarzschild BH is closer, while the gap between the secondary image of the Hayward BH and the innermost stable circular orbit is very small, which is a characteristic different from other types of BHs.When fitting the functional model to the observed flux at various observation inclinations, we find that the exponential function effectively describes the relationship between the observed flux and the magnetic charge parameter.
After conducting our study, we have concluded that distinguishing between regular BHs and Schwarzschild BHs based solely on direct images poses a significant challenging.However, we have identified that the small spacing feature between direct and secondary images can be utilized to differentiate Hayward BHs.Furthermore, we have observed that Bardeen BH images can function as a degenerate projection of Schwarzschild BHs, with the anti-degeneracy effect of mass and magnetic charge closing the gap between them.

Fig 1 .
Fig 1. Effective potential as a function of radius for three types of BHs.The black line corresponds to the Schwarzschild BH, the green line corresponds to the Bardeen BH, and the blue line represents the Hayward BH.The BH mass is set as M = 1 and magnetic charge is g = 0.5.

Fig 2 .
Fig 2. The light trajectories of different BHs in the polar coordinate (b, ϕ).The BHs are shown as the black disks, and the dashed black lines represent the BHs photon ring orbits.Left Panel: Schwarzschild BH.Middle Panel: Hayward BH.Right Panel: Bardeen BH.The BH mass is taken as M = 1 and magnetic charge is g = 0.5.

Fig 3 .
Fig 3.The direct (solid line) and secondary (dashed line) images of three types of BHs accretion disks with different the observation angles.Left Panel: Schwarzschild BH.Middle Panel: Hayward BH.Right Panel: Bardeen BH.The BH mass is taken as M = 3M ⊙ and magnetic charge is g = 1.

F-FFig 4 .
Fig 4. The radiation flux images of three types of BHs accretion disks with different observation angles.Left Panel: Schwarzschild BH.Middle Panel: Hayward BH.Right Panel: Bardeen BH.The BH mass is taken as M = 3M ⊙ and magnetic charge is g = 1.

Fig 5 .
Fig 5.Fitting of Fourier function (black solid line), Gaussian function (bule solid line), and Exponential function (green solid line) to the accurate value of the radiation flux F and magnetic charge g (red dots).The BH mass is taken as M = M ⊙ .

Fig 6 .
Fig 6.Red shift distribution of direct and secondary images of accretion disk at different observation angles.Left Panel: Schwarzschild BH.Middle Panel: Hayward BH.Right Panel: Bardeen BH.

Fig 7 .
Fig 7. The observation images of three types of BHs accretion disks with different the observation angles.Left Panel: Schwarzschild BH.Middle Panel: Hayward BH.Right Panel: Bardeen BH.The BH mass is taken as M = 3M ⊙ and magnetic charge is g = 1.

Fig 8 .
Fig 8. Fitting of Fourier function (black solid line), Gaussian function (bule solid line), and Exponential function (green solid line) to the accurate value of the observed flux obs and magnetic charge g (red dots).Top Panel: Hayward BH.Bottom Panel: Bardeen BH.The BH mass is taken as M = M ⊙ .

Fig 9 .
Fig 9. Shadow diameter of the three types BHs as a function of the magnetic charge.The green and orange shaded regions represent the regions of 1σ and 2σ confidence intervals, respectively.
Figure 10 displays images of three types of BHs corresponding to specific mass and magnetic charge parameters, with an observed inclination angle of θ 0 = 17 • .It is apparent that the structures of the Schwarzschild and Bardeen BHs are essentially the same, featuring similar circle diameters of approximately 40µas.Therefore, we can conclude that the observed Bardeen BH image can be considered a degenerate projection of the Schwarzschild BH image.The anti-degeneracy effect of mass and magnetic charge eliminating the gap between them.

Fig 10 .
Fig 10.The complete apparent images of the Schwarzschild BH with mass 1.732 × 10 6 M ⊙ the Hayward and Bardeen BHs with mass 1.068 × 10 6 M ⊙, and the magnetic charge is g = 1.5.the observation alangle and BH distance are θ 0 = 17 • and 5kpc.