Black Hole Quasinormal Modes and Seiberg-Witten Theory.

We present new analytic results on black hole perturbation theory. Our results are based on a novel relation to four-dimensional N = 2 supersymmetric gauge theories. We propose an exact version of Bohr-Sommerfeld quantization conditions on quasinormal mode frequencies in terms of the Nekrasov partition function in a particular phase of the Ω -background. Our quantization conditions also enable us to find exact expressions of eigenvalues of spin-weighted spheroidal harmonics. We test the validity of our conjecture by comparing against known numerical results for Kerr black holes as well as for Schwarzschild black holes. Some extensions are also discussed.

Introduction

Finding analytic solutions in spectral theory of quantum mechanical operators is hard. Nevertheless it is recently recognized that a geometric perspective of spectral theory [1, 2] often provides us with many useful tools, developed in supersymmetric gauge theories [3-6] and topological string theory [7, 8], to obtain exact solutions for new families of quantum spectral problems.

In this paper we apply such a geometric/gauge theoretical perspective to the spectral problems governing black hole perturbation theory. More precisely we study the so-called quasinormal modes (QNMs). These modes do not correspond to bound states (or normal modes) but rather to resonance states (or dissipative modes) in quantum mechanics.1 Their spectra are discrete and complex. The QNMs are responsible for the damped oscillations appearing, for example, in the ringdown phase of two colliding black holes and have a direct connection to gravitational waves observations [9]. We refer to [10-14] for a review on the subject and a more exhaustive list of references.

We point out in this work that these spectral problems can be “solved” by using four-dimensional supersymmetric gauge theories in a particular phase of the -background [15, 16]. We refer to it as the Nekrasov–Shatashvili (NS) phase [3]. Our first step is to identify the corresponding matter contents on the gauge theory side. We do so by comparing a master wave equation in black hole perturbations with a differential equation originating from Seriberg–Witten (SW) theory [17, 18]. Interestingly, the matter contents encode the dimension of the black hole as well as the type of asymptotic geometries. Higher dimensional black holes are described by four-dimensional gauge theories.

In this work, we deeply look into four-dimensional asymptotically flat Schwarzschild and Kerr black holes. We find that QNM frequencies of these black holes are determined by Bohr-Sommerfeld-type quantization conditions in SU(2) SW theory with three fundamental hypermultiplets (). Moreover, the extremal limit of the Kerr black holes turns out to correspond to the decoupling limit in SW theory, where one of the masses is sent to infinity and we are left with two fundamental hypermultiplets (). This kind of quantization conditions has already been proposed in [19] for Schwarzschild black holes in the complex WKB approach. However, the proposal in [19] takes the form of a formal power series in the Planck parameter. This series is known to have zero-radius of convergence. Therefore one has to truncate the infinite sum at an optimal order as was done in [19] or to resum it by the Borel summation technique. We emphasize that our quantization conditions overcome this difficulty. Our equations still have an infinite sum, but it has finite-radius of convergence, in particular we can expressed them using a new class of special functions: the Nekrasov-Shatashvili functions. The situation is thus crucially different from [19]. One more advantage is that it is easy to extend it to the Kerr black holes. We explicitly show that our quantization conditions reproduce the numerical QNM frequencies correctly. We also present a new analytic result on the spin-weighted spheroidal eigenvalues that are eigenvalues of the angular part of the Teukolsky equation. We find an exact expression in terms of a gauge theoretical function.

We also have to note that there is already a similar attempt to map problems in black hole perturbation theory to those in two-dimensional conformal field theories with central charge [20-25].2 Via the Alday-Gaiotto-Tachikawa (AGT) correspondence [27], they turn out to correspond to the graviphoton phase of the -background. On the contrary our framework corresponds to the NS phase of the -background which, via the AGT correspondence, makes contact with two-dimensional conformal field theories where .3 This approach leads to a simpler and more systematic solution of the problem. For instance, we find a simple closed form expression of the spin-weighted spheroidal eigenvalues in terms of the Nekrasov partition function in the NS phase [see Eq. (4.16)]. In addition, we can perform explicit computations of the quasinormal modes frequencies and compare with the known numerical data.

This paper is structured as follows. In Sect. 2 we present the general idea behind the geometric/gauge theoretic approach to spectral theory. The main building blocks in this setup are the quantum periods which we compute explicitly by using the NS phase of the -background, see Eqs. (2.17) and (2.18).

In Sect. 3 we study the simplest example: the four-dimensional asymptotically flat Schwarzschild black holes. We propose an exact quantization condition for the corresponding QNMs frequency, see Eq. (3.10). We test our proposal against available numerical data.

In Sect. 4 we study the four-dimensional asymptotically flat Kerr black holes. We find an explicit expression of the angular eigenvalues, see Eq. (4.16), and we propose an exact quantization condition for the radial Teukolsky equation, see Eq. (4.18) for the generic situation and Eq. (4.19) for the extremal limit.

In Sect. 5 we briefly discuss higher dimensional black holes. We then conclude by presenting some future directions. In Appendix A we recall the definition of the NS free energy.

A Geometric Approach to Spectral Theory

We begin with a basic review on the geometric approach to spectral theory. The building blocks in this approach are the so-called quantum periods. Once we obtain such periods we can easily determine other spectral quantities, such as the Fredholm determinant or the quantization condition. The key idea of the geometric approach is to relate the quantum periods to gauge theoretical quantities. The quantum periods can be introduced as follows.4 Given an operatorwe associate it with a classical curve and a one formwhere x and p are complex variables. For a given operator (2.1), its classical limit is unique, but inversely a classical curve (2.2) generates an inifinite number of quantum operators. One has to fix a quantization scheme to identify a quantum operator.

In this work we will always deal with operators (2.1) such that the corresponding classical curve (2.2) coincides with a SW curve of a suitable four dimensional gauge theory. In this situation we refer to the operator (2.1) as the quantum SW curve. Moreover, if the curve has genus g, we choose a basis of cyclesand define the classical periods by integrating over such cyclesIn the gauge theoretic framework, (2.4) are identified with the SW periods and they encode the central charges and masses of the BPS particles in the theory [17, 18]. In this context E is usually denoted by u and it parametrises the moduli space of vacua.

Next we define WKB quantum periods by promoting the differential to a quantum differentialThis is formally defined in such a way thatWe then introduce WKB quantum periods byAt the leading order (), it reproduces the classical periods (2.4). When n is large one typically finds thatTherefore (2.7) has zero-radius of convergence and to define them non-perturbativelty we need to find a way to resum their -expansions. The resummation problem is highly non-trivial. This is why supersymmetric gauge theory plays a crucial role in the problem: it allows to resum the -expansions of the WKB quantum periods into well-defined non-perturbative object: the quantum periods [3, 5, 31, 32]. We denote them byThis resummation is usually done either by using thermodynamic Bethe ansatz (TBA) equations [5, 33] or by using instanton counting [3] (one often refers to this construction as Bethe/gauge correspondence). See [29, Sec. 3 and 4] for a more detailed discussion. In this paper we will use the instanton counting approach.

The main examples we consider in this work are operators which arise in the quantization of the SU(2) SW theories with . These were first studied in [34] and they readThe quantities and are parameters that depend on the matter content of the underlying SW theory. We haveas well asIn SW theory are the masses associated to the fundamental hypermultiplets while is related to the gauge coupling/dynamical scale of the theory. Moreoverwhere one can choose without loss of generality. In this work we chose such that is a second-order differential operator. For instance, if we will take5. The classical SW curve behind (2.102.11) has genus one. The corresponding quantum periods are encoded by the Nekrasov-Shatashvili free energyas well as the four-dimensional quantum mirror mapWe refer to Appendix A for a definition of these quantities. More precisely the quantum A period is given bywhile the quantum period B isThe very important fact is that we have a combinatorial formula of the Nekrasov partition function [15]. It directly computes the NS free energy exactly in . The quantum mirror map is also exactly related to the NS free energy. As a consequence, the quantum periods can be exactly reconstructed by only the NS free energy. This is a main reason why the geometric/gauge theoretical approach is so powerful in analyzing spectral theory: supersymmetric gauge theories provide us with a new class of perfectly well-defined special functions, the four dimensional NS functions, which are precisely the functions allowing for an exact solution of such spectral problems. In the same way in which the zeros of the inverse function compute the spectrum of the harmonic oscillator, similarly the NS functions give us the spectrum of four dimensional quantum SW curves.

According to general expectations coming from the Bethe/gauge correspondence, the discrete part of the spectrum of is captured by the following quantization conditionwhere is a numerical constant. This equation has indeed a discrete set of solutions denoted by . Which quantization condition, A or B, we should impose depends on problems or on boundary conditions. Equation (2.19) is regarded as a quantum corrected Bohr–Sommerfeld rule. In fact, at the leading order of , (2.19) reduces to the Bohr–Sommerfeld condition because the quantum periods reduce to the classical periods. Notice that, even though Nekrasov function provides a resummation of the expansion, in principle it could be that such resummation is still missing some non-perturbative contributions of the form which are necessary to get the correct operator spectrum. This is indeed the case if the operator corresponds to a five dimensional quantum Seiberg–Witten curve, see for instance [7]. However for operators corresponding to four dimensional quantum Seiberg–Witten curves (such as the ones studied in this paper) this does not happens. In this case the four-dimensional Nekrasov function contains all the non-perturbative information on the operator spectrum. See for instance [3, 29, 35–40] for proofs/derivations and numerical evidence of this.

Finally, let us stress that the spectral properties of clearly depend on the values of the parameters . On one hand one usually imposes suitable conditions onso that has a real, discrete spectrum. In this case one asks for eigenfunctions of to be in as in [3, 34]. On the other hand if are generic, we can think of the spectral problem in terms of resonances similar to what was done in [37, 41, 42]. One nice aspect of the geometric/gauge theoretic approach is that (2.19) seems to be able to capture the discrete part of the spectrum independently on whether the operator is self-adjoint with normalizable eigenfunctions or not. We refer to [37] for a simple class of unbounded potentials studied within this framework.

Since black hole quasinormal modes are nothing but resonances, we strongly expect that their spectra are computed in the geometric framework. What we will do in the following is to reinterpret the QNM eigenvalue problem geometrically, and connect their defining equations to suitable quantum Seiberg–Witten geometries. For the examples discussed in this paper the relevant quantum curves are (2.102.11) with . We will then impose the quantization condition (2.19) and check that it reproduces the correct numerical QNM frequencies as listed for instance in [43].

In order to make contact with QNMs it is useful to express (2.102.11) in a more convenient form. For instance, by following [44], we can rewritein a formwhere for is given byandLet z beand we redefine the wave function byThen, the new function satisfies the wave equation in a normal form:wherewithThe important observation is that the differential Eq. (2.27) has two regular singular points at and an irregular singular point with Poincaré rank one6 at . Such a differential equation is well-known as the confluent Heun equation [45]. It is also well-known that the master equations in perturbations of the Kerr black holes as well as of the Schwarzschild black holes have the same singularity structure. In the proceeding sections, we will explicitly show the correspondence among the parameters.

In the similar manner, the quantum SW curve for with leads towhereThis is known as the double confluent Heun equation [45]. We will see that the same differential equation appears from the radial part of the Teukolsky equation in the extremal limit.

Quasinormal Modes of Schwarzschild Black Holes

Schwarzschild black holes are static and spherically symmetric solutions to the Einstein equation in the vacuum. The four-dimensional asymptotically flat solution is given bywhereScalar (), electromagnetic () or odd-parity gravitational () linear perturbations of the metric (3.1) are governed by the Regge–Wheeler type equation [46]where is the field encoding the radial part of the perturbation. The potential in (3.3) is:In addition the differential Eq. (3.3) is supplied by the following quasinormal mode boundary conditions [47]:These boundary conditions are satisfied only for special discrete complex values of the frequency .

The computation of the QNM frequencies for the Schwarzschild black holes is already non-trivial. Though there are a lot of numerical ways to compute them, analytic approaches have been less developed. Our goal is to compute such frequencies by using a gauge theoretical approach.

To make contact with the early result (2.27), we rewrite (3.3) byand then we obtain the normal form (see for instance [48-50]):wherewithWe now compare the equations in the form (3.6) to (3.8) with the SU(2) quantum Seiberg–Witten curve for as given in (2.27)–(2.29). Setting , the parameter correspondence is quite simple:Therefore, if we think of the Regge-Wheeler equation from the point of view of the supersymmetric gauge theories it is natural to ask what is the meaning of the quantizaton (2.19) in the context of black holes. We find evidence that, by using the dictionary (3.9), the quantization condition (2.19) for the B-period indeed computes the QNM frequencies. Our conclusion is therefore given by For a given set of quantum numbers , this equation admits a discrete family of complex solutions .

Since the actual computation is intricate, we briefly illustrate it. Using (2.18), the left-hand side in the first equation of (3.10) is expressed by the Nekrasov–Shatashvili free energy. This free energy is computed by Nekrasov’s combinatorial formula (A.7) systematically. For , we have (A.15) with (A.13). The problem is that the NS free energy includes the parameter a that is not directly related to the black hole parameters. To avoid it, we use the Matone Eq. (A.16). This exact relation allows us to express a in terms of E.7 The Matone relation is just the inverse relation of (2.17). Therefore we can finally eliminate a from the NS free energy, and thus we can solve the quantization condition (3.10) with respect to .

Recall that the WKB quantum periods (2.7) are formal power series in . Its radius of convergence is just zero. We cannot plug into it naïvely.8 On the contrary, the (non-perturbatively defined) quantum periods (2.18) are given by the NS free energy which is exact in . All the quantum corrections are already resummed, and we can set without any problems. In this sense, we refer to (3.10) as the exact quantization condition. However, one has to keep in mind that there still remains the sum in the instanton counting parameter . This sum has a finite radius of convergence, and its treatment is easier than the divergent WKB series.

We have performed the procedure above and have checked that it indeed matches the numerical values of the Schwarzschild black hole QNMs as obtained in [10, 51].9 Some examples are given in Tables 1, 2 and 3.

Table 1

The solution to the quantization condition (3.10) for and

Nb	2Mω_0(0,0)
3	0.21453301-0.20342058i
8	0.22088781-0.20978038i
12	0.22090951-0.20979131i
Num	0.22090988-0.20979143i

The matching digits are shown in boldface. We denote by the order at which we truncate the instanton counting series in (A.13). We apply Padé approximants to improve the convergence of the instanton counting series. The numerical value is obtained from [43]

Table 2

Solutions to the quantization condition (3.10) for and

Nb	2Mω_0(1,1)	2Mω_1(1,1)
4	0.493115-0.180881i	0.43066732-0.5887236i
8	0.496470-0.184999i	0.42899228-0.5873530i
12	0.496526-0.184974i	0.42903098-0.5873354i
Num	0.496527-0.184975i	0.42903084-0.5873353i

The matching digits are shown in boldface. We denote by the order at which we truncate the instanton counting series in (A.13). We apply Padé approximants to improve the convergence of the instanton counting series. The numerical values are obtained from [43]

Table 3

Solutions to the quantization condition (3.10) for and

Nb	2Mω_0(2,2)	2Mω_1(2,2)	2Mω_2(2,2)
3	0.7480-0.1985i	0.6947713-0.550331i	0.600036-0.953084i
7	0.7446-0.1890i	0.6933273-0.548018i	0.602154-0.956237i
12	0.7472-0.1777i	0.6934216-0.547829i	0.602101-0.956556i
Num	0.7473-0.1779i	0.6934220-0.547830i	0.602107-0.956554i

The matching digits are shown in boldface. We denote by the order at which we truncate the instanton counting series in (A.13). We apply Padé approximants to improve the convergence of the instanton counting series. The numerical values are obtained from [43]

One issue that we encounter in the computations is that the NS free energy (A.10) is given by the natural series expansion in the parameter . Even though this series converges, the convergence is not very fast. In that perspective it may be useful to compute the quantum periods by using TBA equations as was done in [5, 29, 33] instead of using the NS free energy.10 We leave this issue as future works.

Quasinormal Modes of Kerr Black Holes

Kerr black holes are stationary and axially symmetric solutions to the Einstein equation in the vacuum. The four-dimensional asymptotically flat solution in the Boyer-Lindquist coordinates is:where M is the mass and is the angular momentum. Perturbations of rotating black holes are described by the Teukolsky equation [53, 54]. The Tuekolsky equation is a separable partial differential equation in the Boyer–Lindquist coordinates.

After separation of variables, its angular part reads (see for instance [10, eq. (25)])where and s is the (minus of) spin of a perturbing field. Moreoverwhere for integer spins and for half integer spins. In the black hole perturbation, the parameter c is related to the angular momentum and the frequency byThe eigenfunction is called the spin-weighted spheroidal harmonics in the literature. Its eigenvalue is determined by the regularity condition of at . For general s, l, m and c, no closed form of is known so far. However, for the spheroidal harmonics reduces to the spin-weighted spherical harmonics and one hasThe radial Teukolsky equation is more complicated and reads (see for instance [10, eq. (25)]),where . The potential iswhere . Note that the radial differential Eq. (4.5) has (regular) singular points at corresponding to the Cauchy and event horizons. In addition, (4.5) is supplied by the following boundary conditions (see for instance [10, eq. (80)])whereBoth the angular and the radial parts of the Teukolsky equation have the same singularity structure as the confluent Heun equation, see for instance [48-50].

For the angular part, we change the variable , and define . Then we obtainwhere Q(z) takes the formThe coefficients in Q(z) are computed straightforwardly. Similarly, defining and for the radial part, we obtain the same form as (4.9) and (4.10) with different coefficients. Hence, the Teukolsky equation also has a connection with the quantum Seiberg–Witten geometry with gauge group SU(2) and hypermultiplets. To find the precise dictionary we need to compare (4.9) and (4.10) with (2.27)–(2.29).

For the angular part, we findFor the radial part, we haveWhen , it reproduces the identification in the Schwarzschild case by exchanging . This relabelling comes from the fact that the Teukolsky equation at does not take the form of the Regge-Wheeler equation.

Notice that the extremal limit in Kerr black holes corresponds to . Given the above dictionary, this translates intoFrom the gauge theory point of view this is precisely the decoupling limit under which the theory flows to the .

Hence the radial part of the Teukolsky equation in the extremal limit corresponds to the quantum Seiberg–Witten curve of SQCD with , provided we use the following dictionaryThis result is also obtained by starting with (4.5) in the extremal case . By changing the variable , the radial equation is finally written as the normal form of the double confluent Heun equation which we can compare with (2.30).

The identifications (4.11) and (4.12) allow us to find equations determining both and by using quantities in SW theory. We will demonstrate this in the next two subsections.

Exact Quantization Condition for the Angular Teukolsky Equation

Note that in the angular equation (4.2) or the identification (4.11), the multipole number does not appear explicitly. In the perspective of the quantization conditions, this must appear as a quantum number. It turns out that the angular eigenvalues are exactly determined by the A-period quantization condition:where . The difference between the A-period condition here and the B-period condition in the previous section (and also in the next subsection) is explained as follows. In the angular problem we impose the boundary conditions at the two regular singular points . In the radial problem, we have to impose the conditions at the regular singular point and at the irregular singular point . Therefore we have to consider different period integrals in the WKB approximation.

Now we rewrite the conjecture (4.15) in a more elegant form. The condition (4.15) is nothing but the quantization condition for a (see (2.17)). Also, recall the Matone relation (A.16), which is the inverse of (2.17). Therefore, we finally conclude that Eq. (4.15) is equivalent towhere and is defined in Appendix A. The first few terms readwhere , and we have used defined in (A.13). The identity (4.16) can be compared with the small-c expansion of in [55, 56] up to . Note that the similar consideration is found in [25], but (4.16) looks simpler and more direct. The expression of up to can be found in the Mathematica file attached to the ArXiv version of the manuscript.

Exact Quantization Condition for the Radial Teukolsky Equation

The story for the radial equation is the same as for the Schwarzschild case. The radial equation has a discrete set of complex frequencies , which have been computed numerically for instance in [10, 51, 57]. From a gauge theoretic perspective, and given the identification (4.12), it is natural to conjecture that the frequencies can be obtained by imposing the following B-period quantization conditionwhere is as in (4.16) while is defined in (2.18). We checked that the solutions to (4.18) indeed reproduce the correct QNM’s frequencies as computed numerically in [10, 51] (we took the data from [43]). An example is given in Table 4.

Table 4

Solutions to the quantization condition (4.18) of Kerr BH for with and quantum number

Nb	Mω_0	Mω_1
3	0.1073438-0.1016159i	0.089515-0.330273 i
8	0.1105221-0.1047959 i	0.086036-0.347811 i
11	0.1105330-0.1048013i	0.086216-0.347686i
Num	0.1105331-0.1048015 i	0.086203-0.347664i

The matching digits are shown in boldface. We denote by the order at which we truncate the instanton counting series in (A.13). The numerical values are obtained from [43]

In the extremal limit, we can get a simplified quantization which now involves the quantum period of the theory and readswhere is defined by (2.18). We check that the quantization condition (4.19) reproduces the correct numerical QNM’s frequencies in the extremal case. Some examples are given in Table 5.

Table 5

Solutions to the quantization condition (4.19) of Kerr BH in the extremal limit with for and quantum number

Nb	Mω_0	Mω_1
2	0.10626294-0.088291652 i	0.071401-0.3174727 i
6	0.11024328-0.089430682 i	0.062368-0.3186182i
10	0.11024545-0.089433151 i	0.062353-0.3188423i
12	0.11024548-0.089433184 i	0.062472-0.3188436 i
Num	0.110245-0.089433 i	0.062473-0.318840 i

The matching digits are shown in boldface. We denote by the order at which we truncate the instanton counting series in (A.14). We also use Pade approximant to accelerate the convergence. The numerical values are obtained from [57]. Note that in the extremal limit the numerical values are not as precise as in the general case and it seems that for we already get a few additional digits as compared to [57]

Simple Extensions

In the above sections we focused on four-dimensional black holes with asymptotic flatness. Our gauge theoretic approach is not restricted to these particular examples. Here we briefly illustrate some other examples. More detailed analysis will be reported elsewhere.

A lesson we have learned from the previous examples is that it is important to understand the singularity structure of differential equations. The singularity information tells us the matter contents of gauge theories. Our strategy is the following. We first read off the singularity structure of master wave equations for various black holes. Next we look for gauge theory counterparts by comparing the Riemann sphere with punctures associated with quiver gauge theories. Such punctures describe singularities of quantum SW curves. They were analyzed for regular punctures by Gaiotto [58], and generalized for irregular punctures by several people [59-61]. After identifying a gauge theory, the third task is to find relations between the parameters. This last step is quite complicated.

Four-Dimensional Schwarzschild (A)dS Black Holes

A simple extension is four-dimensional asymptotically (A)dS black holes. It is known that the master equations of these black holes have four regular singular points [62]. This differential equation is well-known as the Heun equation. Here we quickly look at it for the Schwarzschild (A)dS case. The higher-dimensional case is discussed in the next subsection. The radial master equations in the scalar/electromagnetic/odd-parity gravitational perturbations take the same form as (3.3), but the functions are now modified asandwhere is the cosmological constant. After the redefinition of , we get the normal form:whereThe algebraic equation has three roots in general. It is easy to see that these three points as well as are regular singular points of the master equation. The infinity point is subtle. To see the behaviour near , it is convenient to change the variable and , and we findWe look at the behaviour of q in the limit :It is obvious that for , the infinity is not a singular point. For however it is a regular singular point.11 We conclude that for the scalar perturbation the master radial equation has five regular singular points, while for the electromagnetic and the odd-parity gravitational perturbations they have four regular singular points. In both cases, all the singular points are regular, and thus the differential equations are Fuchsian.

For the master equation is equivalent to the Heun equation. It turns out that this case corresponds to the gauge theory with four fundamental matters (). The detailed analysis in this case will be reported soon [63]. The singularity structure for is realized in an quiver gauge theory [58]. We do not explain it anymore in this work.

Higher-Dimensional Extensions

Let us proceed to the higher-dimensional cases. For simplicity, we focus on asymptotically (A)dS Schwarzschild black holes in d-dimension. The metric in this geometry iswhere is the line element of the -dimensional unit sphere . The function f(r) takes the formwhere is proportional to the mass with a non-trivial coefficient. We consider the higher-dimensional analog with the odd-parity gravitational perturbation of Regge and Wheeler. In the terminology of Kodama and Ishibashi [64], it corresponds to the vector-type gravitational perturbation. The master equation in this case is again the same as (3.3) with the potential [64]:For , it actually reduces to the one for in the previous subsection.

We look at singularities of the master equation. If , the equation has roots. These are regular singular points in general. It is also easy to check that is a regular singular point. To see the behavior at , we rewrite the master equation in the form (5.5). The coefficient function has the Laurent expansion:We conclude that the infinity is a regular singular point except for . We observe that for , is actually an apparent singularity. For and , the master equations are thus Fuchsian differential equations with singular points. The quantum SW curve having the same singular structure is an quiver gauge theory [58]. It is interesting to see that the dimensional information of the black hole is reflected in the number of quiver gauge groups.

For the flat case , the equation has roots. These as well as are regular singular points. It turns out that the infinity is an irregular singular point with Poincaré rank one. Therefore in this case, we have to consider quiver theories associated with the Riemann sphere with an irregular puncture [59-61].

Outlook

Inspired by recent developments in the gauge theoretical approach to spectral theory, we analysed the black hole QNMs in this framework. This approach provided us with some new analytic results on the QNM frequencies. We pointed out that their master equations can be written as quantum SW curves and we obtained an exact, analytic expression for their quantization condition. We mostly focused on four-dimensional asymptotically flat Schwarzschild and Kerr black holes and we check our results against available numerical data. We also presented a preliminary analysis for asymptotically (A)dS case, which is connected to , as well as the higher dimensional examples. In this situation it would be interesting to investigate how the existence of unstable BH solutions is reflected on the SW theory side12. A more detailed study will appear elsewhere. Likewise, even though in this work we focused on uncharged black holes, we expect our analysis to carry on for Reissner-Nordström and Kerr-Newman black holes as well13.

In our approach a key role is played by the quantum periods of the underlying SW theory which we computed by using the Nekrasov–Shatashvili free energy (A.10). Even though this quantity is a convergent series in the instanton counting parameter, the convergence is a bit slow. From that perspective it would be good to compute these quantum periods by using the alternative TBA approach as in [5, 29, 33]. Relatedly, it is an important point to understand more in details the radius of convergence of the Nekrasov–Shatashvili free energy and in particular what is its meaning from the gravity point of view. In addition it would be interesting to study the behaviour of the singularities in the Borel plane and see if there is any interpretation of wall crossing phenomena from the black hole viewpoint, in particular in relation to the Hawking-Page transition for the case.

Although we have been focusing on the quantization condition for the QNMs frequencies, the geometric/gauge theoretic approach to spectral theory also allows for the computation of the eigenfunctions. It would be important to pursue this direction in more details and eventually provide a more rigorous derivation of the quantization condition proposed in this paper.

Our result also indicates that it should be possible to compute efficiently the WKB expansion for the QNMs by using the holomorphic anomaly equation for the underlying gauge theory, similar to what was done in [65-67]. It would be interesting to investigate this aspect more in detail as a possible alternative to the existing approaches, see for instance [68-72].

An additional interesting point is the connection with the work of [73-78] on the highly damped QNMs. From our perspective the large n behaviour of the frequencies is encoded in the asymptotics of the quantum periods/Nekrasov–Shatashvili partition function which is accessible analytically. We hope to report on this in the near future.