Warped AdS6 × S2 in Type IIB supergravity II: global solutions and five-brane webs

Motivated by the construction of holographic duals to five-dimensional superconformal quantum field theories, we obtain global solutions to Type IIB supergravity invariant under the superalgebra F (4) on a space-time of the form AdS6 × S2 warped over a two-dimensional Riemann surface Σ. In earlier work, the general local solutions were expressed in terms of two locally holomorphic functions A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{A} $$\end{document}± on Σ and global solutions were sketched when Σ is a disk. In the present paper, the physical regularity conditions on the supergravity fields required for global solutions are implemented on A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{A} $$\end{document}± for arbitrary Σ. Global solutions exist only when Σ has a non-empty boundary ∂Σ. The differentials ∂A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{A} $$\end{document}± are allowed to have poles only on ∂Σ and each pole corresponds to a semi-infinite (p, q) five-brane. The construction for the disk is carried out in detail and the conditions for the existence of global solutions are articulated for surfaces with more than one boundary and higher genus.


JHEP05(2017)131
1 Introduction In a recent paper [1], we solved the BPS equations of Type IIB supergravity for configurations with 16 residual supersymmetries on a space-time of the form AdS 6 × S 2 warped over a Riemann surface Σ. The isometry algebra SO(2, 5) × SO(3) of this space-time is uniquely completed into the exceptional Lie superalgebra F (4) with 16 fermionic generators [2,3]. In the present paper we shall impose on the local solutions of [1] the positivity and regularity conditions needed to obtain physically sensible supergravity fields and we shall construct the corresponding general global solutions. The first part of the present paper provides detailed derivations and discussions of results which were described in a short note [4]. The motivation for this work stems from the study of five-dimensional supersymmetric quantum field theories with the help of holographic methods. Minimal Poincaré supersymmetry in five dimensions has eight supersymmetry generators and SU(2) R-symmetry. In their conformal phase, five-dimensional supersymmetric field theories will have SO (2,5) conformal symmetry, SU(2) R-symmetry and sixteen supersymmetries and therefore be invariant under the F (4) superconformal algebra. Unlike in space-time dimensions three, four, and six, where the maximal superconformal algebra has thirty two supersymmetries, in five dimensions the superconformal algebra F (4) is unique and has sixteen supersymmetries.
The existence of physically consistent interacting five-dimensional gauge theories is far from obvious since standard perturbative methods show them to be non-renormalizable. Arguments for the existence of non-Abelian gauge theories were advanced in [5] and in [6] and drew heavily on the structure imposed by five-dimensional Poincaré supersymmetry and the dynamics on the Coulomb branch. It was further argued that the theory in the nonconformal Coulomb phase can flow to a strongly coupled fixed point with superconformal symmetry. The five-dimensional gauge theories of [5,6] can have gauge groups SU(N ) or USp(N ) for arbitrary N , as well as exceptional gauge groups. Consequently, one should expect that, in the large N limit, the theories with gauge groups SU(N ) and USp(N ) will JHEP05(2017)131 admit weakly coupled AdS 6 gravity duals. The explicit construction of such gravity duals is the goal of [1] as well as of the present paper.
Several brane realizations of five-dimensional gauge theories with eight Poincaré supersymmetries have appeared in the string theory literature. In [5,[7][8][9] Type IIA supergravity solutions were constructed using N D4-branes near a stack of D8/O8 branes. These solutions suffer from singularities associated with the presence of the D8/O8 branes. The Type IIA supergravity solutions can be T-dualized to solutions in Type IIB supergravity, which, however, introduces additional singularities [10,11].
In [12][13][14] a large class of five-dimensional theories was constructed in Type IIB string theory using brane webs. The five-dimensional field theories live on D5 branes which are of finite extent in one dimension and are embedded in a web of semi-infinite (p, q) five-branes as well as optional seven-branes. The brane webs realize the moduli space and the matter content of such theories geometrically, generalizing the construction of three-dimensional theories [15] to five dimensions.
The solutions constructed in [1] are local in the sense that all BPS conditions, Bianchi identities and field equations are satisfied point-wise by the choice of two locally holomorphic functions A ± on a two-dimensional Riemann surface Σ whose global structure is as yet undetermined. 1 In [1], we followed the strategy developed some time ago in [19][20][21], and utilized the integrable structure of the reduced BPS equations [22] along with the holomorphic structure of the Riemann surface Σ to its fullest. These techniques have generated the complete local solution to the BPS equations.
The local solutions of [1] need to be supplemented by positivity and regularity conditions in order to produce genuine global supergravity solutions with physically sensible behavior of all supergravity fields. These positivity and regularity conditions are as follows. First, all supergravity fields must satisfy the reality conditions appropriate for Type IIB, such as for example the space-time metric and the dilaton field must be real. Second, the space-time metric must have the correct ten-dimensional Minkowski signature. Third, the space-time metric must be geodesically complete. Fourth, any singularities that might occur in the solutions, such as curvature singularities or a divergent string coupling, must be associated with physically sensible sources, such as BPS branes already present in Type IIB.
The main new result of the present paper is the complete construction of global solutions which satisfy all the positivity and physical regularity conditions enumerated above. Our solutions evade a recently obtained theorem [23] stating that no supersymmetric warped AdS 6 solutions exist in Type IIB supergravity, since our solutions are warped over a non-compact Riemann surface which violates one of the assumptions of the no-go theorem of [23]. We shall show that the resulting global solutions provide viable candidates for holographic duals to five-dimensional superconformal field theories in the large-N limit.
The organization of this paper is as follows. In section 2, we present a brief review of the local supergravity solutions in terms of holomorphic functions obtained in [1]. We articulate the positivity and regularity conditions in quantitative terms, formulate our JHEP05(2017)131 strategy for the construction of global solutions, and show that no physically sensible solutions exist unless Σ has a non-empty boundary. In section 3, we present a new class of global solutions for which the Riemann surface Σ is the upper half plane and which are parameterized by holomorphic functions with simple poles on the boundary of Σ. For these general solutions we provide a detailed counting of the number of free moduli of the space of solutions. In section 4 we present explicit examples with three and four poles and discuss natural connections to (p, q) 5-brane webs. In section 5 we extend the general construction to the case where Σ is an annulus, and in section 6 to the case where Σ has arbitrary genus and an arbitrary non-zero number of boundary components. In the case of more than one boundary component some global regularity conditions remain to be imposed on the solutions. Finally, in section 7, we conclude with a discussion of our results and present some directions for future research.

Local solutions and strategy for global solutions
In this section, we begin by summarizing from [1] the structure of the local solutions to the BPS equations in Type IIB supergravity with 16 supersymmetries and SO(2, 5) × SO(3) isometry on a space-time of the form AdS 6 ×S 2 warped over a Riemann surface Σ. We then review the reality and positivity conditions required on the physical supergravity fields, and impose mild regularity conditions. Finally, we discuss our strategy for the construction of global solutions starting from an analogy with two-dimensional electrostatics.

Supergravity fields in terms of holomorphic functions
The SO(2, 5)×SO(3) isometry dictates the Ansatz for the bosonic supergravity fields, while the fermionic fields vanish. The five-form field strength vanishes by symmetry, while the metric and the complex three-form flux fields are as follows, Here, f 6 , f 2 , ρ, C, and the dilaton-axion scalar field B are functions of Σ only. The metrics ds 2 AdS 6 and ds 2 S 2 are the SO(2, 5) and SO(3) invariant metrics with unit radius respectively on AdS 6 and S 2 , and vol S 2 is the corresponding volume form on S 2 . Finally, we have introduced local complex coordinates w,w in terms of which the metric on Σ is conformally flat.
The supergravity fields for the local solutions were obtained in [1], and expressed in terms of locally holomorphic functions A ± on the Riemann surface Σ. Of fundamental importance to the structure of the local solutions and, as we shall see shortly also to the structure of the global solutions, are the following real functions, (more properly, one should think of κ 2 as a (1, 1) differential on Σ), where the locally holomorphic function B is defined up to an additive constant by the relation, In Einstein frame, the metric functions f 2 6 , f 2 2 and ρ 2 are given by, where the real positive function R is defined by the following quadratic equation, The axion-dilaton field B is given by, while the flux potential function C for the three-form field strength F (3) takes the form, Residual gauge transformations, which act on C by a constant shift, leave F (3) invariant.

SU(1, 1) symmetry and monodromy
Under the SU(1, 1) symmetry of Type IIB supergravity, the metric in the Einstein frame and the five-form field strength are invariant while the axion-dilaton field B and the threeform field strength transform non-trivially. Upon reduction to the AdS 6 × S 2 × Σ Ansatz, the metric factors f 2 6 , f 2 2 , ρ 2 are invariant under SU(1, 1), the five-form field strength is zero which is SU(1, 1)-invariant, while B and C transform as follows, where the additive constant in C allows for a gauge transformation. Correspondingly, the building blocks κ 2 and G of the metric factors must be invariant under SU(1, 1). The invariance of κ 2 combined with the transformation rule of B dictate the transformation rule for the differentials ∂ w A ± , Integration of the differentials ∂ w A ± and ∂ w B gives the transformation rule of A ± and B, where a ± and b 0 are complex constants. The function G is invariant provided we impose, The combined transformations on A ± form a group G = SU(1, 1) × C, whose action on the differentials ∂ w A ± reduces to the action of the subgroup SU(1, 1). When Σ has a non-trivial fundamental group π 1 (Σ) the functions A ± are allowed to have G-valued monodromy around non-trivial closed loops on Σ. The monodromy on the differentials ∂ w A ± is simplest and given by a representation of π 1 (Σ) → SU(1, 1). Mathematically, this representation could be arbitrary. Physically, for a supergravity embedded into Type IIB string theory, however, the monodromies generated by parabolic elements acting on the axion field by a shift correspond to the presence of D7 branes, and have a well-known interpretation. It is also possible for the differentials ∂ w A ± to have a branch point in the interior of Σ around which axion monodromy is allowed.
In the present paper, we shall consider solutions only for which the dilaton/axion field is single-valued and no D7 branes are present. We plan to investigate solutions with nontrivial monodromies in the future.

Positivity, regularity and boundary conditions
Reality and positivity of the function R in (2.5) require κ 2 and G to have the same sign. Reality of ρ 2 and positivity of the metric factors require κ 2 to have the same sign as 1 − R. Without loss of generality, we may choose the branch 0 ≤ R ≤ 1 in (2.5) resulting in the following positivity conditions in the interior of the Riemann surface Σ, These necessary conditions are sufficient to guarantee that f 2 6 , f 2 2 , ρ 2 are real and positive. The topology of Σ, namely its genus and number of its boundary components, remains to be determined by global considerations. The boundary ∂Σ of the surface Σ does not correspond to a boundary of the physical space-time manifold of the supergravity solution. Instead, ∂Σ corresponds to the shrinking of S 2 as part of the slicing of a regular S 3 cycle, thereby ensuring the geodesic completeness of the full space-time geometry. As such, we may characterize ∂Σ by f 2 2 = 0 and f 2 6 = 0 or, equivalently, by the following conditions, We shall now prove the following Proposition,

1.
There are no regular solutions satisfying (2.12) when Σ is compact without boundary.
A practical implication of this proposition is that it will suffice to establish the property G = 0 on ∂Σ in order to guarantee the condition G > 0 of (2.12) in the interior of Σ.

JHEP05(2017)131
To prove the Proposition we begin by evaluating the mixed derivative of G in the second line of (2.2) in terms of κ 2 with the help of the first line of (2.2), (2.14) To prove item 1 of the Proposition, we assume that Σ is compact without boundary and that κ 2 and G are sufficiently regular. Integrating both sides of (2.14) over Σ, the left side vanishes since Σ has no boundary, yet the right side must be strictly negative since κ 2 > 0 in the interior of Σ. Therefore, (2.14) has no solutions consistent with the requirements of (2.12). This statement is in agreement with results found in a recent paper [23] stating that there are no supersymmetric AdS 6 solutions which are warped over a compact space.
To prove item 2 of the Proposition, we solve equation (2.14) along with the boundary condition on G in (2.13) to obtain the following integral equation, 2 Here, G(w, z) is the scalar Green function on Σ, which is symmetric G(z, w) = G(w, z) and vanishes on the boundary ∂Σ, For any two points w, z in the interior of Σ, the function G(w, z) is strictly positive, as we shall review below, and therefore G > 0 in view of κ 2 > 0 in the interior of Σ.
To prove positivity of the Green function G(w, z) for arbitrary points w, z in the interior of Σ, we fix z and define Σ ε to be the Riemann surface Σ minus the disk of radius ε > 0 centered at the point z. In the interior of Σ ε the function G(w, z) is real harmonic in the variable w. By the minimum-maximum principle for harmonic functions, G(w, z) takes its minimum and maximum values on the boundary of Σ ε . Since G vanishes on ∂Σ by (2.13) and, for sufficiently small ε is positive on the circle of radius ε centered at z, it follows that G must be strictly positive in the interior of Σ ε and thus in the interior of Σ as we take the limit ε → 0.

Strategy for global solutions without monodromy
The construction of global supergravity solutions has now been reduced to solving for two locally holomorphic functions A ± on a Riemann surface Σ with non-empty boundary ∂Σ, subject to the positivity condition κ 2 > 0 in the interior of Σ of (2.12) and the boundary conditions κ 2 = G = 0 on ∂Σ. The positivity condition G > 0 in the interior of Σ in (2.12) is then automatic by the Proposition of the preceding subsection. In the present subsection, we shall present our general strategy for solving these conditions, with the assumption that the differentials ∂ w A ± are meromorphic on Σ, and in the absence of SU(1, 1) monodromy.

JHEP05(2017)131
We fix a Riemann surface Σ with non-empty boundary ∂Σ. We shall use an analogy with two-dimensional electrostatics to solve for the meromorphic differentials ∂ w A ± . We begin by introducing the meromorphic function λ on Σ, Our strategy will consist of the following steps, 1. solving for the function λ using the properties of κ 2 and an electrostatics analogy; 2. constructing the meromorphic differentials ∂ w A ± and their integrals A ± ; 3. enforcing the condition G = 0 on the boundary of Σ.

Constructing the function λ via an electrostatics analogy
We begin by constructing the meromorphic function λ. The further properties of λ follow from its definition (2.17) and the properties of κ 2 in (2.12) and (2.13), and we have, • |λ| < 1 in the interior of Σ in view of κ 2 > 0; • |λ| = 1 on the boundary of Σ in view of κ 2 = 0 on ∂Σ; • λ is holomorphic in Σ, since it is meromorphic and bounded in Σ.
To solve for λ subject to these properties, we use an electrostatics analogy. We view the function − ln |λ| 2 as an electrostatic potential for an array of N point-like charges. This potential is real, and the properties of λ imply that − ln |λ| 2 is locally harmonic, strictly positive in the interior of Σ, and zero on the boundary ∂Σ. The potential − ln |λ| 2 must have singularities, since otherwise it would vanish throughout. Its singularities are located at the zeros of λ in the interior of Σ and are logarithmic. Combining these results, we find, where G is the scalar Green function defined in (2.16). The charges q n are real and strictly positive for all n = 1, · · · , N . Since we have |λ| = 1 on ∂Σ, the points s n must lie in the interior of Σ, away from the boundary. Conversely, when the charges q n are all positive, the potential − ln |λ| 2 constructed in (2.18) is positive throughout the interior of Σ. This result may be established by using the positivity of the scalar Green function G(w, z), already proven at the end of subsection 2.3. To construct a single-valued λ(w) from (2.18) further requires the charges q n to be integers. The holomorphic function λ(w) may be obtained by holomorphically splitting (2.18), and the result is unique up to an arbitrary phase factor. The precise expressions will be given in the sequel when Σ is the upper half plane, an annulus, or a Riemann surface of higher genus.

Constructing the differentials ∂ w A ±
To construct the differentials ∂ w A ± we decompose λ into a ratio of holomorphic functions (or holomorphic differentials of equal weight) λ ± , The functions λ ± are clearly not unique, but we may impose the irreducibility condition that λ − has no zeros in the interior of Σ or on its boundary. In that case, the zeros of λ + coincide with the zeros of λ and are given by the points s n of (2.18). Choosing local complex coordinates w,w near a given boundary component in which the boundary is represented by a segment of the real line, we require λ ± to obey the following conjugation property, which guarantees the relation |λ(w)| = 1 when w is real and on the boundary of Σ. The meromorphic differentials ∂ w A ± are then given as follows, where ϕ(w) is a meromorphic function or differential on Σ. The properties of ∂ w A ± , κ 2 and G narrow down the properties of ϕ. The combinations λ ± ϕ must be single-valued meromorphic 1-forms on Σ. The form ϕ can have neither zeros nor poles in the interior of Σ, since a zero would make κ 2 vanish in contradiction with (2.12), and a pole would produce local monodromy in the supergravity fields. Thus, the zeros and poles of ϕ must all lie on ∂Σ, so that ϕ is an imaginary function satisfying ϕ(w) = −ϕ(w) in local complex coordinates where the local boundary corresponds to real w. Therefore we have the following conjugation property for the meromorphic differentials ∂ w A ± , On an arbitrary Riemann surface Σ, with non-empty boundary ∂Σ, the conjugation condition may be formulated globally in terms of the double surfaceΣ, equipped with an anti-conformal involution I, such that Σ =Σ/I, and each point on the boundary ∂Σ is mapped to itself under I. This formalism will be developed and used to construct ∂ w A ± in section 6.

Enforcing the condition G = 0 on ∂Σ
To analyze the vanishing of G on ∂Σ, we locally integrate the conditions (2.22). By choosing the integration constants, we may impose the following conjugation property, While the differentials ∂ w A ± are meromorphic and have poles on ∂Σ, their integrals A ± will have logarithmic branch points on ∂Σ whose branch cuts must be chosen consistently with the conjugation condition (2.23).

JHEP05(2017)131
The condition (2.22) suffices to guarantee that G is constant along any line segment of an arbitrary boundary component which is free of poles in ϕ. To prove this statement, we consider the derivative tangent to the boundary component of the real function G. Since here we have parametrized the boundary component locally by the real part of w we find, using (2.2) and (2.22), Along a boundary component away from poles in ∂ w A ± , the variable w is real and the differences of the derivative terms vanish identically. Therefore ∂ w G + ∂wG vanishes and G is constant on any such segment of an arbitrary boundary component. To promote the piecewise constancy of G to the vanishing of G throughout the boundary will require restrictions on the free parameters of the functions A ± . We shall establish later that there will be one real condition for each pole. When Σ has several disconnected boundary components, there will be further conditions to ensure that G = 0 along each boundary component separately.

Global solutions for the upper half plane
In this section, we shall construct explicitly the global solutions for which Σ has genus zero and a single boundary component. In this case Σ may be mapped to the unit disk whose boundary is the unit circle, or to the upper half complex plane whose boundary is the real axis. Our construction will follow closely the strategy outlined in section 2.4, and it will be convenient to identify Σ with the upper half plane Σ = H and ∂Σ = R.

The function λ
The starting point is the general formula for − ln |λ| 2 given in (2.18), for an array of N positive integer charges q n located at points s n in the upper half plane H with n = 1, · · · , N . Without loss of generality, we set q n = 1 for all n and obtain higher integer charges by coalescing charges at different points s n if need be. The set-up is represented schematically in figure 1. The scalar Green function G for H is given by, The resulting formula for λ is, Here, λ 0 is an arbitrary phase factor, which may be viewed as effecting an SU(1, 1) transformation A + → λ 0 A + and A − →λ 0 A − on the solution with λ 0 = 1. Following the process leading to (2.19), we express the function λ, which is holomorphic on H and meromorphic JHEP05(2017)131 Figure 1. The upper-half plane Σ = H; its boundary ∂Σ = R; an array of zeros of the function λ at the points s n and poles ats n , with n = 1, · · · , N and Im (s n ) > 0; and an array of poles of ∂ w A ± at the points p ∈ R with = 1, · · · , L.
in C, as the ratio of functions λ ± which are both holomorphic in C, and which are complex conjugate functions of one another, The irreducible solution for λ ± , in the sense of subsection 2.4.2, is unique and given by, We note that λ has a smooth limit as any one of its complex zeros approaches the real line, Im (s n ) → 0, so that the condition Im (s n ) > 0 could actually be smoothly relaxed to Im (s n ) ≥ 0, with the understanding that at least one of the zeros s n must satisfy Im (s n ) > 0.

The differentials ∂ w A ±
The holomorphic functions λ ± and the meromorphic differentials ∂ w A ± satisfy the relation, Eliminating λ ± between this relation and its complex conjugate relation obtained from (3.3), requires the doublet of differentials ∂ w A ± (w) to be proportional to the doublet of differentials ∂ w A ∓ (w) by a meromorphic function we shall denote γ(w). Involution of the complex conjugation operation then imposes a conjugation condition on γ, given as follows,

JHEP05(2017)131
To construct the differentials ∂ w A ± we analyze the zeros on both sides of (3.5). Since λ − has no zeros in the upper half plane, the set of zeros of ∂ w A + is the union of the set of zeros of λ + and the set of zeros of ∂ w A − (counted with multiplicities). But if the differentials ∂ w A ± have a common zero, then κ 2 = 0 must vanish there. Strict positivity of κ 2 in the interior of Σ prohibits such zeros. Thus, any zero common to ∂ w A ± must be on the real line. All poles must be common to ∂ w A ± since their ratio λ is bounded in the upper half plane. A simple pole in the interior of the upper half plane will produce a logarithmic branch cut in A ± . In appendix A we give an argument that no such poles occur and, by complex conjugation, that no poles occur in the lower half plane as well. Thus, all poles in ∂ w A ± must be common and on the real axis. We shall denote the positions of the poles by p ∈ R with = 1, · · · , L, for an as yet undetermined number L.
Putting all together, we find the following general expression for the differentials ∂ w A ± , The zeros s n satisfy Im (s n ) ≥ 0 for all n = 1, · · · , N with at least one zero having strictly positive imaginary part; the poles p are real for all = 1, · · · , L; and the constant ω 0 is complex. Regularity of the differentials ∂ w A ± at w = ∞ imposes the condition ∂ w A ± ∼ w −2 as w → ∞, which provides us with a relation between L and N , given by, The function γ(w) of (3.6) is then a constant phase, given by γ =ω 0 /ω 0 . In summary, the differentials ∂ w A ± given in (3.7) automatically satisfy κ 2 > 0 in the interior of Σ and κ 2 = 0 on ∂Σ = R. It remains to enforce the condition G = 0 on ∂Σ. To compute G, we need A ± and B which we compute in the subsequent subsections.

The functions A ±
The functions A ± may be obtained from the meromorphic differentials ∂ w A ± by integration. Performing this integration is greatly facilitated by expressing the meromorphic differentials of (3.7) in terms of their partial fraction decomposition with residues Z ± , The residues Z ± may be read off from the product formulas (3.7) and are given by,

JHEP05(2017)131
Since λ must have at least one zero in the upper half plane we have N ≥ 1 and the number of poles satisfies L ≥ 3. We now integrate the differentials to get A ± , as follows, where A 0 ± are two complex integration constants. Throughout, we shall choose the branch cuts of the logarithm so that ln(z) is real for z real and positive, and use the following continuation formula by circling the branch point z = 0 in the upper half plane, As a result, the functions A ± (w) have monodromy around each pole p , whose value is given by iπZ ± . In view of the third equation in (3.9), there is no monodromy at w = ∞, consistent with treating infinity as a regular point.

The functions B and G
Using (3.9) and (3.11), we integrate the definition for B given in (2.3), where B 0 is an arbitrary complex integration constant, w 0 is an arbitrary reference point, and we have introduced the following convenient abbreviation, A change in w 0 in (3.13) may be compensated by a change in B 0 , so that the parametrization provided above is subject to a redundancy which will be conveniently fixed shortly. Substituting the expressions for A ± ,Ā ± , B, andB into the formula for G given in (2.2), we obtain the following expression, where G 0 is a real constant given by, A change in w 0 will shift G 0 . The constants w 0 , γ, A 0 ± and B 0 , as well as the positions of the zeros s n and the poles p remain undetermined at this stage. The integrals are related to dilogarithms, and cannot be evaluated in terms of elementary functions.

Necessary conditions for vanishing G on ∂Σ
We shall now enforce the condition G = 0 on the boundary ∂Σ = R, and derive the corresponding conditions on the zeros s n , poles p and other parameters of the solution. Since the integrals on the last line of (3.15) are related to dilogarithms in w, the function G can vanish for all w ∈ R only provided the contribution of the dilogarithms vanishes by itself. This requires γ 2 = 1, given the branch cut choice we have made. Since γ is then real, the contribution of the second line in (3.15) vanishes for w ∈ R since the combination Z [ , ] is anti-symmetric in , , while the logarithmic factor becomes symmetric in , for real w. Cancellation of the terms on the first line of (3.15) for any w ∈ R further requires, The conjugation conditions on the residues then reads as follows for all = 1, · · · , L, and the combinations Z [ , ] are imaginary-valued. We shall require the following relation between the integration constants,Ā 0 Indeed, we may do so without loss of generality since the only supergravity field that involves the combinationsĀ 0 + +A 0 − is the flux field C, and the condition (3.19) simply amounts to a (constant) gauge choice for C. Having imposed (3.19) we see that the above conditions on Z ± precisely imply the conjugation relations of (2.23) advocated on general grounds.

Sufficient conditions for vanishing G on ∂Σ
Due to the presence of branch points in A ± on the real axis, the above conditions are necessary but not sufficient, however, as the monodromy makes the value of G jump across the branch points. To impose the further condition that G has zero discontinuity on the real axis across a pole p , we proceed as follows. Having chosen the branch cuts such that ln(z) is real for z real and positive, we may order the real poles p for = 1, · · · , L as follows, We shall choose the arbitrary reference point w 0 to be real with w 0 > p 1 . We then have G = 0 for w real and w > p 1 provided we set G 0 = 0.
We make the convergence of the combined integrals in (3.15) as w → ∞ explicit by anti-symmetrizing in , , and letting w 0 → ∞. The only dependence of the resulting formula for G on the integration constants A 0 ± is through the combination

JHEP05(2017)131
The resulting expression for G is given by, To derive the discontinuity of G across a pole we evaluate the function G as w approaches any point in the open interval (p k+1 , p k ) from above. To do so, it will be convenient to partition the sums over , according to their position with respect to k. The integrals do not contribute for the partition k < < since the logarithms are manifestly real there.
We evaluate the remaining integrals as a function of the position of w on the real axis with respect to the points p , and p , using the contours of figure 2. The integral evaluated with the help of contour (a) of figure 2 for w in the open interval [p , p ] is given by, The integral evaluated with the help of contour (b) of figure 2 for w in the open interval [−∞, p ] is given by, Substituting the evaluated integrals into the formula for G, we obtain after the straightforward cancellation of all w-dependence, valid for k = 1, · · · , L. Retaining the above condition for k = 1 unchanged, and subtracting from the condition for k ≥ 2 the condition obtained by letting k → k − 1, we obtain a total of L conditions, expressed below for k = 1, · · · , L, Both terms are purely imaginary in view of the conjugation property (3.18) and the definition (3.14). Only L − 1 of these conditions are linearly independent since the sum over all k of the above conditions vanishes by anti-symmetry in , k of the pre-factor Z [ ,k] . We also note that the final form of the conditions does not depend on a particular ordering of the poles on the real axis. Subtracting the k → k − 1 condition localizes the calculation around the pole and, while convenient for the derivation of the conditions, the assumption on the ordering can be safely relaxed in (3.25).
We close this subsection with the explicit expression for G which vanishes on the boundary (away from the poles). Using (3.25) to eliminate A 0 in (3.21) converts the first line in (3.21) to a double sum similar to the remaining terms. Anti-symmetrizing in , yields, where the complex conjugation applies to the entire expression inside the bracket. Note that G is real in view of the fact that Z [ , ] is purely imaginary.
and the functions A ± in (3.11) are invariant under SL(2, R) provided the data transform as follows, The residues Z ± in (3.9) are then invariant, as are the conditions of (3.25). The SL(2, R)invariance therefore allows us to fix the positions of three of the poles at will.

Counting free parameters
We label the families of solutions by a positive integer N which is the number of zeros s n of ∂ w A + in the upper half plane. The independent parameters of the solutions constructed in section 3 are then given by, Note that, even though λ 0 and ω 0 are complex, they satisfy |λ 0 | = 1 andω 0 = −ω 0 which accounts for them representing only one real parameter each. Furthermore, of the JHEP05(2017)131 two complex constants A 0 ± , only the combination 2A 0 = A 0 + −Ā 0 − enters into the metric factors and the dilaton-axion field, while in the flux potential, the combination A 0 + +Ā 0 − produces merely a gauge transformation of C. Therefore, of the constants A 0 ± only A 0 is a genuine physical parameter. Finally, in counting the number of poles p , we have made use of the relation L = N + 2 of (3.8) which relates the number of poles to the number of zeros. The total number of parameters is thus 3N + 6 = 3L.
The parameters of (3.28) must satisfy the L − 1 = N + 1 regularity relations obtained in (3.25), which reduces the number of free parameters by N +1. Furthermore, the SL(2, R) automorphism group of the upper half plane maps to physically equivalent solutions and therefore further reduces the number of free physical parameters by 3. As a result the number of moduli of the global solutions is given by,

Asymptotic behavior near the poles
The global solutions we have constructed above are regular everywhere in the interior of Σ, and on the boundary of Σ away from the poles p . In this subsection, we shall derive the asymptotic behavior of the supergravity fields near the poles p , and then provide a physical interpretation of this behavior in terms of five-brane webs. The original Ansatz and solutions were formulated in Einstein frame metric, where SU(1, 1) symmetry is manifest, and we shall begin by obtaining the asymptotics in this frame. Subsequently to make contact with brane behavior, we shall translate the asymptotics to the string frame.

Asymptotic behavior near the poles in Einstein frame
We will use the expressions for the Einstein frame metric factors f 2 6 , f 2 2 , ρ 2 in terms of the function G, its derivatives, and the functions κ 2 and R given in (2.4). To do so, we derive first the asymptotic behavior of κ 2 , G and R near a pole p m , and parametrize the coordinate w near a pole p m as follows, where 0 ≤ θ ≤ π and 0 < r |p m − p | for all = m. The expansion of ∂ w A ± near p m in terms of the parametrization by r, θ is then given by, which implies that the asymptotic expansion of κ 2 is given by, The coefficients κ m are real in view of the fact that κ 2 > 0 for r > 0 and 0 < θ < π, and may be taken to be positive.

JHEP05(2017)131
To calculate the asymptotic form of G we use (3.26). To exhibit the asymptotics of the integral near the pole p m we use the fact that G vanishes on the real line away from the pole and choose the point u = p m + r where G(u) = 0 to write G(w) = G(w) − G(u) and thus, The leading contributions are from the terms in the sum where either = m or = m.
Since the entire contour for the integral is close to the pole, we can then expand the integrand appropriately, and find, The leading behavior of the metric functions near the pole p m is given by, The dilaton φ and axion χ fields are obtained from (2.6) with To leading order, as r → 0, φ and χ are given by, Clearly, the dilaton diverges at the pole p m provided κ m = 0.

Asymptotic behavior near the poles in string frame
Since the dilaton diverges near the poles, there is a qualitative difference in behavior between the metric in the Einstein frame (where duality is manifest) and the metric in the string frame (where the connection with D-branes is transparent). We shall denote the JHEP05(2017)131 metric factors in the string frame byf 2 6 ,f 2 2 andρ 2 . Their relation with the metric factors f 2 6 , f 2 2 , ρ 2 of the Einstein frame, and their asymptotic behavior near the pole p m are given by, Converting the metric ds 2 in the Einstein frame of the Ansatz of (2.1) to the metric ds 2 in the string frame, and using the relation dwdw = dr 2 + r 2 dθ 2 , the near-pole metric in string frame becomes, The last two terms combine into the round metric of a smooth S 3 sphere without conical singularities provided we fix the value of c 6 , This produces a regular geometry in string frame. The curvature radius of the AdS 6 part diverges as the pole is approached, which turns it into a six-dimensional Minkowski space. 3 Since flat space does not have an intrinsic scale, we can absorb overall constants into a rescaling of the coordinates and find the near-pole string-frame metric, It is important to note that while our supergravity solutions are singular at the poles p , we shall argue in the discussion below that these singularities have a clear and compelling physical interpretation in terms (p, q) five-branes. This is to be contrasted with previously constructed Type IIB solutions [16][17][18].

Asymptotic behavior of the flux form near the poles
Finally, the shift in the flux potential C, as the pole at p m is being circumnavigated in the upper half plane from the right to the left, namely from θ = 0 to θ = π, is determined by examining (2.7). The variables ∂ w A ± , κ 2 , and G do not shift. The shift in ∂ w G cancels, and with A ± → A ± + iπZ m ± we find the shift in C, This expression is exact. Here, and below, we continue to exhibit c 6 in order to make the dependence of C on this parameter explicit, but for solutions with regular (p, q)-brane 3 To see this more explicitly, say we start with the Poincaré patch and ds 2 AdS 6 = dy 2 + e −2ky ηµν dx µ dx ν . Setting y =ỹ/ | ln r| and x µ =x µ / | ln r| yields, up to terms which are subleading in the near-pole limit, | ln r| ds 2 AdS 6 ≈ dỹ 2 + ηµν dx µ dx ν , and we indeed recover flat space.

JHEP05(2017)131
singularities at the poles, one should set c 6 = 1. Evaluating the asymptotic behavior of the flux potential itself we find, where we have used (3.40) to set c 2 6 to one. As r → 0, dr/r is regular and the leading neglected term vanishes. Finally, the factor sin 2 θdθ ∧ e 78 gives the volume form of the S 3 sphere in (3.41). Integrating (3.43) from θ = 0 to θ = π reproduces (3.42).

Identification of the poles with (p, q) 5-branes
With the asymptotic behavior of the solution and the shift in the flux potential around the poles in hand, we can now give a physical interpretation of the poles. To this end we compare the near-pole behavior to the well-known supergravity solution for a (p, q) 5-brane as given in [24]. With the string-frame metric (3.41), axion and dilaton (3.37) and complex three-form field strength (3.43), we find an exact correspondence with the geometry near a (p, q) 5-brane, with matching axion, dilaton, and 3-form field strength. More precisely, matching the 3-form field strength given above with the results presented in [24] and denoting p = q 1 Q and q = q 2 Q yields the identification, The near-pole behavior of our global solutions then matches the near-brane behavior of the (p, q) 5-brane solution not only in scaling and functional form of the various fields and involved geometry, but also in the overall coefficients. 4 This demonstrates that the poles in our solutions correspond to external 5-branes as they typically appear in the brane web diagrams. As we have already advocated briefly in [4], this evidence strongly suggests that the solutions correspond to fully localized intersections of (p, q) 5-branes.

Example solutions and relation to 5-brane webs
In this section we provide explicit realizations and corresponding numerical plots for the solutions on the upper half plane in the case where the number of poles is three and four. These numerical results are used to illustrate the validity of the general regularity arguments made in the previous sections and to provide explicit physically acceptable supergravity solutions. In addition, we confirm the claim, briefly advocated already in [4] and discussed in greater detail in the present paper, that the general L-pole solutions for the upper half plane can be identified precisely with fully localized intersections of (p, q) 5branes. We show that the corresponding features discussed in the previous sections indeed JHEP05(2017)131 admit a natural interpretation as properties of the intersection. In 4.3 we shall extend this discussion and speculate on how the solutions relate to (p, q) 5-brane webs and the corresponding SCFTs.

Solutions with three poles
In this subsection, we discuss the family of minimal solutions with N = 1 zeros and L = 3 poles in the differentials ∂ w A ± . It follows from the counting given in (3.29) that this family of solutions has four real free parameters. Using the SL(2, R) automorphism group of the upper half plane, we may choose the poles to reside at 0, ±1. As argued in the previous section, the single zero s must lie strictly in the upper half plane, The four real parameters are the complex zero s and the overall complex normalization ω 0 λ 0 . The residues are calculated from (3.10), and we obtain, the residues Z − being given by Z − = −Z + for = 1, 2, 3. We then have, and the relations (3.25) are solved by A 0 = ω 0 λ 0 s ln 2. The metric functions, dilaton/axion fields and flux potential all involve polylogarithms and do not admit a simple form. Therefore, we shall directly resort to a numerical evaluation of the supergravity fields of the solution. In figure 3 we illustrate the behavior of the metric functions in the Einstein frame, as well as the axion, dilaton and the two-form potential C for the following choice of the four free parameters, For stable numerical evaluation we note that, since G assumes a maximum in the interior of Σ and ∂ w G = 0 at that point, the expression for ρ 2 in (2.4) has to be treated with care. The expression is regular since R → 0 when ∂ w G → 0, and to make the regularity of ρ 2 explicit we rewrite it in a manifestly regular form. This leads to the following exact alternative expressions for the metric factors, The plots show that the exponentiated dilaton is positive, as required, and goes to zero at the poles, as derived in section 3.9. The metric factors are real and positive, as desired, JHEP05(2017)131

JHEP05(2017)131
and behave at the poles precisely as discussed for the Einstein frame in section 3.9. At the boundary of Σ, f 2 vanishes as required to have a smooth ten-dimensional geometry without boundary, and f 6 stays finite except for at the poles. The two-form potential C is piecewise constant at the boundary, and jumps at the poles in accordance with (3.42). Moreover, it approaches the same value to the left and to the right of all poles, reflecting charge conservation.

Solutions with four poles
Solutions with four poles have six real free parameters. We use SL(2, R) invariance to fix the positions of three of the four poles on the real line. We may choose the remaining free parameters to be the two complex zeros and one complex overall normalization, with A 0 and the position of the fourth pole determined from (3.25). A more physically motivated choice of parameters is as follows. Via the identification in section 3.10, the poles correspond to four stacks of semi-infinite external 5-branes, and the residues are related to their (p, q) charges. Alternatively, we may choose three complex residues at three of the poles as the free parameters, the residue at the fourth pole then being determined by the condition that the sum of the four residues vanishes in (3.9). We will look at the special family of solutions where the residues are related as follows, That is, the stacks of semi-infinite external 5-branes have pairwise opposite charges. Due to the relation Z + = 0 the two relations above are equivalent to one another. Using SL(2, R) to fix the position of three poles as follows, the position of the remaining pole, p 4 remains a free parameter. The conditions (4.6) and (3.25) are solved by, Since the relation between s 1 and s 2 is an SL(2, R) transformation, both s 1 , s 2 will be in the upper half plane provided either one is. The resulting charges are, with Z 3 + and Z 4 + given by (4.6). We will take a closer look at the particular solutions obtained by choosing Plots illustrating the solution are shown in figure 4. The qualitative features are the same as for the 3-pole solution: the metric factors are positive and the fields satisfy the desired regularity conditions throughout, except for the poles where they behave as derived in section 3.9.

Relation to 5-brane webs
We have argued already briefly in [4] that our supergravity solutions can be identified with fully localized intersections of (p, q) 5-branes, and the detailed derivations presented in this paper provide the justification for the arguments used in [4]: by construction, our supergravity solutions have the correct superconformal F (4) symmetry, and by the results of section 3.8 the parameter count precisely matches the parameter count of (p, q) fivebrane webs in the conformal limit. Moreover, as the discussion of section 3.10 shows the solutions also have the correct external states and the parameters can be directly translated to those specifying a 5-brane intersection. Many other features of the solutions admit a natural interpretation in the context of 5-brane intersections as well. The minimal number of poles being three, as discussed in section 3.3, corresponds to the fact that three external branes are needed to produce a codimension-1 intersection. Moreover, there are no solutions with either only D5 charge or only NS5 charge. From (3.14) we have Z [ ,k] = 2iIm (Z k + Z + ). Therefore, if the Z + are either all real or all imaginary, Z [ ,k] vanishes for all , k, which implies that G = 0 everywhere and the solution degenerates. Regular solutions therefore necessarily involve D5 and NS5 charge, and this corresponds to the fact that D5 and NS5 charges are needed to realize a codimension-1 intersection with 5-branes.
The identification of our solutions with 5-brane intersections certainly suggests a holographic relation to the 5d SCFTs obtained by taking the conformal limit of 5-brane webs describing 5d supersymmetric gauge theories [12,13]. Holographic relations usually involve some form of a large-N limit, and we indeed note that the charges of the external 5-branes in our solutions are assumed to be large for the supergravity description to be valid. There is therefore no constraint from charge quantization, similarly to the familiar case of e.g. N = 4 SYM and its AdS 5 × S 5 dual. Generically, the brane intersections described by our solutions therefore involve both, large D5 charge and large NS5 charge. Brane webs related by the SL(2, Z) duality of type IIB string theory describe the same field theory, and in type IIB supergravity this is enhanced to an SL(2, R) symmetry. 5 We thus expect solutions related by SU(1, 1) ∼ SL(2, R) to describe field theories with the same "large-N " limit.  , and ρ 2 , the real and imaginary parts of the two-form potential C and of the axion and dilaton corresponding to the parameters given in (4.10), with N = 2, M = 3. All supergravity functions are regular throughout Σ except for the poles, where their behavior is precisely as discussed in section 3.9.

JHEP05(2017)131
N N N Figure 5. The N -junction, a (p, q) 5-brane intersection with three external 5-brane stacks, obtained by taking N copies of the basic junction of a D5 and an NS5 brane.
We will now use the explicit solutions with three and four poles for a more explicit discussion and to illustrate further points. The solutions with three poles correspond to 5-brane intersections with three external (p, q) 5-branes, with charges given by (3.44) with (3.40). They have a total number of four parameters, corresponding to the choice of charges subject to charge conservation. For the field theory interpretation, SU(1, 1) ∼ SL(2, R) further reduces the number of parameters by 3, which leaves only one free parameter. Indeed, a generic 3-fold intersection of (p, q) 5-branes can be mapped to the form shown in figure 5 by an SU(1, 1) transformation. This "N -junction" is obtained by combining N copies of the basic 5-brane junction with external charges (1, 0), (0, 1) and (−1, −1) (in all in-going convention) and allows for a large-N limit. A choice of parameters for the supergravity solution to realize this form of the web via (4.1), (4.2) is Brane webs of the form shown in figure 5 have been identified in [25] as five-dimensional uplifts of the four-dimensional T [A N −1 ] theories, obtained by wrapping N M5-branes on a sphere with 3 punctures. For the field theory interpretation it is crucial whether and how the external 5-branes defining the intersection end on 7-branes. In [12,13] the external 5-branes were taken to be semi-infinite, but one can also terminate them on 7-branes such that they are of finite extent in the plane in which the brane web is drawn [14]. Taking groups of 5-branes to end on the same 7-brane yields a different field theory than having each 5-brane semi-infinite or terminate on its own 7-brane, and for the intersection of figure 5 various options were discussed in [25]. In our solutions we only see the external 5-brane geometries with no indication for the presence of 7-branes, and at least at the level of disk solutions there are also no obvious moduli corresponding to the choice which 5-branes terminate on which 7-brane. This suggests that the solutions correspond to the original form of the brane webs and intersections, where the external 5-branes are indeed semi-infinite [12,13]. However, one could also argue that our solutions only cover the near-intersection limit, where the external 7-branes may not be directly accessible and more indirect methods will be required to precisely pin down the dual field theory. We will leave a more detailed investigation for the future and leave this question open for the remaining discussion.
The lift of the isolated 4d T [A N −1 ] theories to 5d via intersections of the form in figure 5 offers a chance to find deformations that permit a Lagrangian description as gauge theories in the IR, and such deformations have been constructed in [26]. For the case where each external 5-brane ends on a separate 7-brane, it results in a web describing the quiver That is, a gauge theory with product gauge group SU(N − 1) × · · · × SU(2), with one hypermultiplet in the bi-fundamental representation for each pair of adjacent gauge group factors and in addition N massless hypermultiplets in the fundamental representation of SU(N −1) along with 2 massless hypermultiplets in the fundamental of SU(2). The 5-brane intersection only has one large parameter controlling both D5 and NS5 charge, and we see that in this case it also translates to two large parameters in the field theory defining the UV fixed point: the length of the quiver and the rank of the largest gauge group factor. Moving on to 4 external branes, the 5-brane intersection with the external charges of section 4.2 is shown on the left hand side in figure 6. Once again the field theory interpretation depends on whether and how the external 5-branes end on 7-branes. Taking N = M and all 5-branes within each external 5-brane stack to end on the same 7-brane would yield a 5-brane construction for the USp(N ) theory [8], while the original form of the webs with no 7-branes again leads to long quivers: the configuration is SL(2, R) dual (up to a rescaling of the charges M and N ) to the intersection shown on the right hand side in figure 6. Without the introduction of 7-branes this intersection has been discussed already in [13], and a deformation of the fixed-point SCFT to a gauge theory is described by the quiver N − SU(N ) × · · · × SU(N )

JHEP05(2017)131
This example explicitly exhibits the presence of two independently large parameters, the rank of the gauge group and the length of the quiver, corresponding to the D5 charge and NS5 charge in the brane intersection. Moreover, we note that the number of matter fields is directly linked to the rank of the gauge group factors, and the large-N limit is a form of a Veneziano limit rather than a 't Hooft limit.

The annulus
In this section, we shall investigate the existence of supergravity solutions when Σ has the topology of an annulus, or equivalently of a finite cylinder. This is the next-simplest topology after the case of the disk, since the annulus has genus zero but two boundary components. We shall provide an explicit parametrization of the differentials ∂ w A ± under the assumption that these differentials are single-valued in Σ and no axion monodromy is allowed. We will explicitly construct the functions A ± as well and derive the general physical regularity conditions. We have used a combination of analytical and numerical methods to explore whether the physical regularity conditions on the supergravity fields can all be satisfied simultaneously for the annulus. These explorations are not exhaustive, but the outcome so far has consistently been negative, and no example of a physically regular solution has been found to date. We have no analytical proof that such solutions cannot exist, but the regularity conditions are structurally different when there are multiple boundary components and the negative results so far may well be due to the fact that such solutions do generally not exist. We will leave a more systematic analysis for the future and turn to Riemann surfaces with non-trivial topology and a single boundary component in the next section.

Parametrization of the annulus
Compared to the situation of the disk, the annulus has two new features. First, the annulus has two disconnected boundary components instead of one for the disk; second the annulus has a non-trivial fundamental group π 1 (Σ) = Z.
The function theory on the annulus, which we need to construct the differentials ∂ w A ± and the associated functions A ± , is conveniently obtained in terms of the function theory on the double surface. The annulus Σ may be represented in the complex plane by a rectangle with two opposing edges periodically identified. The double surfaceΣ of the annulus Σ is a torus whose periods may be chosen to be 1 and τ and whose modulus is purely imaginary τ = iτ 2 with τ 2 ∈ R + . The surface Σ and its boundary ∂Σ may then be represented by, both periodically identified under w ≡ w + 1, as represented in figure 7. We chooseΣ by symmetry across the real axis, such that 0 ≤ Re (w) ≤ 1 and |Im (w)| ≤ τ 2 /2. Complex conjugation I is the anti-conformal involution which maps between the components ofΣ in the upper and lower half planes. We may view Σ as the quotient Σ =Σ/I, and the boundary ∂Σ as the fixed set under I.

JHEP05(2017)131
The annulus Σ is represented by the dark blue rectangle with periodic identification w ≡ w + 1; its boundary ∂Σ has two components denoted by thick black lines; the extension to the double surfaceΣ is indicated by the light blue rectangle; four zeros s n of λ in the interior of Σ, and their complex conjugatess n , are indicated by black dots, while four poles p n distributed amongst the two boundary components are indicated by crosses.

The function λ
To investigate the existence of physically regular supergravity solutions, we closely follow the general strategy outlined in subsection 2.4, and begin with the construction of the holomorphic function λ. To construct λ(w|τ ) we use the scalar Green function G(w, z|τ ) = G(z, w|τ ) on Σ which vanishes whenever w is on ∂Σ. To construct G, we make use of the scalar Green function G 0 on the double surfaceΣ = C/(Z + Zτ ), where ϑ 1 is the Jacobi ϑ-function. By construction, G 0 is a real and doubly periodic function, in view of the standard translation properties of the Jacobi ϑ-function, The corresponding Green function G on the annulus is given by, The properties of reality and double periodicity of (5.3) ensure G(w, s|τ ) = 0 whenever w ∈ R + Z τ 2 , so that G vanishes on both boundary components of Σ.

JHEP05(2017)131
The function λ(w|τ ) for the annulus with modulus τ ∈ iR + may be inferred from the electrostatic potential obtained by adding the contributions from an array of positive unit charges placed at points s n ∈ Σ with n = 1, · · · , N , and satisfying s n ∈ R + Z τ 2 . Therefore, the general electrostatic potential used in (2.18) takes the from, To extract the holomorphic function λ, we use the fact that G is real and harmonic in w ∈ Σ away from s to split (5.6) into a sum of holomorphic and anti-holomorphic parts in w, leaving a constant phase factor λ 0 undetermined. This expression is automatically singlevalued under w → w + τ , and is single-valued under w → w + 1 provided we impose the following condition on the points s n , N n=1 (s n −s n ) ∈ Zτ (5.8) We recognize this relation as the divisor condition for meromorphic functions on the torus, applied to the special case where the zeros and poles come in complex conjugate pairs. The condition implies that we must have N ≥ 2, since it is clearly unattainable with N = 1 in view of the requirements that at least one of the zeros s n should lie in the interior of Σ.

The differentials ∂ w A ±
To construct the meromorphic differentials ∂ w A ± , we begin by splitting λ into two complex conjugate functions λ ± , The splitting is not unique, as a common factor of a real function to ∂ w A ± cancels out in λ. The irreducible solution is given as follows, Neither function λ ± is single-valued onΣ. The differentials ∂ w A ± are given in terms of a meromorphic function ϕ onΣ by (2.21) so that,

JHEP05(2017)131
Applying the general conditions on the meromorphic 1-forms ∂ w A ± of subsection 2.4.2 to the case of the annulus allows us to narrow the choices of ϕ as follows.
Since ∂ w A − cannot vanish in the interior of Σ in view of the condition κ 2 > 0, the function ϕ cannot have zeros in the interior of Σ. Furthermore, ∂ w A ± cannot have poles in the interior of Σ since otherwise the functions A ± would have monodromy in the interior of Σ. Therefore, all zeros and poles of ϕ must be on the boundary ∂Σ. Zeros on the boundary may be viewed as corresponding to a degenerate situation of the general case where all zeros a n are in the interior, and we shall therefore assume that no zeros occur on the boundary, so that ϕ has no zeros at all. Since meromorphic differentials on the toruŝ Σ have equal numbers of zeros and poles, ϕ must have precisely N poles onΣ. As a result, the zeros of ∂ w A + in the interior of Σ are precisely the N zeros s n of λ + .
Finally, we shall adopt the condition (2.22) which for the torus reads, With the above assumptions, we parametrize ϕ by its N poles at points p n with n = 1, · · · , N distributed amongst the two boundary components of ∂Σ given in (5.1). Combining the above conditions, we find the following expressions for the differentials ∂ w A ± , where ω 0 is a complex constant while Λ ± are given by, (s n − p n ) (5.14) By construction, and with the help of (5.4), the differentials ∂ w A ± are invariant under w → w + τ , but their invariance under w → w + 1 requires the extra conditions, These conditions on Λ ± combined imply the divisor relation (5.8) derived earlier.
It remains to enforce the conjugation condition of (5.12). In the special case where all the poles p n lie on the real boundary component, the condition amounts to requirinḡ ω 0 = −ω 0 . In the general case when poles are allowed to lie on both boundary components, enforcing (5.12) is more delicate, and we have instead the general relation, where ω 0 is constant and is subject to the relationω 0 = −ω 0 . Satisfying the combination of the conditions (5.15) and (5.16) thus provides single-valued meromorphic differentials ∂ w A ± on the annulus which satisfy the conjugation condition (5.12).

The functions A ±
In this subsection, we shall integrate the differentials ∂ w A ± to obtain the functions A ± . The product representation (5.13) obtained in the preceding subsection is inconvenient to carry out this integration. Instead, we shall derive here an equivalent representation as a sum over meromorphic Abelian differentials which may be easily integrated.
Since the meromorphic differentials ∂ w A ± of (5.13) are single-valued on the double surfaceΣ, the sum of their residues Z ± at the poles p must vanish, and the differentials can be expressed as a sum over meromorphic Abelian differentials ∂ w ln ϑ 1 (w − p |τ ), which are single-valued under w → w + 1, and transform under w → w + τ by a constant shift, The residues Z ± may be obtained from (5.13) by matching poles and are given by, The presence of the constants α ± is due to the fact that the torusΣ has a one-dimensional space of holomorphic Abelian differentials which, in local complex coordinates w, is generated by the differential dw. A meromorphic one-form such as ∂ w A ± (w) will generically have a non-trivial component along dw. The values of α ± may be obtained by using (5.12) and the vanishing of ∂ w A + at any of the zeros s * ∈ {s 1 , · · · , s N }, The integrals of ∂ w A ± are now easily computed, and we have, where A 0 ± are complex integration constants. To secure useful conjugation properties, care is needed in the choice of branch cuts for the logarithm. Given a choice of branch cut in A + , the proper branch cut in A − is actually JHEP05(2017)131 Figure 8. The annulus Σ is represented as in figure 7. The contours of integration along the upper and lower boundary components of ∂Σ are indicated with thick red lines. more properly written as, where the branch cuts chosen in A ± are now the same. With this choice, we may use a constant gauge transformation of the flux field C to set A 0 − +Ā 0 + to zero, just as we did in the case of the upper half plane, upon which we obtain the following conjugation relation, In the next subsection, we obtain the conditions under which G = 0 on the boundary.

Conditions for G = 0 on the boundary
The general arguments of subsection 2.4.2 along with the conjugation conditions (5.12) and (5.22) lead us to conclude that G is constant along any line segment free of poles p on the boundary ∂Σ. To enforce the boundary condition G = 0 on the entire boundary, it will therefore suffice to enforce G = 0 on any one line segment free of poles on each boundary, and then require that the monodromy across every pole on that boundary component vanishes. We begin by requiring the vanishing of the monodromy of G across an arbitrary pole p k along the contour illustrated in figure 8. To keep track of the behavior of the functions A ± near their branch cuts, we introduce ε k defined by, where δ k = 1 for poles on the real line and δ k = −1 for poles on the second boundary component. The resulting ε k has a small positive real part and an imaginary part reflecting the contours in figure 8. To find the jump conditions we start out from, with C k the half circle part of the contour around the pole p k . For the change in A ± we find which leads to Working this out more explicitly yields Using p = p − (1 − δ )τ /2 and the quasi-periodicity of the ϑ 1 functions, the log terms in the round brackets can be evaluated more explicitly. The integrals in ∆ k B + ∆ kB for the jump are restricted to a half circle around the pole, which means the integrand can be expanded for |w − p k | 1. Upon explicit evaluation we find that the contribution from ∆ k B + ∆ kB precisely matches the change in |A + (w|τ )| 2 − |A − (w|τ )| 2 and merely produces an overall factor of 2. Altogether, we find The sum then evaluates to The sum manifestly vanishes if α + = 0 and all poles are on one boundary component. In this case, the conditions (5.29) simplify considerably, and we have, which is analogous to the condition for the disk in (3.25).

Minimal number of poles
For the case of the disk, a minimum of three poles in ∂ w A ± was required in order to have a single zero in the upper half plane, lest the solution be trivial. For the case of the annulus, the number of zeros of ∂ w A ± must equal the number of its poles instead. We shall now analyze the minimal number of zeros and poles required for the annulus case. When all the poles p of ∂ w A ± are on one boundary component the requirement Λ + = N n=1 (s n − p n ) ∈ Zτ implies Λ − ∈ Zτ , and hence both conditions in eq. (5.15). Periodicity of the functions A ± as w → w + 1 requires α ± = 0. This is one complex condition due to α − = −ᾱ + , and is equivalent to the condition, where s is any one of the zeros of ∂ w A ± (cf. (5.19)). The condition Λ + ∈ Zτ has no solutions with two poles on the same boundary, since this would require two zeros in the interior of Σ, whose imaginary parts add up to a non-zero multiple of τ by conditions (5.15). The only possibility would be to have their imaginary parts both equal to τ 2 /2, but this implies in turn that both zeros are on the boundary of Σ leading to an unphysical solution.
The condition Λ + ∈ Zτ can be solved with two poles and two zeros inside the annulus, provided we have one pole on each boundary component. Without loss of generality, we use translation invariance to set p 1 = 0 and p 2 ∈ R + τ 2 . Assuming furthermore that α + = 0, the condition (5.30) with k Z k ± = 0 then implies Re (p 2 ) = 0 and thus p 2 = τ /2. The condition α + = 0, via (5.32) and with k Z k ± = 0, then implies, This equation only has two unacceptable solutions below the real line, and no solutions for s inside the annulus. We thus find that there are no solutions with two poles, and that at least three poles are a necessary condition for a physically regular solution to exist.

Investigating solutions with at least three poles
Starting with 3 poles it becomes considerably more involved to disentangle the collection of conditions developed earlier. While we cannot offer an exhaustive analysis here, we shall present the conclusions we draw from a partial numerical analysis we have carried out. With N poles for N ≥ 3, the divisor conditions (5.15) allow for all poles to be distributed over both boundary components including the case where all N poles are on one boundary component. Our numerical analysis of the conditions α ± = 0, which are required to make A ± single-valued under w → w + 1, in the case of N = 3, 4, 5 poles appears to exclude systematically the cases where all poles of ∂ w A ± are not on the same boundary. Indeed, when poles occur on both boundary components, we find that the condition α ± = 0, for a given distribution of poles on the boundary and for N −2 zeros inside Σ, always forces one of the remaining zeros of ∂ w A + to be outside of Σ. So far, we have no analytical proof of this claim.

JHEP05(2017)131
With N ≥ 3 poles, and all N poles on a single boundary component, it is possible to satisfy the conditions α ± = 0 and have all N zeros of ∂ w A + in the interior of Σ. A simple family of such solutions may be obtained by choosing the distribution of poles and zeros to be invariant under Z N translations, and given by, For given N the zeros are well in the interior of Σ. Therefore, by continuity, we know that there will exist an open set in the space of all solutions which contains the above solution as a point, and has maximal dimension, given by the real coordinates of poles N , the complex coordinates of N zeros, minus the complex condition α + = 0, and minus the real divisor condition, adding up to 3N − 3.
Having established that there are at least the partial solutions exhibited above, we more generally explored configurations with all poles on a single boundary, and showed numerically that the conditions for continuity of the function G across the poles, namely the conditions ∆ n G = 0 for n = 1, . . . , N , may be solved as well. The solutions to all these conditions combined now provide configurations for which all zeros of ∂ w A + are in the interior of Σ, the condition α + = 0 is satisfied, and G is constant on each boundary component.
There now remain two further conditions to be implemented, namely that G = 0 on both boundary components. There is a natural free parameter, namely the integration constant B 0 of the composite holomorphic function B which may be used to set G = 0 on one boundary component. The value of G on the other boundary component is then determined by the parameters of the solution. Our numerical analysis for the case of N = 3, 4, 5 poles appears to show that the condition G = 0 can not be satisfied simultaneously on both boundary components. Rather, the difference between the values of G on the boundary with no poles and on the boundary with all poles appears to consistently be positive. Again, at this time, we have no analytical proof of this claim, but hope to return to this question in the future.
If the indications from the numerical analyses carried out so far are correct, then the annulus may not support physically regular solutions in the absence of axion monodromy. More generally, we may speculate that the obstruction to the existence of solutions resides in the presence of more than one boundary component, in the case of the annulus as well as for surfaces of arbitrary genus.

Riemann surfaces of arbitrary topology
In this last section, we shall construct suitable differentials ∂ w A ± for general Riemann surfaces, and spell out the conditions required for the corresponding supergravity solutions to be physically regular. The final equations are so complicated, however, that so far we have not succeeded in constructing acceptable solutions, even numerically.
Consider an orientable Riemann surface Σ of genus g with ν ≥ 1 boundary components which topologically are circles. Functions and differential forms on Σ may be constructed in terms of functions and differential forms on the double surfaceΣ, equipped with an anti-conformal involution which we shall denote by I (for a standard reference, see [27]).

JHEP05(2017)131
The boundary ∂Σ of Σ is fixed, point by point, under I so that I(∂Σ) = ∂Σ. We may view Σ as the quotient Σ =Σ/I. The genusĝ ofΣ is related to the genus g and the number of boundary components ν of the original surface Σ by the standard relation, This construction is depicted in figure 9 for the case g = 2 and ν = 3. The conventions for the basis of homology cycles indicated in the figure will be given in section 6.2 below.

Generalizing the electrostatics analogy
The electrostatics analogy, used to construct the function λ for the cases of the upper half plane and the annulus, may be generalized without complications to the case of an oriented Riemann surface of arbitrary genus g, and an arbitrary number ν of boundary components the topology of each of which is that of a circle. In the electrostatics analogy, we seek first to construct an electrostatic potential − ln |λ| 2 which vanishes on the boundary ∂Σ and is strictly positive everywhere in the interior of Σ. To achieve such an electrostatic potential, we ground the system to zero potential on every component of ∂Σ and place an arbitrary arrangement of positive charges in the interior of Σ. By the min-max principle for harmonic functions, used already earlier in subsection 2.4, the potential − ln |λ| 2 is then guaranteed to be positive everywhere in the interior of Σ. Mathematically, we start from the original Riemann surface Σ of genus g with ν boundary components and construct its doubleΣ endowed with a conformal involution I. To obtain a potential which is strictly positive in the interior of Σ, and which vanishes on the boundary ∂Σ, we place an arrangement of an arbitrary number N of positive electric charges q n at arbitrary points s n in the interior of Σ, with n = 1, · · · , N , and place mirror charges −q n at the mirror locations I(s n ). Since by construction the potential − ln |λ| 2 is odd under the involution I, it is guaranteed to vanish on the boundary ∂Σ. Since all charges on Σ are positive, and the potential vanishes on ∂Σ, we use again the min-max principle for harmonic functions to argue that − ln |λ| 2 must then be strictly positive in the interior of Σ, just as we had already done in the case of the upper half plane and the annulus. Clearly, such solutions will always exist. To obtain explicit formulas for λ, we shall provide the necessary mathematical set-up in the next subsection.

Double surfaceΣ and involution I
We shall sort the homology cycles on the double surfaceΣ, and their dual holomorphic 1forms, according to their parity under I. Our labeling of the cycles generalizes to arbitrary g and ν the labeling indicated in figure 9 for g = 2 and ν = 3, and we have [27]  There is one further boundary component, A g+ν which is homologically trivial onΣ. Throughout, the indices I, I , i will run over the ranges given in (6.2) and we shall conveniently use the composite indexÎ = (I, I , i). We define the involution matrix S by, where I g and I ν−1 are the identity matrices respectively in dimensions g and ν − 1. The cycles may be arranged to have definite parity under I, which we choose as follows, The normalization of the holomorphic 1-forms on the A-cycles and their corresponding behavior under the pull-back I * of the involution I to 1-forms, are given as follows, whereωĴ is the complex conjugate of ωĴ . The period matrix ΩÎĴ ofΣ is defined by, BÎ ωĴ = ΩÎĴ (6.6)

JHEP05(2017)131
Its expression in terms of the prime form and Abelian integrals simplifies as follows, ωĴ (6.14) In order to produce a function λ without branch cuts the positive charges should be chosen to be all equal to unity, giving the following general solution for the electrostatic potential, To split the holomorphic dependence in w from the anti-holomorphic dependence we begin by splitting the line integral using (6.7), Splitting off the holomorphic dependence in w now clearly amounts to splitting off the dependence in w from the dependence in I(w). Using furthermore the fact that S 2 = I and that S commutes with Y , and exponentiating the holomorphic part to obtain λ, we find, While |λ(w|Ω)| was single-valued onΣ, the splitting into λ(w|Ω) and its complex conjugate generally introduces monodromy by phase factors. To obtain a single-valued λ will necessitate divisor relations on the zeros s n , which we now obtain. The monodromy properties of the function λ(w|Ω) as w is moved around the homology cycles AÎ and BÎ onΣ are readily evaluated using (6.11). We state here, without detailed derivation, the conditions for the absence of monodromies around all homology cycles, where the column matrices m, n have integer entries, restricted in the following way, where the integers m I , n I for I = 1, · · · , g and n j for j = g + 1, · · · , g + ν − 1 are arbitrary. These conditions state that the zeros {s 1 , · · · , s N } and poles {s 1 , · · · ,s N } indeed satisfy the customary divisor condition of a meromorphic function λ on a Riemann surfaceΣ which possesses an involution I (see for example [27,28]).

The differentials ∂ w A ±
To construct the differentials ∂ w A ± we split the meromorphic function λ onΣ into the ratio of two conjugate multiple-valued holomorphic differentials onΣ as follows, The differential form λ + has its zeros at s n , while λ − has its zeros at I(s n ). The splitting of the prime form factors in (6.17) is manifest. To split the exponential factor, we make use of (6.16) but with s n replacing w. The splitting of λ is now manifest, and we find, where the point w 0 lies on the boundary of Σ, so that we have I(w 0 ) = w 0 . As constructed above, the differential forms λ ± have weight (−N/2, 0) in w. Next, we construct the holomorphic differentials ∂ w A ± to be well-defined and singlevalued onΣ. The basic equations for ∂ w A ± are as follows, The differential form ϕ must have weight (1 + N/2, 0) and has neither zeros nor poles in the interior of Σ, so that its poles and zeros are on ∂Σ. As in the cases of the upper half plane and the annulus, we shall view any real zeros as degenerations of pairs of conjugate zeros (s n , I(s n )) to the boundary. We shall denote the poles by p with = 1, · · · , L. To have well-defined and single-valued meromorphic differentials ∂ w A ± onΣ, the number of their poles and zeros are related as follows, Forĝ ≥ 2, one may view the positions of the L poles and of N −ĝ + 1 of the zeros as arbitrary, while the remainingĝ − 1 zeros are determined by the divisor condition for forms of weight (1, 0). Taking these considerations into account, we obtain the following expressions for ∂ w A ± , where Λ ± J are given by, The holomorphic form σ(w) has weight (ĝ/2, 0) and has neither zeros nor poles. Its role in (6.24) is to guarantee that the forms ∂ w A ± have the weight (1, 0) in w, and are singlevalued. The monodromy of the form σ(w) around A-cycles onΣ is trivial; its monodromy around B-cycles and its expression in terms of the prime form are given below, Here, ∆Î (w 0 ) is the Riemann vector on the surfaceΣ. Choosing the base point w 0 to be invariant under the involution I, the Riemann vector is real in the Jacobian Cĝ/(Zĝ + ΩZĝ). The conjugation property of σ under conjugation follows that of the prime form and we have, I * (σ(w)) = σ(I(w)) = σ(w) (6.27) Using it, we readily establish the conjugation relations for the differentials, ∂ w A ± (w) = −∂ I(w) A ∓ (I(w)) = −∂ w A ∓ (I(w)) (6.28) The transformation laws given above for the prime form, the exponential in (6.24), and σ allow us to compute the monodromies of ∂ w A + around AÎ and BÎ cycles. Requiring the monodromies to vanish around all cycles imposes the following conditions on Λ ± where m ± and n ± are column matrices whose entries are integers which satisfy, where the integers m ± I , n ± I for I = 1, · · · , g and n ± j for j = g + 1, · · · , g + ν − 1 are arbitrary. These conditions imply the relations for single-valuedness of λ of (6.18) with, The total number in (6.29) is 2ĝ real conditions. JHEP05(2017)131

The functions A ±
To integrate ∂ w A ± we decompose these differentials onto the following Abelian differentials of the third kind, with simple poles at w = z and w = w 0 with unit residues of opposite signs. In view of the cancellation of the sum of the residues Z ± at the poles p , all dependence on w 0 cancels out, and we have the following decomposition, where αÎ ± are the coefficients of the holomorphic one-forms ωÎ onΣ. The residues satisfy, and their values may be read off from the product representations in (6.24), The coefficients αÎ ± may be obtained by evaluating (6.33) atĝ of the zeros, denoted by s n , of ∂ w A + andĝ of the zeros I(s n ) of ∂ w A − , and then solving eachĝ ×ĝ linear system in turn. The functions A ± are obtained by integrating (6.33), The functions A ± do not necessarily have to be single-valued onΣ. Indeed, monodromy of A ± by constant shifts, as given in (2.10) subject to the condition (2.11), is allowed as it gives rise to single-valued supergravity fields. We shall now investigate the conditions this

JHEP05(2017)131
kind of allowed monodromy imposes on the parameters of the holomorphic functions A ± . First, by comparing the complex conjugate of the second line in (6.36) with the first line of (6.36), we obtain the relation, or simplyᾱ − = −Sα + in matrix notation. When the indexÎ = i corresponds to a boundary cycle, this relation generalizes the one obtained for the annulus in (5.19). Second, the monodromies around cycles AÎ and BÎ are given by, There are no conditions on the monodromy around B i cycles, as such cycles extend beyond the surface Σ. Thus, the conditions for allowed monodromy around these cycles are given by, or simplyᾱ − = α + and (I − S)(β − − β + ) = 0 in matrix notation. Combining these relations, we obtain the following constraints on the parameters of the solution, expressed here in matrix notation,ᾱ while the second set of equations gives α I ± in terms of the residues Z k ± and the poles p k . We see that, by contrast with the case of the annulus, for higher genus the parameters α ± need not be forced to vanish by the monodromy conditions alone.

The conditions required for G = 0 on the boundary
The general arguments of subsection 2.4.2 guarantee that G is piecewise constant on any line segment of the boundary which is free of poles. The condition for continuity of the function G across a pole p is given by the same arguments as in the case of the upper half plane and the annulus. Here, we shall specialize to the case where we have a single boundary component and hence ν = 1, since the expressions considerably simplify in this case. We shall also assume that the condition for single-valuedness around AÎ cycles has been imposed, so that αÎ ± = 0. The conditions are local around p , and are given as follows,

JHEP05(2017)131
When only a single boundary component is present, the vanishing of G on ∂Σ may be achieved by solving for the integration constant B 0 of the holomorphic function B. When ∂Σ has ν ≥ 2 boundary components, there are ν − 1 further real conditions required to guarantee the vanishing of G on each boundary component. This completes the construction of the ansatz for general Riemann surfaces and the derivation of the regularity conditions. We leave a detailed investigation of the existence of solutions with non-trivial topology for the future and close by noting that, based on the numerical investigations for the annulus described in subsection 5.7, we may speculate that no solutions will exist if Σ has multiple boundary components.

Discussion
In this paper we have constructed a large family of Type IIB supergravity solutions which are candidates for holographic duals to five-dimensional superconformal field theories. The space-time of these solutions consists of the manifold AdS 6 × S 2 warped over a twodimensional Riemann surface Σ with boundary. The supergravity solutions preserve sixteen supersymmetries and give a holographic realization of the F (4) superconformal algebra. Unlike corresponding solutions in dimensions 3, 4 and 6, there is no asymptotic enhancement of the supersymmetry, in accord with the fact that there exist no five-dimensional superconformal theories with thirty two supersymmetries.
The solutions become singular at isolated points on the boundary of the Riemann surface, where the poles in the differentials ∂ w A ± which characterize the solutions are located. The singularities have a clear physical interpretation, as the supergravity fields near the poles take the form of the singular "near horizon" limit for a (p, q) five-brane. This gives a natural identification of the poles as remnants of the external five-branes which are used to construct brane web realizations of five-dimensional field theories. The SO(2, 5) symmetry of our solutions implies that we are describing the conformal fixed points where the web collapses to a five-brane intersection with all branes intersecting at a single point. This offers a clear path to identifying the dual SCFTs, as we have discussed for example solutions with three and four poles corresponding to intersections with three or four external branes.
For the large class of explicit solutions presented in this work we have assumed the topology of the Riemann surface to be that of a disk. For higher topologies, we have presented the general set-up for the construction of physically regular supergravity solutions, but we have not shown that such solutions exist. For the case of the annulus, our numerical investigations point in the direction that physically regular solutions without axial monodromy may not exist, though we have no analytical proof. From the analysis and numerical results for the annulus, it is natural to speculate that physically regular solutions will exist only when there is precisely one boundary component. Whether solutions with higher genus exist, and what their physical interpretation in terms of five-brane webs may be, is at this time an open question, to which we hope to return in the future.
Many aspects of the supergravity solutions that we have constructed and the proposed dual superconformal field theories to these solutions remain to be explored as well. We end JHEP05(2017)131 the paper with a discussion of some of the open questions and avenues for future research.
One might wonder whether the presence of the singularities at the poles of the differentials ∂ w A ± , and the fact that the supergravity solution near the singularities approaches the near horizon limit of a five-brane, can invalidate the interpretation of the solution as a dual of a five-dimensional theory. It is a natural next step to study field theory observables holographically in order to address this question.
In the AdS/CFT correspondence the spectrum of operators and correlation functions are the most direct observables which can be obtained by considering the KK-spectrum of linearized perturbations and Witten diagrams, respectively. However, due to the fact that the AdS 6 × S 2 space is warped over the Riemann surface Σ this analysis is considerably more complicated than in more familiar cases such as the AdS 5 × S 5 solution of Type IIB. It would be interesting to investigate the operator spectrum for our warped AdS solutions, as little is known about the operator spectrum of the five dimensional SCFTs.
Other CFT observables such as the entanglement entropy can be calculated more straightforwardly using the Ryu-Takyanagi prescription [29,30]. Such a calculation for a spherical region has been performed for the singular Type IIA AdS 6 supergravity solutions in [31], and we are currently investigating the holographic entanglement entropy for the solutions presented in this paper [32]. Preliminary results suggest that the presence of the poles does not obstruct the holographic computations of either the entanglement entropies or of the free energy directly.
Another interesting prospect for the future is to further generalize our solutions. Sduality acts naturally on the holomorphic functions which parameterize them, but the solutions presented in this paper all exhibit vanishing monodromy under the SU(1, 1) Sduality group. This means that as one moves along a closed loop on the Riemann surface Σ all supergravity fields come back to the same value. It is well known, however, that seven-branes induce nontrivial monodromies on some supergravity fields. For example, a D7-brane induces a shift of the axion along a closed contour enclosing the D7-brane in the two transverse directions. We note that seven-branes can naturally be incorporated into five-brane webs [14] and they are indeed needed for the most general construction of five dimensional field theories and their moduli spaces. Consequently, it is an interesting question whether it is possible to incorporate seven-branes and their monodromies into the framework of our solutions. We plan to address this question in future work.

JHEP05(2017)131
A Absence of poles for ∂ w A ± in the interior of Σ In the main part of the paper we have assumed that all poles of ∂ w A ± are located on the boundary of Σ. In this appendix we shall rule out the presence of poles in ∂ w A ± in the interior of Σ. The arguments are local. For concreteness, we shall present them here for the explicit solutions on the upper half plane, but their locality guarantees that their validity extends to surfaces of arbitrary topology.
The starting point for the construction of the physically regular supergravity solutions is the function λ which, for the upper half plane, is given by, To construct the differentials ∂ w A ± from λ, we had originally assumed that all zeros of λ become zeros of ∂ w A + , while all poles of λ become zeros of ∂ w A − . But this is not the most general possibility available for given λ. Indeed, some of the zeros of λ can actually become poles of ∂ w A − and the accompanying complex conjugate poles of λ can become poles of ∂ w A + . We illustrate this by singling out the pair of zeros and poles of λ, (s 1 ,s 1 ), The locality of the argument will guarantee that it applies equally well when more than one pair of zeros and poles are reversed. We now construct ∂ w A ± accordingly, Note that we have L = N poles to ensure regularity at infinity. We automatically have κ 2 > 0 in the interior of Σ, since this condition was already guaranteed by the form of λ alone. Recall that the pre-factor of i is present (previously denoted ω 0 ) to enforce the conjugation property ∂wA + (w) = −∂ w A − , which is essential for the piece-wise vanishing of G along the real line. Decomposing ∂ w A ± into partial fractions, we have, and integrating,

JHEP05(2017)131
Next, G is computed by first evaluating its derivative, The dominant behavior near the pole w = s 1 in the upper half plane is carried by A − which diverges logarithmically there, while A + is regular. Hence near w = s 1 we have, This implies that the leading behavior of G is given by Hence, near w = s 1 , the function G is negative and divergent. Hence G fails to be positive everywhere in the interior of Σ and the solution fails to be regular.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.