Quasinormal Modes in Noncommutative Schwarzschild Black Holes: A Spectral Analysis

We present a comprehensive analysis of quasinormal modes (QNMs) for noncommutative geometry-inspired Schwarzschild black holes, encompassing both non-extreme and extreme cases. By employing a spectral method, we calculate the QNMs in the context of scalar, electromagnetic, and gravitational perturbations. Our findings not only challenge previous claims in the literature regarding the instability of these black holes but also reveal remarkable stability for both non-extreme and extreme Schwarzschild black holes under various perturbations.


I. INTRODUCTION
What exactly is a quasinormal mode (QNM)?In essence, we can think of it as the real-world counterpart to the idealized normal mode.To illustrate this concept, let us consider the example of a gently perturbed wine glass.When we tap the glass lightly, it produces a sound that gradually fades away as energy dissipates into the surroundings.If this were a normal mode, the ringing would persist indefinitely.In our current study, we venture into a realm of objects far more intriguing than an ordinary wine glass, namely, black holes.When we delve into more complex scenarios, accounting for possible dissipative phenomena, it becomes natural to expect that an expansion in normal modes is no longer feasible.In such cases, it will be quite natural for the solutions to exhibit decay over time instead of oscillating indefinitely.For this reason, the frequencies will be represented by complex numbers.Let us consider the possibility that the nature of the equilibrium point of a physical system around which we perform an expansion is not known a priori.The QNMs we determine may be characterized by frequencies whose imaginary part does not describe the damping of solutions around the aforementioned point but instead signifies an exponential growth over time.It is easy to deduce that, in this case, the equilibrium point is, in fact, unstable.The significance of this concept is profound: often, the expansion in QNMs is employed to ascertain the stability of highly complex systems, and it was precisely this idea that prompted Regge and Wheeler in their seminal paper [1] to initiate the study of black holes within the perturbative framework.
Over the years, various techniques have been developed to study the QNMs, each with its own peculiarities.One of the most important methods for determining them is the WKB (Wentzel-Kramers-Brillouin) method [2].This particular approach extensively exploits the analogy between the physical phenomenon we are studying here and the study of energy eigenvalues for wave functions using the Schrödinger equation in Quantum Mechanics.Another technique useful for the determination of QNMs makes use of group theory.The proposed technique, known as the inverted potential method [3], is specifically applicable only in the case of a Schwarzschild black hole, not a Kerr black hole.It is based on the similarity between the effective and reversed Pöschl-Teller potentials.Furthermore, the continued fractions method analysis was first introduced by [4,5] while the Laplace transform method was exploited in [6].In addition to the aforementioned techniques, it is worth mentioning that the Asymptotic Iterative Method (AIM) has also been applied to determine the QNMs for Schwarzschild, Reissner-Nordström, and Kerr black holes in the case of spin 0, 1/2, and 2 perturbations [7][8][9][10][11][12][13]. Recently, in a study by [14], QNMs of a Schwarzschild black hole were revisited using an innovative approach called the spectral method.Spectral methods, though less common in the general relativity (GR) quasinormal mode (QNM) community, offer substantial advantages both analytically and numerically.These methods harness spectral expansions derived from constructive approximation theory.The coefficients within these expansions are computed numerically, resulting in an approximate analytical representation of the solution, effectively blending the strengths of both analytical and numerical approaches.
In contrast to widely used methods such as WKB, inverse potential, and AIM methods, spectral techniques do not rely on local approximations of the potential near its maximum nor require constraints on the angular momentum parameter.The WKB technique also struggles with potentials whose second derivative at the maximum is zero, thus leading to divergent behaviours in formulae like (4.1) in [15].Finally, for a thorough analysis of the validity of the WKB approximation in calculating the asymptotic quasinormal modes of black holes we refer to [16].Asymptotic iteration methods, while useful, are highly sensitive to the choice of initial values in the parameter space.An incorrect choice during the initialization process can affect the stability of the method and the number of iterations required to compute eigenvalues or quasinormal modes.A recent attempt to find the optimal selection can be found in [17].
On the other hand, spectral methods utilize the full analytical form of the potential without any need for local approximations or heuristic modifications.This comprehensive approach not only provides high accuracy but also ensures exponentially fast convergence for the regular part of the solution of a radial equation once its singular behaviour has been completely determined.This feature is particularly advantageous in the search of QNMs of astrophysical objects.For the aforementioned reasons, in the current paper, we will employ an improved version of the spectral technique to explore the QNMs of the Schwarzschild black hole inspired by noncommutative geometry in the case of a massless particle of spin s = 0 (scalar perturbation), s = 1 (electromagnetic perturbation), and s = 2 (vector-type gravitational perturbation) [18].
While a substantial body of literature exists dedicated to the investigation of QNMs in the context of the classic Schwarzschild spacetime [10,[19][20][21][22][23][24][25][26], the situation is considerably less comprehensive for the noncommutative Schwarzschild manifold.In the latter case, the literature contains only partial and occasionally contradictory findings, which we will briefly address in the following discussion.For example, [27] analyzed the asymptotic QNMs of the nonextreme black hole inspired by noncommutative geometry under the assumption that the black hole mass significantly exceeds its extreme mass.As previously highlighted by [28], in this regime, the QNMs already align with those of its classical counterpart.Hence, it comes as no surprise that [27] concluded that the asymptotic QNMs remain proportional to ln 3, thus mirroring the corresponding behaviour in the case of the classic Schwarzschild solution.Recently, [29,30] conducted a study on QNMs of massless scalar field perturbations in the nonextreme regime utilizing the third-order WKB approximation.The author specifically focused on cases with ℓ ∈ {1, 2, 3}.In this work, the quasinormal frequencies have been computed for values of the noncommutative parameter θ across the interval [0.01, 0.2758].If we recall that √ θ represents a quantum of length, it is reasonable to assume that √ θ is on the order of the Planck length, L P .Consequently, the interval chosen in [29,30] corresponds to a selection of black hole masses within the range 3.6M P < M < 19M P , where M P represents the Planck mass.In this specific range and employing the third-order WKB approximation, [29,30] observed stable results for the QNMs.It is noteworthy that, in contrast, [28] identified instabilities using WKB approximations from third sixth order not only in the narrower range of 1.9M P < M < 2.3897M P but also in the extreme case µ = 1.9MP .[28] emphasized that extending the WKB calculations to the sixth order reveals a lack of convergence precisely when the presumed instabilities occur.To gain insight into the nature of these instabilities, [28] resorted to alternative methods.The inverted potential and Asymptotic Iterative Methods (AIM) failed to detect any instability.However, it is worth noting that both of these methods entail certain approximations, and the AIM method, in particular, is known to encounter convergence issues.For a more in-depth discussion of this aspect, we refer to [17].Last but certainly not least, [31] conducted an examination of the stability within the interior of the noncommutative Schwarzschild black hole.This analysis involved the study of massless scalar field propagation between the two horizons.Their findings demonstrated that spacetime fuzziness, influenced by higher momenta of the field, effectively mitigates the classical exponential blue shift divergence.This mitigation suppresses the emergence of infinite energy density in regions proximate to the Cauchy horizon.Conversely, [32] investigated the surface gravity of both inner and outer horizons and posited that the Cauchy horizon might become unstable when subjected to perturbations from infalling matter.In light of this ongoing discussion, there is a significant interest in investigating the (in)stabilities of nonextreme and extreme noncommutative geometry-inspired black holes by analyzing their respective QNMs.In the present study, we tackle this challenge by employing the spectral method to analyze quasinormal frequencies associated with scalar, electromagnetic, and vector-type gravitational perturbations.Our results unequivocally demonstrate the absence of instabilities in both the extreme and nonextreme cases.Notably, in the extreme case, where the event horizon essentially serves as a de facto Cauchy horizon, our findings align with the conclusions derived in [31], further bolstering their validity.
The paper is organised as follows.In Section II, we introduce a suitable rescaling of the noncommutative Schwarzschild manifold and derive the effective potential for a massless scalar field in the aforementioned geometry.Moreover, we also establish the corresponding QNMs boundary conditions of a massless scalar field in the presence of nonextreme and extreme noncommutative Schwarzschild black holes.Section III describes the numerical method used in our study, providing the basis for our computational approach.In Section IV, we validate our numerical method by applying it to large-mass noncommutative Schwarzschild black holes.This choice is strategic as, in this high-mass regime, these black holes converge towards the characteristics of classical Schwarzschild black holes, providing a robust benchmark for validation.Here, we also delve into the numerical findings for both nonextreme and extreme noncommutative black holes.Finally, Section V offers conclusions and an outlook, reflecting on the implications of our findings and suggesting directions for future research.

II. EQUATIONS OF MOTION
We consider a massless scalar field ψ immersed in the noncommutative geometry inspired Schwarzschild background whose line element in units where c = G N = 1 reads [18] with the mass function where M is the total mass of a gravitational object, θ is a parameter encoding noncommutativity and having the dimension of a length squared.At the same time, γ(•, •) is the incomplete lower gamma function.Considering the relationship between the upper and lower incomplete gamma functions as expressed by the formula [33] we can represent the g 00 metric coefficient in (1) as a combination of the typical Schwarzschild term augmented by a noncommutativity-induced perturbation, namely Note that when r/ √ θ → +∞, the term Γ(3/2, r 2 /4θ) tends to zero, yielding the classic Schwarzschild metric.If we introduce the following rescaling coupled with the identities [33] we can cast (4) in the form where erf(•) represents the error function.The picture emerging from Figure 1 where f has been plotted against x, is the following: an extreme noncommutative geometry inspired Schwarzschild black hole appears when the rescaled mass takes the critical value µ e = 1.904119076 . ... In this regime, the Cauchy horizon x c and the event horizon x h coincide, and their values correspond to x e = x c = x h = 0.7936575898 . . .Whenever µ < µ e , the line element (1) represents a naked gravitational droplet regular at x = 0 [18].Finally, two distinct horizons exist for µ > µ e .Differently as in the case of the Schwarzschild metric, if we probe into shorter distances, i.e. r ≪ √ θ, it is possible to show the absence of a curvature central singularity which is instead replaced by a regular deSitter core [18].In the present work, we focus our analysis on the QNMs for the nonextreme and extreme cases.
Considering this, it is advantageous to implement the rescaling z = x/x h , as this approach effectively maps the event horizon to 1 in both the extreme and nonextreme cases.Hence, the g 00 metric coefficient, as expressed by (7), becomes FIG. 1: Plot of the metric coefficient g00 = f (x) defined in (7).An extreme black hole occurs when µe = 1.904119076 . . .(solid line) and the horizon is located at xe = 0.7936575898 . . .A nonextreme black hole with two horizons is present for µ > µe (see the dotted line corresponding to the choice µ = 2.3).A naked gravitational droplet occurs when 0 < µ < µe (see the dashed line for µ = 1.3).
while in the presence of a noncommutative geometry-inspired Schwarzschild black hole, the equation governing a massless Klein-Gordon field (assumed to have a time dependence of the form e −iωt and an angular component described by spherical harmonics) is as follows [28] f (r) with ℓ = 0, 1, 2, . . .and ǫ = −1 (massless scalar perturbation s = 0), ǫ = 0 (electromagnetic perturbation s = 1), and ǫ = 3 (vector-type gravitational perturbation s = 2).A comprehensive analysis of the effective potential, along with its graphical representations for both configurations, has been extensively covered in [28].By means of the substitution r = 2M x h z, the above equation can be recast in the equivalent form with f (z) as given in (8).As a consistency check, note that in the limit corresponding to the classic Schwarzschild case, i.e. x h → 1 and r/ √ θ → ∞, the above equation correctly reproduces equation ( 2) in [4], where the term ω 3 should be corrected to ω 2 .In the following analysis, we focus on computing the QNMs for the spectral problem stated in (10).For this purpose, we represent Ω as Ω = Ω R + iΩ I , where Ω I < 0 ensures that the perturbation is dumped in time.The boundary conditions are set so that the radial field exhibits inward radiation at the event horizon and outward radiation at spatial infinity.This necessitates a detailed examination of the solution asymptotic behaviour in (10), both near the event horizon (z → 1 + ) and at large spatial distances (z → +∞).Moreover, to compute the QNMs using the spectral method, we must recast the differential equation in (10) and the appropriate boundary conditions over the compact interval [−1, 1].This adjustment is essential as the method expands the regular part of the eigenfunctions by means of Chebyshev polynomials.

A. The non-extreme case
This scenario focuses on mass parameters µ > µ e .The metric coefficient g 00 is given as defined in (8).It is important to note that the Cauchy and event horizons are distinct in the present case, and therefore, g 00 exhibits a simple zero at z = 1.To establish the QNM boundary conditions at the event horizon and at infinity, we first need to determine the asymptotic behaviour of the radial solution ψ Ωℓǫ as z → 1 + and as z → +∞.We can then extract the QNM boundary conditions from this asymptotic data.As a consistency check, we will also verify that in the limit r/ √ θ, these conditions correctly reproduce the corresponding ones for the classic Schwarzschild case specified in equation ( 4) of [4].We split our analysis by examining the behaviour of the radial field in two different regions.
1. Asymptotic behaviour as z → 1 + : Given that z = 1 is a simple zero of f (z), we can represent the latter in the form f (z) = (z − 1)g(z) where g is an analytic function at z = 1 with the property that g(1) = f ′ (1) = 0. Here, the prime symbol stands for differentiation with respect to z.This representation enables us to reformulate (10) in the form Since p and q have poles of order one and two at z = 1, respectively, this point is classified as a regular singular point of (11), according to Frobenius theory [34].Hence, we can construct solutions of the form The leading behavior at z = 1 is represented by the term (z − 1) ρ where ρ is determined by the indicial equation with The roots of ( 14) are ρ ± = ±2ix h Ω/f ′ (1) and the correct QNM boundary condition at z = 1 reads .
We can easily obtain a compact formula for f ′ (1) if we observe that can be further simplified by taking into account the condition f (1) = 0. Solving the latter with respect to the error function leads to which in turns allows to express (17) as follows At this point, a remark is in order.We note that as √ θ → 0, the rescaled mass parameter µ = M/ √ θ → +∞.This implies that f ′ (1) → 1 aligns with the expected behaviour in the classic Schwarzschild case.Applying this observation to (16) and considering that in this limit x h → 1, it becomes evident that the prescription (16) accurately reflects the QNM boundary condition at the event horizon, as delineated in equation ( 4) of [4] 2. Asymptotic behaviour as z → +∞: We first rewrite (10) as with V ǫ (z) defined in (10).In order to classify the point at infinity, it is convenient to introduce the transformation τ = 1/z.Then, the above differential equation becomes Since for k = 1 we have with D 0 = 0, according to [35], the point at infinity is an irregular singular point of rank one.Consequently, the asymptotic behaviour of the solutions to equation ( 20) can be efficiently deduced using the method outlined in [36].For this purpose, we start by observing that Given that at least one of the coefficients f 0 , g 0 , g 1 is nonzero, a formal solution to ( 20) is represented by [36] ψ where λ 1 , λ 2 , µ 1 and µ 2 are the roots of the characteristic equations A straightforward computation shows that λ ± = ±2ix h Ω and µ ± = ±2iΩ.As a result, the QNM boundary condition at space-like infinity can be expressed, as It is gratifying to observe that it reproduces correctly the corresponding condition for the classic Schwarzschild case.
At this point, we can proceed as in [4], that is, we transform the radial function ψ Ωℓǫ (z) into a new radial function Φ Ωℓǫ (z) such that the QNM boundary conditions are automatically implemented and Φ σℓǫ (z) is regular at z = 1 and at space-like infinity.To this aim, we consider the transformation Note that in the limit of a vanishing noncommutative parameter, the above ansatz reduces to the corresponding one represented by equation ( 5) in [4] for the classic Schwarzschild case.If we substitute (28) into (10), we end up with the following ordinary differential equation for the radial eigenfunctions, namely with Let us now introduce the transformation z = 2/(1 − y) mapping the point at infinity and the event horizon to y = 1 and y = −1, respectively.Furthermore, a dot denotes differentiation with respect to the new variable y.Then, equation ( 29) becomes TABLE I: Classification of the points y = ±1 for the relevant functions defined by ( 42), ( 36), (37) and (38).The abbreviations z ord n and p ord m stand for zero of order n and pole of order m, respectively.y f (y) Vǫ(y) S2(y) S1(y) S0(y) −1 z ord 1 z ord 1 z ord 4 z ord 3 z ord 3 +1 +1 z ord 2 +1 p ord 2 p ord 2 where with Notice that we must also require that Φ Ωℓǫ (y) is regular at y = ±1.As a result of the transformation introduced above, we have Table I shows that the coefficients of the differential equation ( 35) share a common zero of order 3 at y = −1 while y = 1 is a pole of order 2 for the coefficients S 1 and S 0 .Hence, in order to apply the spectral method, we need to multiply (35) by (1 − y) 2 /(1 + y) 3 .As a result, we end up with the following differential equation where with It can be easily verified with Maple that (i, j) lim where and In the final step leading to the application of the spectral method, we recast the differential equation ( 43) into the following form with Moreover, in Table II, we have summarized the L ij appearing in (60)-(62) and their limiting values at y = ±1.

B. The extreme case
In this scenario characterized by µ = µ e , we still introduce the rescaling x = r/(2M ).However, since the Cauchy and the event horizon coincide, i.e. x e = x c = x h = 0.7936575898 . .., it is useful to introduce a second rescaling given by ξ = x/x e .Hence, the g 00 metric coefficient, as expressed by (7), becomes while the radial part of a massless Klein-Gordon field can be expressed in the equivalent form with f e (ξ) as given in (63).The fact that f e (ξ) exhibits a zero of order two at ξ = 1, implies that f e (1) = 0 = f ′ (1) where the prime denotes differentiation with respect to ξ.This observation allows us to derive the following functional relations which, in turn, plays an important role in simplifying the forthcoming computations.In order to derive the QNM boundary conditions at the event horizon and at infinity, we first need to determine the asymptotic behaviour of the radial solution ψ Ωℓǫ as ξ → 1 + and as ξ → +∞.We can then accurately extract the QNM boundary conditions from this asymptotic data.
1. Asymptotic behaviour as ξ → 1 + : Taking into account that ξ = 1 is a double zero of f e (ξ), we can represent the latter in the form f e (ξ) = (ξ − 1) 2 h(ξ) where h is an analytic function at ξ = 1 with the property that h(1) = f ′′ e (1)/2 = µ 2 e x 2 e − 1 ≈ 1.2837.This representation enables us to reformulate (64) in the form Because q has a fourth-order pole at ξ = 1, it follows that ξ = 1 qualifies as an irregular singularity, rendering Frobenius's theory inapplicable in this scenario.On the other hand, since for k = 1 we have with d 0 = 0, then, according to [35], ξ = 1 is an irregular singular point of rank one.Consequently, the leading behaviour of the solutions to equation (66) in a neighbourhood of the event horizon can be efficiently deduced using the method outlined in [36].To this purpose, we start by observing that by means of the transformation τ = (ξ − 1) −1 mapping the event horizon at infinity and infinity to zero, (66) becomes Since at least one of the coefficients c 0 , d 0 , d 1 is nonzero, a formal solution to (66) is given by [36] ψ where λ ± , and µ ± are the roots of the characteristic equations A straightforward computation shows that At this point, it is important to note that a radial field exhibiting solely inward radiation near the event horizon (ξ → 1 + ) corresponds, under the transformation τ = (ξ − 1) −1 , to an outward radiating field as we approach spatial infinity (τ → +∞ − ).This indicates that we should choose the plus sign in the formulas above.Consequently, the correct QNM boundary condition at ξ = 1 is as follows with µ + and λ + defined in (73).
2. Asymptotic behaviour as ξ → +∞: By means of the transformation η = 1/ξ, it is not difficult to verify that the point at infinity is again an irregular singular point of rank one.Therefore, the asymptotic behaviour of the solutions to equation (66) can be derived according to the method outlined in [36].To this purpose, we observe that With the help of ( 26), we immediately find that λ ± = ±2ix e Ω and µ ± = ±2iΩ.Hence, the QNM boundary condition at space-like infinity can be expressed as Let us transform the radial function ψ Ωℓǫ (ξ) into a new radial function Φ Ωℓǫ (ξ) such that the QNM boundary conditions are automatically implemented and Φ σℓs (ξ) is regular at ξ = 1 and at space-like infinity.To this aim, we consider the transformation If we rewrite it in a more compact form, namely and we replace it into (64), we end up with the differential equation with Let us introduce the transformation ξ = 2/(1 − y) mapping the point at infinity and the event horizon to y = 1 and y = −1, respectively.Furthermore, a dot denotes differentiation with respect to the new variable y.Then, equation (80) becomes where S 0e (y) = Ω 2 Σ 2e (y) + iΩΣ 1e (y) + Σ 0e (y) (89) and the requirement that Φ Ωℓǫ (y) is regular at y = ±1.As a result of the transformation introduced above, we have and Table III shows that the coefficients of the differential equation ( 86) share a common zero of order 4 at y = −1 while y = 1 is a pole of order 2 for the coefficients S 1e (y) and S 0e (y).Hence, in order to apply the spectral method, we need to multiply (86) by (1 − y) 2 /(1 + y) 4 .As a result, we end up with the following differential equation where with It can be easily checked with Maple that Finally, in order to apply the spectral method, we rewrite the differential equation ( 95) into the following form

III. NUMERICAL METHOD
In order to solve the differential eigenvalue problem (106) to determine the QNMs along with the corresponding frequencies Ω, we have to discretize the differential operators L (e) j [•] with j = 1, 2, 3 defined in (107)-(109).Since our problem is posed on the finite interval [−1, 1] without any boundary conditions, more precisely, we only require that the QNMs be regular functions at y = ±1, then, it is natural to choose a Tchebyshev-type spectral method [37,38].Namely, we are going to expand the function y → Φ Ωℓǫ (y) in the form of a truncated Tchebyshev series where N ∈ N is kept as a numerical parameter, {a k } N k=0 ⊆ R, and {T k (y)} N k=0 are the Tchebyshev polynomials of the first kind After substituting expansion (110) into the differential equation ( 106), we obtain an eigenvalue problem with polynomial coefficients.In order to translate it into the realm of numerical linear algebra, we employ the collocation method [38].Specifically, rather than insisting that the polynomial function in y is identically zero (a condition equivalent to having polynomial solutions for the differential problem as per equation ( 106)), we impose a weaker requirement.This involves ensuring that the polynomial vanishes at N + 1 strategically selected points.The number N + 1 coincides exactly with the number of unknown coefficients {a k } N k=0 .For the collocation points, we implemented the Tchebyshev roots grid [39] In our numerical codes, we also implemented the second option of the Tchebyshev extrema grid The users are free to choose their favourite collocation points.Notice that we used the roots grid in our computation, and in any case, the theoretical performance of the two available options is known to be absolutely comparable [38,39].Upon implementing the collocation method, we derive a classical matrix-based quadratic eigenvalue problem, as detailed in [40] In this formulation, the square real matrices M j , each of size (N + 1) × (N + 1) for j = 0, 1, 2, represent the spectral discretizations of the operators L (e) j [•], respectively.The problem (114) is solved numerically with the polyeig function from Matlab.This polynomial eigenvalue problem yields 2(N + 1) potential values for the parameter Ω.To discern the physical values of Ω that correspond to the black hole's QNMs, we first overlap the root plots for various values of N in equation ( 110), such as N = 100, 150, 200.We then identify the consistent roots whose positions remain stable across these different N values.
In order to reduce the rounding and other floating point errors, we performed all our computations with multiple precision arithmetic that is built in Maple and which is brought into Matlab by the Advanpix toolbox [41].All numerical computations reported in this study have been performed with 200 decimal digits accuracy.

IV. NUMERICAL RESULTS
In this section, we delve into the analysis of the QNMs for both non-extreme and extreme noncommutative geometryinspired Schwarzschild black holes.The central revelation from our numerical computations is the stability of both black hole regimes under scalar, electromagnetic, and gravitational perturbations.
The fact that the noncommutative geometry-inspired Schwarzschild black hole reproduces its classical counterpart for large values of the rescaled mass parameter µ provides a basis for validating our numerical method.By selecting a sufficiently large µ-value, computing the QNMs, and demonstrating their agreement with the results obtained by [4,14,15] for the classic Schwarzschild black hole in the s = 0 case, we can substantiate the accuracy of our approach.To this end, we refer to Table V, where the numerical values listed in the last column were derived using our high-precision computational approach, employing 200 polynomials with an accuracy of 200 digits.
As an additional note, we would like to highlight that in the extreme case, the transition to the Schwarzschild limit is more complex due to potential non-analyticities in the behavior of the metric and associated physical quantities as θ tends towards zero.The extreme configuration in noncommutative geometry, characterized by a specific critical mass µ e at which the inner and outer horizons converge, does not straightforwardly reduce to the Schwarzschild solution by merely setting θ to zero.This non-analyticity arises because the extreme condition itself alters the underlying spacetime structure in a way that does not directly map onto the classical Schwarzschild case without noncommutativity.
It is interesting to observe that [28] utilized a sixth-order WKB approximation to show that the WKB method does not converge in critical cases marked by third-order instabilities.These occur for specific choices of the rescaled mass parameters close to µ = µ e for a nonextreme, noncommutative geometry-inspired Schwarzschild black hole.Instabilities were also detected by [28] in the case of an extreme black hole.This indicates that conventional methods have exhibited inconsistencies, particularly when applied to non-commutative geometry settings.There are several reasons for this.For instance, the WKB approximation relies on the local properties of the potential near its peak.However, it fails in cases where the second derivative of the potential at the peak is zero.As for the Inverse Potential and Asymptotic Iteration Methods, the former depends on a rough approximation of the potential at the maximum, while the latter depends on the choice of initial parameters or specific local conditions that may not be optimally TABLE V: This table presents the quasinormal frequencies for scalar perturbations (s = 0) of the classic Schwarzschild black hole alongside those for a noncommutative geometry-inspired Schwarzschild black hole with a large mass parameter (refer to the last column).The third column showcases numerical values determined by [4] using the continued fraction method.The fourth column includes third-order WKB results reported by [15], and the fifth column presents numerical values from [14], derived via the spectral method employing 40 polynomials.Here, Ω represents the dimensionless frequency as defined in equation (10).The notation 'N/A' indicates data not available.representative or easily determinable in complex geometries like those modified by the non-commutative parameter.These dependencies can lead to instabilities or divergences in computations, as indicated by the lack of convergence observed in higher-order WKB calculations, suggesting that these methods may reach their analytical limitations in such scenarios.Unlike local methods, the spectral method utilizes the global properties of the potential, thus avoiding inaccuracies due to local approximations.Moreover, it provides high accuracy and exponential convergence rates for the solutions of differential equations once the irregular part of the solution has been analytically captured.This is particularly beneficial in the quest for QNMs of astrophysical objects.Furthermore, the method is less sensitive to the choice of initial parameters (in contrast with the Asymptotic Iterative Method), making it more robust and reliable for theoretical explorations in gravitational physics where exact solutions are often challenging to compute.For noncommutative geometry-inspired black holes, the spectral method allows us to capture the modified dynamics without the constraints and approximations imposed by traditional techniques.This capability is crucial in revealing that previously reported instabilities are likely artifacts of the limitations of those methods rather than inherent physical instabilities in the black hole models themselves.
Consequently, by harnessing the global analytical capabilities of the spectral method, we have circumvented the shortcomings of traditional approaches like the WKB method.This has allowed us to demonstrate that the instabilities previously reported in non-commutative geometry-inspired black holes are not inherent to the black holes themselves but are rather artifacts resulting from the limitations of local approximation methods.
Therefore, our application of the spectral method provides a clearer, more accurate depiction of the stability characteristics of these black holes.For a comprehensive presentation of our numerical results on massless scalar, electromagnetic, and gravitational perturbations, we invite readers to consult Tables VI, VII, and VIII.These tables detail the calculated QNMs of a nonextreme noncommutative black hole, illustrating variations across different values of the mass parameter.In contrast, Tables IX, X, and XI focus on the extreme case, providing analogous data for the same spectrum of perturbations.
Our findings suggest a remarkable stability of the nonextreme noncommutative Schwarzschild black hole across various perturbations.This observation is particularly noteworthy given the existing literature, such as the study by [32], which posited potential instability in similar scenarios.Our results, however, indicate that even in extreme scenarios where the Cauchy horizon coincides with the event horizon, the black hole exhibits stability under the class of perturbations considered.This result extends the work of [31], where the stability of the nonextreme noncommutative Schwarzschild black hole's interior was established by investigating the propagation of a massless scalar field between the two horizons.

V. CONCLUSIONS AND OUTLOOK
In this study, we have presented a comprehensive analysis of QNMs for noncommutative geometry-inspired Schwarzschild black holes, encompassing both non-extreme and extreme cases.Our investigation, employing the spectral method, covered scalar, electromagnetic, and gravitational perturbations.The significance of our findings TABLE VI: QNMs for scalar perturbations (s = 0) of the noncommutative geometry inspired Schwarzschild metric (nonextreme case) for different values of the mass parameter µ.The third and fifth columns represent the numerical values found by [28] by means of a sixth-order WKB approximation.The fourth and sixth columns report the corresponding results obtained by means of the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation (10).The 'SM' abbreviation denotes the 'Spectral Method'.   is the remarkable stability exhibited by both non-extreme and extreme versions of these black holes across various perturbations.Contrary to previous claims in the literature, our results suggest that the noncommutative geometryinspired Schwarzschild black holes are stable, a finding that is crucial for the understanding of their physical properties.This is particularly relevant in light of the existing debate over the stability of such black holes.Our study, therefore, provides a new perspective on the dynamics of these exotic objects and contributes to a more nuanced understanding of their behaviour under perturbations.[28] using the WKB approximation and the inverted potential method (PT), respectively.The fifth column includes the results obtained via the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation (64).The 'SM' abbreviation denotes the 'Spectral Method'.
ℓ n Ω (WKB) [28] Ω (PT) [28] Ω (SM) 0 0 0.0395 + 0.1367i 0.1003 − 0.1232i 0.1027 − 0.0928i 1 0.1524 + 0.4582i 0.1003 − 0.3697i 0.0299 − 0.3374i 1 0 0.2699 − 0.0744i 0.2983 − 0.0997i 0.2882 − 0.   Our future endeavours will focus on expanding this study in several key directions.Firstly, we plan to investigate the QNMs in the context of massive perturbations.This extension will allow us to explore the influence of mass on the stability and dynamics of these black holes, offering further insights into their complex nature.Additionally, we aim to apply our methods to the study of QNMs of noncommutative geometry-inspired wormholes [42,43].This will not only broaden our understanding of wormhole physics but also provide a comparative framework to evaluate the similarities and differences between black holes and wormholes within the context of noncommutative geometry.Another intriguing direction for future research involves the investigation of the so-called dirty black holes TABLE XI: This table presents the quasinormal frequencies for vector-type gravitational perturbations (s = 2) of the extreme noncommutative geometry-inspired Schwarzschild black hole via the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation (64).as proposed by [43].These black holes, characterized by additional matter fields surrounding them, present a more realistic and complex scenario.Analyzing their QNMs will shed light on how external matter influences the stability and quasinormal spectra of black holes, further enriching our understanding of these objects.In summary, our research opens new avenues for exploring the intricate dynamics of noncommutative geometry-inspired black holes and wormholes.By extending our analysis to include massive perturbations and dirty black holes, we anticipate uncovering new aspects of these entities, contributing to the evolving landscape of gravitational and astrophysical research.

TABLE II :
Definitions of the coefficients Lij and their corresponding behaviours at the endpoints of the interval −1 ≤ y ≤ 1.The symbols appearing in this table have been defined in (44)-(56).

TABLE III :
Classification of the points y = ±1 for the relevant functions entering in (87), (88) and (89).The abbreviations z ord n and p ord m stand for zero of order n and pole of order m, respectively.

TABLE IV :
Definitions of the coefficients L ij and their corresponding behaviours at the endpoints of the interval −1 ≤ y ≤ 1.The symbols appearing in this table have been defined in (96)-(105).

TABLE VII :
(10) for electromagnetic perturbations (s = 1) of the noncommutative geometry inspired Schwarzschild metric (nonextreme case) for different values of the mass parameter µ.The corresponding numerical results are obtained by means of the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation(10).

TABLE VIII :
(10) for vector-type gravitational perturbations (s = 2) of the noncommutative geometry inspired Schwarzschild metric (nonextreme case) for different values of the mass parameter µ.The corresponding numerical results are obtained by means of the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation(10).

TABLE IX :
This table presents the quasinormal frequencies for scalar perturbations (s = 0) of the extreme noncommutative geometry-inspired Schwarzschild black hole.The third and fourth columns showcase numerical values determined by

TABLE X :
This table presents the quasinormal frequencies for electromagnetic perturbations (s = 1) of the extreme noncommutative geometry-inspired Schwarzschild black hole via the spectral method employing 150 polynomials with an accuracy of 200 digits.Here, Ω represents the dimensionless frequency as defined in equation (64).