Constraints on Quasinormal modes from Black Hole Shadows in regular non-minimal Einstein Yang-Mills Gravity

This work deals with the scalar quasinormal modes using higher order WKB method and black hole shadow in non-minimal Einstein Yang-Mills theory. To validate the results of quasinormal modes, time domain profiles are also investigated. We found that with an increase in the magnetic charge of the black hole, the ring-down gravitational wave increases non-linearly and damping rate decreases non-linearly. The presence of magnetic charge also results in a decrease in the black hole shadow non-linearly. It is found that for large values of the coupling parameter, the black hole changes to a solitonic solution and the corresponding ring-down gravitational wave frequency increases slowly with a decrease in the damping rate. For the solitonic solutions, the shadow is also smaller. The constraints on the model parameters calculated using shadow observations of M87* and Sgr A* and an approximate analytic relation between quasinormal modes and shadow at the eikonal limit is discussed.


I. INTRODUCTION
Black holes are one of the most fundamental objects in the universe.They originate from collapsing stars and are governed by Einstein's theory of general relativity (GR).Regardless of complex interior structure of dying stars, once formed, black holes become one of the simplest celestial objects.According to the no-hair theorem, black holes are uniquely described by three physical parameters i.e, mass, angular momentum and electric charge.However, there are several scenarios where the no-hair theorem can be evaded.For instance, in the Einstein-Yang-Mills (EYM) theory and its variants, several numerical black holes with Yang-Mills hair are studied [1].When asymptotic structure of spacetime is modified to boxlike boundary, various hairy black holes are found [2].
Moreover, there has been much interest in gravitational theories coupled with nonlinear electrodynamics sources.Such marriage leads to a novel type of black hole with event horizon but possesses no essential singularity or regular black hole [3].Regular black holes are also studied in modified theories of gravity, for instance, f (R) gravity [4,5], Gauss-Bonnet gravity [5][6][7].An exact spherically symmetric Wu-Yang monopole solution is derived and studied [8][9][10].This is a non-minimally extended of Einstein-Yang-Mills equation with SU (2) Wu-Yang ansatz.Later, this regular magnetic black hole is generalized to include a cosmological constant [11].
As predicted by GR, any accelerated massive objects create disturbances in spacetime around them.The disturbances propagate through spacetime in the form of gravitational waves (GWs).These waves contain crucial information about their astrophysical sources.When two black holes collide and merge into a single massive black hole, during this process GW is emitted continuously.After the merger, a single black hole emits GWs signal which is prominently dominated by a ringdown mode.The black hole ringdown is generally described by quasinormal modes (QNMs) [12][13][14][15].An associated complex frequency is determined by black hole's mass and angular momentum.The real component of the frequency is the emission frequency while the imaginary component denotes exponential decay.Thus, studying QNMs of black holes is crucial since they allow us to attain various characteristics of black holes.
Remarkably, the detection of GWs [16] in 2015 has opened a new window to testing theories of gravity.Vast number of works have been devoted to study the QNMs of various black holes.QNMs of nonminimally coupled scalar field in three dimensional spacetime are explored in [17].In addition, QNMs of nonminimally coupled scalar field in a five dimensional Einstein-Power-Maxwell background are studied via the Wentzel-Kramers-Brillouin (WKB) method and pseudospectral Chebyshev method [18].Optical properties and QNMs of black holes dependence on the Yang-Mill charges are investigated in Einstein-Power-Yang-Mills theory [19].QNMs of black holes within the framework of generalized uncertainty principle (GUP) are analyzed [20][21][22][23].Due to variety of black hole solutions in modified theories of gravity, there are many works considering QNMs of black holes in various models e.g., f (Q) gravity [24], f (R, T ) gravity [25], Horndeski gravity [26,27], Rastall gravity [28][29][30] and de Rham-Gabadadze-Tolley (dRGT) massive gravity [31][32][33][34][35]. Very recently, QNMs of black hole in emergent gravity framework are computed via the WKB and asymptotic iteration methods [36].
In recent times, the visual representation of black holes has garnered significant attention, driven by observations of celestial objects like M87 ⋆ [37][38][39][40][41][42] and Sgr A ⋆ [43].These observations have sparked interest in the realms of both GR arXiv:2402.06186v1[gr-qc] 9 Feb 2024 and modified gravity (MOG), offering an intriguing opportunity to depict and analyze hypothetical images associated with these observed phenomena.This avenue allows us to engage in a geometric and topological exploration, comparing the theoretical depiction of black hole images with the actual images of M87 ⋆ and Sgr A ⋆ .Consequently, the visualization of black holes assumes a captivating optical dimension, particularly as we endeavour to uncover the relationships between optical features, such as the shadows they cast, through the application of thermodynamics [44,45] and quasi-normal modes [46,47].Additionally, when considering the shape data, the shadow images can take the form of either distorted geometries on the image plane [48][49][50][51][52] or a pattern of concentric circles, contingent on the presence of a rotation parameter within the spacetime.To gain a deeper understanding of these shadow behaviours, the Hamilton-Jacobi formalism [53] emerges as a compatible analytical framework.It is based on the concept that massless photons in the vicinity of a black hole generate specific orbits around the region, a phenomenon referred to as the geodesic null background.In the realm of black hole imaging, ongoing research aims to provide a comparative analysis of observed images, such as those of M87 ⋆ and Sgr A ⋆ , in terms of their size and shape.
Motivated by these studies, in this work, we aim to investigate the QNMs and black hole shadows in the framework of regular non-minimally EYM theory in four-dimensional spacetime.We also introduce constraints on the model parameters from black hole shadow observations and try to relate the shadow behaviour with the QNMs of the black hole spacetime.
This paper is organized as follows.In Section II, we discuss the EYM theory and the associated black hole solution in brief.In Section III, massless scalar perturbations on the black hole spacetime have been investigated and the QNMs spectrum is analysed.Section IV deals with the evolution of scalar perturbation.In Section V, we investigate the black hole shadow.In Section VI, we discuss how QNMs are linked with shadows of the black hole.Finally in Section VII, we provide a concluding remark of our work.

II. REGULAR NON-MINIMAL EINSTEIN YANG-MILLS BLACK HOLE
In this section, we shall discuss the theory and the black hole solution in brief.The associated action of the 4-D regular non-minimally EYM theory is [8,11] where R stands for the Ricci curvature scalar.The indices labeled with Greek letters range from 0 to 3, while those labeled with Latin letters range from 1 to 3. Additionally, the notation F where ∇µ represents the covariant derivative and the symbols f (b) (c) (a) denote the real structure constants of the threeparameter YM gauge group SU (2).The tensor Rαβµν is defined as [11] Here, R αβ and R αβµν represent the Ricci and Riemann tensors, respectively, and ξ i (i = 1, 2, 3) denote the non-minimal coupling parameters between the YM field and the gravitational field.Assuming the gauge field is described by the Wu-Yang ansatz and taking ξ 1 = ξ, ξ 2 = 4ξ, ξ 3 = −6ξ with ξ > 0, a regular, static, and spherically symmetric black hole solution was discovered [8,11].This black hole solution is described by the metric [8,11], where the metric function has the following form: with M representing the mass of the black hole and Q being the magnetic charge.
The electromagnetic field four-potential for the regular nonminimal magnetic black hole is given by: It's evident that when Q = 0, the metric (4) simplifies to the standard Schwarzschild black hole metric, while when ξ = 0, it corresponds to the Reissner-Nordström black hole metric but with a magnetic charge instead of an electric charge.The quantity of event horizons of the non-minimal EYM black hole depends on both the non-minimal parameter and the magnetic charge, potentially resulting in multiple horizons.In this work, we shall use the black hole line element (4) to study the massless scalar QNMs and the black hole shadow.

III. MASSLESS SCALAR QUASINORMAL MODES
In this section, we shall investigate the QNMs associated with massless scalar perturbation.In this case, we shall assume that the test field has negligible impact or influence on the black hole spacetime defined in the previous section and the associated scalar field is massless in nature.We obtain a Klein-Gordon-type equation associated with the QNMs while considering the pertinent conservation relations associated with the spacetime under consideration.We shall implement Padé averaged 6th order WKB approximation method to obtain the QNMs numerically.
Taking into consideration of axial perturbation only, it is possible to express the perturbed metric in the following form [24,54]: In this context, the parameters p 1 , p 2 , and p 3 are intricately linked to the perturbation occurring in the spacetime surrounding the black hole.These parameters play a crucial role in characterizing the modifications induced by the perturbation.Specifically, they contribute to the altered dynamics and structure of the black hole spacetime.Now, for the massless scalar field, as it is considered that the effect of the field on the spacetime is minimal, the perturbed metric Eq. ( 7) takes the following form: The Klein-Gordon equation in curved spacetime in this case is written as Eq. ( 9) describes the QNMs associated with the scalar perturbation in the considered black hole spacetime.
Utilising spherical harmonics, it is possible to express the wave function of the scalar field in the following decomposed form: In the above equation, Y lm (θ, ϕ) is the spherical harmonics part of the wave function and the parameters l and m are the associated parameters of spherical harmonics.Now, one can use this decomposed form of the wave function in the Eq.(10) to obtained a reduced wave equation given by This is a standard differential equation where the differentiation is done with respect to the tortoise coordinate r * .The tortoise coordinate is connected with the normal coordinate r by the following definition: In the Eq. ( 11), parameter V s (r) is the effective potential for the massless scalar perturbation and it is given by: In this context, the symbol l is employed to signify the multipole moment associated with the QNMs of the black hole.
The variation of the scalar potential V s (r) with respect to the multipole moment l is shown in Fig. 1.The overall effective potential increases with ℓ.
The other model parameters Q and ξ also have noticeable impacts on the potential V s (r) which are depicted in Fig. 2. It is seen that with an increase in the black hole charge Q, the peak value of the potential increases non-linearly and the corresponding r shifts towards the black hole event horizon.Similarly, an increase in the value of ξ also increases the peak value of the potential.However, the variation is comparatively smaller at the peak position and the peak value seems to change uniformly with respect to ξ.This suggests that the model parameter Q might have a non-linear impact on the spectrum of QNMs and ξ might have an almost linear impact on QNMs.The details will be investigated in the next part of the study.

A. The Padé averaged WKB approximation method
In this study, we implement a well-known method called WKB approximation method to calculate the QNMs numerically for the massless scalar perturbation in the abovementioned black hole spacetime.We consider Padé averaged 6th-order corrections to the WKB method [55][56][57][58] to obtain QNMs with higher accuracy.
In higher-order WKB approximation method, the oscillation frequency ω of ring-down GWs is expressed as: Here, the parameter n, representing the overtone number, takes on discrete values starting from zero (n = 0, 1, 2, . ..).
Additionally, V 0 represents the value of the potential function V (r) at a specific radial coordinate r max , while V ′′ 0 signifies FIG.2: Variation scalar potential with respect to r using M = 1 and l = 3.On the first panel, we have used ξ = 0.1 and on the second Q = 0.7.For ξ = 0.9 and ξ = 1.3, these configurations are not black holes but represent solitonic solutions.
the second derivative of V (r) with respect to r evaluated at the same position.The pivotal location r max corresponds to where the potential function achieves its maximum value, rendering this radial coordinate particularly significant.
Integral to our computational framework is the correction terms Λk , which play a crucial role in refining the precision of our calculations.These correction terms are instrumental in capturing subtle nuances in the determination of the oscillation frequency, contributing to the overall accuracy of the framework.
To elaborate on the procedural aspects of our approach, we draw upon established references in the field.Specifically, Ref.s [55][56][57][58] provide comprehensive insights into the mathematical formulations of the correction terms and elaborate on the intricacies of the Padé averaging method applied in our study.
By incorporating these computational techniques, our research aims to uncover the nuanced behaviour of ring-down GWs, making significant contributions to the broader understanding of astrophysical phenomena in such types of theories of gravity.Within this tabular representation, the term ∆ rms assumes significance as it quantifies the root mean square error asso-ciated with the 6th-order Padé averaged WKB approximation.
Moreover, the term ∆ 6 is introduced to gauge the error between two adjacent order approximations.This error term is computed by assessing the absolute difference between ω 7 , representing QNMs computed using the Padé averaged 7th order WKB approximation method, and ω 5 , denoting QNMs obtained from the Padé averaged 5th order WKB approxima-tion method.Mathematically, ∆ 6 is expressed as follows: A notable trend within the table is the observed reduction in error associated with QNMs as the multipole moment l increases.This behaviour aligns with the characteristic limitations of the WKB approximation method, particularly evident when the overtone number n surpasses the multipole moment l.The diminishing error with increasing multipole moment is a distinctive feature of the WKB method, signifying its challenges in providing accurate results when overtone numbers exceed the corresponding multipole moment values [21,23,24,[59][60][61][62].By keeping this fact in mind, in this investigation, we have considered n < l to calculate QNMs.
We have shown the variation of the massless scalar QNMs in Fig. 3 and 4 with respect to Q and ξ, respectively.In Fig. 3, one can see that with an increase in the charge of the black hole Q, the real QNMs increase non-linearly.On the other hand, the damping rate or decay rate of ring-down GWs increases very slowly initially up to around Q = 0.6 and beyond this the damping rate decreases drastically showing a highly non-linear pattern.
In Fig. 4, the initial panel intricately illustrates the nuanced shift in real quasinormal frequencies concerning the parameter ξ, while the second panel meticulously delineates the corresponding alteration in the damping rate of ring-down GWs with respect to ξ. Evidently, as the parameter ξ experiences an incremental surge, the real quasinormal frequencies manifest a discernible linear escalation.In contrast, a rise in the value of ξ induces a linear reduction in the damping rate.It is noteworthy, however, that this variation is comparatively modest when juxtaposed with the scenario involving the parameter Q.
Consequently, when contemplated through the prism of QNMs, the influences exerted by Q and ξ exhibit notable disparities.Looking forward, as we anticipate a wealth of data from state-of-the-art GW detectors such as Laser Interferometer Space Antenna (LISA), there emerges the exciting prospect of establishing stringent constraints on the Yang-Mills field by harnessing insights derived from the meticulous examination of QNMs.

IV. TIME EVOLUTION OF SCALAR PERTURBATIONS
In the previous section, we conducted numerical computations to unravel the characteristics of QNMs and scrutinized how they hinge on the model parameters.Shifting our attention to the subsequent section, our exploration now centers on the temporal behaviours of massless scalar perturbations.To elucidate these profiles, we adopt the time domain integration framework delineated by Gundlach et al. [63,64] To progress further, we denote ψ(r * , t) = ψ(i∆r * , j∆t) = ψ i,j , and V (r(r * )) = V i .Here r * represents the tortoise coordinate and t denotes time.Subsequently, we express the scalar wave equation ( 9) as: The initial conditions are specified as ψ(r and ψ(r * , t)| t<0 = 0, where k 1 and σ represent the median and width of the initial wave-packet, respectively.Employing these initial conditions, we compute the time evolution associated with the scalar perturbation using the iterative scheme: We implement this iterative scheme to obtain the time domain profiles by keeping a fixed value of ∆t ∆r * .However, it is important to mention that one must consider ∆t ∆r * < 1 to satisfy the Von Neumann stability condition throughout the numerical procedure.Throughout our numerical calculations, we have considered this ratio around 0.7 to obtain the time evolution of the massless scalar perturbations on the black hole spacetime.We have set k 1 = 0 and σ = 5 for the numerical calculations of time domain profiles.
In Fig. 5, we have shown the time domain profiles in the units of M .One may observe that the variation of time domain profiles with respect to multiple moment l stands in agreement with our Table I.From the second and third panel, it is evident that the impact of the coupling parameter ξ on the QNMs spectrum is very small in comparison to the charge parameter Q as predicted earlier.These results are also in agreement with the results obtained from Padé averaged WKB approximation method.

V. OPTICAL BEHAVIOUR OF THE BLACK HOLE: SHADOW
The investigation into supermassive black holes (SMBHs) holds significant importance in astrophysics, given their mysterious nature and their central positions in galaxies, where they exert exceptional gravitational forces with event horizons preventing light escape.Recent advances in observational techniques, exemplified by the Event Horizon Tele-   In this part, we explore the mysterious concept of black hole shadows, uncovering the hidden areas within the complex structure of space and time that cover the edge of these cosmic objects.Black holes have incredibly strong gravity that forces even light to be pulled in when it crosses the boundary called the event horizon [65,66].Consequently, the outer boundary of a black hole manifests as an indistinct yet solemn circle, projecting its inscrutable silhouette against the surrounding cosmic tableau of matter or luminous radiance.The dimensions and intricate contours of this umbral figure furnish profound insights into the fundamental characteristics of black holes and the essence of gravity itself.Recent advancements in both astronomical observation and technological capabilities have conferred upon us the capacity to capture and visually scrutinize these elusive black hole shadows, representing a paradigmatic leap in our pursuit to comprehend these enigmatic celestial phenomena [67].As we peer into these cosmic shadows, we not only unveil the visual facets of black holes but also unlock new strata of scientific comprehension, shedding light on the intricacies veiled within the expansive cosmic tapestry.These advancements herald a gateway to a deeper understanding of the universe's intricacies, underscoring the remarkable synergy between cutting-edge technology and our relentless scientific inquisitiveness about the cosmos.Such a study will extensively contribute towards our understanding of different modified gravity theories and the validity of GR in light of current observational aspects.
The Euler-Lagrange equation is given by where the Lagrangian is expressed as For static and spherically symmetric black hole, the Lagrangian becomes Here the dot over the variables denotes the derivative with respect to the proper time τ .Choosing the equatorial plane where θ = π/2, the conserved energy E and angular momentum L can be computed by utilizing the killing vectors ∂/∂τ and ∂/∂ϕ as given by In the case of photon, we can write the geodesic equation as, Utilizing Eq. ( 22) in conjunction with the conserved quantities, namely, E and L, the orbital equation for a photon can be derived as.[67] dr dϕ where V ef f is defined as Expressing Eq. ( 23) in radial form yields, In this expression, ζ represents the impact parameter and is defined as ζ = L/E.Additionally, the reduced potential V r (r) can be expressed as: To scrutinize the shadows of black holes, our focus is directed towards a specific point along the trajectory denoted as r ph .This particular point corresponds to the turning point, signifying the location of the light ring that encircles the black hole.It can also be interpreted as the radius of the photon sphere, holding significant relevance in the analysis of the black hole's shadow characteristics.At this pivotal turning point of a black hole [65,[67][68][69][70][71][72], It is feasible to determine the impact parameter ζ at the turning point by Therefore, the radius of the photon sphere, denoted as r ph , can be determined using [67]: This equation can be reformulated as here we have A(r) = h(r)/f (r) with h(r) = r 2 .Now, to derive the shadow of the black hole, we express Eq. ( 23) using Eq. ( 28) in terms of the function A(r) as Utilizing Eq. ( 31), one can calculate the shadow radius of the black hole.In the scenario where a static observer is positioned at a distance r 0 from the black hole, we can ascertain the angle α between the light rays originating from the observer and the radial direction of the photon sphere.This angle can be computed as [67] cot α = 1 In conjunction with Eq. ( 31), the aforementioned equation can be represented as Once more, the above equation can be reformulated using the relation sin 2 α = 1/(1 + cot 2 α) as By substituting the expression for A(r ph ) from Eq. ( 28) and using A(r 0 ) = r 2 0 /f (r 0 ), the shadow radius of the black hole for a static observer at r 0 can be approximated as [67] Once more, as r 0 → ∞, i.e., for a static observer at a large distance, f (r 0 ) → 1.Consequently, for such an observer, the shadow radius R s becomes, Finally, through the stereographic projection of the shadow from the black hole's plane to the observer's image plane with coordinates (X, Y ), the apparent form of the shadow can be determined.These coordinates are defined as [22,67,73] In this context, θ 0 denotes the angular orientation of the observer in relation to the plane of the black hole.The investigation presents stereographic projections of the black hole shadow for the specific black hole under examination.Fig. 6 illustrates these projections, offering visual depictions of the shadow's characteristics corresponding to various model parameters.In the initial panel of Fig. 6, the stereographic projections of the black hole are displayed for different values of the model parameters Q and ξ.It is evident that an increase in the model parameter Q results in a decrease in the black hole shadow.The second panel illustrates the black hole shadow for various values of the parameter ξ, revealing a gradual decrease in the shadow with an increase in ξ.For enhanced visualization, Fig. 7 presents the square of the black hole shadow radius concerning the model parameters Q and ξ.This representation indicates that ξ has a minimal impact on the black hole shadow.Therefore, the investigation demonstrates that both model parameters exert distinct influences on the black hole shadow.
A notable feature of SMBHs is the presence of a shadow cast by their event horizons, a phenomenon arising from the gravitational bending of light.The size and shape of this shadow furnish valuable information about the SMBH and its surrounding environment.To comprehend and quantify the shadow's dimensions, the spatial separation D between the SMBH and the galactic center must be considered.
A conventional method for determining the classical diameter of the black hole shadow involves the application of the arclength equation: where d sh denotes the shadow's diameter, θ sh represents the angular size of the shadow, and M corresponds to a specific unit of measurement.
For M87*, one of the extensively studied SMBHs, the measurement of the shadow diameter has been realized through EHT observations [37,43].The reported diameter of the M87* shadow is d M87* sh = (11 ± 1.5)M , serving as a crucial data point for scrutinizing the properties of this SMBH.
However, it is imperative to consider the uncertainties associated with these measurements for a more accurate comprehension of the physical parameters.The determination of such uncertainties requires meticulous analysis, and in the case of the M87* shadow diameter, Ref.s [74,75] offer detailed insights into the employed methodology.
This study aims to delve deeper into the implications of the measured M87* and Sgr A* shadow diameters by incorporating uncertainties from the aforementioned references.Our specific focus is on constraining the potential values of the black hole charge Q by calculating the 1σ limits based on results obtained for M87* and Sgr A*.By considering these constraints, we aspire to augment our understanding of the model behaviour and the model parameters from the physical properties of M87* and Sgr A* and contribute to the broader comprehension of SMBHs as well as possible constraints on QNMs from them.In this section, we shall briefly discuss the approximate analytical connection of QNMs with the black hole shadow radius.One may note that from the numerical results of QNMs and black hole shadow, there is a possible correspondence between them.
In this scenario, we shall consider only 3 rd order WKB expansion given by One may note that even we consider higher order expansion such as 6th order or 12th order expansion, in the eikonal limit, the first significant terms are identical with those obtained by the 3rd order WKB approximation method.The term V i stands for the i-th derivative of the potential V .We also considered ν = n + 1 2 , where n is the overtone number of QNMs.By expanding (40), at the eikonal regime, one can obtain [19,76] Comparing the imaginary part and the real part Further exploring the constraints on the visual representation of a black hole i.e. shadow can be achieved through the examination of observational data concerning its QNMs.It is essential to emphasize the direct relationship between the features of the black hole's silhouette and the specific values of its QNMs as this will help us to uncover different properties of the black hole hairs.Consequently, the availability of significant observational data from space-based gravitational wave detectors like LISA presents an avenue for imposing stricter constraints on the parameters of models used to understand black holes.Moreover, this process enables the assessment of consistency between two distinct approaches-shadow analysis and QNMs-in testing the underlying theoretical framework.
From the approximated relations ( 42) and (43), at the eikonal limit, we observe that the black hole shadow has an impact on the behaviour of the QNMs spectrum.Most importantly, an observational constraint on the black hole shadow will also put a constraint on the real and imaginary parts of QNMs of the black hole at the eikonal limit.However, as evident from the table II, such a constraint on the QNMs will be very weak at the current stage.In the near future, observational results from LISA might put a very strong constraint on QNMs which may be useful along with the observational results of black hole shadow to check for the consistencies between them and to constrain the theory more stringently.

VII. CONCLUDING REMARKS
In this work, we considered a charged black hole solution in non-minimally coupled EYM theory and studied the scalar perturbation on the black hole spacetime and the associated QNMs.From the behaviour of the potential, we found that for higher values of the coupling parameter ξ, the line element represents solitonic solutions instead of a black hole solution.
In order to obtain the QNMs, we used Padé averaged WKB approximation method.The analysis showed that black hole charge Q has a more significant impact on the quasinormal mode spectrum.With an increase in charge, the oscillation frequency of ring-down GWs increases non-linearly.In the case of the damping rate of GWs, we found that initially, with an increase in Q, the damping or decay rate increases slowly and beyond some threshold value around Q = 0.65, it starts to decrease rapidly as Q approaches 1.One may note that some previous Ref.s [19,29] also found similar variation trends of QNMs with respect to black hole charge Q.However, due to the presence of coupling term and non-linear charge distribution, in this case, we observe slightly different behaviour of QNMs towards higher values of Q.
Unlike black hole charge Q, the coupling parameter ξ has a linear effect on the QNMs spectrum.With an increase in the value of ξ, the oscillation frequency of QNMs increases while the damping rate decreases slowly.
We also investigated time domain profiles for different model parameters and found that the profiles depict similar results as those found by using the WKB method.
In the next part of our work, we considered the shadow of the black hole and found that Q has a noticeable impact on the shadow behaviour of the black hole.For a charged black hole, the shadow is expected to be smaller in size.A similar behaviour is observed in the case of ξ also.However, the variation of shadow with respect to ξ is small.Thus the impact of the coupling parameter as well as the presence of solitonic solutions may not be well differentiated by utilising shadow observations.As a result, observational data on black hole shadows, do not have a stringent constraint on ξ.
We also studied the relationship between the QNMs and the shadow of the black hole and found that the shadow observations can provide a constraint on the QNMs spectrum.In the near future, observational data of QNMs from LISA may provide us with scope to constrain the theory more efficiently and to check whether the shadow observations are well consistent with it or not.
Our study has uncovered that in the standard non-minimally coupled EYM theory, both the charge parameter Q and the coupling term ξ exert notable influences on the quasinormal spectrum and shadow of the black hole.The effects of these parameters on the quasinormal mode spectrum differ significantly in their nature.
µν represents the Yang-Mills (YM) tensor and is connected to the YM potential A (a) µ through the equation:

FIG. 5 :
FIG. 5: Time domain profiles of scalar perturbation with M = 1.On the first panel, Q = 0.3 and ξ = 0.5; on the second panel l = 1 and Q = 0.3 and on the third panel l = 1 and ξ = 0.5.

TABLE II :
Variation of shadow with 1σ ranges of Q and ξ based on the shadow radii of Sgr.A* and M87* as depicted in Figs.8 and 9.
VI. CONNECTION BETWEEN QUASINORMAL MODES AND SHADOW OF THE BLACK HOLE