A classification of supersymmetric Kaluza-Klein black holes with a single axial symmetry

We extend the recent classification of five-dimensional, supersymmetric asymptotically flat black holes with only a single axial symmetry to black holes with Kaluza-Klein asymptotics. This includes a similar class of solutions for which the supersymmetric Killing field is generically timelike, and the corresponding base (orbit space of the supersymmetric Killing field) is of multi-centred Gibbons-Hawking type. These solutions are determined by four harmonic functions on $\mathbb{R}^3$ with simple poles at the centres corresponding to connected components of the horizon, and fixed points of the axial symmetry. The allowed horizon topologies are $S^3$, $S^2\times S^1$, and lens space $L(p, 1)$, and the domain of outer communication may have non-trivial topology with non-contractible 2-cycles. The classification also reveals a novel class of supersymmetric (multi-)black rings for which the supersymmetric Killing field is globally null. These solutions are determined by two harmonic functions on $\mathbb{R}^3$ with simple poles at centres corresponding to horizon components. We determine the subclass of Kaluza-Klein black holes that can be dimensionally reduced to obtain smooth, supersymmetric, four-dimensional multi-black holes. This gives a classification of four-dimensional asymptotically flat supersymmetric multi-black holes first described by Denef et al.


Introduction
Black holes are in the focus of gravitational research. In four-dimensional vacuum gravity, or Einstein-Maxwell theory, asymptotically flat black holes have a surprisingly simple moduli space due to the well-known uniqueness theorems (see e.g. [1,2]). In contrast, higher dimensional general relativity has a much richer structure, and black hole uniqueness does not hold even in the asymptotically flat vacuum case (for review see e.g. [3]), which became clear with the discovery of rotating vacuum black holes with S 2 × S 1 horizon topology, known as black rings [4]. Rather surprisingly, for a range of asymptotic charges, black rings coexist with the spherical Myers-Perry black holes [5], providing an explicit example of non-uniqueness.
Much is known in general about higher dimensional stationary black holes. The topology of the horizon is restricted to be of positive Yamabe type, i.e. they admit metrics with positive scalar curvature [6], which becomes less restrictive as we go to higher dimensions. Further restrictions have been derived in the literature for black holes with an axial symmetry [7]. This assumption was motivated by the rigidity theorem, which states that analytic solutions with rotating black holes must also admit an axial U(1) symmetry [8][9][10]. The topology of the domain of outer communication (DOC) is also restricted by topological censorship. In the asymptotically flat case it is required to be simply connected [11]. This has been generalised to asymptotically Kaluza-Klein spacetimes [12], where the quotient space of the DOC by the symmetry group corresponding to translations in the compact dimensions must be simply connected. In cases when the spacetime admits biaxial U(1) 2 symmetry, a uniqueness theorem for black holes has also been established [13][14][15], however it has been conjectured (for the vacuum case) that solutions with fewer symmetries must exist [16]. Evidence for the existence of such solutions has long been gathering in the literature [17,18], but the first explicit examples have been constructed only recently in five-dimensional supergravity [19].
Black hole non-uniqueness is also present in higher dimensional supergravity theories, even among supersymmetric solutions. For five-dimensional minimal supergravity the first known black hole was the BMPV solution [20] with a spherical horizon. Later a black ring solution [21], and concentric black ring/black hole solutions [22] were also found. Asymptotic charges of the latter one can overlap with those of the BMPV solution, and even more surprisingly, they can have greater total horizon area than the corresponding BMPV solution. This has since been shown for even single-black hole solutions of this theory [23,24]. This finding is particularly puzzling given that the microscopic derivation of Bekenstein-Hawking entropy in string theory for given charges matches the entropy of the BMPV solution [20,25]. This 'entropy enigma' is yet to be resolved, but microscopic counting of entropy provides further motivation to determine the moduli space of supersymmetric black holes.
Although the full moduli space is yet to be explored, this theory is known to admit various black hole solutions. On top of the aforementioned spherical and ring solutions, black holes with lens space horizons L(p, 1) have been constructed with flat asymptotics [24,[26][27][28][29], but later also with Kaluza-Klein asymptotics [30]. An important feature that contributes to the richness of solutions is that this theory admits solitonic solutions, termed 'bubbling spacetimes', which are smooth, horizonless solutions admitting non-trivial topology in the form of non-trivial 2-cycles, supported by magnetic fluxes of the Maxwell field [31,32]. It is also possible to construct black hole solutions with 'bubbling' domain of outer communication, which include examples that contribute to the single black hole entropy enigma mentioned above [23,24,33].
The general local form of a supersymmetric solution of minimal five-dimensional super-gravity is known [34,35]. Using Killing spinor bilinears, one can show that such a solution admits a globally defined causal Killing field, a scalar function, and three 2-forms. When the supersymmetric Killing field is timelike, the solution takes the form of a timelike fibration over a hyper-Kähler base manifold for which the three supersymmetric 2-forms are the complex structures. Furthermore, if the solution admits a triholomorphic U(1) isometry (i.e. an isometry that preserves the complex structures of the base) the base manifold takes the form of a multi-centred Gibbons-Hawking space [36][37][38], and the solution is locally fully determined by four harmonic functions on R 3 [34]. This is the case if the spacetime has biaxial symmetry (U(1) 2 ) [29], or a single axial symmetry which preserves the Killing spinor [19].
Using the aforementioned local results, a number of classification theorems have been proven for supersymmetric black holes. The near-horizon geometries have been determined in [16] (assuming that the supersymmetric Killing field becomes timelike outside the horizon), and the possible geometries are locally isometric to 1 S 3 or S 2 × S 1 . In the asymptotically flat case global results are known as well. First it was shown that a locally spherical black hole with a supersymmetric Killing field that is timelike on the domain of outer communication must be isometric to the BMPV solution [16]. This assumption, however, is restrictive and excludes the majority of the moduli space of black holes (black lenses and black holes in bubbling spacetimes). A common feature of all near-horizon geometries, and in fact all the solutions mentioned so far, is biaxial (U(1) 2 ) isometry. In [29] a classification of asymptotically flat black hole solutions with such symmetry has been achieved. These solutions have Gibbons-Hawking base with the associated harmonic functions having simple poles at collinear centres on R 3 , corresponding to horizon components or fixed points of the triholomorphic U(1) symmetry.
Recently, these results have been generalised to asymptotically flat solutions admitting only a single axial symmetry that preserves the Killing spinor [19]. Similarly to the biaxial case, these solutions are of multi-centred type, with harmonic functions having simple poles at generic (not necessarily collinear) points. This provided the first explicit construction of a higher dimensional asymptotically flat black hole with just a single axial symmetry, confirming the conjecture of [16] for supersymmetric black holes.
In this paper we will focus on asymptotically Kaluza-Klein solutions, which are asymptotically diffeomorphic to a circle fibration over flat Minkowski space. This includes trivial circle products, and also non-trivial fibrations like the Kaluza-Klein monopole. Many explicit Kaluza-Klein black hole solutions of five-dimensional supergravity are known in the literature [30,[39][40][41][42][43][44][45][46][47][48][49][50], furthermore, a uniqueness theorem for non-supersymmetric, biaxisymmetric, spherical black holes of this theory has also been proven [51]. A classification of more general or supersymmetric solutions is not yet available, however. The main motivation of studying such solutions is that they can be dimensionally reduced to four dimensions, which might have more physical relevance for our universe. For five-dimensional minimal supergravity, the bosonic sector of the reduced theory contains gravity, two Maxwell fields, and two scalars, a dilaton and an axion. Previously, an interesting connection has been unveiled between five-dimensional Kaluza-Klein and four-dimensional asymptotically flat black holes of supergravity theories [39,40,52,53]. Five-dimensional supersymmetric (multi-)black holes/rings with Taub-NUT base space correspond to (multi-)black hole solutions of the four-dimensional theory first derived by Denef et al. [54][55][56].
The purpose of this paper is to generalise the results of [19] to Kaluza-Klein asymptotics, and classify all supersymmetric black hole or soliton solutions with an axial symmetry that 1 In [16] T 3 has also been derived as a possible geometry, however it is excluded by [6].
'commutes' with the supersymmetry, i.e. it preserves the Killing spinor. The latter assumption is the supersymmetric generalisation of the usual requirement that the axial symmetry commutes with the stationary symmetry. Furthermore, when considering reductions of the solutions to 4D, as a result of this assumption we obtain supersymmetric solutions in the lower dimensional theory.
The main difference to the asymptotically flat case is that the supersymmetric Killing field is no longer naturally identified with the stationary Killing field. Instead, we assume that the stationary Killing field is a constant linear combination of the axial and supersymmetric Killing fields. This leads to -from the five-dimensional perspective -two qualitatively different classes of solutions depending on whether the supersymmetric Killing field is timelike on a dense submanifold or globally null. The former one is similar to the one found in the asymptotically flat case [19], and its classification is provided by the following theorem (for the full statement see Theorem 2).
Theorem A. A supersymmetric, asymptotically Kaluza-Klein (in the sense of Definition 1) black hole or soliton solution of D = 5 minimal supergravity with an axial symmetry that preserves the Killing spinor and a supersymmetric Killing field that is not globally null must have a Gibbons-Hawking base (on a dense submanifold), and is globally determined by four associated harmonic functions on R 3 which are of 'multi-centred' form, with parameters satisfying a set of algebraic constraints. The centres either correspond to fixed points of the axial Killing field, or connected components of the horizon, each of which has topology S 3 , L(p, 1) or S 2 × S 1 .
Despite the similarity to the asymptotically flat classification, there are technical differences in its proof. The main complication comes from the fact that the DOC is not simply connected in general, hence certain closed 1-forms do not necessarily define the global functions that play a key role in the asymptotically flat proof. In order to overcome this, first we show that the axial Killing field must be tangent to the Kaluza-Klein direction, then after excluding the possibility of orbits with a discrete isotropy group (exceptional orbits), we apply topological censorship to deduce simple connectedness of the orbit space of the axial Killing field [12]. This allows us to define certain U(1)-invariant functions globally. From this point, the proof is almost identical to that of the asymptotically flat case.
The lack of exceptional orbits follows purely from the axial Killing field preserving three linearly independent two-forms (Killing spinor bilinears, which are the complex structures of the hyper-Kähler base in the timelike case). This result, together with the orientability of the spacetime, restricts the possible horizon topologies to S 3 , L(p, 1) or S 2 × S 1 . These are the only allowed 2 Seifert three-manifolds (with orientable fibres and base) that do not contain any orbifold points. Such orbifold points would require the presence of exceptional orbits in the DOC [7], which are ruled out as mentioned above.
Even in the case when the supersymmetric Killing field is timelike on a dense submanifold, the limit of its norm at infinity can be zero. We call this case asymptotically null, in contrast to the asymptotically timelike case, when we normalise the supersymmetric Killing field to have unit norm at infinity. These should not be confused with the globally null case, which we detail next.
The other main class of solutions contains spacetimes on which the supersymmetric Killing field is globally null, for which we have the following result (the full statement can be found in Theorem 3). 2 Only orbifolds with positive Euler characteristic are allowed [7]. Theorem B. An asymptotically Kaluza-Klein (in the sense of Definition 1), supersymmetric black hole or soliton solution of D = 5 minimal supergravity with an axial Killing field W preserving the Killing spinor and for which the supersymmetric Killing field V is globally null has a metric of the form where W = ∂ u , V = ∂ v , and G, Q are harmonic functions on R 3 with simple poles at centres corresponding to connected components of the horizon, each with topology S 2 × S 1 .
To our knowledge these 'null' Kaluza-Klein black ring solutions have not been previously described. As the supersymmetric Killing field does not become timelike outside the horizon, the proof of the near-horizon classification in [16] is no longer valid. Therefore, we derive the near-horizon classification (with no symmetry assumptions) in the null case, and we find that the near-horizon geometry in the null case agrees with the null limit of the timelike case. Similar 'null' solutions have been previously found in [57], which describe static black strings, those, however, cannot be compactified to obtain smooth, asymptotically Kaluza-Klein black holes, hence they are not part of our classification.
It turns out that, if one relaxes a condition on the harmonic functions of Theorem A so that they allow the supersymmetric Killing field to be globally null (which is a priori not clear that one is allowed to do), then one obtains precisely the solutions of Theorem B. Presumably, this is because they have common higher dimensional origin, as both timelike and null solutions can be uplifted -at least locally -to obtain supersymmetric solutions of six dimensional minimal supergravity [58]. A general feature of the six-dimensional solutions is that the supersymmetric Killing field is null everywhere, similarly to the solutions in Theorem B. Interestingly, the cartesian coordinates in both Theorem A and B (those of the Gibbons-Hawking base, and x i in (1), respectively) originate from a Gibbons-Hawking base that is used to construct the six-dimensional solutions. It would be interesting to investigate the classification from a six-dimensional perspective.
In all cases we find that the axial Killing field is tangent to the Kaluza-Klein direction, thus it is natural to consider the reduction of these solutions to four dimensions, which are also supersymmetric. We prove that the obtained four-dimensional solution is smooth on and outside the horizon provided that there are no fixed points of the U(1) Killing field. This is automatically true in the null class, but restricts the timelike class. Conversely, given a four-dimensional supersymmetric, asymptotically flat black hole solution of a certain four-dimensional supergravity with two scalar fields and Maxwell fields, it uplifts to a supersymmetric, asymptotically Kaluza-Klein solution of minimal supergravity in D = 5. A condition for a smooth uplift is that one of the Maxwell fields corresponds to a principal U(1)-bundle, and hence its magnetic charges are quantised. Thus, with this assumption, we obtain a classification theorem of four-dimensional supersymmetric black holes of the theory (for the detailed statement see Theorem 6).
Theorem C. Consider a supersymmetric, asymptotically flat black hole solution of D = 4 N = 2 supergravity coupled to a vector multiplet in which one of the Maxwell fields is the curvature of a connection on a principal U(1)-bundle over the spacetime. Furthermore, assume that the supersymmetric Killing field is timelike on the DOC. Then the solution must belong to the class of multi-black holes derived in [54][55][56].
The structure of this paper and the proof of theorems is as follows. In Section 2 we briefly summarise the general form of such solutions based on [34], and present the local solutions in the timelike and null case separately. In Section 3 we derive that U(1) Killing field must be tangent to the Kaluza-Klein direction, by using that it preserves the Killing spinor bilinears, and showing that any Killing fields must approach that of a flat space with a compact direction. In Section 4 we perform a near-horizon analysis of the null case (following and completing that of [16]). In Section 5 we analyse the structure of the orbit space. In particular, we show that any exceptional orbits are excluded, then using topological censorship [12] we show that the cartesian coordinates defined by the axial Killing field and the two-form bilinears define a global chart on the orbit space. This also restricts the general form of the solution to those with associated harmonic functions having simple poles in both the null and timelike case. In Section 6 we derive the sufficient conditions for smoothness of such multi-centred solutions. In Section 7 we state and prove the main classification theorems (Theorem 2-3). In Section 8 we perform a Kaluza-Klein reduction to four dimensions, and show that the fourdimensional solution is smooth on and outside the horizon provided there are no fixed points of the axial symmetry in 5D. Finally, we use this correspondence to prove a classification theorem for four-dimensional asymptotically flat black hole solutions to minimal supergravity coupled to a vector multiplet (Theorem 6).

Supersymmetric solutions in five dimensions with axisymmetry
The theory in consideration is the bosonic sector of five-dimensional minimal supergravity, given by the action where F = dA. We work in the conventions of [16], so the signature of the metric is 'mostly plus'. The most general local form of a supersymmetric solution has been derived in [34] using Killing spinor bilinears which globally define a function f , a vector field V , and three 2-forms X (i) , i = 1, 2, 3. These satisfy the following algebraic equations.
where ǫ ijk is the Levi-Civita symbol with ǫ 123 = 1. Using the Killing spinor equation, one further obtains In particular, V is Killing, X (i) are closed, and V preserves X (i) and the Maxwell field. Further specification of the local solution depends on whether V is timelike or null in some open region, which we will describe in detail in the following sections.
In this paper we are considering asymptotically Kaluza-Klein spacetimes, for which we use the following definition. to R × Σ 0 , with Σ 0 being a circle fibration over R 3 \ B 3 , where B 3 denotes a 3-ball, and the R factor corresponds to orbits of a timelike Killing field (stationary Killing field), (ii) the metric on this end can be written as g µν =g µν + O(r −τ ) for some decay rate τ > 0 whereL is a positive constant, u 0 and (u i ) 3 i=1 are the pull-back of the cartesian coordinates on R × R 3 , i.e. ∂ 0 is the stationary Killing field, andr := u i u j δ ij ,ψ is a 4π-periodic coordinate on the fibres, and in these coordinates (iii) the components of the Ricci tensor in these coordinates fall off as R µν = O(r −τ −2 ) and its R 3 derivatives as

Remarks.
1. Definition 1 only requires that the spacetime is a circle fibration at infinity. This includes trivial fibrations (standard Kaluza-Klein asymptotics), as well as non-trivial fibrations such as the Kaluza-Klein monopole. The nature of the fibration is determined by subleading terms (in coordinates of Definition 1) in the metric (see later in Section 6.1.1).
2. The fall-off of the components of the Ricci tensor is in fact equivalent (through the Einstein equations) to requiring that the Maxwell field falls off as F ∼ O(r −τ /2−1 ) at spatial infinity. This can be seen by looking at T 00 , which is a positive definite quadratic in the components of the Maxwell field. The fall-off for the third derivative of the metric and the derivative of the Ricci tensor is a technical assumption, to ensure that components of the Riemann tensor fall off as R µ νλκ = O(r −τ −2 ) as shown in Appendix A. Alternatively, one can assume the fall-off of the Riemann tensor directly.
We now list our assumptions. These are, except for asymptotics, equivalent to Assumptions 1 and 2 of [19]. We assume that (M, g, F ) is a solution of (2) such that (i) it admits a globally defined Killing spinor ǫ (i.e. supersymmetric), (ii) the DOC is globally hyperbolic, that is, it admits a Cauchy surface Σ, From global hyperbolicity of the DOC and completeness of V follows that the DOC has topology ⟪M⟫ ≃ R × Σ, where Σ is a smooth manifold, which we can identify with the orbit space of V 3 .
We emphasise that we do not assume that the stationary Killing field coincides with the supersymmetric Killing field, in contrast to the asymptotically flat solutions considered in [19]. Neither do we assume that the axial Killing field is tangent to the Kaluza-Klein direction at infinity. In Section 3 we show that the latter is indeed true, so the U(1) Killing field (and the length of its orbits) is bounded in the asymptotic region. This means that its linear combination with the causal vector field V can be timelike, so assumption (viii) is natural to make. The solution in the coordinates adapted to these Killing fields generally will not be in the rest frame. Indeed, this is the case for the asymptotically Kaluza-Klein black holes considered in [30,39].
Another important difference to the asymptotically flat case is that in general the DOC is not simply connected. Therefore, we cannot a priori assume that closed 1-forms globally define scalar functions. However, in Section 5 we will show that this is indeed the case for the U(1)-invariant one-forms considered during the proof.
An important consequence of assumption (vii) is that all the Killing spinor bilinears are preserved by W , that is The assumption that W preserves F is redundant. We will see that either the timelike or null region is dense in the DOC, and in both cases F can be expressed by U(1)-invariant quantities (17) and (37). Thus, by continuity of L W F , the Maxwell field is U(1)-invariant everywhere.

Timelike case with axial symmetry
Let us define M ⊂ ⟪M⟫ as the region where f = 0, that is, where V is timelike. Around each point of M one can define a local chart in which the metric is given by where V = ∂ t , ω and h is a 1-form and metric on the four-dimensional riemannian base manifold B := M/R V ⊂ Σ. Note that h is invariantly defined on M, and ω is defined by ι V ω = 0 and dω = −d(f −2 V ) up to a gradient. Since L V X (i) = ι V X (i) = 0, X (i) can be regarded as 2-forms on B. Using properties of the Killing spinor, one can show that (B, h, X (i) ) is hyper-Kähler, i.e. 4 where ∇ (h) is the Levi-Civita connection of h, and ⋆ h denotes the Hodge star on the base with orientation η defined by the spacetime orientation f (dt + ω) ∧ η. On M the Maxwell-field is given by where G + = 1 2 (1 + ⋆ h )f dω. From our assumptions, using Lemma 1 of [19] it follows that the two Killing fields commute, and W defines a triholomorphic U(1) action on B (i.e. L W X (i) = 0). We will use the gauge L W t = 0, thus W also preserves ω. It is well-known that a hyper-Kähler four-manifold with triholomorphic U(1) action can locally be written in Gibbons-Hawking form [38] where x i , i = 1, 2, 3 are cartesian coordinates on R 3 , W = ∂ ψ , H is a harmonic function on R 3 , and the 1-form χ satisfies Here ⋆ 3 denotes the usual Hodge star operator on R 3 with respect to the euclidean metric. In this chart the cartesian coordinates are related to X (i) via Note, that (20) can be regarded as an equation on B ⊂ Σ, as both X (i) and W can be regarded as tensors on Σ, the orbit space of V . An important difference to the asymptotically flat case is that (20) does not define the functions x i globally on the DOC. A useful result of [34] is that if the solution admits a U(1) Killing field commuting with V that is also triholomorphic on B, then the whole solution is locally determined by H and three further harmonic functions, K, L, M on R 3 as follows. Let us define a function and two 1-forms (up to a gradient) on R 3 by Then f and ω can be written as while the Maxwell field takes the form We would like to emphasise again that this is a fully local result on M. Given a local solution, the corresponding set of harmonic functions H, K, L, M is not unique. Indeed, one can check that yield the same solution for any c ∈ R.
Following [29], let us define a key spacetime invariant as where the last three equalities are valid on M when N > 0. Note that N is preserved by both Killing fields, and its zeros in ⟪M⟫ exactly coincide with F = {p ∈ ⟪M⟫|W p = 0} by our assumption that the span of Killing fields is timelike on ⟪M⟫ [19]. It is also worth noting that from (20) and (6) we have irrespective of whether V is timelike. We now use the facts above to prove the following result for a general (not necessarily timelike) solution.
Proof. Let N be the set on which V is null in the DOC, i.e. N := f −1 ({0}) ∩ ⟪M⟫. This is closed in the DOC by the continuity of f , hence N = int N ∪ ∂N . If int N = ∅, then M is dense in ⟪M⟫. If int N = ∅, then assume for contradiction that there exists some p ∈ ∂ int N , and look at a simply connected neighbourhood U of p in ⟪M⟫. Since U is simply connected, we can integrate (20) to obtain x = (x 1 , x 2 , x 3 ) : U → R 3 . Since f (p) = 0, p is not a fixed point of W (otherwise the span of Killing fields would be null contradicting assumption (ix)). It follows that N(p) > 0, and thus by continuity it is also positive on (a possibly smaller) U, and by (29), x is a submersion, therefore open. We extend the definition of H to U by H := f /N, which is smooth on U, and harmonic on the open set x(U) ⊂ R 3 . In particular, it is zero on the open set x(U ∩ int N ). Since H is harmonic and therefore analytic in x i , it follows that H ≡ 0 on U, but then f ≡ 0 on U, which is a contradiction. We conclude that ∂ int N = ∅, and hence it must be that N = ⟪M⟫, so V is null on ⟪M⟫.
In the context of classifying global solutions, we will refer to the first case in Lemma 1 as timelike case, and to the latter one as null case.

Null case with axial symmetry
In this section we consider the case where V is null on the DOC. The metric locally takes the form [34] g where and We now assume the existence of an axial Killing field that preserves V and X (i) , and deduce the following lemma. Proof. We will use properties of W together with the gauge freedoms that preserve the form of the solution [34].
so W u only depends on u.
Then, in the new coordinates (omitting primes) preserves the form of the solution (after redefining Q and a), and the choice dv Finally, the remaining gauge freedom allows us to change the v = const surfaces by x) the axial Killing field becomes W = ∂ u as claimed. In this coordinate system all metric functions are independent of u.
It is useful to note that the remaining gauge freedom that preserves the form of the metric, the two-forms X (i) , and the Killing fields is This, however, changes Q and a as [34] hence these quantities are not gauge-invariant. The lack of dependence on u simplifies the analysis of [34], and the full solution can be obtained as follows. a i is determined up to a gradient term (corresponding to the gauge freedom (34)(35) where ⋆ 3 is the Hodge-star on flat R 3 , and K(x i ) is another harmonic function on R 3 . Choosing positive orientation to be given by dv ∧ du ∧ dx 1 ∧ dx 2 ∧ dx 3 , (9-10) yields 5 Defining D a := a i ∂ i and W ij : Note, that (38) only determines Q up to a harmonic function Q 0 , hence the local solution is determined by three harmonic functions G, K, Q 0 . From Lemma 2 follows that locally dx i = ι W X (i) (as in the timelike case (20)). Our assumption that the span of Killing fields must be timelike excludes any fixed points of W in the DOC for the null case, and therefore N > 0 on the DOC. It follows that G is globally defined by We can flip the sign of G by redefining u → −u and Q → −Q, so without loss of generality we will take on the DOC. On the horizon W must be orthogonal to the generators of the horizon (since every Killing field is tangent to the horizon), which must be proportional to V , thus N = 0 and G must diverge. Note that from (37) and Lemma 2 it follows that The left-hand side is invariantly defined, but since the DOC is not necessarily simply connected, this does not define K globally. Still, in each local (sufficiently small) patch we find that K/G is bounded. In particular, near the horizon K can diverge at most as G does. Later we will show that (41) indeed globally defines K on ⟪M⟫.

Asymptotics
In this section we determine the asymptotic behaviour of the U(1) Killing field W , and the cartesian coordinates x i using Definition 1. As discussed in the previous section, generally we take the stationary Killing field, ∂ 0 in asymptotic coordinates (11), to be a linear combination of the other two, i.e.
with γ, v H constants. For this we first need to look at the asymptotic form of Killing fields in an asymptotically Kaluza-Klein spacetime.

Asymptotic form of Killing fields
We first narrow down the possible form of the axial Killing field near spatial infinity. Proposition 2.1 of [59] for asymptotically flat spacetime states that the Killing fields asymptotically approach those of Minkowski. The statement carries over to asymptotically Kaluza-Klein spacetimes.
Lemma 3. Let (M, g) be stationary and asymptotically Kaluza-Klein as in Definition 1, and K a Killing field that commutes with ∂ 0 . Then there exist constants Λ ij = −Λ ji such that If all Λ vanish then there exist constants A µ such that If all Λ = A = 0 then K = 0.
Proof. The proof is identical to the one in Appendix C of [60], but for completeness we outline it here in a bit more detail. It is well-known that for a Killing field K Using this, we can write This means that for ther-derivatives we obtain where we used ∂r =r −1 u i ∂ i . In Appendix A we derive from Definition 1 that where the explicit form of the matrix C AB (u i ,ψ) can be obtained from (49) and C ′ > 0 is a constant. For a uniform bound on C AB (second relation in (51)) we used that due to (50) C AB (r,θ,φ,ψ) ≤ B(θ,φ,ψ)r 0 ≤ C ′ , whereθ,φ (together denoted by θ A ) are angular coordinates on S 2 and B is some function. Therefore, by integrating (51), there exists a β . Let us assume that β ≥ 1 + τ . Using our estimates for K, ∇K in the right-hand side of (49) and then (48), we obtain an improved estimate with β → β − τ . We iterate this procedure until 1 ≤ β < 1 + τ . Then (47) yields From the Lemma 6 of Appendix A in [61] it follows that ∇ µ K ν − Λ µν = O(r β−τ −1 ) with some Λ µν constants, which we first assume that are not all zero. This substituted back into (47) improves the estimates to β → 1 and K µ = O(r). Finally, we obtain This should hold for ∀k ∈ Z, which implies that Thus, we obtained (43).
. By a similar procedure we can improve this estimate by −τ at each iteration until

. Then again by Lemma of Appendix A in [61] and the arguments above, there is a constant
If Λ ij = 0 and A µ = 0, by the iterative process W µ = O(r −κ ) for any κ > 0. Then integrating (51) as we obtainr 2C ′ 0 X(r 0 ) ≤r 2C ′ X(r) → 0 asr → ∞, thus X(r 0 ) = 0, which means that the Killing field and its first derivative is zero at a point, hence it is zero everywhere.

Chart at spatial infinity from supersymmetry
In this section we construct a chart at spatial infinity defined from Killing spinor bilinears as described in Section 2. We also show that W is tangent to the Kaluza-Klein direction.
The proof is quite different depending on whether v H = 0 or v H = 0 in (42), and for the former, we first need to derive the asymptotic form of the hyper-Kähler structure. In this case the supersymmetric Killing field V coincides with the stationary one, which means V is timelike in the asymptotic region (thus we are in the timelike case since V is not globally null). Furthermore, its limit at infinity is also timelike (i.e. limr →∞ g(V, V ) < 0), hence we are in the asymptotically timelike case, as opposed to the (timelike or null) case when limr →∞ g(V, V ) = 0, which we call asymptotically null.
For the case v H = 0, let us normalise the Killing spinor such that γ = 1, i.e. V = ∂ 0 in the asymptotic coordinates (11). Definition 1 implies that on the asymptotically Kaluza-Klein end we have 7 and the asymptotically Kaluza-Klein end Σ 0 ⊂ B is a circle fibration over R 3 \ B 3 . For the asymptotic form of the hyper-Kähler structure we prove the following lemma.

Lemma 4.
Assuming v H = 0, on the asymptotically Kaluza-Klein end the complex structures of (B, h) can be written as 8 where Ω Proof. The proof is analogous to the asymptotically flat case (Lemma 4 of [19]). From the quaternion algebra X (i) · X (i) = −4 (no sum over i), hence to leading order X which after integration yields X ab are constants. Let us now defineX (i) (58)). Then using we deduce by anti self-duality of X (i) that the constantX (15) implies thatX (i) obeys the quaternion algebra with respect to h 0 . Since (60) forms a basis of ASD 2-forms with respect to h 0 also satisfying the quaternion algebra, we can always perform a global SO(3) rotation of X Next we use Lemma 3 and triholomorphicity to deduce the asymptotic form of W . The following result holds for both timelike and null cases.

Lemma 5.
For any values of v H , on the asymptotically Kaluza-Klein end we can choose coordinates such that the metric is of the form (11) with the stationary Killing field ∂ 0 , and the U(1) Killing field is given by Proof. By Lemma 3 the leading order behaviour of W is determined by constants Λ ij and A µ (using the notation of the Lemma). Since W has closed orbits, A 0 = 0 and W 0 must be subleading inr. First we will consider the general case v H = 0, and assume that Λ ij are not all zero for W . Then W generates a rotation on R 3 (possibly simultaneously with a rotation in the Kaluza-Klein direction), therefore, without loss of generality we can write where (r, θ, φ) are the usual spherical coordinates on R 3 and Λ is a constant. Then the norm It follows that Λ = 0, therefore Λ ij = 0, and by Lemma 3 Since W has closed orbits, where the normalisation has been chosen such that W has 4π-periodic orbits.
In the spacetime the integral curves of W wind around the Kaluza-Klein direction, hence we can adapt coordinates To see this, note that W commutes with the stationary Killing field, thus ∂ 0 λ µ = 0, and the stationary Killing field is unchanged by the coordinate transformation, i.e.
, the metric only receives O(r −τ ) corrections, so g ′ has the same form as (11). Thus, (after omitting primes) we get the claimed result.
For the special case v H = 0 we work on the hyper-Kähler base and use the triholomorphic property of W . L W V = 0 and L W f = 0, thus the projection of π * (W 0 , W a ) = W a defines a Killing vector on the base space. To leading order W is equal to π * W (thus in the following we do not distinguish between the two). Again, assuming that Λ ij are not all zero for W , without loss of generality we can write it as (64). Using the asymptotic form of the complex structures from Lemma 4 one can check that ∂ φ preserves only one of them, which would contradict triholomorphicity. It follows again that Λ = 0, and by the same arguments as for the general case, W must have the claimed form.
Next, we construct the asymptotic charts for the cases when V is asymptotically timelike or null separately. For the asymptotically timelike case we have the following result.
Lemma 6. Assume that V is asymptotically timelike. Then the base of the asymptotically Kaluza-Klein end is covered by a single chart in which (together with the vertical coordinate ψ) the spacetime metric is of Gibbons-Hawking form ((13) with (18)). The Gibbons-Hawking coordinates are related to the asymptotic coordinates by where (with appropriate normalisation of V ) and the Killing fields are given by V = ∂ t and W = ∂ ψ , and the cartesian coordinates x i provide a surjection to R 3 \ B 3 R for some R > 0, Proof. The supersymmetric Killing field is given by V = γ∂ 0 − γv HL −1 ∂ψ. Since it is asymptotically timelike, we can normalise it such that which in terms of the constants γ, v H yields (69) with |v H | ≤ 1. Let us define so that W = ∂ ψ and V = ∂ t . In these coordinates the base metric becomes This has the same form as the base metric in Lemma 4 withL 2 →L 2 γ 2 . One can use similar arguments to deduce that the hyper-Kähler two-forms are Hence, for the cartesian one-forms we get Since ι W (ι W X (i) ) = 0 = L W (ι W X (i) ), and W is tangent to the circle fibres, the closed 1-form ι W X (i) descends to the base of the fibration R 3 \B 3 , which is simply connected. Hence, by (20) the cartesian coordinates x i can be globally integrated on (and uplifted to) the asymptotic end to get (68). Since (u 0 , u i ) form a single chart of the base of the asymptotically Kaluza-Klein end, so do the Gibbons-Hawking coordinates (t, x i ). On this chart N =L 2 γ 2 + O(r −τ ) > 0, hence by (28) the metric is invertible.
Finally, we consider the case when V is asymptotically (or globally) null, i.e. limr →∞ g(V, V ) = 0, and thus without loss of generality we can take We then have the following result.
Lemma 7. Assume V is asymptotically (possibly globally) null as in (75). Then the base of the asymptotically Kaluza-Klein end is covered by a single coordinate chart of (t, x i ) with Let parametrise the fibres. In such a chart V = ∂ t and W = ∂ ψ .
Proof. From (4) it follows that the most general form of X (i) is jk . From the 00 component of (6) using (11) we obtain hence T (11), (75) and (3) follows f = O(r −τ /2 ), and using (5), one can check that R where we used that the relevant Christoffel symbols decay as O(r −τ −1 ) (Appendix A) and that F = O(r −τ /2−1 ) (see Remark after Definition 1). After integration of (80) we get for some constants T (i) j . Furthermore, these constants satisfy (79), hence T ∈ O(3). By a global orthogonal transformation of u i , we can arrange that and thus By the same argument as in Lemma 6 x i can be integrated to get and (t = u 0 , x i ) is a global chart on the base of the asymptotically Kaluza-Klein end.
The notation of the coordinates (t, ψ, x i ) matches that of the timelike (but here asymptotically null) case ( (13) and (18)), but the result is valid for the globally null case as well. Then the coordinates of (30) are given by v = u 0 and u = ψ instead.
We have thus in each case constructed a chart on the asymptotically Kaluza-Klein end that is adapted to the Killing spinor bilinears. In particular, a suitable time coordinate t and the cartesian coordinates x i , related to the bilinears X (i) by (20), provide a single chart of the base of the fibration at infinity.

Classification of null supersymmetric near-horizon geometries
The near-horizon analysis and classification of near-horizon geometries in [16] assumes that the supersymmetric Killing field becomes timelike in some neighbourhood of the horizon, which is no longer true in the globally null case. In this section we show that the same results apply in the null case, and we deduce some technical results for later use.
Following [16] we introduce Gaussian null coordinates near a connected component of the horizon. In these coordinates the horizon is at λ = 0, V = ∂ v , and the metric takes the form where γ and h are a family of metrics and 1-forms on the 3-manifolds H v,λ , which are the v = const , λ = const hypersurfaces. Note that since V = ∂ v is Killing, the metric components do not depend on v. Let H = H v,0 denote a spatial cross-section of the horizon. Equations (4)(5) imply that we can choose a coframe Equation (9) implies that V is hypersurface orthogonal, which is equivalent tô whered denotes the exterior derivative projected onto H v,λ (i.e. it does not include λ derivatives). Closedness of X (i) is equivalent to (87) together witĥ To determine the Maxwell field, following [34], it is convenient to work in the spacetime coframe Equation (7) implies that F = F +i e + ∧ e i + 1 2 F ij e i ∧ e j . F ij is determined by (9), while + + j component of (10) determines F +i . Choosing e − ∧ e + ∧ e 1 ∧ e 2 ∧ e 3 to be positively oriented, after some algebra we get where ⋆ γ is the Hodge star of γ with orientation e 1 ∧ e 2 ∧ e 3 , and we used that X (i) = e + ∧ e i . Equation (90) agrees with (3.35) of [16] in the null limit, but its derivation does not rely on the assumption that V becomes timelike outside the horizon. The analysis of [16] determines the leading order behaviour of the metric quantities h, Z (i) , γ, which applies without modification. For completeness, we sketch it here as well. (87) and Bianchi identity for (90) impliesd respectively, from which follows where∇ and Ric are the Levi-Civita connection of γ and its Ricci tensor, and · is with respect to γ −1 . The j + k component of (10) at λ = 0 yields an expression for∇Z (i) on H, which after taking another derivative and antisymmetrising yields Then considering the integral I = H |∇h| 2 γ dvol, after integration by parts and using (91-93), one can show that∇ h = 0 on H.
Equation (88) implies that Z (i) are hypersurface orthogonal, and without loss of generality there exist coordinates z i and a function K(z) such that Equation (94) imposes a condition on K, which has two solutions, one corresponding to T 3 horizon geometry, which is not allowed [6], the other one is where ψ = 1 2 log(z i z i ) and K 0 is some constant. This corresponds to S 2 ×S 1 horizon geometry with where we introducedx i = z i exp(−ψ) that satisfyx ixi = 1 and dx i dx i is the round metric on S 2 .

Imposing axial symmetry
Building on the above results of [16], it has been shown in [19] that a U(1) Killing field that preserves V and X (i) must be of the form W = W ψ ∂ ψ + W v ∂ v + O(λ) for a constant non-zero W ψ and a function W v on H. In fact, since W is spacelike on the horizon 10 , one can choose v = const surfaces such that W is tangent to H, i.e. W v = 0, in which case one can show that exactly in some neighbourhood of the horizon (see Remark after Lemma 6 of [19]). In the following we will work in such a coordinate system. Another result of [19] is that ι W X (i) ∝ dλ 11 , hence if one can integrate (20) to obtain functions x i , these are constant on each connected component of the horizon (i.e. each connected component of the horizon is mapped to a point in R 3 ). In the null case the horizon topology is S 2 × S 1 and W is tangent to the circle fibres, thus in a neighbourhood of the horizon one can integrate (20) to obtain (up to an additive constant which for simplicity we set to zero) [19] x and the radial distance on R 3 from a horizon component is For later reference, we now look at the next order in λ and prove the following lemma. Proof. From the leading order metric and (98) follows that g(W, W ) = (W ψ K 0 ) 2 + O(r) thus α 0 = (W ψ K 0 ) 2 is trivially constant. In the second-to-leading order we introduce where Z 1 andĥ 1 are one-forms on S 2 . (Note that h and Z (i) cannot depend on ψ, since ∂ ψ is Killing that also preserves X (i) .) (87) implies that and from (88) it follows that for some functions ζ i on S 2 . The norm of W is given by We now show that h 1,ψ and ζ ixi are constants. Bianchi identity for (90) in second-toleading order yields that where ⋆ 2 and , are the Hodge star operator and the inverse metric on the unit S 2 , where the orientation is related to the horizon orientation by ǫ γ = K 3 0 dψ ∧ ǫ S 2 . From (106) and the compactness of S 2 follows that h 1,ψ + ζ ixi is constant. (107) can be written as Acting with ⋆ 2 d⋆ 2 on (108) yields that ζ ixi is harmonic on S 2 , hence it is constant. Thus, both h 1,ψ and ζ ixi are constants, and the norm of W has the claimed form.
The above proof solely relies on the near-horizon geometry, and does not use the form of the solution in the DOC. If one uses (30), (37), and Lemma 2, there is an alternative way of proving that h 1,ψ and ζ ixi are constants. Integrating (20) to second-to-leading order yields Using this, the inner product of Killing fields becomes On the other hand, from (30) g(V, W ) = −G −1 , so G diverges as ∼ 1/r and harmonicity on for some constant C. The two expressions for the Maxwell field (37) and (90) evaluated on the horizon using (31), (86), and (99) yields respectively where G =: where we used the boundedness of K/G and the harmonicity of K. Thus, by comparison one obtains where for (115) we have applied ⋆ 2 on the one-form in the second line of (113). Taking ⋆ 2 d of (115) and substituting into (114) yields and therefore the first term of (112) vanishes, which, to leading order, corresponds to the first term of (90). Thus, the Maxwell field must be of the form where we used (111). The dλ ∧ dψ ∧ dx terms of the Bianchi identity for (117) imply for some function ζ(x), and by (115) ζ = ζ ixi is a constant, and by (111) so is h 1,ψ .

Corollary 1.
There exists a gauge choice of (34-35) such that around a horizon component Proof. The invariant g(W, W ) = −Q/G + G 2 |a| 2 . By (36) and using that K = O(r −1 ) and G = O(r −1 ), one can see that (with an appropriate gauge choice of (34-35)) we have G 2 |a| 2 = O(r 2 ). Thus, from harmonicity of G and Lemma 8 follows the claim.

Orbit space and the general global solution
The next step is to determine the structure and topology of the orbit space. We will show that even if the DOC is not simply connected, it is still possible to define invariants (scalar functions) from certain closed 1-forms that are left invariant by V, W and thus descend to the three-dimensional orbit spaceΣ := ⟪M⟫/[R V × U(1) W ]. Using these functions, one can 'almost globally' define the harmonic functions, and prove that the solution is globally defined by a set of harmonic functions of multi-centred type, i.e. with isolated singularities where they diverge as 1/r on R 3 .

Orbit space analysis and invariants
Let us now look at the structure of the orbit spaceΣ, which has been analysed in detail in [7]. It is a topological manifold [62,63] i are non-isolated fixed points corresponding to 'bolts' [64], andĤ j are quotients of horizon components. Generally, the orbit space can be written asΣ =F ∪Ê ∪L, whereF ,Ê,L denotes fixed points (U(1) isotropy), exceptional orbits (discrete isotropy), and regular orbits (trivial isotropy), respectively.L is open inΣ, and has the structure of a smooth manifold.Ê consists of curves inΣ that are either closed, end on an isolated fixed point, or on a horizon componentĤ i .F consists of isolated fixed points or the previously mentioned 2-spheres S 2 i . A crucial observation is that the orbit space of the asymptotic regionΣ 0 := Σ 0 /U(1) W ≃ R 3 \B 3 is simply connected, so we can use topological censorship for Kaluza-Klein asymptotics (Theorem 5.5 of [12]) to conclude thatΣ = Σ/U(1) W is simply connected. This, however, is not sufficient to define potentials for closed 1-forms, as the orbit space, in general, fails to be a smooth manifold at fixed points and exceptional orbits. To rule the existence of the latter out, we have the following lemma.

Lemma 9.
There are no exceptional orbits in ⟪M⟫.
Proof. Let us assume for contradiction that there exists an exceptional orbit E e ⊂ ⟪M⟫ with isotropy group Z p , and let e ∈ E e be a point on that orbit, i.e. e i2π/p · e = e, where e iα ∈ U(1) with α ∼ α + 2π. First, we will construct a local chart around E e as follows.
By assumption S e := span{W e , V e } ⊂ T e M is timelike and non-degenerate, hence S ⊥ e is spacelike. Since L W V = L W W = 0 and W is Killing, e i2π/p * is an isometry that preserves S e and S ⊥ e , so R := e i2π/p * | S ⊥ e is a three-dimensional orientation-preserving 12 rotation. Furthermore, R p = Id S ⊥ e , thus R must be a rotation by 2πn/p for some n ∈ Z. Let Z e ∈ S ⊥ e be the unit vector preserved by R, and X e , Y e ∈ S ⊥ e such that {Z e , X e , Y e } is an orthonormal frame in S ⊥ e and R X e Y e = cos(2πn/p) sin(2πn/p) − sin(2πn/p) cos(2πn/p) Now let us extend Z e , X e , Y e along E e by (The factor of 2 only appears because W is normalised such that it is 4π-periodic.) Simply Lie-dragging X, Y along W would rotate them according to (119), which is cancelled by the right-hand side of (120) (recall that the exceptional orbit is 2π/p-periodic), so extending X and Y as in (120) is necessary in order for them to be single-valued along E e . Finally, let us extend X, Y, Z along the integral curves of V to some 2-surface E by This is possible since [V, W ] = 0. Thus, we have constructed an orthonormal frame along E in the normal bundle NE. Let λ parametrise the integral curve of 2p −1 W along E, so that λ ∼ λ + 2π on E, and τ be the affine parameter distance from E e along an integral curve of V . By the Tubular Neighbourhood Theorem (see e.g. [65]) we can introduce coordinates {x, y, z} by exponentiating linear combinations of {X, Y, Z} such that E is at x = y = z = 0, which together with λ, τ form a local chart in some neighbourhood U of E e . Since W, V are Killing, they map geodesics to geodesics, hence V = ∂ τ , and components of W can be deduced from its action on the orthonormal frame in NE, that is (cf. (2.20) of [7]) Now we look at X (i) on E. From (4) it follows that X (i) have no τ leg. The two-forms are preserved by W , thus for a = λ, z. Taking a λ derivative of these equations, and substituting the originals back yields that and similarly for X (i) ay . Equation (125) admits solutions which are 2πp/n-periodic in λ. The exceptional orbits, however, are 2π-periodic in λ, which means n/p ∈ Z. In that case the neighbouring orbits (for x, y > 0) would have the same isotropy as the exceptional one, which cannot happen. Hence, the only solution to (125) is X (i) ax = X (i) ay = 0 for a = λ, z. Thus, at each point of the exceptional orbit X (i) ∈ span{dλ ∧ dz, dx ∧ dy}, which is two-dimensional, contradicting linear independence of the three two-forms following from (6).

Remark.
The key assumptions for Lemma 9 is that the axial Killing field preserves three linearly independent 2-forms possibly having legs in four directions (recall that they have no dv legs). This follows from supersymmetry in both timelike and null case, but it is a more general result.
From the assumption that the span of Killing fields is timelike (assumption (ix)), it follows that F ⊂ int B, so F is a set of fixed points of a triholomorphic Killing field, which must be isolated (see [38] and the proof of Lemma 8 in [19]). This rules out any 'bolts' [64] in the orbit space. Thus, we have the following result (cf. Lemma 8 of [19]). Corollary 2. The orbit spaceΣ =L ∪F , whereL corresponds to regular orbits, andF is a finite set of points corresponding to isolated fixed points of the U(1) action. Its boundary is Corollary 3.L is simply connected.
Proof. Corollary 2 is an immediate consequence of Lemma 9. Corollary 3 follows from the fact that a simply connected 3-manifold stays simply connected after the removal of finite many points.
Lemma 10. The equations define smooth functions x i , Ψ globally on ⟪M⟫ ∪ H (up to an additive constant).
Proof. From their definition, we have and similarly for dΨ using (7), so they descend as smooth 1-forms onL. SinceL is a simply connected smooth manifold by Corollary 3, they define smooth functions x i and Ψ onL. In fact, since W | H = 0 and the orbits of connected horizon componentsĤ i = H i /(R × U(1)) ≃ S 2 are simply connected [7], the functions can be defined on the horizon as well. Then we can uplift these functions to (⟪M⟫ ∪ H \ F )/R by L W x i = L W Ψ = 0. We then extend the functions to the isolated fixed points continuously. This is possible, since the functions can be integrated on some small 4-ball around a fixed point with a constant chosen such that the function agrees with the one on (⟪M⟫ ∪ H \ F )/R. Finally, we uplift x i and Ψ to the spacetime by L V x i = L V Ψ = 0.

General global solution
We are now ready to deduce the global form of the general solution that satisfies our assumptions, heavily relying on results of [16,19,29].
where h, h i are constants and a i are the positions of centres corresponding to connected components of the horizon or fixed points of the axial symmetry (the latter only possible when f ≡ 0). In the timelike case (f ≡ 0) the metric is of Gibbons-Hawking form determined by (21)(22)(23)(24), and the Maxwell field is determined by (26). In particular, (128) and (26) smoothly extends in coordinates (t, ψ x i ) to regions where f = 0.
Proof. The functions x i are globally defined on ⟪M⟫ ∪ H by Lemma 10. By Lemma 1 V is either globally null, or V is timelike on some dense subset of ⟪M⟫. The near-horizon analysis of Section 4 for the null case, and that of [16,19] for the timelike case implies that each connected component of the horizon is mapped to a single point of R 3 by x (Lemma 6 of [19]). Furthermore, there are no exceptional orbits of W (Lemma 9). Thus, we can apply Lemma 9 of [19] to deduce that x :L → R 3 \ x(Ĥ ∪F ) is a global diffeomorphism. Hence, on the three-dimensional orbit space with each horizon component added as a single point, Σ ∪ i {H i } is in bijection with R 3 . We next analyse the null and timelike cases in turn.
Null case (f ≡ 0): We have seen that the solution must have the form (30)(31)(32)(33)(34)(35)(36)(37). a i must correspond to horizon components, as fixed points are ruled out (see discussion after Lemma 2). G is globally defined and non-zero on ⟪M⟫ by (39). This means that G is a non-zero harmonic function on R 3 \ ∪ i {a i }, which diverges at each a i as ∼ 1/|x − a i |. This follows from the near-horizon analysis (see proof of Lemma 8), or by Bôcher's Theorem (see e.g. [66]). This determines the singular structure of G, and its regular part is a harmonic function on R 3 which approaches G →L −1 by (11), (39) and (75), therefore it is constant.
From (41) and (126) follows that (up to an additive constant) which defines K on ⟪M⟫. It also follows that K has (at most) simple poles at horizon components. Since dΨ ∼ F ∼ O(|x| −τ /2−1 ) (see Remark after Definition 1), Ψ → const at infinity, which we are free to choose to be zero. (This corresponds to the 'gauge' freedom K → K + cG for some constant c.) Thus, the regular part of K on R 3 is constant (in fact zero in this 'gauge'). Finally, let us look at Q. The norm of the axial Killing field g(W, W ) = − Q G + G 2 |a| 2 . From the asymptotic behaviour of G and K, by (36) we deduce that |a| = O(|x| −1 ), and thus Q = −L + O(|x| −1 ). By Corollary 1, near a horizon component there exists a gauge such that with constants q i −1 and q i 0 . It follows that we can impose the following boundary conditions on the inhomogeneous part Q I : we require that it vanishes at infinity, and it is bounded at each horizon component. This fixes the inhomogeneous part 13 , and thus the homogeneous at the i th centre and approaches −L at infinity. By the same arguments as for G, Q 0 has the claimed form.
Timelike case (f ≡ 0): We have already established that the base has Gibbons-Hawking form on M \ F . Now let us focus on the associated harmonic functions. Since N > 0 on ⟪M⟫ \ F , by Lemma 1 of [29] the metric (128) is smooth, invertible with smooth inverse, the Maxwell field given by (26) is smooth, and the harmonic functions are smooth, and they are invariantly defined by even on the set f = 0, where N must be non-zero by our assumption that the Killing fields span a timelike vector space. Here we used that Ψ is globally defined due to Lemma 10 up to an additive constant, changing of which corresponds to shifting the harmonic functions as in (27). Now that we have established that the harmonic functions are smooth at generic points of the spacetime, let us look at their behaviour at fixed points of W and the horizon. For fixed points, note that the expression for H in (132) implies that the zeros of H and f must coincide on ⟪M⟫. Also, f (p) = 0 at a fixed point p and (by continuity) on some neighbourhood, thus H is non-zero in some neighbourhood of p. As N(p) = 0, by (132) H must diverge at p and by Bôcher's theorem H must have a simple pole at x(p) = a i for some i. From (132) it follows that all harmonic functions have the same type of singularity at a i . For the horizon components, using the near-horizon analysis of [16], triholomorphicity implies that the harmonic functions have (at most) simple poles at the horizon components (Lemma 9 of [19]).
This fixes the form of the harmonic functions up to a globally defined harmonic function on R 3 , which, as in the null case, is determined by the asymptotic conditions. By calculating the invariants f, Ψ, g(V, W ), g(W, W ) at infinity and using (132) we find that the regular part of the harmonic functions must be a constant. Depending on whether f → 1 or f → 0 at infinity (asymptotically timelike and null cases respectively), the values of these constants are different.
In the asymptotically null case we obtain by a similar calculation that where for simplicity we set Ψ 0 = 0.
Remark. Note that for all cases harmonicity on R 3 sets τ = 1 in Definition 1 for the fall-off at infinity.

Regularity and asymptotic conditions
Theorem 1 necessarily includes the most general global solution under the stated assumptions, however it does not guarantee that all such solutions are smooth, asymptotically Kaluza-Klein black hole solutions. In order to establish the sufficient criteria for this, we need to check the asymptotics and regularity of these solutions. This analysis is different in the timelike (f ≡ 0) and globally null (f ≡ 0) case, so we consider them in turn.

Timelike case
By Theorem 1 the solution is globally determined by harmonic functions and the asymptotic values of the harmonic functions depend on whether f → 1 or f → 0 as |x| → ∞, as in (133-134), that is respectively and First we analyse the asymptotic geometry, then we determine the regularity conditions near the horizon and around fixed points.

Geometry at spatial infinity
At spatial infinity (|x| → ∞), we can expand the harmonic functions as and similarly for K, L, M. As we will shortly see, in the timelike case the geometry of the 'sphere' at infinity is governed by the second-to-leading order terms in H. After integration, (19) yields where we definedh 0 := N i=1 h i , on R 3 we use standard spherical coordinates (r, θ, φ), andχ 0 is a constant of integration. We can set the latter to an arbitrary value by a coordinate change (ψ, φ) → (ψ + cφ, φ) which shiftsχ 0 →χ 0 − c. We will shortly see that the bundle structure requires thath 0 is an integer, so for convenience we work in a gauge in whichχ 0 ≡h 0 mod 2.
Since we have setχ 0 ∓h 0 to be even, and ψ is independently 4π-periodic, this implies that One can check that the geometry of the 'sphere' at infinity is L(|h 0 |, 1) or S 3 if |h 0 | = 1.

Regularity
Regularity has to be established at (i) generic points, (ii) the horizon, (iii) at fixed points of W . Recall that assumption (ix) implies that N > 0 on ⟪M⟫ \ F . For generic points, smoothness of the solution is established by Lemma 1 of [29], which states that if and H, K, L, M are smooth, then (g, F ) is smooth, and g is invertible with smooth inverse. Therefore, (144) must be imposed on the harmonic functions everywhere on their domain. Much like in the asymptotically flat case, it is currently an open problem to reformulate (144) as an explicit condition on the parameters h i , k i , l i , a i . The smoothness of a solution determined by harmonic functions of the form (135) has been analysed for the asymptotically flat case, and the sufficient conditions have been determined at the horizon and fixed points in [29] for axisymmetric, and in [19] for general (non-symmetric) harmonic functions. Ref. [19] (a) uses a local expansion of harmonic functions in terms of spherical harmonics, (b) assumes that the periodicity of the angular coordinates are determined by asymptotic conditions, and they are ψ ∼ ψ + 4π and φ ∼ φ + 2π independently, where φ is the azimuthal angle in a spherical coordinate system of the R 3 base of (128). For (a), around each centre (which we take to be the origin) we expand the harmonic functions as and similarly for the other harmonic functions. This has the same form irrespective of asymptotics, thus the same sufficient conditions hold for asymptotically Kaluza-Klein spacetimes. For (b) we have seen in Section 6.1.1, that the angle coordinates ψ, φ admit the same identifications as in the asymptotically flat case. This means that the regularity analysis is identical. We here only present the results, details can be found in [19]. If the centre corresponds to a horizon component, existence of a coordinate transformation to Gaussian null coordinates (and also positivity of the horizon area) requires that Integrating (19) and (22) yields where χ 0 , ω 0 are constants of integration, andχ,ω contain higher order terms in r. Smoothness at the axes θ = 0, π is equivalent to The horizon topology is determined by For h −1 = 0 the topology is S 2 × S 1 , for h −1 = ±1 it is S 3 , otherwise L(|h −1 |, 1). Correct identification of the angles requires that If the centre is a fixed point of W , then there is a curvature singularity at the centre unless The Killing fields having timelike span implies that f = 0 at the centre, which is equivalent to and the spacetime has the correct signature around a fixed point if and only if Smoothness of ω is equivalent to and the correct identification of the angle coordinates, similarly to the horizon, requires that

Null case
By Theorem 1 the solution is determined by the harmonic functions on Our assumption that there exists a timelike linear combination of the Killing fields in the DOC (assumption (ix)) excludes any fixed point of W , hence all centres must correspond to connected components of the horizon. We start our analysis at the horizon.
Let us focus on a single horizon component, around which we will use standard spherical coordinates (r, θ, φ) on R 3 , so the horizon is at the origin. Then we can expand locally the harmonic functions and, using Corollary 1, Q as with some constants γ −1 , γ 0 , q −1 , q 0 , κ −1 , κ 0 . The quantities with tilde are of O(r), and in the case ofG,K harmonic 15 . Recall that G is positive in the DOC and must diverge at the horizon, which means that In Section 4 we have seen that the axial Killing field W is spacelike on the horizon. Recall that 0 < g(W, W ) = −Q/G + G 2 |a| 2 and by (36), G 2 |a| 2 = O(r 2 ) (with a suitable gauge choice of (34-35)), which implies that q −1 < 0.
where c 0 is a constant of integration, Y m l are the spherical harmonics, ⋆ 2 is the Hodge star on the two-sphere, and the coefficients are defined by Note thatc is analytic in r and smooth on the two-sphere. The metric (30) has a coordinate singularity at the horizon. The coordinate change that removes this singularity is of the form 15Q contains the inhomogeneous part of Q, hence not harmonic.
with A 0 , A 1 , B constants. In these coordinates the metric becomes Bcdr Gr where the norm is with respect to the flat metric on R 3 . Setting the 1/r 2 and 1/r terms in g rr and the 1/r terms in g ur to zero is equivalent to The near-horizon geometry ((r, v) → (εr, v/ε), in the limit ε → 0) is given by which corresponds to a black ring, in agreement with the near-horizon analysis in Section 4. Regularity of (166) at the axes θ = 0, π requires that c is a smooth 1-form on S 2 . Setting c 0 = ∓1 makes it smooth at θ = 0, π, respectively, however the required coordinate change between the two charts (covering the northern and southern hemisphere) by (34)(35) is v ′′ = v ′ −2φ, which is well-defined only if v is periodic, but that cannot happen. Hence, smoothness at the horizon is equivalent to Note that the same constraint as (170) for K and G has been derived in the near-horizon analysis in (116). The regularity condition (170) has the following consequence.
Proof. By Theorem 1 G(x) = N i γ i /|x−a i |+L −1 and K(x) = N i κ i /|x−a i |, and all centres correspond to horizon components. In terms of the parameters of the harmonic functions, (170) at each centre (horizon component) is equivalent to where a ij := |a i −a j |+δ ij 16 . We look at this as a system of N linear equations for N unknown κ i of the form A ij κ j = 0, where the matrix A is given by SinceL −1 and γ i are positive, A is a strictly row diagonally dominant matrix, that is As a consequence, A is invertible (see e.g. (5.6.17) of [67]), and thus κ i = 0 for all i, so K ≡ 0.
Proof. As a consequence of Lemma 11 and (169), using (36), a is pure gauge, so it can be set to zero by a transformation of the form (34)(35). Thus, the right-hand side of (38) vanishes, hence the inhomogeneous part Q I is harmonic, and by our boundary conditions in Theorem 1, Q I = 0.
As a consequence, the Maxwell field is simply given by In (160) γ 0 +G is smooth at r = 0, thus so is their contribution in (174). Also, which is smooth on S 2 , so the Maxwell field is smooth at the horizon. We have thus seen that given (161-162), (170), the solution is smooth at the horizon.
In the DOC G > 0 and Q are smooth, hence all metric components are smooth. The metric is invertible, its inverse is given by which is explicitly smooth. The Maxwell field is trivially smooth on the DOC by (174). This concludes the regularity analysis.

Classification theorem for supersymmetric Kaluza-Klein black holes
We are now ready to present our main result, the classification theorem for supersymmetric, axisymmetric, asymptotically Kaluza-Klein black holes (cf. Theorem 3 of [19]). Recall that by Lemma 1 the supersymmetric Killing field is either generically timelike or globally null. We present the two cases in turn.
Theorem 2. An asymptotically Kaluza-Klein (in the sense of Definition 1), supersymmetric black hole or soliton solution (M, g, F ) of D = 5 minimal supergravity with an axial symmetry satisfying assumptions (i)-(x) of Section 2, with supersymmetric Killing field which is not globally null, must have a Gibbons-Hawking base (wherever f = 0), and is globally determined by four associated harmonic functions, which are of 'multi-centred' form, i.e.
where r i := |x − a i |, a i = (x i , y i , z i ) ∈ R 3 , and the parameters are given by in the asymptotically timelike (f → 1) and asymptotically null case (f → 0), respectively, with constantsL > 0, |v H | < 1 and γ = (1 − v 2 H ) −1/2 . The centres either correspond to fixed points of the axial Killing field, or connected components of the horizon. The 1-forms can be written as where χ i 0 are integers such that and The parameters h i , k i , l i , m i must satisfy for each centre i = 1, . . . , N, and for the asymptotically null case Moreover, if a i corresponds to a fixed point of the axial Killing field, satisfying whereas if a i corresponds to a horizon component h i ∈ Z and The horizon topology is S 1 × S 2 if h i = 0, S 3 if h i = ±1 and a lens space L(|h i |, 1) otherwise. Finally, for all x ∈ R 3 \ {a 1 , . . . , a N }, the harmonic functions must satisfy The topology of the 'sphere' at infinity (t = const , |x| → ∞) is The functional form of the harmonic functions (177-179) is required by . By Lemma 1 f = 0 on a dense submanifold, and since zeros of f and H coincide by (28) and assumption (ix), H cannot be identically zero, which implies (186). Smoothness requires that h i = ±1 at a fixed point (153), and h i ∈ Z at a horizon component (151), which also determines the horizon topology. The 1-forms χ, ξ then obtained by simple integration of (19,23), where we introduced constants of integration χ i 0 such that χ i 0 + h i ∈ 2Z for all i. The latter requirement follows from the correct identification of the angles around a horizon component (152) and fixed point (158). For the integration of (22), we follow [27,32,68], and introduce 1-forms β ij as a solution to with r ij = |a i − a j |. One can easily check that β ij as given in (184) is a smooth 1-form away from the centres on R 3 , in particular, it is free of string singularities. (22) is then solved bŷ where Notice that for a local expansion around a centre as in (145) the first two coefficients for the harmonic functions are given by Thus, we see that at each centre (157) and (150) is equivalent to ω i −1 = 0, which yields (185) and (181). With the lack of string singularities inω all conditions of (157) and (150) are satisfied. The remaining smoothness conditions at fixed points (153-156) give (187-188) and at a horizon component (146) yields (189). (190) is the necessary and sufficient condition to smoothness of the solution at generic points (144). Up until this point the proof is in essence identical to that of the asymptotically flat classification in [19].
All that remains to be checked is that the solution is asymptotically Kaluza-Klein. In the asymptotically timelike case, using (21) and (24) we have f = 1 + O(|x| −1 ), ω ψ =Lγv H + O(|x| −1 ). In the asymptotically null case it is simpler to evaluate the metric components directly by using (10) of [29] to obtain In both cases from (191) and (181)ω = O(|x| −2 )dx i , and χ is given by It is straightforward to check that the metric is asymptotically Kaluza-Klein with the coordinates defined in (67)(68) for the asymptotically timelike case, and (76-77) for the asymptotically null case. From (182) it follows thath 0 +χ 0 ∈ 2Z, so the given solution is compatible with our initial assumption made in Section 6.1.1, and the geometry of the 'sphere' at infinity is indeed that of S 3 (forh 0 = ±1), S 2 × S 1 (forh 0 = 0), or L(|h 0 |, 1).
Theorem 3. An asymptotically Kaluza-Klein (in the sense of Definition 1), supersymmetric black hole or soliton solution (M, g, F ) of D = 5 minimal supergravity with an axial symmetry satisfying assumptions (i)-(x) of Section 2, for which the supersymmetric Killing field V is globally null, must be of the form with with constantsL and a i ∈ R 3 correspond to connected components of the horizon with topology S 2 × S 1 . The topology at infinity is S 2 × S 1 .
Proof. Theorem 1 says that the solution is globally determined by three multi-centred harmonic functions through (30) and (36)(37)(38), where each centre corresponds to a connected component of the horizon. The constant terms in (201) are fixed by asymptotic behaviour of the metric. Regularity of the horizon requires that K = 0 by Lemma 11, and thus by Corollary 4 the metric and the Maxwell field are of the claimed form with Q = Q 0 . The constantL is positive by Definition 1, and (161-162) for each centre translate to (202). In Section 6.2 we have seen that this is sufficient for the solution to be smooth on and outside the horizon. One can also easily check that the metric asymptotes to (11) with coordinates defined by One can check that the geometry at infinity (v = const , r → ∞) is S 2 × S 1 .

Remarks.
1. If one removes the condition (186) from Theorem 2, which violates the assumption that the solution is in the timelike class, one exactly obtains the solutions in Theorem 3 with The proof of Theorem 2 heavily relies on the timelike Gibbons-Hawking ansatz, so it is not obvious a priori that one can relax (186). A possible explanation is that these solutions have a common six-dimensional origin, from which one can obtain the timelike and null class by a Kaluza-Klein reduction along different directions [58].
2. The null solutions in Theorem 3 can also be obtained as a limit of certain asymptotically null solutions in Theorem 2, where we definẽ and we take ǫ → 0 while keepingx,ã i , γ i and q i fixed. Thenx becomes the cartesian coordinate of the null solutions andã i are the positions of the centres. One can check that the parameter constraints of Theorem 2 are consistent with those of Theorem 3.
3. The constants have the following physical meaning. v H is velocity of the horizon in the Kaluza-Klein direction with respect to the asymptotic observer [39], and γ is the corresponding relativistic factor. This is apparent from V = γ(∂ 0 − v HW ), whereW is the unit vector in the Kaluza-Klein direction, and V is tangent to the generators of the horizon.L sets the length of the Kaluza-Klein direction at infinity.

4.
As with flat asymptotics, in the timelike case it is not known whether (185-190) guarantees that the DOC is globally hyperbolic. In fact, it is not clear what the sufficient conditions are for it to be stably causal (g tt < 0), which is a consequence of global hyperbolicity. In [69] it has been conjectured that positivity of N −1 = K 2 + HL (which is necessary for smoothness at generic points) implies the lack of closed timelike curves for soliton solutions, which has been supported by numerical evidence. In line with this conjecture, in [19,24] numerical tests found no smooth asymptotically flat black holes with positive ADM mass that violated stable causality. In contrast, in the null case by (176) g vv = GQ < 0, so the spacetime is stably causal automatically with no further constraints on the parameters.
5. All solutions of Theorem 3, and those of Theorem 2 for which K ≡ 0 (which includes all static solutions) have been argued to be exact string backgrounds [58].

A classification of four-dimensional supersymmetric black holes from Kaluza-Klein reduction
In this section we consider the dimensional reduction of the five-dimensional solutions classified in Theorem 2-3, and determine the subclass for which the reduced solutions are smooth on and outside the horizon. We perform the Kaluza-Klein reduction along the direction of W in coordinates adapted to it 17 , so that W = ∂ ψ . For the dimensionally reduced theory, we will follow the field definitions of [39], which are given by Here we used that since L W F = 0, we are free to work in a gauge in which L W A = 0. g (4) is the four-dimensional metric, A (4) and A are one-form potentials, and Φ, ρ are scalar fields. The action (2) then reduces to where G 4 = G/4π, and R (4) denotes the Ricci scalar of g (4) . It is important to establish which fields are physical, as we will require the smoothness of those only. Physical fields must be invariant under five-dimensional coordinate changes of the form ψ ′ = ψ + µ(t, x i ) and gauge transformations which preserve the gauge condition L W A = 0. Since ι W dA is invariant, and 0 = L W (A ′ − A) = d(ι W A ′ − ι W A), the allowed gauge transformations must be of the form A ′ = A + dλ(t, x i ) + cdψ with some constant c. Under such transformations the fields transform as g (4)′ = g (4) , hence the physical fields are those in (210).

Remarks.
1. Even though the last term of (208) contains gauge-dependent fields, the theory is gauge invariant since the five-dimensional theory is. Indeed, one can check that the equations of motion of (208) contain only the physical fields (210).
2. Ref. [40] uses an alternative definition for the 2-form field F (4) , that is thus it is closed (as opposed to F (4) ). It is evident from its definition that smoothness ofF (4) is equivalent to the smoothness of F (4) (assuming the other fields in (210) are smooth), since for solutions with a simply connected DOC (which are the relevant ones here due to topological censorship), ρ is globally defined by dρ up to an additive constant, so dρ is smooth if and only if ρ is smooth. Now we establish the subclass of solutions classified in Theorem 2 and 3 that reduce to a four-dimensional solution that is smooth on and outside the horizon.
Theorem 4. A solution to five-dimensional minimal supergravity as in Theorem 2 defines a four-dimensional, asymptotically flat black hole solution of (208) if and only if all centres correspond to horizon components (i.e. there are no fixed points of the axial Killing field) and for all x ∈ R 3 \ {a 1 , . . . , a N }. Then the four-dimensional solution is given by where N −1 = K 2 + HL, c is an arbitrary constant, H, K, L, M are given by (177-179), and 1-forms χ,ω, ξ are given by (180-184).
Theorem 5. A solution to five-dimensional minimal supergravity as in Theorem 3 defines a four-dimensional, asymptotically flat black hole solution of (208), given by where G and Q are given by (201) with (202).
Proof of Theorem 4 and 5. It is easy to see that the four-dimensional fields are related to five-dimensional smooth invariants by  (10) of [29] yields the right-hand sides of (213-216).
Finally, using the asymptotic explicit form of the harmonic functions, one can easily check that the four-metric approaches the Minkowski metric in coordinates 19 and in the asymptotically timelike and null case 20 , respectively. Remarks.
Definition 2. A four-dimensional spacetime is asymptotically flat if it has an end diffeomorphic to R×(R 3 \B 3 ), and on this end the metric g (4) = −du 0 du 0 +δ ij du i du j +O(R −α )du a du b for some α > 0, where (u 0 , u i ), i = 1, 2, 3 are the pull-back of the cartesian coordinates on R × R 3 and R := δ ij u i u j , and the k th derivatives of the metric fall off as O(R −α−k ) for k = 1, 2 in these coordinates.
We assume that (M 4 , g (4) , Φ, ρ, F (4) , G (4) ) is a solution of (208) such that  Proof. By assumption (ix) we can uplift the solution to five dimensions, identifying M with the total space of the U(1)-bundle, on which we define the five-dimensional metric g and Maxwell field F as where π : M → M 4 is the bundle projection. The Killing spinor ǫ (4) lifts to a fivedimensional Killing spinor ǫ, invariant under the U(1)-symmetry [72] (also see Appendix B). Thus, (M, g, F, ǫ) is a supersymmetric solution of (2) (as we have just undone the dimensional reduction). We will now show that it satisfies the assumptions of Theorem 2 or 3, hence (M 4 , g (4) , Φ, ρ, F (4) , G (4) ) must belong to the class of Theorem 4 or 5. Note that we include the latter solutions by not assuming (186) for the harmonic function H (see also Remark 1 after Theorem 5). We first prove that the five-dimensional DOC, ⟪M⟫ is globally hyperbolic by showing that Σ = π −1 (Σ 4 ) is a Cauchy-surface. Let p ∈ ⟪M⟫, γ an inextendible causal curve through p in (M, g), and U its tangent vector. From causality of γ thus π * U defines a causal curve γ 4 in M 4 . By assumption (ii) γ 4 goes through Σ 4 , hence γ goes through Σ. Acausality of Σ follows from a similar argument, therefore it is a Cauchy surface, and ⟪M⟫ is globally hyperbolic. Next we show that (M, g) is asymptotically Kaluza-Klein according to Definition 1. By compactness of the fibres and assumption (vi), Definition 1(i) is satisfied. The Dirac-currents of ǫ and ǫ (4) define the supersymmetric Killing fields V and V (4) on M and M 4 , respectively. Let W be the generator of the U(1) action normalised such that its integral curves are 4πperiodic 22 . [W, V ] = 0 since ǫ is U(1)-invariant, and π * V = V (4) (for details see Appendix B). Let us adapt local coordinates to the vertical vector field so that W = ∂ ψ and ψ ∼ ψ + 4π. It is obvious that W is a Killing field of g (that also preserves F ). In such a chart the connection is given by η = dψ + A. V preserves g and W , so it must also preserve η = W ♭ /g(W, W ). We may partially fix the gauge by requiring L V A = 0, then we have thus V ψ = c for some constant c. By assumption (iv) on the asymptotic end π * V = V (4) = ∂ 0 , hence V = ∂ 0 + c∂ ψ .
It follows from assumption (viii) that Our gauge condition (224) and (225) imply with some constantc. Equation (226) implies that the rest of A can be written in a gauge without changing the form of the Killing fields such that The leading order behaviour of the five-dimensional metric then becomes for some τ = min{α, β}, and η ab denoting the 4D Minkowski metric. By defining ψ ′ = ψ +cu 0 and rescaling the coordinates by constants, we get a metric of the form (11), with Φ 0 determining the asymptotic length scaleL of the Kaluza-Klein direction. By assumption (viii) and because there is no dependence on ψ, the first two derivatives of the metric have the fall-off as in Definition 1(ii), and the components of the Riemann tensor fall off as O(R −τ −2 ). We also see from (225) and from the final coordinate change that ∂ 0 is a constant linear combination of V and W . Finally, we need to check if the span of the supersymmetric and U(1) Killing field is timelike. For this, by (28) and (222) which is positive on ⟪M⟫ by assumption (iv). That is, at each point of ⟪M⟫ the determinant of the inner product matrix of Killing fields is negative, hence there exists a timelike linear combination of V and W .

Remarks.
1. Let us emphasise that we have not assumed any isometry apart from stationarity, which is guaranteed by supersymmetry (for a class of D = 4 supergravities see [73]). Indeed, generically, solutions of Theorem 6 only have a single Killing field.
2. Assumption (ix) quantises the magnetic charges of the black holes associated with G (4) , which in terms of the harmonic functions means that h i ∈ Z for H in (177). Omitting this requirement leads to a more general class of black holes in four dimensions, however those cannot be uplifted to get a smooth black hole solution in five dimensions.
3. The requirement that the supersymmetric Killing field is timelike in assumption (iv) was also required for uniqueness of four-dimensional supersymmetric, asymptotically flat black holes in minimal supergravity [74], which shows that the general solution belongs to the Majumdar-Papapetrou (MP) class. Alternatively, in minimal supergravity, one can assume the existence of a maximal hypersurface with a finite number of asymptotically flat or weakly cylindrical ends, and prove that the Killing field is strictly timelike and static (see Theorem 1.2 of [74]). This reasoning would not work for the black holes considered in the current work, as these solutions are not static in general.
4. The black holes of Theorem 5 are a generalisation of the aforementioned MP black holes. They are static solutions that depend on two harmonic function on R 3 , and they are, in general, solutions to the full theory (208). If Q = −G, we obtain solutions to fourdimensional minimal supergravity 23 (Einstein-Maxwell), which describe magnetically charged MP black holes. In detail, In contrast, the electrically charged and dyonic MP black holes (with quantised charges) belong to the asymptotically timelike class of Theorem 4 with harmonic functions K = 0, L = γ 2 H, M = γ 3 v H H, and parameterL = 1. For the detailed derivation see the end of Section 3.7 in [34], with the only difference that we used 'gauge' freedom (27) to set K ≡ 0 in agreement with our previous choice in Theorem 1 for the asymptotic values of the harmonic functions. The solution is given by Thus, for 0 < |v H | < 1 we obtain dyonic MP black holes, while v H = 0 yields electrostatic MP black holes. 5. In Remark 4 all static solutions of minimal four-dimensional supergravity (with quantised charges) have been obtained. It is an interesting task to determine all static solutions of the full theory (208) (again, with quantised charges). This includes all solutions in Theorem 5, but also requires deriving the set of harmonic functions that yield dω = 0 and satisfy the constraints of Theorem 2 and 4. Staticity in the timelike case implies 23 Consistent truncation of (208) to 4D Einstein-Maxwell is achieved by Φ = 0, dρ = 0, F (4) = √ 3 that f −1 is harmonic [34], thus by (24) K = cH for some constant c (so we are in the asymptotically timelike case by (178-179)). By (22) dω = 0 requires M + 3cL/3 = kH for some other constant k. Again, changing the 'gauge' of (27) such that K ≡ 0 yields that M = c ′ H for a constant c ′ , and L is unconstrained. H, L, c ′ must be such that (189)(190) and (212) are satisfied. Asymptotic constants (178) fix c ′ =L 3 γ 3 v H . Inequality (190) yields (together with (178)) that H > 0 and L > 0 on their domain, i.e. h i > 0 and l i > 0. Then (212) (which implies (189)) is satisfied if and only if l i > |c ′ | 2/3 h i at each centre, and does not impose any further constraint on the parameters v H andL.
In summary, as in the null case (Theorem 5), the static solutions of the timelike case are determined by two harmonic functions, with parameters satisfying l i >Lγ|v H | 2/3 h i > 0 and h i ∈ Z.
Using the algebraic Bianchi identity, we obtain R cψ ab = −R c bψa + R c aψb = O(r −τ −2 ), and by the fall-off of the Ricci tensor Rψψψ a = −R dψ da + O(r −τ −2 ) = O(r −τ −2 ). Finally, we have When f = 0 and the base is well-defined as a Riemannian manifold, it is Ricci flat. The above argument works with just replacing abc... indices with ijk....

B Killing spinor equation in four and five dimensions
In this section we reduce the five-dimensional Killing spinor equation to four dimensions, and derive the four-dimensional Killing spinor. We also prove that its Dirac-current is the four-dimensional projection of the five-dimensional supersymmetric Killing field.
The Killing spinor equation (KSE) of five-dimensional minimal supergravity is given by [34] ∇ a + 1 4 √ 3 F bc γ a bc + 4δ b a γ c ǫ = 0 , where η is the U(1)-connection, and {E A } 3 A=0 is a co-frame of the four dimensional metric g (4) , and the U(1) Killing field in this frame is W = e −Φ/ √ 3 e 5 . The lie derivative of the Killing spinor with respect to W (which we assume to vanish) is defined as [75] where ω ab is the spin-connection. One can check that in this frame (dW ♭ ) ab = −2ι W ω ab , hence (240) simplifies to L W ǫ = e −Φ/ √ 3 e 5 (ǫ) = 0 .
We now derive the four-dimensional Killing spinor equation from the five-dimensional one. Let us write the following field components in the four-dimensional co-frame: In this frame (238) is given by where in the penultimate step we used that π * e 5 = 0 since e 5 is vertical. The supercovariant connection of (249) satisfies (247), so the four-dimensional Dirac current V (4) is Killing, and coincides with the projection of the five-dimensional Killing vector field, as claimed.