Excitations of the Myers-Perry Black Holes

We demonstrate separability of the dynamical equations for all p-form fluxes in the Myers-Perry-(A)dS geometry, extending the earlier results for electromagnetic field. In the physically important cases of p=(1,2,3,4), we explicitly write the ODEs governing the dynamics of separable solutions.


Introduction
Black holes are important laboratories for studying classical and quantum gravity. We are entering a new era in exploration of these objects that was launched by the first direct observations of electromagnetic [1] and gravitational [2] waves produced by them. Reconstruction of images from these measurements required detailed understanding of classical fields in the vicinity of black holes, and the relevant calculations are usually done numerically [3,4]. Nevertheless, in the idealized case of an isolated black hole, all excitations can be studied analytically [5], and these classic results have played an important role in our understanding of classical and quantum properties of the holes. In four dimensions, light excitations of black holes contain four types of fields (scalars, spinors, vectors, and gravitons), and all of them are covered by the Teukolsky's equations [5]. In higher dimensions, one encounters additional excitations, and exploration of their dynamics is the main goal of this article. The motivation for our investigation comes from the desire to understand classical scattering from black holes, as well as quantum effects encoded in the Hawking radiation [6].
Hawking radiation is one of the most important quantum phenomena in gravitational physics. While the essence of this process is the same for all fields [6], the detailed wavefunctions of the emitted particles depend on the properties of the underlying dynamical equations. Furthermore, since black holes radiate all light fields at comparable rates, quantitative understanding of quantum dynamics requires analysis of all available excitations. In four dimensions, wavefunctions for all relevant fields have been computed in [5]. Although eventually we are interested in explaining four-dimensional phenomena, in the last three decades great insights into physics of black holes have been gained from analyzing these objects in higher dimensions within the framework of string theory [7,8,9]. In this case light bosonic excitations are not exhausted by scalars, vectors, and symmetric rank-two tensors, but they also contain higher forms predicted by string theory. Detailed study of such fields, which would shed new light on classical and quantum properties of black holes, is the main objective of this article.
Hawking radiation originates from dynamics of quantum fields on nontrivial gravitational backgrounds. For the static black holes, the original calculations for radiation rates and scattering amplitudes have been extended to various fields in higher dimensions [10], and important lessons extracted from such investigations inspired formulation of the AdS/CFT correspondence [11]. Success of these studies relied on rich SO(d − 1) × U(1) isometries of static black holes in d dimensions, which give rise to (d − 1) commuting conserved charges and guarantee full separation of variables in equations for all dynamical fields. Unfortunately, rotating black holes, described by the Myers-Perry geometry [12], have only [U (1)] [(d+1) /2] isometries, which give rise to [ d+1 2 ] conserved quantities, and they are not sufficient for ensuring full integrability of the dynamical equations. Nevertheless, equations of motion for probe particles and scalar fields in the Myers-Perry geometry turn out to be fully separable due to additional conserved quantities associated with rank-two Killing tensors. The hidden symmetries responsible for this separation were first discovered for the Kerr black hole in [13], and in the last two decades they have been extended to rotating black holes in arbitrary dimensions in [14,15,16,17,18], including some geometries carrying charges [19,20]. Remarkably, separation of variables in the wave equation persists for the GLPP geometries [21] 1 , which generalize the Myers-Perry black holes to solutions with non-zero values of the cosmological constant, so we will briefly discuss excitations of these backgrounds in section 7. In contrast to the static black holes, where decomposition into spherical harmonics persists for fields of arbitrary spin 2 , for rotating spacetimes separation of variables for various fields has been worked out only on a case-by-case basis. In the four-dimensional Kerr geometry, separability of equations for electromagnetic and gravitational waves was proven in the classic work by Teukolsky [5]. Recently separability of Maxwell's equations has been demonstrated for rotating Myers-Perry and GLPP black holes in all dimensions [25] 3 , and the goal of our article is to extend these results to dynamics of higher forms 4 which play an important role in string theory. A better understanding of the black hole excitations will give new handles on probing these fascinating objects. This paper has the following organization. In section 2 we review the known results pertaining to the Myers-Perry geometry and its excitations. Specifically, we recall the special frames [15,16,17], which have played crucial role in separating Maxwell's equations in [25], and review the procedure for separating variables in the equations for the scalar and vector fields. The notation established in section 2 is used throughout this article.
In section 3 we demonstrate separability of the dynamical equation for the two-form, which can be used to describe either NS-NS or Ramond-Ramond fluxes in string theory. We derive the most general separable ansatz in the Myers-Perry geometry and verify that the resulting ODEs contain the expected number of separation constants. In section 4 these results are extended to the three-form potentials. While equations for the higher forms follow the same pattern, the resulting ODEs become rather complicated, so we just present the general structure in section 6 and write the explicit equation for the four-form in the physically interesting ten-dimensional case in section 5. All our results are extended to the GLPP geometry in section 7. Some technical details are presented in the appendices.

Summary of the known results
This article is dedicated to analyzing dynamics of various fields in the backgrounds of the Myers-Perry black holes [12], so we begin with reviewing some general properties of these geometries. In particular, section 2.1 establishes the notation used throughout this paper. The main conclusion of this article is that the procedures for solving equations for all p-forms follow similar patterns, so in sections 2.2 and 2.3 we review the known results for the scalar fields [14,15,16,17,18] and for the one-form potentials [25]. These constructions will be extended to higher forms in the remaining parts of this article. Readers familiar with the earlier work can go directly to section 3.

Myers-Perry black hole and its symmetries
Let us review some well-known properties of the d-dimensional Myers-Perry geometry [12]. The form of metric differs between even and odd values of d, so we begin with quoting the solution in d = 2n + 2 dimensions [12,28]: Here variables (µ i , α) are subject to a constraint (r 2 + a 2 k ). (2.3) To study equations for various fields in the background (2.1), we need to review the symmetries of this geometry, which are encoded the Killing-Yano tensors associated with it [15,16,17,20,18].
Metric (2.1) has (n + 1) Killing vectors corresponding to constant shifts of t and φ i , but these isometries are not sufficient to ensure full separation of the wave equation or dynamical equations for tensor fields. For the wave equation such separation is guaranteed by a family of Killing tensors [15,20], while already for the Maxwell field separation requires a more rigid structure associated with Killing-Yano tensors (KYT) [25]. The general equation for an anti-symmetric KYT [29] is ∇ µ Y ν 1 ...νp + ∇ ν 1 Y µν 2 ...νp = 0, (2.4) and for the Myers-Perry black holes (2.1), all Killing-Yano tensors are given by a very elegant formula [15,16,17] Y 2(n−k) = ⋆ ∧h k . (2.5) Here h is a special two-form that will be written below. As demonstrated in [25], the eigenvectors of the antisymmetric tensor h play the central role in separation of Maxwell's equations, and we expect them to be important in the study of higher forms as well.
Using the eigenvectors of a symmetric tensor t µν = h µα h ν α as frames, one finds very simple expressions for the two-form h and for the metric [15,16,17]: (2.6) The frames were constructed in [15,16,17], and we will write only the components with upper indices, e µ A , using the notation introduced in [20]: Here we defined convenient expressions The map between coordinates {µ i } and {x k } can be found in [20], we just recall that in the variables (r, . (2.9) Note that, apart from the common overall factors, the components (e µ t , e µ r ) of the frames depend only on r, while the components (e µ i , e µ x i ) depend only on x i . This crucial fact is responsible for separation of variables in all equations discussed in this article.
As found in [25] and reviewed in section 2.3, separation of the Maxwell's equations occurs not in components (A t , A r , A φ i , A x i ), but in projections of the gauge fields to particular combinations of the frames (2.7). Specifically, it is convenient to define m Then, as demonstrated in [25], it is the set of projections [m (I) ± ] µ A µ that separates. Note that the components m (j)µ ± depend only on x j , and m (0)µ ± depend only on r. The metric and the generator h of the Killing-Yano tensors (2.5)-(2.6) can be written as The frames (2.10) will play an important role throughout this article.
Let us now discuss the Myers-Perry black holes in odd dimensions (d = 2n + 1). The metric is given by [12,28] (2.12) and coordinates µ i are subject to a constraint The counterparts of the special frames (2.6) are given by [15,20] Most relations (2.8), (2.9) still hold, with two exceptions: 14) The Killing-Yano tensors still have the form (2.5)-(2.6), although the metric acquires an extra term (e ψ ) 2 , and we refer to [20] for the detailed discussion. Finally, the light-cone frames m (I) ± are given by counterparts of (2.10), and there is an additional frame n µ : To study dynamical equations, we will need the expression for the inverse metric in terms of these frames: Note that the components [m ± ] µ depend only on r, and n µ are constants. This feature is crucial for ensuring separation of variables in the Klein-Gordon equation and in dynamical equations for p-form potentials. In this article we are focusing on massless fields, so we begin with reviewing separation of the wave equation.

Separation of the wave equation
As we will see in the later sections, dynamics of a p-form in the Myers-Perry geometry is governed by a separable scalar function that satisfies a system of ordinary differential equations. To get an inspiration for the resulting ODEs, it is useful to review separation of variables in the wave equation. Our discussion of the scalar field in the Myers-Perry geometry will be very brief, and we refer to [25] for the details.
Let us consider the wave equation 2.17) in the background (2.1). Imposing a separable ansatz 18) we conclude that the functions (Φ, X i ) satisfy a system of differential equations [25]: Here P n−1 [y] is an arbitrary polynomial of degree (n − 1), and it is crucial that all equations (2.19) contain the same function P n−1 . Note that the separable ansatz (2.18) reduces one PDE (2.17) to a system of d = 2(n + 1) ordinary differential equations for function Ψ: (n + 1) of them come from (2.19) upon replacement (Φ, X i ) → Ψ, and the remaining (n + 1) relations are (2.20) Thus equations (2.18), (2.19) guarantee full separability of the wave equation, and as expected, the resulting solution depends on d − 1 = 2n + 1 integration constants: (i) (n + 1) parameters (ω, n i ); (ii) n free coefficients of the polynomial P n−1 .
The separation (2.18) also extends to rotating black holes with cosmological constant [15,16,17], and we refer to [25] for the counterparts of equations (2.19) for that case.
Remarkably, at least in some special situations, solutions of the resulting ODEs can be written in terms of Painleve transcendental functions [30]. It would be interesting to see whether a similar reduction to Painleve transcendentals occurs for the ODEs encountered in sections 3-6 of our article 5 .
We conclude this short subsection by writing the counterpart of equations (2.19) for the geometry (2.12) in d = 2n + 1 dimensions [25]: As before, the relations (2.18), (2.21) guarantee full separability of the wave equation, but the counting of d − 1 = 2n free parameters is slightly different: (i) (n + 1) parameters (ω, n i ); (ii) (n − 1) free coefficients of the polynomial P n−2 .
In the remainder of this article we will extend the separation (2.18), (2.19), (2.21) to the dynamical equations for all p-form potentials. We begin with reviewing the p = 1 case that has been solved in [25].

Maxwell's equations in the Myers-Perry geometry
The first extension of (2.18), (2.19), (2.21) beyond the scalar field has been accomplished in [25], and the ansatz found in that article was inspired by the classic analysis of the Maxwell field in the four-dimensional Kerr geometry [5]. We refer to [25] for the detailed discussion of the relation between the higher-dimensional ansatz for the gauge field and the four-dimensional variables introduced by Teukolsky [5]. Here we just mention that, in contract to a gauge-invariant formulation used in [5], separation of variables in higher dimensions occurred in certain projections of the gauge potential, so it worked only in a specific gauge. As we will see in later sections, similar preferred gauges persist for all p-forms, so it is useful to recall the construction of [25]. It is convenient to separate the discussions of odd and even dimensions.

Even dimensions
As demonstrated in [25], the most general separable ansatz for the Maxwell field in even dimensions has the form [m 22) where the scalar function Ψ is given by (2.18). Relations (2.22) have a free parameter µ that will play a role of one of the separation constants. Equations (2.22) can be easily solved for the gauge potential, but we find the expression (2.22) to be more compact and more useful for extending to higher forms. Substitution of the ansatz (2.22) into Maxwell's equations leads to a system of ODEs [25]: Various functions appearing in these equations are defined by Note that the wave equation (2.19) is also covered by a minor modification of the system (2.24): scalar : D r = D j = 1, ∀Λ, P n−2 → P n−1 . (2.26) This implies that in the process of going from scalar to vector systems (i.e., from (2.19) to (2.24)), one coefficient in the polynomial P n−1 is traded for a separation constant µ.
Thus the vector counterpart of the counting presented on page 9 is (i) (n + 1) parameters (ω, n i ); (ii) (n − 1) free coefficients of the polynomial P n−2 ; (iii) constant µ.
As expected, the total number of free parameters is still d − 1 = 2n + 1. Furthermore, as demonstrated in [25], the ansatz (2.22) covers all d−2 polarizations of the electromagnetic field: they correspond to different ranges of µ. Note that the ansatz (2.22) fully fixes the gauge since it is not preserved under a gauge transformation. However, after equations (2.22), (2.24) were derived in [25] as the system describing the most general separable solution of the Maxwell's equations, it was observed in [26] that the resulting potential happened to be in the Lorentz gauge: This observation allowed the authors of [26] to develop an elegant shortcut in deriving the system (2.24) by imposing the gauge condition (2.27) from the beginning. Unfortunately a similar strategy seems to be failing for higher forms 6 , so the original derivation presented in [25] appears to be more suitable for the extension to p-forms.

Odd dimensions
The counterpart of the ansatz (2.22) for odd dimensions was also derived in [25], where it was shown that the most general separable solution of the Maxwell's equations has the form [m 28) with Ψ given by (2.18). The resulting ordinary differential equations are [25] D Here functions (W j , W r ) are still given by (2.25), while the constants (A , Ω) are defined as As in even dimensions, one can show that the system (2.28)-(2.29) describes the most general separable solution of the Maxwell's equations, and it depends on d − 1 separation parameters and covers d − 2 independent polarizations [25]. In the rest of this article we will extend the systems (2.22), (2.24) and (2.28)- (2.29) to equations governing the dynamics of p-form potentials.

Two-form potential in the Myers-Perry geometry
In the last section we reviewed separation of variables in equations for scalar and vector fields in the Myers-Perry geometry. While in four dimensions such excitations, along with spinors and gravitons, exhaust all interesting modes, for d > 4 one also encounters higher forms, such as Ramond-Ramond potentials predicted by string theory. Since the Myers-Perry black hole is a solution of the type II supergravity, it is interesting to study SUGRA excitations of this geometry, and they include higher forms corresponding to either Ramond-Ramond potentials or the Kalb-Ramond field. In this section we will extend the success in separating equations for the scalar and vector perturbations to equations of motion for a two-form field, and dynamics of higher forms will be discussed in sections 4-6. To emphasize the analogy with the Maxwell field and with higher forms, we will denote the dynamical variables by A, and our analysis is equally applicable to the Kalb-Ramond field B and to the Ramond-Ramond potential C (2) . As in the case of electromagnetism, the discussion naturally splits into even-and odd-dimensional cases, which are covered in subsections 3.1, 3.2.

Even dimensions
Let us consider a two-form potential A that obeys a system of linear differential equations Here Hodge dual is taken with respect to the Myers-Perry geometry (2.1). Since the metric has (n + 1) isometries corresponding to translations along the (t, φ i ) coordinates, we can look for solutions in the form As in the case of electromagnetism, we expect that separable solutions are governed by one "master function" Ψ, and now we will determine the relation between the gauge potential and such function. We begin with focusing on a special case Since the metric has a block form 4) and there is no dependence on φ a = {t, φ i }, we conclude that the field equations for configurations (3.3) can be divided into three decoupled groups that involve A (1) , A (2) , or A (3) : 3.5) Thus in the special case (3.3) one would have three independent "master functions" (Ψ (1) , Ψ (2) , Ψ (3) ), but if a configuration (3.5) is viewed as a limit of a general solution with nontrivial (ω, n i ), then all Ψ (a) must be the same, so the dynamics is governed by a single scalar Ψ: It is clear that we are interested only in excitations of this type. For separable solutions of equations (3.1), function Ψ must have the form All three cases (3.5) are analyzed in the Appendix A.2, where it is shown that the most general separable solutions in each class are given by 7 j , Q j ) introduced here have the form and a priori there are independent sets of coefficients (a (a) , b (a) , c (a) ) for three values of the superscript. Configurations (3.8) solve the field equations (3.1) if and only if each of the functions (Ψ (1) , Ψ (2) , Ψ (3) ) obeys a system of ODEs is a polynomial of degree (n−3) in its argument. Imposing the restriction (3.6), we conclude that functions Q (a) j must be the same for all three values of a, and P (a) n−3 must have the same property.
Combining relations (3.8), we find the expression for various projections of the most general gauge potential with a separable function Ψ: Here F is a polynomial that is at most quadratic in each of its arguments, and motivated by the form of (2.22) and (3.8), we require function F to be separable. This leads to identification of (p, q) with (p,q) in (3.8) and to the final expression for the most general separable ansatz in the special case (3.3) 8 : Q I = ax 4 I + bx 2 I + c, V I,α = px 2 I + q + αx I . Here Ψ is a function with factorized dependence on the coordinates (3.7), that satisfies ordinary differential equations 13) with the same polynomial P n−3 for all values of j. Let us now extend the solution (3.12) -(3.13) beyond the special case (3.3).
Using the separable ansatz (2.22) for the Maxwell's equations as a guide, we assume that addition of nontrivial angular dependence does not modify relations (3.12), but results in a new expression for function Ψ: In other words, we assume that coefficients (a, b, c, p, q) in (3.12) define a separable ansatz for all values of (ω, n i ), and the quantum numbers (ω, n i ) affect only the final ODEs. Specifically, using the scalar and vector cases as guides (see equations (2.19), (2.24)) we impose differential equations 15) and assume that (ω, m i ) dependence appears only in functions f j [x 2 j ]. For nontrivial angular dependence, such restrictive ansatz leads to inconsistency of the field equations (3.1) unless coefficients (a, b, c, p, q) are constrained by relations 9 : These restrictions imply factorization of the quartic polynomials Q I : The easiest way to see this is to look at (ij) components of (3.1), but since the relevant expressions are rather unwieldy, we do not quote them here. A detailed derivation of the counterpart of (3.16) for Maxwell's equations was presented in [25], and extension to higher forms follows the same logic.
Once the ansatz (3.12), (3.14) with factorization (3.17) is imposed, all field equations (3.1) reduce to a system of ODEs, just as in the case of electromagnetism.
To present the final answer in a compact form, we begin with simplifying the ansatz (3.12) in the presence of the constraint (3.17) 10 : Substitution of this configuration with separable function Ψ (3.14) into the field equations (3.1) in d = 2(n + 1) dimensions leads to a system of ODEs: Here we introduced convenient notation that will be used throughout this article: We also defined functions D j and D r that are specific to the two-form: Equations (3.18), (3.14), (3.19), (3.20) summarize our main result for the two-form gauge potential in the even number of dimensions. As expected, there are d − 1 = 2n + 1 free separation constants: (i) (n + 1) parameters (ω, n i ); (ii) (n − 2) coefficients of the polynomial P n−3 ; (iii) two ratios (e 2 /e 1 , e 3 /e 1 ) appearing in the expressions (3.14) for the functions h I ± .
Although one can set e 1 = 1 by rescaling the gauge potential, we kept all three parameters (e 1 , e 2 , e 3 ) to make the expressions more symmetric.
The counting presented above works for n ≥ 3, but it clearly fails in four dimensions. Formally, for n = 2 equations (3.19)-(3.20) contain a polynomial P −2 with (−1) degree of freedom 11 , and this clearly indicates the failure of the description (3.19)-(3.20). It turns out that in the degenerate four-dimensional case, the ratio e 2 /e 1 is not an independent parameter, but rather it is determined in terms of (ω, n). A direct calculation for the Kerr black holes gives The differential equations are Four-dimensional system has another peculiarity: a harmonic two-form A can be dualized into a scalarΨ using the standard relation dΨ = ⋆dA. (3.25) A direct calculation shows that a two-form satisfying the ansatz (3.18), (3.14), (3.23) and differential equations (3.24) maps into a separable scalar: As a consistency check, one can verify that functions (X,Φ) indeed satisfy the ordinary differential equations (2.21), which follow from the wave equation for the separable scalar (3.26). This concludes our discussion of a two-form field in even dimensions, and in the next subsection our results will be extended to the odd-dimensional case.

Odd dimensions
The analysis presented in the last subsection can also be carried out for the two-form potential in the odd-dimensional Myers-Perry geometry (2.12). Instead of repeating the algebraic manipulations, we just outline the logical steps: (1) As in section 3.1 we begin with a special case (3.3) of the general two-form (3.2). Then manipulations presented in the Appendix A.2 lead to the most general separable ansatz (3.8), but now there is a new set of projections: n µ A µν = 0. Recall that the frames for the metric (2.12) were introduced in (2.15).
(2) Once the condition (3.3) is relaxed, three types of components (3.5) become mixed in the equations of motion (3.1), so consistency of the separable ansatz (3.8) requires the identification (3.6). Furthermore, by combining various components into counterparts of the relation (3.11), and by requiring separability in the resulting projections, we conclude that the limit (3.3) of the general separable solution must have the form (3.12): This ansatz leads to an odd-dimensional counterpart of equations (3.13): (3) Extension of the solution (3.28) beyond the special case (3.3) imposes consistency conditions (3.16) that lead to factorization (3.17). In the odd-dimensional case, introduction of angular dependence also leads to specific nontrivial expressions for n µ A µν . The final results for the unique separable ansatz and for the ODEs governing the dynamics of the two-form are written below.
The three steps outlined above lead to the full set of components of the gauge potential Various ingredients of function Ψ satisfy the system of ordinary differential equations Functions appearing in these expressions are given by (3.21)- (3.22), and the constants (A , Ω) are defined in (2.30). As expected, solutions of (3.32) contain d − 1 = 2n free separation constants: (i) (n + 1) parameters (ω, n i ); (ii) (n − 3) coefficients of the polynomial P n−4 ; (iii) two ratios (e 2 /e 1 , e 3 /e 1 ) appearing in the expressions (3.14) for the functions h I ± . Furthermore, as in the case of the Maxwell's equations, different ranges of the ratios (e 2 /e 1 , e 3 /e 1 ) lead to the correct number of independent polarizations of the two-form.
The counting presented above fails for d = 5, which corresponds to n = 2. In the last subsection we have encountered a similar situation in the four-dimensional case, where the ratios (e 2 /e 1 , e 3 /e 1 ) were not independent, but rather they were determined by (ω, n). For the five-dimensional black hole there are two options: (a) The general ansatz (3.30)- (3.31) can be kept, as long as three parameters (ω, n 1 , n 2 ) are constrained by Then the differential equations are given by a simpler version of (3.32): The second lines in both equations (3.32) disappear due to the constraint (3.33) and the explicit form of the product H i W i in five dimensions. As expected, there are four separation constants: two ratios (e 2 /e 1 , e 3 /e 1 ) and two angular momenta (n 1 , n 2 ). The frequency ω is fixed by the constraint (3.33).
(b) If the constraint (3.33) is not satisfied, then the ansatz (3.31) becomes inconsistent unless e 2 = 0: The differential equations (3.32) are replaced by Here we defined Again, there are four separation constants: three unconstrained parameters (ω, n 1 , n 2 ) and one ratio e 1 /e 3 .
In five dimensions, a two-form A (2) can be dualized into a vectorÃ (1) using a standard relation The resultingÃ (1) turns out to have a separable form (2.28), and the expression for its "master function"Ψ in terms of Ψ is similar to (3.26)- (3.27).
To summarize, in this section we have demonstrated separability of equations for the twoform potential in all dimensions and found the resulting systems of ODEs that govern the dynamics. In even dimensions, the ansatz and the ODEs are given by (3.18) and (3.19), while the odd-dimensional results are (3.30) and (3.32). In the degenerate cases of four and five dimensions, the separable configurations are given by (3.23)-(3.24) and (3.34)- (3.35). The differential equations contain the correct number of separation constants, and all polarizations are covered by various ranges of the these constants.

3-form potential in the Myers-Perry geometry
In this section we will analyze a three-form potential and construct the most general separable solution of the equations of motion Since the procedure is very similar to the one implemented in the last section, we will only outline the logic and present the results for the most crucial steps. As before, we will have separate discussions for even and odd dimensions.

Even dimensions
In this subsection we construct the most general separable solution of equations (4.2) in the Myers-Perry geometry (2.1). Using the results for the Maxwell field and for the two-form as an inspiration, we impose an ansatz with some unknown functions F and h I ± (x I ). Indices (α, β, γ) in (4.3) are not summed over, and they take values (+, −). The frames m Iµ α are given by (2.10). For separable solutions, function Ψ must have the form As in the case of the two-form, we begin with the special case ω = n i = 0 (4.5) to determine the functions F [x I , x J , x K ] and h I α , then we switch on the angular dependence to find the differential equations for (Φ, X j ). Specifically, we perform the following steps: (1) In the special case (4.5), the ansatz (4.3)-(4.4) and its counterparts for higher forms are analyzed in Appendix A.3, and the result for the three-form potential reads Here polynomials Q I are defined by and coefficients (p k , b k ) are the same for all values of I. Function Ψ satisfies a system of ordinary differential equations (A.71) where P n−4 [y] is an arbitrary polynomial of degree (n − 4) in y. The system (4.6)- (4.8) gives the most general separable solution (4.3) of equations (4.2) in the special case (4.5).
(2) Introduction of angular dependence (i.e., relaxation of the requirements (4.5)) leads to severe constraints on the coefficients (p k , b k ). In particular, the roots of two equations must also be the roots of Q I = 0. A similar situation has already been encountered for the electromagnetic field and for the two-form potential, and in the present case the constraints lead to the final expression for h I ± : . (4.9) This relation is the counterpart of the expression for h I ± from (3.18), although now we have rescaled the gauge potential to set e 1 = 1.
(3) In the case of the three-form potential, there are further constraints on parameter q, and the separable solution (4.3)-(4.4) exists if and only if q satisfies an algebraic equation This constraint is automatically satisfied in the special case (4.5), but any nontrivial angular dependence leads to n solutions 12 for the parameter q. Once the constraint (4.10) is imposed, equations (4.2) reduce to a system of ODEs (4.12) for the function Ψ. Note that the constraint (4.10) is a new feature of the three-forms that has not been encountered in p = 1, 2 cases, and in sections 5, 6 we will show that the counterparts of the restriction (4.10) have to be imposed for the p-forms with p > 3 as well.
To summarize, we have found that the most general separable solution of equations (4.2) has the form (4.11) where functions h I ± are given by (4.9), and parameter q satisfies the constraint (4.10). The differential equations for functions (X j , Φ) are 13 Here we used various functions defined in (3.21), as well as some new ingredients: This concludes our discussion of three-form potentials in even dimensions.

Odd dimensions
In odd dimensions, the metric and the frames are given by (2.12), (2.15). Repeating the analysis presented in the last subsection, we arrive at the most general separable ansatz for the three-form satisfying equations (4.2): Functions h I α are still given by (4.9), where parameter q satisfies an algebraic equation (4.10). For function Ψ that has the form (4.4), the dynamical equations (4.2) are equivalent to a system of ODEs: Here we used functions (H x , H r , ∆, W x , W r , G x , G r , W ⋆ , H ⋆ , G ⋆ ) defined in previous sections, (D x , D r ) given by (4.13), and two new factors (N x , N r ): Unlike the explicit equations (3.19), (3.32), (4.12) encountered earlier, the system (4.15) contains some unspecified function F , and unfortunately we were not able to find a closed form expression for F that holds in all dimensions. Let us list the properties of F which are sufficient for writing equations (4.15) in all d ≤ 11 that are of interest for supergravity.
(a) Function F maps polynomials in x into polynomials.
(b) On the space of polynomials, function F is linear: Here we defined convenient building blocks P k : .
Note that analogues of function F have appeared in equations (3.19), (3.32), (4.12), but we found closed-form expressions for them. For example, the first equation in (4.12) can be rewritten as where functionF is defined bỹ (4.20) ThisF satisfies the properties (a)-(b) listed above. Unfortunately we were not able to find a simple counterpart of (4.20) for the function F appearing in (4.15), so we have to rely on the case-by-case construction listed in the item (c). Relations (4.18)-(4.19) are sufficient for writing the ODEs (4.15) in d ≤ 13 dimensions. Although equations (4.14) give the most compact form of the ODEs corresponding to a separable three-form in odd dimensions, it is convenient to rearrange some terms in (4.14) to inspire the extension to higher forms discussed in section 6. It is sufficient to look only at the alternative form of the equation for X: Here G is defined a linear map between polynomials of x: and properties of G are similar to those of F . In section 6 we will propose a conjecture for an extension of equation (4.21) to all p-forms in odd dimensions.
To summarize, in this section we have demonstrated separability of equations for the three-form potential. In even dimensions, the most general separable ansatz is given by (4.11), and it leads to the system of ODEs (4.12). In odd dimensions the results are (4.14) and (4.15).

Four-form potential in ten dimensions
The main motivation for studying p-form potentials comes from string theory, which contains dynamical fields with p ≤ 4. In previous sections we have analyzed equations of motion for p = (1,2,3), and we will now study the last important case: the four-form potential. Unfortunately, differential equations describing separable configurations of this field in arbitrary dimensions are rather complicated, so in this section we will focus only on the physically relevant situation of ten dimensions and find the full answer for that case. The structure of equations for the four-potential in arbitrary dimensions will be discussed in the next section, which will also cover all p-forms with p > 4.
A four-form potential appears in the ten-dimensional type IIB supergravity [32], so our analysis focuses on space-times of the form where MP d is a d-dimensional Myers-Perry geometry, and T 10−d is a torus. The Ramond-Ramond potential C (4) has four types of components, C µνλσ , C µνλa , C µνab , C µabc , C abcd , where Greek indices correspond to the directions on MP d , and Latin indices cover the torus. Only the first ingredient in (5.2) describes a genuine 4-form on MP d , while other sectors reduce to lower forms. Furthermore, unless d = 10, the field C (4) µνλσ can be dualized to a four-potential that has at least one index on the torus, i.e., to one of the last three sectors in (5.2). Therefore, equations for the four-potential become interesting only for d = 10, and here we will focus only on that case.
In the ten-dimensional type IIB supergravity [32], the five-form field strength F 5 satisfies the self-duality relation and the Bianchi identity 15 To construct separable solutions of these equations using the techniques developed in the previous sections, we begin with the Maxwell-type equation for C (4) ≡ A (4) and build the five-form by adding an exact form dA (4) and its dual 16 : (4) . (5.5) In this section we will focus on solving the first equation.
Using the constructions (2.22), (3.18), (4.11) as inspirations, we impose a separable ansatz where function Ψ has the form In the special case (4.5), the resulting equations (5.4) are analyzed in the Appendix A.3, where it is shown that, for consistency, function F IJKL must have the form 8) and that factors h I ± must be given by Then the Maxwell-type equations (5.4) reduce to a system of ODEs: As in the case of the lower forms, extension of the ansatz (5.6)-(5.9) beyond the special configurations (4.5) leads to certain restrictions on coefficients b k in (5.9). Specifically, we find that factors h I ± must have the form 17 where (q 1 , q 2 ) are two different roots of the algebraic equation (5.12) Note that the same equation (4.10) has been encountered in our discussion of threeforms. Once the ansatz (5.6) with functions (5.8) and (5.11) Here functions (G x , G r , G ⋆ ) are given by (3.21), and functions (F x , D x , F r , D r ) are defined by Equations (5.6)-(5.8), (5.11), (5.13) give the full separable solution of the equation (5.4) for the four-form in ten dimensions.
Although solving equation d ⋆ dA (4) = 0 in arbitrary dimensions is beyond the scope of this paper, the experience with lower forms gained in previous sections suggests that the ansatz (5.6), (5.8), (5.11) should hold for all d ≥ 10. Furthermore, by comparing equations (5.13) with the structure of the systems (2.24), (3.19), and (4.12), we arrive at a reasonable guess for the extension of (5.13) to all even dimensions d = 2(n + 1): 17 To simplify the formulas below, we used the freedom in rescaling A (4) to require h I ± [0] = 1. An analogous choice in (3.18) would have set e 1 = 1. 18 As in previous sections, we write X(x) instead of X j (x j ) to avoid unnecessary clutter.
Function G is a linear operator on the space of polynomials 19 of x, and it will be discussed in more detail in the next section. Equations (5.13) clearly fit the pattern (5.15) with a specific function G , and finding the expression for G in arbitrary dimensions is an interesting open problem.
To summarize, in this section we have demonstrated separability of the equation for the four-form potential and derived the most general separable ansatz (5.6), (5.8), (5.11).
In the ten-dimensional case, we also found the system of ODEs (5.13) that governs the dynamics, and we conjectured an extension of these equations to all dimensions (5.15). In the next section we will discuss additional conjectures for the ODEs describing p-forms with p > 4.

Separability of equations for higher forms
Although supergravity contains only p-forms with p ≤ 4, the construction presented in this article can be extended to cover dynamics of all p-forms in all dimensions. In the Appendix A.3 we demonstrate that equations for such fields in the Myers-Perry geometry are separable, and in this section we summarize the resulting ansatz. We also propose a conjecture for the system of ordinary differential equations governing the dynamics of p-form in all dimensions.
Let us consider a dynamical equation for a p-form field A (p) , d ⋆ dA (p) = 0, (6.1) in the even-dimensional Myers-Perry geometry (2.1). Inspired by the results of the previous sections, we impose an ansatz In the Appendix A.3 we analyze the resulting equations (6.1) in the special case ω = n i = 0. (6.4) and demonstrate that the ansatz (6.2)-(6.3) is consistent if and only if 20 With this choice, the dynamical equations (6.1) reduce to a system of ODEs Here P n−p−1 [y] is an arbitrary polynomial of degree (n − p − 1) in its argument, and it must be the same in all equations (6.6). Relaxation of the condition (6.4) leads to a complicated set of equations, and a full analysis of this system is beyond the scope of this article. However, by comparing the answers (2.22), (3.18), (4.9), (5.11) for p = 1, 2, 3, 4, we conjecture that the ansatz (6.2)-(6.3) still works, but relations (6.5) must be replaced by more restrictive requirements 21 Here {q k } are (p − 2) different roots of the algebraic equation Comparison of the systems (2.24), (3.19), (4.12) leads to the conjecture for the ODEs governing the dynamics (6.1) of the p-form: Here we used the expressions (3.21) and introduced factors (D x , D r ): Equations (6.9) contain a function G p , which depends on p, but not on n, and has the following properties: (a) Function G p maps polynomials in x into polynomials. 21 As in (4.9) and (5.11), we normalized the p-form potential to ensure that h I ± [0] = 1. This avoids unnecessary complications in differential equations (6.9).
(b) On the space of polynomials, function G is linear: (c) Coefficients of the polynomial G p [x m ] depend only on parameters (e 2 , q 1 , . . . q p−2 ) and integers (p, m).
In previous sections we have found explicit expressions for G p in all dimensions for p = (1, 2, 3) (see equations (2.24), (3.19), (4.12)). For the four-form, we looked only at the ten-dimensional case (5.13) that gave limited information about G 4 : for all values of parameters (p, m), which would amount to finding the ordinary differential equations (6.9) for all forms in all even dimensions, is beyond the scope of this paper.
We conclude this section with a brief discussion of the odd-dimensional case. Our results for p = (1, 2, 3) forms suggest that the ansatz (6.2) should modified as Functions F I 1 ...Ip , h I ± , and Ψ are still given by (6.7), (6.3), andF I 1 ...I p−1 is defined aŝ Substitution of the ansatz (6.12) into the dynamical equation (6.1) leads to a system of ODEs that generalizes (2.29), (3.32), (4.15), (4.21): Functions G have the properties (a)-(c) listed on the previous page, and the explicit expressions for G 1 , G 2 , G 3 can be inferred from sections 2.3, 3.2, 4.2.
To summarize, in this section we proposed the separable ansatze (6.2), (6.12) for arbitrary p-forms and conjectured the structures (6.9), (6.14) of the resulting ordinary differential equations. For the special configurations (4.5), the derivation of these results is presented in the Appendix A.3, but full justification of equations (6.9), (6.14) in the general case for p > 4 is an interesting open problem.

Extension to the Kerr-(A)dS geometry
The main motivation for this article comes from the desire to understand the dynamics of p-form fluxes on nontrivial stringy backgrounds. The Myers-Perry geometry, which solves equations of motion of supergravity, is a natural example of such a background, so until now we have focused on excitations of this space. The equations for scalar and vector fields on this geometry have been solved in the past, and interestingly, success in separation of variables for such excitations extends to rotating black holes with cosmological constant [15,25]. This turns out to be true for the higher forms as well, but so far we have focused only on vanishing cosmological constant since, as we will explain below, this seems to be the only physically interesting case. For completeness, in this section all results are extended to black holes with cosmological constant, even though physical applications of equation (6.1) on such spaces are not clear.
In certain limits, string theory reduces to supergravities in ten and eleven dimensions, and light excitations of such theories contain various p-forms. In particular, the D = 10, 11 supergravities admit solutions of the form 1) and in this article we have studied the p-form excitations of such spaces. Since the geometry (7.1) is a pure metric, linearized equations for all fluxes have the form (6.1), so all results obtained in previous sections are directly applicable to string theoretic excitations of the geometry (7.1). In addition to backgrounds (7.1), supergravities also admit solutions of the form 22 One can try to combine the structures appearing in (7.1) and (7.2) by replacing the AdS factor in the last equation by a black hole with a negative cosmological constant. Such "Myers-Perry-AdS" solutions were constructed in the article [21], and we will refer to them as GLPP geometries, whose explicit metrics will be written below. Therefore in string theory it is very natural to study excitations of spacetimes Note that, unlike the pure metric (7.1), the geometry (7.3) necessarily contains fluxes, which modify the equation (6.1) for the p-forms, making separation of (6.1) with p > 1 irrelevant for studying supergravity excitations of (6.1). Let us explain this in more detail.
The dynamics of scalars and vectors on the GLPP geometry was analyzed in [15,25], where it was shown that the relevant equations separate, just as in the case of the Myers-Perry black hole. Furthermore, a wave equation for a scalar on the full space (7.3) reduces to a Klein-Gordon equation on GLPP, and similar reduction occurs for Maxwell's equations, so the results of [15,25] lead to full integrability of dynamical equations for scalar and vector fields on (7.3). Unfortunately, equations of motion for higher forms depend on the structure of background fluxes, and they do not have the universal form (6.1). For example, there are two two-form excitations (B (2) , C (2) ) of the space 7.4) and linearized equations of supergravity mix them [32,34]: (2) , d ⋆ dC (2) = F 5 ∧ dB (2) . (7.5) This system is more complicated than a single equation d ⋆ dA (2) = 0, which appears to be irrelevant for analyzing supergravity excitations of (7.4). Nevertheless in this section we will study a formal equation d ⋆ dA (2) = 0, as well as its generalization (6.1), in the GLPP geometry and demonstrate full separation of variables, extending the results from the previous section. The analysis of the physical equations, such as (7.5) on (7.4), is an interesting open problem.
To demonstrate separability of equations (6.1) in the GLPP geometry, we should begin with giving the explicit form of the metrics. As in the Myers-Perry case, the odd and even dimensions must be treated separately. The geometries were constructed in [21], and the special frames that generalize (2.7) and (2.13) were found in [15,16,17]. In even dimension one gets The notation is the same as in (2.7), with two new ingredients: Parameter L introduced in [21] is related to the value of the cosmological constant, and expressions (2.7) are recovered for L = 0. The counterparts of the null vectors (2.10) are defined by To solve the equation (6.1) in the even-dimensional GLPP geometry, we introduce the ansatz (6.2), (6.3), (6.7)- (6.8). This results in the generalization of the ODEs (6.9), Note that the cosmological parameter L appears only in the terms (7.10) while all other ingredients of equations (6.9) remain the same as in the Myers-Perry case. Modifications (7.10) can also be easily implemented in the more explicit equations (2.24), (3.19), (4.12), (5.13) for p = (1,2,3,4).
In odd dimensions, the GLPP metric [21] can be written in terms of separable frames [15,16,17,25] 7.11) which are obtained by inserting the factors (Q r , Q i ) given by (7.7) into (2.13). The separable ansatz for the p-forms ((6.12), (6.13), (6.7), (6.3)) is built using the null vectors (7.8) 23 , as well as n = e ψ . Substitution of the resulting form into the dynamical equation (6.1) leads to a system of ODEs This is a very simple modification of equations (6.14), and similar insertions of the Qfactors should be made in the special cases (2.29), (3.32), and (4.15).
To summarize, in this section all results obtained in our article have been extended to the GLPP geometries, which generalize the Myers-Perry solutions to black holes in the presence of cosmological constant. If would be very interesting to find further extensions to equations like (7.5), which describe the dynamics of p-forms in the spaces with effective cosmological constants generated by fluxes. Such extensions will provide valuable tools for studying excitations of black holes in string theory.

Discussion
In this article we have demonstrated separability of the dynamical equations describing all p-forms in the backgrounds of the Myers-Perry and GLPP black holes, generalizing the results known for the scalar and vector fields. The most general separable ansatze are given by ((6.2), (6.7)) in even and by (6.12)-(6.13) in odd dimensions. The structures of the resulting ordinary differential equations are summarized in (6.9) and (6.14), and in the special cases of p = (1,2,3,4), the explicit ODEs are given by (2.24), (2.29), (3.19), (3.32), (4.12), (4.15), (5.13). Our procedure relied on existence of the Killing-Yano tensors, and very restrictive uniqueness theorems for such objects [33] imply that it would be hard to extend the methods presented here beyond the Myers-Perry and the GLPP geometries, if one focuses only on vacuum solutions of Einstein's equations. However, one should be able to incorporate charged black holes since they also admit Killing-Yano tensors [20]. Unfortunately, equations for interesting excitations of such spacetimes mix various fluxes (see the discussion around formula (7.5)), so we leave exploration of these modes to future work. Another interesting direction involves the study of integrable excitations which are not covered by separable ansatze. Although separability of dynamical equations guarantees integrability, the inverse relationship is not true. In particular, there are several examples of gravitational backgrounds where the equation for the scalar field is fully integrable, while there is no separation of variables [35]. It would be interesting to see whether a similar phenomenon occurs for the higher forms. Finally, it would also be very interesting to extend the framework introduced this article to gravitational waves.

Acknowledgments
This work was supported in part by the DOE grant DE -SC0017962.

A Derivation of the separable ansatz
This appendix is dedicated to justifying the ansatze (3.18), (4.11), (6.2), (6.7). Specifically, we focus on configurations which do not depend on the cyclic coordinates and determine the relations between various components of the gauge potential and the "master" scalar function. The extensions of these relations to configurations with general angular dependence are straightforward, but the intermediate formulas in the derivations are rather cumbersome, so we just write the final answers in the appropriate sections 24 .
We begin with reviewing the results for electromagnetism by presenting a modification of the analysis performed in [25] that can be easily extended to higher forms. Then in section A.2 we derive the ansatz (3.18) for the two-form, and in section A.3 the results are extended to all p-form potentials.

A.1 Maxwell's equations
In this subsection we will justify the ansatz (2.22) for special configurations with Although the result (2.22) has been already derived in [25] using a slightly different method, it is the logic presented here that can be easily extended to higher forms, so it is instructive to begin with applying the new method to the Maxwell field before extending it to higher forms in sections A.2 and A.3.
We begin with studying Maxwell's equations in even dimensions, and some peculiarities of the odd-dimensional case will be discussed after equation (A.38). Since the frame [m (I) ± ] µ defined by (2.10) is associated with coordinate x I , a reasonable requirement for separation of variables is 25 with some unknown functions (S r,± , S i,± ). For Ψ that does not depend on (t, φ k ), the relation above simplifies In the absence of (t, φ k )-dependence the Maxwell's equations for (A r , A x i ) decouple from the ones for (A t , A φ i ), so it is convenient to analyze these sectors separately. Similar decoupling will take place for the higher forms discusses in sections A.2, A.3. For the Maxwell field, we will analyze the (A r , A x i ) components in subsection A.1.1 and focus on (A t , A φ i ) components in section A.1.2.

Let us consider the Maxwell's equations involving
A r and A a ≡ A xa , a = 1 . . . n.
Recalling the frames (2.10), we find The cyclic components of the gauge potential, (A t , A φ k ), are more complicated, but they do not mix with (A.4) in the Maxwell's equations. In four dimensions, n = 1 and configurations (A.4) describe a pure gauge (see [25] for details), therefore here we focus on n ≥ 2.
Using indices (i, j, . . . ) to cover r and x a and defining a function we find the expressions for the gauge potential and the field strength in the sector (A.4) Furthermore, evaluating the determinant of the Myers-Perry metric (2.1) in coordinates (r, x a , t, φ a ), we write the i-th component of Maxwell's equations as Simplifications lead to the final system of equations j =i Separable ansatz in d = 2(n + 1) dimensions implies that function Ψ can be written as then the set of equations (A.8) reduces to a system 26 , and x i -dependence of this function is still undetermined. Substitution of (A.10) back into (A.8), leads to a system of algebraic relations for P (i,j) n−2 : j =i The left-hand side of the last relation is a meromorphic function of all x j for j = i, then analysis of the asymptotic behavior and residues at various poles demonstrates that polynomials P (j,i) n−2 must be the same for all values of j: Once the relation (A.12) is satisfied, equation (A.11) becomes an identity. To summarize, for the separable ansatz (A.9), the Maxwell's equations (A.8) are equivalent to the system of ODEs (A.10) with restrictions (A.12). Furthermore, as in the scalar case (2.19), the polynomials P (i) n−2 contain some number of separation constants, in particular, setting these constants to zero, one can get a special case of the system (A.10) with P (i) n−2 = 0: The current discussion is applicable to n > 2, then taking j = 0 and combining equations (A.13) with i = 1, 2, we find In particular, this implies that .14) This relation and its counterparts with other values of (i, j) in (A.13) can be summarized as a system of equations with some function g(x). The analogues of the functional equations (A.15) will play the crucial role in derivations of separable ansatze for higher forms, which will be presented in sections A.2, A.3. Explicit evaluation of the derivatives in (A.14) gives The most general solution of this differential equation for function S is with arbitrary constants (α, β, γ, δ). Although the result (A.17) was derived for the special choice of separation constants corresponding to reduction of (A.10) to (A.13), it is clear that function S cannot depend on the separation constants, so relations (A.17) also hold in the general case (A.10). Substitution of (A.17) into (A.10) with (A.12) gives .18) This equation implies that where P n−2 (x 2 j ) is a polynomial with constant coefficients. To summarize, we have demonstrated that the polarization (A.4) is consistent with the separable ansatz (A.9) if and only if (S r , S a ) are given by (A.17), and functions X j satisfy ordinary differential equations .20) Let us now discuss the equations for (A t , A φ i ).

A.1.2 Angular components of the gauge field
The angular components of the gauge field are determined from combinations of ( and substituting the ansatz (A.21) into Maxwell's equations, we find (n + 1) independent equations: It is convenient to look at projections of M µ on various basis vectors, and straightforward calculations give We will focus on n scalar equations, + ] µ M µ can be analyzed in the same way. To use the differential equation (A.20), it is convenient to rewrite the derivatives appearing in the last term of (A.24) as

Then equation (A.24) becomes
[m Here we defined .

(A.26)
Substituting the separable ansatz (A.9) into equation (A.25), one arrives at several consistency conditions, in particular, all W C,i must vanish, and all ∂ i derivatives must disappear from the last line of (A.24). The latter condition leads to differential equations for functions X i : (A.27) while equations W C,i = 0 allow one to determine the functionsS i . Specifically, defining Treating this relation as an algebraic equation for Q 1 , we find that .29) for i = 1. Equation W 2,1 = 0 leads to a similar relation for Q 2 , and substituting the result into (A.28), we conclude that coefficients (c 1 , c 2 ) can be arbitrary, but they must be the same for all values of i. This leads to the expressions forS i and to a more explicit form of equation (A.27): The coefficient in front of ∂ C Ψ vanishes if and only if Let us define a new set of functionsP Rewriting equation (A.32) in terms of these functions and recalling the expression (2.8) for d k , we conclude thatP (k) n−1 must be a polynomial of degree (n − 1) in its argument. Furthermore, taking the limit x i → x C in (A.32), we find thatP n−1 must be the same for all values of k, then Substitution of the last formula into equation (A.32) leads to a relatioñ which is identically satisfied for an arbitrary polynomial 27P n−1 . Thus functions P (k) must have the form (A.34), and relation (A.31) reduces to To avoid pole at x 2 C = −c 2 /c 1 , we must require thatP n−1 [−c 2 /c 1 ] = 0, then relation (A.34) can be rewritten as As expected, this "angular" polarization and its "radial" counterpart (A.20) lead to the same system of ODEs for X k = (R, X a ). Although the derivation of (A.38) turned out to be somewhat tedious, we note that the main results pertaining to the Maxwell field, equations (A.20) and (A.17), were obtained using very transparent calculations, and it is this new derivation that will be easily generalized to higher forms in the remaining part of this appendix.
We conclude the discussion of the Maxwell field by making comments about odd dimensions. Substituting the ansatz (A.4) into Maxwell's equations, we arrive at a counterpart of (A.8): Additional factors of x j are caused by different expressions for the frames and the determinant of the metric. The rest of the derivation proceeds as before, leading to (A.17) and the counterpart of (A.20): For the angular polarization we again find (A.21) and the expression forS i from (A.38). Furthermore, the static configurations (A.1) must have n µ A µ = 0.
Although all results of this subsection have already been obtained in [25], the new derivation presented here can be easily extended to higher forms, and such generalizations are carried out in the next two subsections.

A.2 Two-form potential
Let us now extend the construction presented in the last subsection to the two-form potential. As in the case of Maxwell's field, we will look for separable solutions of linear equations d ⋆ dA = 0 in the absence of angular dependence. We will mostly focus on the x-components of the gauge field, but as we explicitly saw in the last subsection, the resulting equations would apply to other components as well. The verification of this fact is an straightforward extension of the arguments presented in the Appendix A.1.2, although the intermediate formulas become rather complicated.
We begin with looking at x-components of a two-form field in even dimensions and imposing a natural generalization of the ansatz (A.4) 28 Evaluating the field strength, (A.42) and substituting the result into equation d ⋆ F = 0, we arrive at a counterpart of equation (A.8): 28 Recall that in this Appendix indices (i, j, k) cover n coordinates x a , as well as x 0 = −ir.
Simplifications of the last relation leads to a system of equations labeled by indices i and k: Imposing a separable ansatz for function Ψ, (A.44) one arrives at a system of differential equations The last relation should become a second-order ODE for X j (x j ), which means that, upon multiplication by an overall function, the coefficients in front of (∂ 2 j X j , ∂ j X j , X j ) must become functions of x j only. In particular, this implies the following separation of variables Recalling the definition of S ijk given by (A.42), the last relation can be written in a more symmetric form Although a priori we should have used g (i,j,k) (x) instead of g(x) in the last equation, the final ODE (A.47) ensures that all g (i,j,k) (x) must be the same, so we avoided unnecessary complications already in (A.45). The functional equation (A.45), a counterpart of (A.15), is the main consistency condition of the separable ansatz for the two-form potential. Substitution of (A.45) into the field equations (A.43) gives Imposing the separable ansatz (A.44) and recalling that d j is a polynomial of degree 2n, we arrive at the system of ordinary differential equations where P n−3 is a polynomial of degree (n − 3). 29 In particular, for n = 2 equation (A.47) does not contain the last term. Furthermore, for n = 1, which corresponds to four 29 To see this, we observe that, upon multiplication of (A.46) by m ) the terms that do not contain ∂ m become polynomials of (n − 3) in x 2 m .
dimensions, the logic outlined here breaks down, and this degenerate case is discussed in the end of section 3.1.
Let us now solve the functional equations (A.45) to determine functions g(x) and S ij (x i , x j ). Introducing new variables y k = x 2 k and definingg k =g(y k ) ≡ g( √ y k ), we can write some integrability conditions of (A.45) as (no summations) Treating this relation as a first-order differential equation forg i , we find g(y) = 1 with some constants (a, b, c). Substitution of this relation back into (A.48) yields an identity. Let us now look at another integrability condition of (A.45): Integration of this relation gives Substitution of this formula into equation (A.45) produces restrictions on (A ij , B ij ): This leads to the final formula for the most general S ij : where constants C ij satisfy a system of homogeneous equations With this assignment the functional equation (A.45) becomes an identity.
To summarize, we have found the most general solution of the functional equation (A.45): Here the set of constants {C ij } satisfies the relations (A.53). Note that these constants, as well as A(y), don't contribute to S ijk and to the field strength, so the physical degrees of freedom are described by a simpler version of equation (A.54): To turn on (ω, n i ), we would need a slight generalization of the last relation. Using (A.55) as an inspiration, we require a multiplicative separation in S ij : Then consistency with (A.54), leads to differential equations for function f , and the most general solution is 30 We conclude that the most general separable expression for (g, S ij ) is a special case of (A.54): In section 3 we demonstrate that this ansatz works for nontrivial (ω, n i ) as well. To summarize, we have shown that the most general "radial" polarization (A.41) with separable function Ψ leads to a system of ODEs (A.47). Furthermore, S ij and g must be given by (A.59).
One can also analyze the "angular" polarization by modifying the relevant parts of section A.1.2. This logic of such extension from the Maxwell field to the two-form potential is rather straightforward, but the intermediate formulas are somewhat messy, so here we just quote the result. Starting from an ansatz for the cyclic components 31 of the gauge field A ab , 30 Strictly speaking, relation (A.57) imposes an additional linear constraint on (p, q). However, this constraint can be eliminated by a rescaling of function g and polynomial P n−3 that introduces a multiplicative constant in the right-hand side of equation (A.45) without affecting the ODE (A.47). Thus we can keep four free parameters (a, b, c, p, q) in (A.58). 31 Recall that for cyclic components A ab , indices (a, b) take values (t, φ i ).
with a separable function Ψ, (A.44), and substituting this A ab into the field equations, we find the system of ODEs (A.47) with Although a priori polynomials Q j and Q (2) Function Ψ is given by (A.44), and its ingredients satisfy a system of differential equations (A.47) with This ansatz is the starting point for our discussion of separable two-form potentials in section 3. The analysis of odd dimensions follows the same route, and the final result is presented in section 3.2.

A.3 Higher forms
Let us now extend the results obtained in the last subsection to p-form potentials with p > 2. We will mostly focus on the "radial polarization" (A.41) in even dimensions, and comment on the counterparts of (A.60) and (A.63) in the end of this subsection. The extension to the odd-dimensional case is straightforward, and it is discussed in sections 4.2 and 6. The generalization of the ansatz (A.41), (A.44) is Here S i 1 ...ip is antisymmetric its indices, and we can also introduce an antisymmetric tensor S i 1 ...i p+1 by generalizing the definition (A.42): Then the field equations, d ⋆ dA = 0, lead to a counterpart of the system (A.43): This gives one equation for every set of indices (i 1 , . . . , i p ). Imposing the separable ansatz (A.66) for Ψ, and repeating manipulations that led to (A.45) and (A.46), we can rewrite the system (A.68) as where tensor S i 1 ...i p+1 and function g are subject to constraints This functional equation for S i 1 ...ip and g is a generalization of the constraints (A.15) and (A.45) encountered for the one-and two-forms. As before, equations (A.69) with separable function Ψ lead to a system of ODEs and all these equations contain the same P n−p−1 , a polynomial of degree (n − p − 1) in its argument. Equations (A.71) apply only for n ≥ (p + 1), and degenerate cases of small n can be analyzed separately. leading to an integrability condition for the system (A.70): Treating the last relation as a first order ODE for the function g(x i 1 ), we find the most general solution of (A.73): Substituting (A.74) into the functional equation (A.70) and solving for S i 1 ...ip , we find 32 whereŜ i 1 ...ip satisfies the homogeneous equation Taking various derivatives of this relation and using symmetries of the tensorŜ, we find the most general solution of (A.76): Here C i 1 ...ip are arbitrary constants satisfying the homogeneous equation Using the separable ansatz (A.56) as an inspiration, we focus on a special case of (A.75): We have already found an example that corresponds toŜ i 1 ...ip = 0: and now we will determine the most general function f that solves equation (A.70). Since expressions S i 1 ...ip must have the form (A.75) with (A.77), we conclude that the ansatz (A.79) solving equation (A.70) must satisfy an integrability condition Viewing the left-hand side of this relation as a holomorphic function of x i 1 , we conclude that the poles of f (x i 1 ) must coincide with zeroes of Q i 1 . Furthermore, function f must have the form To simplify the discussion, we assume that all roots of Q are distinct, but the extension to degenerate roots is straightforward. In the vicinity of y ip = c p , the leading contribution to the left-hand side of (A.83) is proportional to (y ip − c p ) −2 , and equation (A.83) can hold only if the residue at this pole vanishes: Coefficients h I k α k are rational functions of x J : As demonstrated in sections 4, 5 and 6, the ansatz (A.91) can be easily extended to configurations with nontrivial (ω, n i ), although we managed to derive the resulting ODEs only on a case-by-case basis.
To summarize, in this Appendix we have derived the most general separable ansatze for arbitrary p-form fields in the Myers-Perry geometry. Our final result is summarized by (A.91)-(A.92), and differential equations produced by this ansatz are discussed in the main body of this article.

B Divergence of the two-form potential
In the Appendix A we derived the separable ansatze for the p-forms and outlined the derivation of the resulting ordinary differential equations for the "master field" Ψ. In the presence of angular dependence, the analysis of the equation d ⋆ dA (p) = 0 (B.1) was somewhat messy even in the simplest case of the Maxwell field, i.e., for p = 1. After such a derivation had been performed in [25], the authors of [26] observed that the resulting field A (1) satisfied the Lorentz condition Then by imposing the ansatz (2.22) introduced in [25] along with the Lorentz gauge (B.2), the authors of [26] found an elegant shortcut for arriving at (2.24), (2.29) which simplified the original derivation of [25]. In this Appendix we explain the origin of this shortcut and demonstrate that, unfortunately, it cannot be directly extended to higher forms.
Let us begin with the discussion of the one-form A ≡ A (1) in even dimensions. The ansatz (2.22), is inconsistent with the dynamical equations (B.1) for the separable configurations (6.2)-(6.3) with p > 1. To avoid unnecessary complications, we will focus on a special class of solutions although the conclusion about satisfaction or failure of the Lorentz gauge holds in the presence of (t, φ i ) dependence as well.
In the special case (B.5), the non-cyclic components of the Maxwell field (B.3) become 34 Substitution into the left-hand side of (B.2) gives Here we used the expressions for g ij from (2.11), as well as equation (A.7) for the determinant of the metric. As we have shown in the Appendix A, the Maxwell's equations for the separable configurations (B.3) are equivalent to the system of ODEs (2.24), then in the present case ω = n i = 0 we find (no summation) As we saw in the Appendix A.1.1, the right-hand side of the last equation vanishes identically, so we have shown that any separable solution (B.3), (B.5) of the Maxwell's equations satisfies the Lorentz gauge ∇ µ A µ = 0. 34 As in the Appendix A, here we use the notation A j = (A r , A xa ).
The conclusion holds even for a more general function E given in (2.18), but the proof requires longer algebraic manipulations, and we will not present them here.
Let us now look at the separable two-form potential (3.18): Once again we will focus on the special configurations whose "master function" has the form (B.5). The non-cyclic components of the gauge field are given by Following the steps used to evaluate ∇ µ A µ , we find the non-cyclic components of the divergence of the two-form: Here we used the ODEs (3.19) for function Ψ, as well as identities (A.8)-(A.10) for the polynomial P n−2 . The last equation can be summarized as Although this expression is very compact, it is clear that the two-form does not satisfy the Lorentz condition (B.4). The cyclic components of the divergence can be computed in a similar fashion, and the result reads Once again, this expression does not vanish. The expressions (B.10), (B.11) can be extended to configurations with nontrivial (t, φ i ) dependence and generalized to higher forms, and although the explicit formulas are more complicated, the conclusion remains the same: the separable configurations (6.2) satisfying the equations of motion (6.1) do not belong to the Lorentz gauge (B.4) for any p > 1. Thus existence of an elegant shortcut for deriving the ODEs (2.24) proposed in [26] seems to be a peculiarity associated with Maxwell's equations, and for higher forms one has to follow the longer route discussed in the Appendix A.