Shadow, lensing, quasinormal modes, greybody bounds and neutrino propagation by dyonic ModMax black holes

Motivated by recent work on the Modified Maxwell (ModMax) black holes [Phys.Lett.B 10.1016/j.physletb.2020.136011], which are invariant in duality rotations and conformal transformations founded in [ Phys.Rev.D 10.1103/PhysRevD.102.121703], we probe its effects on the shadow cast, weak field gravitational lensing, and neutrino propagation in its vicinity. Using the EHT data for the shadow diameter of Sgr. A* and M87*, and LIGO/VIRGO experiments for the dyonic ModMax black hole perturbations, we find constraints for ModMax parameters such as $Q_\text{m}$ and the screening factor $\gamma$. We also analyze how the shadow radius behaves as perceived by a static observer and one that is comoving with the cosmic expansion. The effect of the ModMax parameters is constant for a static observer, and we found That it varies when the observer is comoving with cosmic expansion. We also analyzed its effect on the weak deflection angle by exploiting the Gauss-Bonnet theorem and its application to Einstein ring formation. We also consider the finite distance effect and massive particle deflection. Our results indicate that the far approximation of massive particle gives the largest deflection angle and amplifies the effect of $Q_\text{m}$ and $\gamma$. Then we also calculate the quasinormal modes and greybody bounds which encode unique characteristic features of the dyonic ModMax black hole. With the advent of improving space technology, we reported that it is possible to detect the deviation caused through the shadow cast, Einstein rings, quasinormal modes, and neutrino oscillations.

.121703], we probe its effects on the shadow cast, weak field gravitational lensing, and neutrino propagation in its vicinity. Using the EHT data for the shadow diameter of Sgr. A* and M87*, and LIGO/VIRGO experiments for the dyonic ModMax black hole perturbations, we find constraints for ModMax parameters such as Qm and the screening factor γ. We also analyze how the shadow radius behaves as perceived by a static observer and one that is comoving with the cosmic expansion. The effect of the ModMax parameters is constant for a static observer, and we found That it varies when the observer is comoving with cosmic expansion. We also analyzed its effect on the weak deflection angle by exploiting the Gauss-Bonnet theorem and its application to Einstein ring formation. We also consider the finite distance effect and massive particle deflection. Our results indicate that the far approximation of massive particle gives the largest deflection angle and amplifies the effect of Qm and γ. Then we also calculate the quasinormal modes and greybody bounds which encode unique characteristic features of the dyonic ModMax black hole. With the advent of improving space technology, we reported that it is possible to detect the deviation caused through the shadow cast, Einstein rings, quasinormal modes, and neutrino oscillations. holes and try to probe the ModMax BH parameters by studying gravitational waves (GWs) at the ringdown stage.
In 1974 Hawking showed that there are 'grey holes' and not black holes because of their radiation, which is impossible classically for ingoing particles [120]. The quantum field theory calculations near the black hole's horizon provide the emission of quantum radiation due to the creation and annihilation of particles. The Hawking radiation is moving on a curved spacetime ( as a potential barrier) so that some radiations are reflected into the black hole and the rest travel to spatial infinity [121]. These deviations from the blackbody radiation spectrum, as seen by an asymptotic observer, are known as the grey-body factor, which is a synonym for transmission probability [110,[122][123][124][125][126][127][128], and there are various methods to calculate it. One can use the matching technique [129], the WKB approximation [130] or the rigorous bound [131,132]. On the other hand, in 1968, Matzner studied the absorption and scattering of a massive scalar field from hitting a black hole [133]; hence he showed that the total cross-section for the absorption process vanishes. In 1978 Sanchez calculated the absorption spectrum of the Schwarzschild black hole and obtained the total absorption cross section in the Hawking formula [134]. Sanchez showed that the absorption spectrum as a function of the frequency makes clear oscillations characteristic of a diffraction pattern, where the oscillations occur around the constant value of the geometrical optics with decreasing amplitude and constant period. The absorption cross-section for black holes has been studied in the literature for various black hole spacetimes [46,108,[135][136][137][138][139][140][141][142][143][144][145][146], especially for the low-frequency behaviour of the absorption and for the zero-frequency limit (it becomes equal to the surface area of the black hole horizon) [136,147,148], as well as the high-frequency behaviour of the scalar absorption via Sinc approximation [134,149,150].
In this paper, we also investigate the propagation of neutrinos in dyonic ModMax BH. We focus, in particular, on the neutrino flavour oscillation and spin-flip of neutrinos when they scatter off BH. Clear evidence that neutrinos may oscillate in different flavours comes from experiments on neutrinos [151][152][153]. These experiments indicate that neutrinos are massive particles and, therefore, that physics beyond the Standard Model is required to incorporate the neutrino particle properly. An important aspect related to the properties of neutrinos is their interaction with external fields. These interactions can be either electromagnetic or gravitational. In these cases, the formulas of the oscillation probabilities of mixed neutrinos (different flavours) are affected by the external field, compared to the formulas computed in vacuum [154]. For example, the interactions of neutrinos with external electromagnetic fields may induce the so-called spin oscillation and/or spin-flavour oscillations, that is a helicity transition of neutrinos with different helicities [155]. It is also well known the gravitational interaction or propagation of neutrinos in the curved background can affect the neutrino oscillations probability, or induce the change of the polarization of a spinning particle . In this paper we study the neutrino oscillation probability for neutrinos propagating in dyonic ModMax BH. We also study helicity transitions of neutrinos, hence the transition ν f L → ν f R , where f = e, µτ , corresponding to the case in which neutrino flavour remains fixed (the generic case of neutrino spin-flavour oscillations has been studied in [177]). This analysis is important because neutrinos are produced with a fixed left-handed polarization in the Standard Model so that if a transition to right-handed polarization occurs, these neutrinos become sterile (they interact only gravitationally). As a consequence, a detector would register a different neutrino flux, providing in such a way a signature of the coupling of neutrinos with the gravitational background. Here we shall consider only the gravitational field described by ModMax BH.
Our work's layout is as follows: Sect. III is devoted to investigating the shadow cast behaviour and the effect of a co-moving observer. In Sect. IV, we analyzed the behaviour of the weak deflection angle of both massive and null particles. Then, in Sect. V -VII, we study the effect of the screening parameter on the spherical infalling accretion, quasinormal modes, and bounds of the greybody factors. Finally, in Sect. VIII, we investigate its effect on neutrino oscillations. In particular, we focus on the neutrino flavour and spin transition in the background described by a dyonic ModMax BH. We shortly analyzed also the effect on nucleosynthesis processes. Conclusions are discussed in Section IX. Throughout this paper, we used geometrized units G = c = 1, and the metric signature (−, +, +, +).

II. DYONIC MODMAX BLACK HOLES
In this section, we study the spherically symmetric metric of Dyonic ModMax black hole solutions derived in [9] as follows; where dΩ := dθ 2 + sin 2 θdφ 2 and (t, r, θ, φ) with the metric function: We plot the above metric function shown in Fig. 1. Here, we can see theoretically that only a single horizon is formed for the Schwarzschild case, while there are inner and outer horizons for the RN case when Q e = 0.75M . As the plot indicates, we can observe how the screening parameter affects these null boundaries. The event horizon of the ModMax black hole is located at: Note that the mass must always be larger than the charge; however, it is not valid for the case of the ModMax black hole because, in the extremal case of the ModMax black hole, it occurs at r + = r − : so that γ > 0, then M ext < |Q m |; similarly with non-linear electrodynamics [178,179].

III. SHADOW CAST
In this section, we first analyze the behaviour of the ModMax black hole's shadow as the screening parameter γ varies. Since spacetime has spherical symmetry, we can analyze null geodesics θ = π/2, along the equatorial plane. To derive the null particle's equation of motion, consider the Lagrangian where we wrote A(r) = f (r), B(r) = 1/f (r) and C(r) = r 2 . With the least action principle, two constants of motion can be derived: where we can define a useful constant called the impact parameter: Going back to the metric, ds 2 = 0 for null geodesics allows us to derive how the r changes with φ. In other words, where we have defined With the above function, one can solve or locate the photonsphere radii via h (r) = 0 where the prime denotes the derivative with respect to r. For our case, we obtained A non-rotating black hole, such as in this study, only produces a circular shadow silhouette. The shadow is caused by photons escaping the photonsphere due to tiny perturbations. In the Schwarzschild case, in particular, r ph = 3M , while the shadow radius is R sh = 3 √ 3M . It is important to note that this shadow radius coincides with the critical impact parameter b crit . In some situations, if the cosmological constant is considered, this is not the case since the photons are affected by the astrophysical environment as it travels toward the receiver, usually a remote and static observer [62].
Considering the static observer and following the formalism found in Ref. [62], a careful inspection of the line element and situation allows us to define or, written in another way, With the orbit equation and elementary trigonometry, it is easy to see that where b crit is associated with the photonsphere radius. It can be derived by satisfying the condition d 2 r/dφ 2 = 0 and imposing r → r ph : where in our case yields Using Eqs. (9), (10), and (15) to Eq. (13), we obtain an exact formula for the shadow radius as Let us first discuss the observational constraint of the screening parameter using the obtained data from M87* [31] and Sgr. A* [32]. First and foremost, these astrophysical black holes are indeed rotating. However, as pointed out in Ref. [180], the constraints for the spin parameter a on Sgr. A* is not clear, but recent constraints place the spin parameter at a 0.1M due to the impact of frame-dragging precession on the orbit of S-stars [181]. Similar uncertainties for the spin of M87 are also present: Using accretion physics, recent estimates have shown that 0.20 < a < 0.50 [182]. Other authors [183] reported that through the measurements of the average power of the relativistic jet and using the upper limit of mass accretion rate for black holes in general relativistic magnetohydrodynamic models of jet formation, the spin parameter of M87* would likely fall from a ≥ 0.40M for the prograde case, and a ≥ 0.50M for the retrograde case. It was also shown in Ref. [180] that slow spin does not affect the shadow radius dramatically, and a ∼ 0 can be safely assumed for constraining parameters coming from different black hole models. Finally, astrophysical black holes such as these are more likely to contain no charge due to the neutralization process of the accreting ionized plasma. In Ref. [184], constraints indicate that Q m /M can be very close to zero. Nonetheless, in the following analysis, we will assume the role played by Q e and Q m , as well as the parameter γ. We will analyze the Modmax black hole for low and average values of Q e . Finally, we need to remark that we will only analyze the classical shadow silhouette, where its boundary is dominated by the glowing accretion disk causing it to become invisible. The dark spot in the EHT image of M87* and Sgr. A* is the imprint of the event horizon's shadow. It is caused by highly redshifted escaping photons coming from accreting matter near the event horizon. As a result, this shadow is smaller than the classical shadow [65]. According to the seminal EHT papers [31,32], the reported classical shadow angular diameter is θ M87* = 42 ± 3 µas. Other parameters are the distance of the M87* from the Earth, which is D = 16.8 Mpc, and the mass of the M87* is M M87* = 6.5 ± 0.90x10 9 M . For Sgr. A* the angular shadow diameter is θ Sgr. A* = 48.7 ± 7 µas (EHT), the distance of the Sgr. A* from the Earth is D = 8277 ± 33 pc and the mass of the black hole is M Sgr. A* = 4.3 ± 0.013x10 6 M (VLTI). Let us now calculate the diameter of the shadow size in units of the black hole mass using which is just the standard arc-length formula. Meanwhile, the theoretical shadow diameter can be obtained via d theo sh = 2R sh . Therefore, by using Eq. (17), we get the diameter of the shadow image of M87* and Sgr. A* as d M87* sh = (11 ± 1.5)M , and d Sgr. A* sh = (9.5 ± 1.4)M respectively. In Fig. 2, we plotted the allowed values of Q m while maintaining the screening parameter γ fixed. Such a plot will give how the shadow diameter behaves at our location from Sgr. A* and M87* as we vary Q e and Q m . We included a very small value for Q e to represent the notion that these astrophysical black holes may contain a negligible charge. We see that if γ = 1, it amplifies the effect of Q e and Q m , in contrast with γ = −1. Nevertheless, we observe how γ affects the curve. An important observation of this plot is that as γ decreases, the line becomes more curved, allowing smaller values of ±Q m within the 1σ uncertainty. It is also worth noting that if Q e is close to zero, the value of γ is irrelevant when Q m is also close to zero. While it happens above the mean for Sgr. A*, it happens below the mean for M87*. As a result, the range for Q m at 1σ level is lesser than Sgr. A*.
Let us now plot how a static observer, at different points r obs , will perceive the shadow size of Sgr. A* and M87*. We plotted this in Fig. 3, along with the Schwarzschild and RN cases, where Q = 0.10M to see its effect. Our current position from these SMBHs is represented in the right inset plot, where we zoomed in on the effects of the screening parameter γ and Q m , manifest in visuals in Fig. 2. The pattern is the same for M87*, except that we are much farther. Such a similarity is expected because the ModMax parameters have no dependence on r. Finally, it is interesting how the shadow behaves near the black hole, where we can see an abrupt increase in the shadow radius until r obs /M ∼ 316. The case where γ < 0 shows that the shadow forms early or nearer the event horizon, compared to the case where γ > 0. However, as r obs increases, γ > 0 gives a larger shadow. The point where the curves intersect is also interesting, for it represents that θ sh = π/2, indicating that in such a location, half of the sky is dark [62].
A. Behavior of the shadow due to a co-moving observer In reality, observers are not static. Thus, this study aims to investigate the effect of the screening parameter on the shadow perceived by an observer co-moving with the cosmological expansion of the Universe. Here, we apply the recent formalism in Ref. [60], first focusing on the dark energy-dominated Universe. Then, we also extend the analysis to matter and radiation-dominated Universes. An observer in a co-moving frame is better understood using the McVittie metric. According to Ref. [185], we can find that where dΩ 2 = sin 2 ϑdϕ 2 + dϑ 2 , and a(t) = e H0t is the scale factor. In addition, H 0 is the present value of the Hubble constant. The relation between the cosmological constant Λ and H 0 is At present, Λ = 1.1056x10 −52 m −2 , which makes the distance of the cosmological horizon to be d cosmo = 9.51x10 25 m. Note the time dependence of Eq. (18), but such dependence is negligible if one is near the black hole. Being near means comparing the distance of the black hole to the observer and the cosmological horizon. Indeed, even if we are so far to M87* or Sgr. A*, our distances from these black holes are so small compared to the scale of the cosmological horizon. Then by [60], t 0 − t << H −1 0 . Thus, as a(t) ∼ a(t 0 ) =constant and using x = a(t 0 )r, we can recast Eq. (18) as which is the RN black hole within isotropic coordinates. Introducing where in the weak field limit (M, Q e , Q m ∼ 0), R ∼ x, we recover a form similar to the RN metric: Let the present time be t 0 where the observer is at some radial position r in , who observes the black hole shadow in the strong field regime. Here, the subscript "in" indicates the inner region. Then we should have x in = a(t 0 )r in and using these, Eq. (21) can be recast as Then using Eq. (22), an observer co-moving with the spacetime in Eq. (20) will then observe the shadow radius R in as where r ph is still given by Eq. (10). We note that Eq. (24) is the "inner" solution to the shadow radius where the effect of the cosmological expansion is considered negligible. Next, we consider the solution to the outer region in the weak field regime of the black hole while considering a stronger effect of the cosmological expansion. It is easy to see how Eq. (18) will reduce to an FRW metric in such a case: In this spacetime, it is well-known that the effective linear shadow radius L sh , now applied to the shadow, is given in terms of the angular size of the black hole shadow Ψ cosmo as where and Here, Ω mat , Ω rad , Ω Λ are present dimensionless density parameters for matter, radiation and dark energy, respectively. The angle Ψ cosmo is considered to be so tiny that the relation sin(Ψ cosmo ) ∼ Ψ cosmo remains valid. Due to the dependence of Eq. (26) z, we can assume that x out ∼ R out due to the enormous distance from the black hole. It enables one to connect z to R out via by [60] Between these inner and outer regions, there is a location where the influence of the black hole and cosmological expansion begins to lose and gain, respectively. This overlap region, denoted by R o , is still remote from the black hole. Thus, the approximation z << 1 is valid since the scale factor is a(t 0 ) at the location of such an observer at time t 0 . Thus, we have which is still equal to the weak field approximation of Eq. (24): Then, we can match L o and L sh [60]. Using R out instead of R o , the effective shadow radius R cosmo in the outer region of rapid expansion is obtained: Finally, the perceived shadow radius R approx by an observer co-moving with the cosmic expansion is approximated via composite solution [60]: Note that in Ref. [60], the authors analyzed the situation using the angular radius, whereas, in this study, we slightly extended the formulation to investigate the shadow radius behaviour. Eq. (33) enables us to calculate the approximated shadow radius using different models for our Universe. In particular, let us take, for example, a Universe dominated by dark energy, where, after evaluating Eqs. (27) and (29), we find Now for brevity, let Then, Eq. (33) implies that the approximate shadow radius seen by a co-moving observer is Consider a Universe dominated by matter where and we find Lastly, for the radiation dominated Universe, and The shadow radius behaviour seen by an observer in these models of the Universe can be visualized by plotting Eqs. (36), (38), (40) numerically for an immediate comparison. In Fig. (4), we note that the Hubble constant is scaled to see the cosmological expansion's overall effect at short distances. However, it is understood that without this scaling and if one uses the experimental value of H 0 , the effect can only be noticeable at distances near the cosmological horizon. Note that at low values of R obs /M , we can see that the difference between these models is almost negligible. However, for large values of R obs /M , we can observe the differences between these models of the Universe. Worth noting is the radiation-dominated Universe since there is some point where a peak is reached, and the shadow rapidly decreases in size for a very small change in R obs /M . For dark energy and matter-dominated Universe, the co-moving observer sees no limit on large the shadow radius will increase. As a final remark, it can also be observed in the plot that the effect of the screening parameter γ is affected by the co-moving state of motion and the type of Universe the observer lives in.

IV. DEFLECTION ANGLE USING GAUSS-BONNET THEOREM IN WEAK FIELD LIMITS
In this section, we aim to investigate the weak deflection angle of both massive and null particles by exploiting the Gauss-Bonnet theorem (GBT), which is originally stated as [77,186], In this equation, the Gaussian curvature K that describes the domain D can be oriented in a 2D curved surface S with infinitesimal area element dS. The boundaries of D are given by ∂D a (a=1, 2, .., N ), and the geodesic curvature κ g is integrated over the path d in a positive sense. Also, the jump angle is denoted by θ a , wherein χ(D) is the Euler characteristic equal to 1 due to the non-singularity of D.
In Ref. [85], the GBT is applied to accommodate black hole metrics with non-asymptotic flatness. Although our metric is asymptotically flat, we can still use this formalism as a special case. With the photonsphere radius r co as being part of the quadrilateral, which serves as the domains of integration, the weak deflection angle can be derived viâ where S and R are the radial positions of the source and receiver, respectively. Furthermore, the infinitesimal curve surface dS is given by and φ RS is the coordinate position angle between the source and the receiver defined as φ RS = φ R − φ S , and g is the determinant of the Jacobi metric. line element of a static, spherically symmetric (SSS) spacetime where dΩ 2 = dθ 2 + sin 2 θdφ 2 , the Jacobi metric is defined as Here, E is the energy of the massive particle defined by where v is the particle's velocity. Due to spherical symmetry, we can analyze along the equatorial plane without losing generality. The determinant of the Jacobi metric can then be sought off as We need the orbit equation to determine the expressions for φ RS . For time-like particles, this can be derived through g µν dx µ dx ν = −1 which results to where we have used the substitution r = 1/u and the angular momentum of the massive particle and b is the impact parameter. With the metric coefficients, we find Solving the above iteratively, we find The Gaussian curvature K, in terms of affine connections and determinant g, is defined as since Γ φ rr = 0 for Eq. (45). Then with the analytical solution to r co , thus, The prime denotes differentiation with respect to r. The weak deflection angle is then [85], Using Eq. (51) in Eq. (54), we find Next, we find the expression for φ. To do this, we use Eq. (51) and solve for φ. For the source and receiver, we find respectively. Careful observation of these equations will allow us to write φ RS = π − 2φ S . Now, we take note of the following relations: The last property cancels the sine terms in Eq. (56). We find cos(φ S ) as In these plots, the dot represents the critical impact parameters in every case. and cot(φ S ) as Using Eq. (46) and by plugging Eqs. (57)- (61) in Eq. (55), we finally obtain which also involves the finite distance u S and u R . The above expression can still be further approximated as b 2 u 2 ∼ 0: For the case of photons where v = 1, we find The above agrees with the result in Ref. [187]. How will the screening parameter γ affect the weak deflection angle, in combination with the finite distance effects? These can be seen in Fig. 5, where we also compared theα caused by both massive and null particles. Our results indicate thatα is increased as both the source and the receiver are approximated at r → ∞. Interestingly, massive particles also give a larger deflection angle than photons. When it comes to the effect of γ, it does not discriminate whether it is massive or null particles. That is, γ > 0 has a larger value forα as compared to γ < 0.
We also remark that γ < 0 gives a much lower bound to the critical impact parameter. It means that near (but greater) than the critical impact parameter,α has more range to exist. A useful astronomical application of the weak deflection angle involves the Einstein ring. Let D SL and D LR be the position of the source and receiver from the black hole (lensing object L), respectively. The thin lens approximation implies that D RS = D SL + D LR , and the position of the weak field images is given by An Einstein ring will be formed when β = 0, and the above equation simplifies to Finally, using the relation b = D LR sin θ ∼ D LR θ we can obtain where the parameter = Q 2 e + Q 2 m e −γ /b 2 . We plot Eq. (67) as shown in Fig. 6. In this plot, we compared four cases, and the solid lines represent the effect of the screening parameter on the Einstein ring. As γ decreases, it can be gleaned that the Einstein rings angular radius decreases for a given impact parameter. Furthermore, how fast these angular radii change with b/M is less than cases with no screening parameter (dashed line). As b/M → ∞, the Einstein rings approach the Schwarzschild case. Note that with these parameters, the angular radius of the formed Einstein ring can be experimentally detected since the lowest value we can get in the plot is ∼ 10000µas. Furthermore, if the deviation from Q m = 0.25M and Q m = 0.50M at γ = −1 is examined, one needs a device capable of detailing a difference of 570µas at b/M = 10. We see that as b/M → ∞, the effect of the dyonic Modmax parameters lessens.

V. SPHERICALLY INFALLING ACCRETION
In this section, we study the spherically free-falling accretion model on the BH from infinity using the method defined in [188,189]. This method gives us realistic visualization of the shadow cast with the accretion disk. In reality, the actual image of the BH can not be seen as an apparent boundary in the universe. Moreover, it is not realistic to use a static accretion disk model because there is a moving accretion disk around a black hole, and also it provides a synchrotron emission from the accretion. To do so, we first study the specific intensity observed at the observed photon frequency ν obs solving this integral along the light ray: It is noted that b γ is the impact parameter, j(ν e ) is the emissivity/volume, dl prop is the infinitesimal proper length and ν e stands for the photon frequency of the emitter. Here we define the redshift factor for the infalling accretion as follows: where the 4-velocity of the photon is k µ =ẋ µ and 4-velocity of the distant observer is Using the relation of k α k α = 0, one can derive k r and k t which is a constant of motion for the photons: Note that sign of ± shows that the photon gets close to (away from) the black hole. Then the redshift factor g and proper distance dl γ can be written as follows and We then consider only the monochromatic emission for the specific emissivity with rest-frame frequency ν * as follows: Afterwards, the intensity equation given in (68) becomes We investigate the shadow cast with the thin-accretion disk of the black hole in ModMax. First of all, we solve the above equation numerically using the Mathematica notebook package [146], (also used in [42,46,190,191]) and this integration of the flux show the effects of the parameters of the ModMax γ on the specific intensity seen by a distant observer for an infalling accretion in Figs. (7, 8 and 9). These plots in Figs. (7, 8 and 9) show the specific intensities for various values of the parameter γ versus b observed by the distant observer.
We observe that increasing the value of b, increases the intensity first. Afterwards, intensity reaches the peak value sharply where the photons are captured by a black hole quickly ( at the photon sphere). It is seen that after the peak value, intensity slowly decreases. Moreover, we show the shadow cast of the black hole in the 2-dimensional image with a photon sphere by a distant observer in (X, Y) plane where the dark centre of it the event horizon is located, and it is circled by a bright ring with a strongest luminosity (photon sphere). It can be seen that brightness decreases gradually after the maximum region. Hence, we show the effect of the screening parameter γ on the black hole luminosity of the shadow cast where the intensity decreases with increasing the value of the screening parameter γ as seen in Fig. 9.

VI. QUASINORMAL MODES OF DYONIC MODMAX BLACK HOLES
Since the first detection of gravitational waves (GWs) from the coalescence of two stellar-mass black holes in 2015 by the LIGO/VIRGO collaborations [192], gravitational wave physics has begun. Then data from gravitational waves are analyzed to test alternative theories of gravity and different models of compact objects. Perturbative analysis of black hole spacetimes dominated by 'quasinormal ringing' is used to do this. Quasinormal modes (QNMs) are oscillations with complex frequencies with energy dissipation. The complex frequencies of QNMs have the characteristic properties of the BHs, such as mass, charge, and angular momentum, independent of the initial perturbations.
In this section, we use the method of WKB approximation to study quasinormal modes (QNMs) of the Dyonic ModMax BHs. This method is known as an effective way to derive QNMs, which is firstly used in [193] by Schutz and Will, then Iyer and Will extended to the third order of WKB approximation in [194]. Recently, WKB approximation has been extended to sixth order and higher-order cases in [130,195]. To study QNMs of the Dyonic ModMax BHs, first, we use a massless scalar field perturbation in the background of the black hole 81. The Klein-Gordon equation, which is used for the scalar field, can be written as follows: where g is the determinant of the metric g µν . Then to solve the above equation, we first use the separation of variables: It is noted that Y l (θ, φ) are spherical harmonics with the multipole number l = 0, 1, 2, . . . . We then obtain the Schrodinger-like wave equation (Regge-Wheeler-Zerilli equations) [117,196]: where the following relaton defines the tortoise coordinate r * : with the effective potential for the scalar field V (r) = A(r) l(l + 1) where the line element is written as: Note that A(r) = B(r) = f (r) in Eq. 2.
One can think of the effective potentials as the potential barriers which provide decaying on the event horizon and at infinity. Apply the following ansatz where e −iωt is the time evolution of the scalar field, then we obtain the time-independent wave equation where ω is the complex QNM frequency (or eigenvalues of the above wavelike equation) written in the form ω = Re(ω) + i Im(ω) to be determined after the appropriate boundary conditions are applied (ψ(r) corresponds to the purely outgoing wave at spatial infinity; ψ(r) behaves as the purely ingoing wave at the event horizon): Here we use the numerical method of the sixth-order WKB formula to calculate QNMS, which is based on the WKB expansion of the wave function at both the event horizon and spatial infinity. It is then matched with the Taylor expansion is near the peak of its effective potential. There are two turning points/monotonic decay for this potential. To do so, we should solve the following equation [194] i where Λ j (n) (N indicates the order of the WKB method) are the WKB correction terms,r * is the location of the maximum of the QNM potential V (r * ) in the tortoise coordinate. Note that Λ 2,3 are defined in [194] 1 and Λ 4,5,6 are given in [195]. The dependence of the QNM frequencies on the screening parameter γ is qualitatively different for lowest and higher multipoles, as seen in Table I. We can see in Table I and II that both real and imaginary parts of the ω decrease when the parameter γ is increased. The scalar field QNMs for different values of l is plotted in 10, as well as for l = 0, 100 for different values of γ are presented in Figs. 11, 12 and 13. Next, we study the eikonal quasinormal modes for dyonic ModMax black holes. This method is also known as the geometric optics method due to its relation with the parameters on the null geodesics [197,198]. The imaginary part of the  quasinormal mode frequency (Im ω=-ω I ), which is responsible for the temporal, exponential decay, can be calculated in the large-l limit (l → ∞) (only the g tt component is relevant) [199,200] i.e., as the angular momentum number describing our mode solution becomes very large, as follows:    with the angular velocity Ω ph : and Lyapunov exponent λ L : where n is the overtone number and take values n = 0, 1, 2, .... Note that the eikonal limit is independent of the spin of the perturbation, so that black holes' scalar, electromagnetic, and gravitational perturbations give the same behaviour in the eikonal limit [201]. Table III shows that the real parts decrease, but the imaginary part of the QNMs increase with the increasing parameter γ. We can conclude that these modes are stable because the imaginary parts of the QNMs frequencies are negative. The decay rates (imaginary part) of QNMs frequencies increase as the screening parameter γ increases. To study time-domain profiles of the scalar field perturbation of the BHs, we use the initial disturbance as a Gaussian wave packet: It is noted that we assume that σ = 2,r * = −40, and A = 10. Then apply the boundary conditions such that the wavefunction vanishes at (r * = −200, r * = 250) to observe the differences between the Schwarzschild black hole and ModMax black hole spacetimes. Afterwards, the time-domain profiles of the scalar field perturbation for l = 0 and l = 2 cases are plotted numerically in Figs. 14 and 15. The ModMax black hole is given in blue, while the Schwarzschild case is in black. The logarithmic plot shows that the ModMax frequency is slightly lower. For the l = 0 case, the plot shows the relatively short period of quasinormal ringing, which makes it hard to extract values of frequencies with good accuracy.
To examine the gravitational signal emitted by an oscillating black hole within the bandwidth of the LIGO/VIRGO interferometers, its frequency must be in the range of around 10 − 40 Hz, or for the LISA interferometer, about 10 −4 − 10 −1 Hz [202]. We assume M = αM , M = 1.48 × 10 5 cm , for the frequencies and damping times obtained by [202]   After calculating the fundamental frequencies, we show the allowed range of masses of oscillating ModMax black holes which LIGO and LISA can detect: for the fundamental frequency range ν ∈ [12 Hz, 1.2kHz], and for LISA the detectable mass range is Note that it corresponds to frequencies ν ∈ 10 −4 , 10 −1 Hz. Moreover, in the future, LISA will be able to detect signals emitted by the oscillations of the massive black hole at the centre of our Galaxy SGR A*, the mass of which is M = 4.1 × 10 6 M [32].

VII. BOUNDS OF GREYBODY FACTORS AND HIGH-ENERGY ABSORPTION CROSS-SECTION
Greybody factor (GF) is a quantity related to the quantum nature of a black hole, and its high value gives a high probability that Hawking radiation can reach infinity, so the GF of test fields is important for estimating the intensity of Hawking radiation. Here, we study the rigorous bound on the GF of the dyonic ModMax black hole to probe the effect of γ on the bound. Rigorous bound of the GF was first proposed in [131,132] which gives qualitative description of a black hole as follows: Note that h (r * ) is a function satisfying the condition h(−∞) = h(∞) = w. [131]. Selecting h = w, and substituting the tortoise coordinate r * , we write Here V is the effective potential for the massless scalar field given in (80). We can calculate the bound as follows: The bound reduces to the Schwarzschild case for (γ, Q e , Q m ) → 0, as T Sch ≥ sech 2 2l(l+1)+1 8mw . We illustrate the effect of the screening parameter γ on the greybody bound for a scalar field in the ModMax black hole in Fig. 16. It is seen that when the value of γ parameter increases, the greybody bound T b grows larger.
Afterwards, the high-energy absorption cross-section is studied by applying the Sinc approximation, which Sanchez first studied for the Schwarzschild black hole. It is shown that increasing the frequency for the ordinary material sphere monotonically increases the absorption cross-section oscillated around the constant geometric-optics value for the black hole (related to the photon sphere) [134] which shows the relation between the impact parameter and cross-section of the photon sphere at the critical value.
On the other hand, it was concluded that at low energy scales, the characteristic properties of BHs and the cross-section of BHs equal the BH area [147]. However, at high energies, one should use the complex angular momentum technique to study the geometrical cross-section of the photon sphere [150]. Decanini et al. use the Regge pole techniques to show the relationship between the oscillatory pattern of the high-energy absorption cross-section and the Sinc(x) function (with Sinc(x) denote the sine cardinal sinc(x) ≡ sin x x ) with the photon sphere. This method extensively discussed in the literature [203][204][205][206][207][208][209].
In the limit of high energy, the oscillatory part of the absorption cross-section can be calculated by using [150]: and Lyapunov exponent λ L is given in Eq. ( (88)). The absorption cross-section can also be written as where Ω ph is the angular velocity with the radius of the photon sphere r ph . The total high energy formula for the absorption cross-section is equal to σ abs ≈ σ osc abs + σ geo [149,150]. In Fig. 17 the total absorption cross section for various values of γ is plotted. The numerical analysis shows that the greater values of screening parameter γ, the total absorption cross section becomes more and more unstable for higher energies and exhibits the largest amplitude. Moreover, there is a regular oscillatory behaviour around the high-frequency limit. Fig. 17 presents the absorption spectrum as a function of the frequency, making clear oscillations characteristic of a diffraction pattern, where the oscillations occur around its constant value of the geometrical optics with decreasing amplitude and constant period.

A. Flavour oscillation
In this Section we will study the effects of the metric in Eq. (2) on neutrinos, following the treatment in Ref. [210][211][212][213]. Even if in the Standard Model neutrinos are assumed massless, in the last fifty years there have been various detection of the neutrino flavour oscillation phenomenon. The common explanation for the latter is to take into account the existence of small neutrino masses and that neutrino mass and flavour eigenstates are not coincident. Therefore, neutrinos will be treated as massive particles in this Section. Firstly, let us define the oscillation length: where dl 2 pr = −g ij dx i dx j is the infinitesimal proper distance, with i, j that runs over the spatial coordinates, Φ αβ = Φ α − Φ β and: where L is the angular neutrino momentum. Considering E m 2 α and L = 0, one obtains Using Eq. (101), we show in Fig. 18 the behaviour of the oscillation length in r for the different parameters of the metric in Eq. (2). As it is possible to see, the main differences in General relativity (GR) occur at low values or r, up to a 10% of changes in the neutrino oscillation length. Moreover, we can discuss the effects of gravitational redshift of the considered metric on the energy spectra of ν e andν e Type II supernova explosion. The gravitational effects affect the r-process nucleosynthesis in astrophysical environments, and we might theoretically estimate that with the Y e parameter [212][213][214][215]: with the local neutron to proton ratio R 0 n/p and Γ are defined as: Lν e Eν e L νe E νe , Γ = g 00 (rν e ) g 00 (r νe ) where Lν e and L νe = are the antineutrinos and neutrinos luminosity respectively, Eν e = 25 MeV, E νe = 10 MeV and rν e = 5 km, r νe = 0.9rν e are the two neutrino-sphere radius forν e and ν e respectively. In Fig. 19, we show the behaviour of Y e for the different ratio of the ν e ,ν e luminosity with the metric in Eq. (2). It is possible to see that, compared to GR (red curve), we obtain a difference up to 5%, which could be important in the supernovae' evolution.

B. Spin oscillation
In this subsection, we treat the neutrino spin oscillation problem in the metric in Eq. (2), following the Ref. [213,216,217]. The evolution of the spin vector s a is given by where ζ and s a are defined as with u the four-velocity in the local Minkowski space. In the metric of our interest, Ω g = (0, Ω 2 , 0) and therefore we can use the following representation for ζ = (ζ 1 , 0, ζ 3 ) = (cos α, 0, sin α). Focusing on the neutrino spin oscillation, we define the helicity h = ζ · u/|u|. The initial condition is Moreover, we can write that u +∞ = + √ E 2 − m 2 /m, 0, 0 and therefore h +∞ = cos α. From that, we write the helicity states as where a 2 + + a 2 − = 1 due to the normalization and a 2 + − a 2 − = cos α = h +∞ . Therefore, one finds that where dt/dτ = U 0 and dr/dτ = U 1 . Finally, the angle α +∞ reads Using Eq. (116) it is possible to obtain the spin-flip probability for metric in Eq. (2), shown in Fig. 20. The important feature is that, with the used metric, the spin-flip phenomenon is suppressed. The analysis with the inclusion of a magnetic field will be faced elsewhere.

IX. CONCLUSIONS
In this work, we probed how the dyonic ModMax black hole parameters would affect the shadow, lensing, and neutrino propagation in its vicinity. To do so, we initially constrained the values of Q m under the effect of γ = 1 and γ = −1 charge screening parameters using the EHT data on the shadow diameter of Sgr. A* and M87*. Focusing on Q e = 0.01M − 0.10M , since astrophysical black holes may contain a nearly zero electric charge, we see that the value of γ is irrelevant near Q m = 0. For Sgr. A* occurs above the observed shadow diameter, while it occurs below for M87*. With these parameters, the We have also shown the GR spin-flip probability (red line) for comparison. We have considered neutrinos with γ = E/m = 10.
deviation to the observed shadow radius is still within 2σ uncertainty, while for M87*, lower bounds for Q m can be observed. See Fig. 2.
We also analyzed the shadow cast behaviour for a static observer and one that is comoving with cosmic expansion for a different model of Universes (i.e., dark energy, matter, and radiation-dominated Universes). With the parameters chosen from the constraints, we saw that the shadow radius decreases due to Q m > 0, and such an effect is amplified as γ decreases. As r obs → ∞, the rate of change in the shadow radius is zero. However, when the observer is comoving with the cosmic expansion, such change is no longer zero. See Fig. 4. While empirical data favours that our Universe is a dark energy-dominated one, it is interesting how we can also confirm this by using the deviation observed from the shadow cast. According to Fig. 4, such a deviation is amplified if one observes the shadow of a black hole near the cosmological horizon.
Next in our analysis is how the Modmax parameters are manifested through the deviations on the weak deflection angle, which is so sensitive even for large b/M . Our results revealed that the effect is amplified when u → 0. Furthermore, it was shown that relativistic massive particles give a larger deflection angle than photons. As γ decreases, so doesα, which is also manifested in Einstein ring formation. With the chosen parameters, these deviations can be detected by current sophisticated experiments. For instance, the EHT is capable of achieving an angular resolution of 10µas, which is used in mapping the stellar distributions near Sgr. A*. Furthermore, the ESA's GAIA mission can provide a sensitivity from 20µas to 7µas, which still depends on certain stellar magnitudes, used to map the Milky Way galaxy [218]. The smaller the angular resolution, the more the device can probe deflection at larger impact parameters. Note that the eikonal limit is independent of the spin of the perturbation, so that black holes' scalar, electromagnetic, and gravitational perturbations give the same behaviour in the eikonal limit [201]. Table III shows that the real parts decrease, but the imaginary part of the QNMs increase with the increasing parameter γ. We can conclude that these modes are stable cause the imaginary parts of the QNMs frequencies are negative. The decay rates (imaginary part) of QNMs frequencies increase with the increase of the screening parameter γ.
In Sect. V, the spherically free-falling accretion model on the ModMax BH from infinity has been investigated to provide realistic visualization of the shadow cast with the accretion disk. To do so, we first calculate the flux numerically to illustrate the effects of the parameters of the ModMax γ on the specific intensity seen by a distant observer for an infalling accretion in Figs. (7, 8 and 9) where show the specific intensities for various values of the parameter γ versus b observed by the distant observer. We have concluded that increasing the value of b, increases the intensity first. Afterwards, intensity reaches the peak value sharply, where the photons are captured by a black hole quickly. There is a peak value, and then intensity slowly decreases. Furthermore, it is shown that the shadow cast of the black hole in the 2-dimensional image with a photon sphere by a distant observer in (X, Y) plane where the dark centre of it the event horizon is located, and it is circled by a bright ring with a strongest luminosity (photon sphere). It can be seen that brightness decreases gradually after the maximum region. As a result, we have concluded that the effect of the screening parameter γ on the black hole luminosity of the shadow cast where the intensity decreases with increasing the value of the screening parameter γ as seen in Fig. 9.
In Sect. VI, we calculate the quasinormal modes (QNMs) of the Dyonic ModMax BHs by applying the method of WKB approximation. To study QNMs of the ModMax BHs, first, we use a massless scalar field perturbation in the background of the black hole 81. The dependence of the QNM frequencies on the screening parameter γ is qualitatively different for lowest and higher multipoles, as seen in Table I. We can see in Table I and II that both real and imaginary parts of the ω decrease when the parameter γ is increased. The scalar field QNMs for different values of l is plotted in 10, as well as for l = 0, 100 for different values of γ are presented in Figs. 11, 12 and 13. Then, we study the eikonal quasinormal modes for dyonic ModMax black holes. This method is also known as the geometric optics method due to its relation with the parameters of the null geodesics. Furthermore, the time-domain profiles of the scalar field perturbation for l = 0 and l = 2 cases have been plotted numerically in Figs. 14 and 15. The ModMax black hole is given in blue; on the other hand, the Schwarzschild case is black. The logarithmic plot shows that the ModMax frequency is slightly lower. For the l = 0 case, the plot shows the relatively short period of quasinormal ringing, making it hard to extract values of frequencies with good accuracy. After calculating the fundamental frequencies, we have provided constraints for the allowed range of masses of oscillating ModMax black holes which LIGO and LISA can detect.
Moreover, in Sect. VII, first we investigated the greybody factor of the ModMax black hole by using the rigorous bound. We show the effect of the screening parameter γ on the greybody bound for the scalar field of the ModMax black hole in Fig. 16 where it is seen that when the value of γ parameter increases, the greybody bound T b grows larger. Second, the high-energy absorption cross-section has been studied by applying the Sinc approximation for the ModMax black hole. We have shown in Fig. 17 that the total absorption cross section for various values of γ where shows that the greater values of screening parameter γ, the total absorption cross section becomes more and more unstable for higher energies and exhibits the largest amplitude. Furthermore, there is a regular oscillatory behaviour around the high-frequency limit. Fig. 17 presents the absorption spectrum as a function of the frequency, making clear oscillations characteristic of a diffraction pattern, where the oscillations occur around its constant value of the geometrical optics with decreasing amplitude and constant period.
We have analyzed the neutrino oscillations and spin-flip phenomena in the dyonic ModMax black holes field. This analysis can be relevant in relation to the recent observations of the event horizon silhouette of a supermassive BH [31,32]. The accretion disk around BHs is a source of photons and neutrinos [219]. The neutrino flavour oscillations and spin oscillations may affect the neutrino flux expected in a neutrino telescope. We have seen that neutrino oscillation lengths get modified of ∼ 10% in dyonic ModMax geometry, as compared to GR spacetime. Spin flip oscillations, instead, are suppressed, as in GR. This is a consequence of the fact that one has to consider relativistic neutrinos. Moreover, we have only considered the effects of the gravitational field, while a complete analysis should include the magnetic field [216]. We have also discussed the influence of the dyonic ModMax geometry on the nucleosynthesis processes, getting a difference ∼ 5%, with respect to GR, which could be relevant in the Supernovae evolution. We conclude by noting that a follow-up of these studies will be: the propagation of neutrino in both gravitational and electromagnetic fields; the coupling of neutrino magnetic momentum in non-linear electrodynamics; the neutrino spin oscillations (in gravitational fields) can be potentially observed in core-collapsing SN (in our Galaxy) since a huge amount of neutrinos are emitted. These effects are expected to be probed with future neutrino telescopes. All these topics will be faced in future works.