BMN Vacua, Superstars and Non-Abelian T-duality

Acting with non-Abelian T-duality on the $S^3$ inside the $AdS_5$ subspace of $AdS_5\times S^5$ with $N$ units of flux, we generate a new half-BPS solution with $SU(2|4)$ symmetry that belongs to the Lin-Lunin-Maldacena class of geometries. The analysis of the asymptotics, quantised charges and probe branes in this geometry suggests an interpretation as the gravity dual to the Berenstein-Maldacena-Nastase Plane Wave Matrix Model, in a particular vacuum associated to a partition of $N$, in which the multiplicity of each $SU(2)$ irreducible representation is equal to its dimension. This vacuum is interpreted in M-theory in terms of giant gravitons backreacting in the maximally supersymmetric pp-wave geometry. Consistently with this, we show that the non-Abelian T-dual solution exactly agrees with the Penrose limit of the superstar solution in $AdS_7\times S^4$. This suggests an interesting global completion of the non-Abelian T-dual solution in terms of an M5-brane geometry.


Introduction
Among the different dualities in Quantum Field Theory [1] and in String Theory [2], non-Abelian T-duality is probably among the least understood ones. The idea of extending Buscher's procedure [3,4] to the case of non-Abelian isometries, proposed by de la Ossa and Quevedo in [5], motivated thorough studies. Some stringy features (mostly the sigmamodel aspects of the transformation) were spelled out in different publications during the mid-nineties-see for example [6,7]. These papers made the point that global aspects of the duality in [5] are ill-understood. Indeed, for world-sheets of genus bigger than zero, it is presently unclear how to make the procedure consistent. It was then suggested that the non-Abelian generalisation of the Buscher procedure and the transformations in [5] were actually not a duality, but a transformation between different world-sheet conformal field theories [7]. All these discussions concerned only the NS-NS sector of the string theory. These studies were followed by a period of low-activity in the subject.
It was the 2010 paper of Sfetsos and Thompson [8], that gave a prescription to consistently transform the Ramond-Ramond fields, which reignited the interest in the field. Indeed, the synergy between non-Abelian T-duality and the AdS/CFT correspondence [9], lead to various lines of research: • The use of non-Abelian T-duality as a solution generating technique resulted in new backgrounds that in some cases avoid previously studied classifications [10][11][12][13][14][15][16][17][18].
• Defining a string sigma-model in these new backgrounds new relations with integrable models have been worked out [25].
Some of the papers quoted above combine more than one of these lines of research.
Since it happens that some of the backgrounds generated by non-Abelian T-duality are singular, the study of the properties of the field theory dual to those singular spaces, is by definition incomplete. Indeed, in such cases one cannot use the string background to learn about non-perturbative aspects of the dual field theory. This is the case for some paradigmatic examples, like the non-Abelian T-dual of AdS 5 × S 5 -see for example [8,15,17]. For such examples (that characteristically have large isometry group and number of SUSYs), the idea proposed in [22] is to use the field theory (that in cases with large global symmetry groups is sharply determined) to amend or complete the background, curing its singular behaviour. In other words, the quantum field theory dictates the geometry. This idea was made explicit in the papers [22] and [23], where an N = 2 super-conformal field theory in four dimensions and an N = 4 super-conformal field theory in three dimensions, respectively, dictate the way to resolve the singular behaviour of the backgrounds obtained using non-Abelian T-duality.
In this paper, a particularly interesting realisation of these ideas is studied. The general outcome of the papers [22] and [23], that backgrounds obtained via non-Abelian T-duality are a small patch of a more general geometry, will work-out explicitly in our example. The precise way in which such patch will be selected is via a Penrose limit. Indeed, as we will see, the new non-Abelian T-dual background studied in this paper is the Penrose (pp-wave) limit of the so-called Superstar solution [26]. This is the first explicit realisation of a loose but intuitive conjecture that can be found scattered in different papers-see for example [8,27].
The precise set-up in which we are able to find this explicit realisation is a R t × SO(3) × SO (6) symmetric solution that we construct acting with non-Abelian T-duality on the AdS 5 subspace of AdS 5 × S 5 . This solution is dual to a particular vacuum of the large class of 1/2 BPS vacua of the Berenstein, Maldacena, Nastase (BMN) Plane Wave Matrix Model [28]. These vacua are interpreted in the dual supergravity solution as D0-branes polarised into D2-branes at weak coupling, or into NS five branes at strong coupling, very similarly to the N = 1 * Super-Yang-Mills vacua structure [29]. Using that the BMN matrix model is a massive deformation of the BFSS matrix model [30], Lin applied the method in [29] to derive a set of equations that the gravity dual solutions to BMN vacua should satisfy [31]. The paper of Lin and Maldacena [32] finally identified the precise electrostatic problem that must be solved to construct explicitly these solutions within the formalism of [33].
Let us be more detailed about the contents and plan of this paper: • We start our paper by going through the main results of some of the previously mentioned papers, of relevance for our work. These are summarised in Section 2.
• In Section 3, we construct the new solution with SU (2|4) isometry supergroup (this is R t ×SO(3)×SO(6) and 16 SUSYs). We study its asymptotics, show that it satisfies the set of equations posed by Lin in [31] and identify the precise way in which it enters in the classification of [32]. We calculate conserved charges and interesting probe branes that characterise our new background and allow to interpret it as dual to a particular BMN vacuum.
• In Section 4, we present a (0 + 1)-dimensional quantum field theory that captures the dynamics of the non-Abelian T-dual solution and dictates a completion of the background. This allows to interpret our solution as a patch of a more general geometry dual to a relevant deformation of the BMN Plane-Wave Matrix Model [28].
Finally, some conclusions and possible future research directions are collected in Section 7, which closes our paper.

The BMN Matrix Model and its gravity dual
In this section, we summarise different results pertaining to the Plane Wave Matrix Model (PWMM) introduced by Berenstein, Maldacena and Nastase in [28]. These results have been derived in a large collection of papers, of which we single out [31]- [41]. The aim of this section is to organise them in a way useful for our treatment in Section 3. We start by considering N = 4 SuperYang-Mills (SYM) defined on R t × S 3 . The field theory has global bosonic symmetries SO(2, 4) × SO(6) and preserves 32 SUSYs. Since the CFT is defined on a three-sphere, it is natural to distinguish inside SO(2, 4) a subgroup SO(4) ∼ SU (2) ×ŜU (2). This field theory has a single vacuum state and is dual to the Type IIB string theory on AdS 5 × S 5 , with the AdS 5 written in global coordinates, making explicit the SO(4) isometry of the gravity solution.
We consider modding out this field theory by subgroups ofŜU (2). By this, we mean erasing or truncating-away from the Lagrangian of N = 4 SYM on R t × S 3 all the fields that transform (are not singlets) under the given subgroup. The subgroups we can choose are Z k , U (1) and the wholeŜU (2). The field theories that result after the modding out are N = 4 SYM on R t × S 3 /Z k , N = 8 SYM on R t × S 2 and the BMN Plane Wave Matrix Model, respectively. All these quantum field theories have global symmetries SO(6) × SO(3) and preserve 16 SUSYs. The three theories have a very interesting vacuum structure with a mass gap.
In the following and for the purposes of this paper, we will focus on the Plane Wave Matrix Model.

The BMN Plane Wave Matrix Model
As discussed above, this quantum mechanical model is obtained by modding out byŜU (2) the CFT with 32 SUSYs defined on R t × S 3 . After a suitable rescaling of fields, the Hamiltonian reads [28], Here the labels take values A = 1, ...., 9; i = 1, 2, 3 and a = 4, 5, 6, 7, 8, 9. The scalar fields X A , the conjugate momenta P A and the omitted 16-component fermions are hermitian N ×N matrices. The parameters characterising the PWMM are g, that is continuous, and N , that is discrete. The model has a gauge symmetry that requires observable operators M , to transform as M → U −1 M U . The rotation among the six X a fields realises the SO(6) global symmetry and the rotation among the three X i does the same for the SO(3) symmetry. The model preserves 16 SUSYs. All these global symmetries are compactly expressed by the supergroup SU (2|4). Notice that the Hamiltonian in eq.(2.1) is a massive deformation of the dynamics of D0-branes. Short wave-length excitations are then described by the BFSS matrix model [30]. The classical vacua of the PWMM are found by choosing X a = 0. The remaining F-term equation is, Its solution is 3gX i = J i , being J i an N -dimensional representation of SU (2). In general, the vacuum state of the PWMM is then in one-to-one correspondence with partitions of N . We can have n k copies of the k-dimensional representation of SU (2)-of dimension k = (2j k + 1), The vacuum then breaks spontaneously the gauge symmetry according to, It is worth stressing that all these vacua are uplifted to the same vacuum in N = 4 SYM, since they give rise to zero energy configurations that differ from the ordinary vacuum by large gauge transformations. It is interesting that the BMN action can also be obtained by quantising the supermembrane action or the system of N D0-branes [34], in the maximally supersymmetric 11d pp-wave background. In the last description the F-term in the vacuum equations arises as a Myers dielectric term in the action of the N D0-branes. The interpretation of the non-trivial vacua is then in terms of N D0-branes expanded into D2-branes with the topology of fuzzy spheres.
This vacua structure is very similar to that of the N = 1 * SYM theory studied by Polchinski and Strassler [29], which also contains fuzzy sphere vacua, associated in this case to D3-branes polarising into D5. In the same vein, the D0 branes represented by eq.(2.1) can also polarise into NS5-branes at strong coupling [35].
Let us now discuss the Type IIA backgrounds that are dual to the different vacua of the BMN matrix model.

The Type IIA String Duals
In the interest of space, we will not summarise here the results of the remarkable paper by Lin-Lunin-Maldacena [33] (LLM), that started the classification of backgrounds preserving 16 SUSYs with a given set of isometries. The formalism of LLM has applications to a wide class of examples, being the PWMM one of them.
It is clear that backgrounds preserving SO(6) × SO(3) × R t in M-theory can only depend on three coordinates (x 1 , x 2 , y). All the warp factors in the metric and the tensor-structure of the four-form flux respecting the isometries and SUSY preservation can be written-as shown by LLM [33], in terms of a function D(x 1 , x 2 , y) satisfying a Toda equation. Further imposing invariance under x 1 translations, we can reduce to Type IIA and after an implicit change of variables (x 2 , y) → (σ, η), we are left with a Laplace problem in three dimensions with an x 1 -axial symmetry for a potential function V (σ, η)-see [33].
Reducing in the x 1 -direction, the Type IIA background and fluxes are completely determined in terms of V (σ, η) and its derivativesV = σ∂ σ V , The most general Type IIA configuration with these isometries and preserving 16 SUSYs reads [33], with the dilaton, NS-NS and RR potentials given by, (2.5) As we mentioned, the function V solves a Laplace equation (with axial symmetry) whose solutions are subject to precise boundary conditions. These were discussed in the paper by Lin and Maldacena [32]. We can write the potential as a superposition of one part associated with the background electric field and one associated with the particular electrostatic configuration imposed by the boundary conditions, V (σ, η) = V back (σ, η)+φ(σ, η).
For the particular case of the PWMM system, the boundary conditions can be summarised as follows: • The background electric field must grow at infinity. Imposing that the background asymptotes to the 11d maximally supersymmetric plane wave geometry in eleven dimensions one finds that • There is an infinite conducting plane at η = 0.
• There is a set of conducting disks at positions η = η i = π 2 N (i) 5 , set by the NS5-branes in the background, with radii R i and charge Q i = π 2 8 N (i) 2 , associated to the D2-brane charges. These NS5 and D2 brane charges are obtained tensoring the S 2 and the S 5 of the background by different lines in the (σ, η) plane. Imposing that the electric field on the rim of each disk is finite forces a relation between R i and Q i .
As we pointed out before, the PWMM describes D0-branes plus some relevant deformations. Therefore, the corresponding solutions must asymptote to the metric of D0-branes. This can be seen studying the leading large distance asymptotic form of the solution, which is obtained with a potential [32], which includes a 'dipole' sub-leading correction, produced by the disks and their images with respect to the conducting plane at η = 0. Plugging eq.(2.8) into the metric given by eq.(2.4), the asymptotics at (η, σ) → ∞ becomes which, changing coordinates r 2 = 16(η 2 +σ 2 ) and considering the limit of very large distances, r → ∞, becomes the standard D0-brane solution (2.10) The number of D0-branes is thus proportional to P , the dipole moment produced by the disks, 2 for a distribution of i disks. Comparing this with eq.(2.2) we see that N In summary [32], in order to find the gravity duals of the different vacua of the PWMM, we need to search for a solution to the Laplace equation in the (σ, η)-space with a set of finite size R i conducting plates positioned at distances η i = π 2 N (i) 5 and one infinite conducting disk at η = 0. The charges of the finite plates satisfy Q i = π 2 8 N (i) 2 . These charges are carefully chosen to produce a smooth geometry. Such solution of the Laplace problem-when replaced in eqs.(2.4)-(2.5), results in a Type IIA background that is conjectured to be dual [32] to the vacuum of the PWMM described by N Notice that in principle, we could have tunneling between different vacua 1 , but in the 't Hooft limit with large coupling this is suppressed (as e − N λ t , as it is mediated by instantons) and it is justified to discuss solutions referring to a particular vacuum, as above.
We will now sketch different methods developed in the literature to solve the electrostatic problem described above.

Solutions to the Laplace Problem
Different papers dealt with the challenging problem of solving the Laplace equation in the (σ, η) plane with the boundary conditions imposed by a configuration of one infinite and many finite size conducting disks. As we explained above, any solution of such problem captures the strong coupling dynamics of the PWMM in a particular vacuum.
Considering the situation depicted in Figure 2, in which the vacuum is described by N 2 copies of representations of dimension N 5 , the dual electrostatic distribution of conductors consists on a plate placed at η = η 0 = π 2 N 5 , that extends in the σ-direction up to a distance σ = R (the potential on such plate is taken to be V = ∆), plus an infinitely long conducting plate at η = 0, where by convention, we choose to set V = 0. Writing the electrostatic potential function as V = V back + φ(σ, η), with V back given by (2.7) and φ(σ, η) satisfying the Laplace equation, it must satisfy continuity on the disk A solution for the function φ(σ, η) satisfying these conditions can be obtained by separating variables. This leads to solutions in terms of exponentials and a Bessel function, basically a Fourier-Bessel transform for each of the two regions depicted in Figure 2 [36] and [37], for an efficient way of imposing the conditions above and after lengthy algebra, a pair of coupled integral equations and relations among the functions A l (u), B l (u) can be obtained.
The resolution of these integral equations is quite challenging. For the case at hand, they were solved implicitly in [36] and [37] in terms of nested integrals. The problem was generalised to the case of many disks and solved-again in an implicit fashion and in terms of nested integrals-by the authors of [38].
The main message is that finding analytic solutions of the precise Laplace problem dictated by the Lin-Maldacena formalism, is quite involved. In the following, we will study an approximated version of the Lin-Maldacena problem for which it is easier to find solutions, and we will discuss the physical implications of such approximation.
This approximation was originally devised by the authors of the paper [39] (written before the work of Lin and Maldacena!). Whilst they analise exact solutions to the Laplace equation (2.6), the 'approximation' takes place in the boundary conditions. Instead of a discrete distribution of plates with finely-tuned charges Q i , placed at precise positions η i and with radii R i , the solutions in [39] describe the situation in which on the η-axis, we have a continuous distribution Q(η) of very short charged conductors. As explained in [39], for the S 2 and the S 5 to shrink smoothly it is necessary that Q(η) is odd and positive, but this is not sufficient to avoid a null singularity, that appears in the geometry at the point where the radial component of the electric field vanishes, and is rather similar to that of extremal Dp-branes. Subsequently, the authors of [40,41], confirmed this and gave an interpretation in terms of a coarse-grained (or Thermodynamical) description around a 'typical state' (or more probable arrangement of conducting plates). Interestingly, the presence of the singularity in these backgrounds is a consequence of this 'average' or coarse-graining of the Lin-Maldacena distribution of conducting plates, as it happens for microstates of black holes, discussed in [43].
Following [39], the electrostatic problem is then reduced to finding a solution of the Laplace equation, subject to the condition that an infinite conducting plane lays at η = 0 and a continuous charge distribution Q(η) is located on the η-axis. Using for instance the method of images one gets that the solution for such electrostatic problem is -see [39][40][41], Given that in the solution ∞ 0 Q(η)dη counts the number of D2-branes and ∞ 0 ηQ(η)dη the number of D0-branes, the proposal of [39] is to make a correspondence between this function and the partitions of N -that we called n k in eq.(2.2), such that In other words, the function Q(η) is giving a 'density' of SU (2) representations of dimension η. Somewhat as it happens for the case of four dimensional N = 2 SCFTs, the η-direction in the dual gravity problem becomes the 'theory space' direction, see for example [22]. This interpretation was confirmed by other arguments in [40,41]. Let us close this section with a brief summary. The method in [39] allows to find a solution to the Laplace equation (2.6). This solution does not satisfy the boundary conditions that a Lin-Maldacena problem demands, hence this leads to singularities in the generated background. Indeed, instead of a discrete set of plates, with precisely defined positions and charge-to-radius relations, a continuous distribution of charge is considered. From the perspective of the field theory, one loses the precise dimensions of each entry in the vacuum state matrices, to just keep a 'density' of SU (2) representations represented by Q(η). This is the method we will use to understand and make sense of the backgrounds generated by non-Abelian T-duality. We now move on to discussing this new background.

A new non-Abelian T-dual to AdS 5 × S 5
In this section, we study the background generated by applying a non-Abelian T-duality transformation to the AdS subspace of AdS 5 ×S 5 . This will be a new solution of the Type IIA equations of motion that can be understood in the context of dual backgrounds to PWMM vacua. This connection will provide us with a precise way of understanding and curing problematic features of our new solution. Let us start by considering the background AdS 5 × S 5 , with the AdS 5 -subspace written in global coordinates, with time and radial-coordinates (t, r). Let us use the left invariant forms of SU (2) to write the three-sphere as dΩ 2 In terms of these, the background reads, The parameter L is the radius of AdS 5 (also the size of the five sphere) and obeys the condition , that quantises the size of the space. Let us perform a non-Abelian T-duality transformation using one SU (2) inside the SO(4) isometry of the previous three sphere. The resulting configuration is a new Type IIA background that reads, This solution is singular at the point r = 0, due to the presence of NS five branes. Indeed, close to r ∼ 0 the metric and dilaton become with ν = r 2 , which is the metric of smeared NS five branes, see [23]. In Section 4, we will make sense of this singular behaviour using the field theory information reviewed in Section 2.
Note that the previous singularity could also be avoided in the presence of a black hole event horizon. Indeed, we could have started with the black hole in global AdS 5 . In this case we have the same RR flux plus a metric that reads, where the function h(r) is the blackening factor and c 1 , c 2 are constants. The event horizon is localised at the point where h(r h ) = 0. One can see that the blackening factors h(r) and h −1 (r) appear in the SU (2) NATD solution associated to this background in front of the g tt and g rr metric components, while the RR and NS-NS fields remain the same 2 . Further study of this solution could be useful in order to understand finite temperature effects. In this paper we will focus however on the zero temperature background described by eq. (3.2) without horizons. Let us then consider the solution in eq.(3.2). This is a half-BPS background (it preserves 16 supercharges) of Type IIA supergravity with R t × SO(3) × SO (6) isometry. Therefore it should belong to the general class of solutions described by eqs.(2.4)-(2.6). Indeed, one can check that with the change of coordinates, we can derive the solution in eq.(3.2) from a potential 3 , where V back is the background potential associated to the PWMM. Indeed, in these coordinates, the background metric reads, which can be checked to be that of eq.(2.4) for the potential V NATD above. Similar expressions can be calculated for the fluxes and dilaton in (σ, η)-coordinates, using eq.(3.2), the change in eq.(3.5) and comparing them to those in eq.(2.5) for V NATD .
In light of what we discussed in Section 2, the singularity of the non-Abelian T-dual solution at σ = 1 should then be a consequence of the fact that our background does not satisfy the boundary conditions itemised below eq.(2.6). We will come back to this issue in the following subsections.
Let us now proceed with a detailed study of our new Type IIA solution.

Asymptotics
As we have pointed out, the PWMM describes D0 branes plus some relevant deformations. Therefore, the corresponding solutions must asymptote to the metric of D0 branes. Let us analyse whether this is the case for the non-Abelian T-dual solution.
The logarithmic term in the potential associated to the non-Abelian T-dual solution, given by eq. (3.6), does not vanish at infinity, unlike the sub-leading dipolar deformation of eq.(2.8). This suggests that the non-Abelian T-dual solution will not have D0-brane asymptotics. This is indeed the case. Consider the limit σ → ∞, η = fixed, in the non-Abelian T-dual background metric. The metric in eq.(3.7) becomes, In order to understand what brane configuration gives this asymptotics, consider D0 branes defined by the metric in eq. (2.10). If we smear the D0 branes on the (η, Ω 2 ) directions, which implies a continuous superposition of these D0 branes [44,45], we find where h s,l denote the warp factor of the smeared and localised D0 brane solutions respectively. Then, the smeared D0 brane metric behaves like, which is precisely the asymptotic form of the NATD solution in eq. (3.8). Therefore, the largeσ asymptotics of the NATD solution is that of smeared D0 branes on R 3 . Had we considered the situation in which the D0 branes had been smeared on the finite range η ∈ [0, ∆ 7/3 ], we would have found a smeared warp factor, h s (σ) ∼ 4π , which has the same large distance asymptotics as the D0 brane metric in eq.(2.10). The conclusion is that the non-Abelian T-dual solution of eq. (3.8) can be interpreted as an infinite superposition of D0 branes. Its different asymptotic behaviour can then be associated to the presence of an irrelevant operator in the matrix model, whose effect is to change the dimensionality of the field theory. Given that an irrelevant operator is never welcomed in a field theory, we will use the ideas summarised in Section 2 to 'cure' this problem.
Similarly to the analysis above, we can consider the limit in which σ is kept fixed, but the coordinate η → ∞. We can see that in this case the NATD background asymptotes to the Abelian T-dual (in the U (1) ψ -direction) of the solution in eq.(3.1). This connection between non-Abelian T-duality and its Abelian counterpart was observed in [17] and later refined in [22]. Here again, the metric represents a dual to a quantum field theory in dimension higher than (0 + 1).

Realisation as a smeared Lin solution
In the previous subsection, we have shown that our new Type IIA background obtained by applying non-Abelian T-duality is a singular solution, and that it can be understood as describing D0 branes that have been smeared along R 3 .
In this section we will make contact with a nice work by Hai Lin [31], where, applying the method used by Polchinski and Strassler [29] to construct the Type IIB duals to the vacua of the N = 1 * field theory, he considered backgrounds of Type IIA supergravity obtained by perturbing the near horizon geometry of N D0-branes by transverse RR 6-form and NS-NS 3-form fluxes. These fluxes are the Type IIA duals to the deformation of the BFSS matrix model by a mass and a dielectric term, thus giving rise to the Type IIA realisation of the BMN matrix model [28]. Lin showed that in the presence of these fluxes the D0-branes polarise to concentric D2-branes with radii proportional to their D0-charge. When lifted to eleven dimensions, these solutions describe giant gravitons that expand into M2-branes on the eleven-dimensional pp-wave background, thus making contact with the ideas proposed by [39], subsequently refined by the authors of [40,41]. We will show that our NATD solution solves the equations proposed by [31], provided that the branes are smeared. This, together with the material discussed in Section 3.1, reinforces the interpretation of the NATD solution in terms of smeared D0 branes.
Let us start by briefly summarising the main results of the work [31]. The Type IIA solutions constructed in [31] are perturbations in the large-r limit around the near horizon geometry of N D0-branes, by magnetic RR 6-form and NS-NS 3-form fluxes. The background metric, dilaton and F 2 flux characterising D0-branes, are perturbed up-to-second order in the fluctuations of the fields H 3 , F 6 . One just needs to solve the equations of motion for the fluxes, which take the simple form [31], where in these expressions * 9 is the Hodge dual in the transverse nine-dimensional space with respect to a flat 9d metric.
Let us now come back to the σ → ∞ asymptotic behaviour of the non-Abelian T-dual solution. As we showed below eq.(3.8), the behaviour of the metric is that of smeared D0branes on R 3 . The same smearing happens for the dilaton and A 1 RR field (they are also written in terms of the smeared warp factor h s ) that read, Thus, as a result of the smearing, the fields H 3 , F 6 and their * 9 -duals are given to leading order in σ as, Using these fluxes, one can show that the first two lines of eq.(3.10) are straightforwardly satisfied. Considering that the warp factor is a function of σ only, h(σ), one can show that the fourth equation in (3.10) is satisfied, while the third line of eq.(3.10) imposes that h −1 ∼ σ 4 . This is, precisely, the same result as the smearing of the D0 branes on R 3 .
In summary, we have shown that the non-Abelian T-dual configuration provides a full solution (beyond the perturbative expansion devised by Lin) for a background of smeared D0 branes deformed by asymptotically sub-leading F 6 and H 3 fluxes. There remain the unpleasant features of the smearing along a non-compact space (which obscures the field theoretical interpretation of the non-Abelian T-dual solution), and the fact that the background is singular at σ = 1 (which can be physically interpreted following [39], [40] and [41]). We will discuss this in Section 4. Before that and to get more clues for the field theoretical understanding of the non-Abelian T-dual solution, we will calculate its conserved charges.

Quantised charges
The study of the quantised Page charges was instrumental in developing the field theory dual of the only backgrounds obtained via non-Abelian T-duality to which a concrete CFT dual has been proposed [22,23]. In these examples, the quiver associated with the field theory was encoded in a brane set-up compatible with the Page charges of the solution. In this section, we study the conserved charges of our Type IIA background.
The presence of F 2 and F 4 fluxes suggests that the solution should be associated with a combination of D0 and D2-branes. The Hodge duals of the fluxes in eq.(3.2) are given by, The relevant "Page flux"F 8 = F 8 − B 2 ∧ F 6 reads in turn, In analogy with other cases studied in the literature-see for example [13,15,19,22,23]the B 2 -field should be defined up-to large gauge transformations, associated to the non-trivial 2-cycle that sits at the singularity at r = 0, where the cone spanned by (r, Ω 2 ) becomes singular. Indeed, B 2 should be defined as, 13) and this naturally divides the ρ-coordinate in intervals of size π, with a large gauge transformation as that in eq.(3.13), performed every time the interval ρ ∈ [kπ, (k + 1)π] is crossed. This affectsF 8 in eq.(3.12) as, (3.14) Using that 2κ 2 10 T Dp = (2π) 7−p g s α ′(7−p)/2 and setting g s = 1, the D2 and D0 charges associated to F 6 andF 8 read, for ρ ∈ [kπ, (k + 1)π], We have also NS five brane charge, The NS five branes are located at ρ = π, 2π, . . . kπ, and r = 0 (note that we integrate over the S 2 at r = 0). There is thus no brane set-up in which color branes are stretched between pairs of NS five branes. This represents an essential difference with respect to previous backgrounds constructed through non-Abelian T-duality, such as the examples in [22,23]. Our discussion in Section 3.2 suggests that the correct way to interpret the D0 and D2 brane charges is instead in terms of D0-branes expanded onto spherical D2-branes. Indeed, the second term in the integral in equation (3.16) originates from and is thus equal to, 19) with N 2 the D2-brane charge given by eq.(3.15). Thus, this contribution to the D0-brane charge comes from the D0-branes that are dissolved in the D2-branes. Conversely, we can think of the D2-branes as originating from D0-branes expanded onto fuzzy 2-spheres through Myers' dielectric effect [46]. Indeed, D2-branes in the non-Abelian T-dual background only carry D0-brane dipole charge, induced by the large gauge transformations of the B 2 field (see the next section). We have used the notation N  is just standard (monopole) D0-brane charge.
The picture that arises is that of N 2 spherical D2-branes with k D0-brane charge dissolved. We will show in the next section that ρ can in fact be interpreted as the radius of smeared D0-branes expanded into the S 2 spanned by the coordinates (χ, ξ) in eq.(3.2). The brane set-up is thus a concentric distribution of expanded, smeared, D0-branes on 2-spheres with radii ρ = kπ. This is depicted in Figure 3. On top of this there are as well standard point-like D0-branes at ρ = kπ, with charges given by eq.(3.20), that do not have a direct interpretation in terms of BMN vacua (see below).
Taking a [0, pπ] interval, the total D0-brane charge distributed among the concentric D2-branes at ρ = π, 2π, . . . , pπ is then given by, This can be seen as a partition of the total number of D0-branes of exactly the same form of the partition in eq. (2.2), that defines the different BMN vacua. Here k is the dimension of the SU (2) representation, k = (2j k + 1), and N 2 its number of copies, n k . As in [32], the dimension of the representations is associated to the number of NS five branes, while their multiplicities are associated to the number of D2-branes of the gravity dual. Let us discuss in more detail the situation at the [0, π] interval. In this interval we have one NS five brane located at ρ = π and N 2 = L 4 16α ′2 point-like D2-branes. This configuration should then be associated to the trivial BMN vacuum. It is interesting that this appears to be an explicit realisation of the conjecture in [35] that the trivial vacuum should be associated to a single NS five brane. It would be interesting to understand this connection more precisely. Accordingly, we have to shift all charges by L 4 16α ′2 in order to fulfil the condition that V = 0 at η = 0 and that the D2-brane charge distribution is odd along the η axis. Doing this shift  we have N 2 = L 4 8α ′2 k at each [kπ, (k + 1)π] interval. Through eq. (3.21), the total D0-brane charge in [0, pπ] is then Finally, following [35], the same number of D0-branes given by eq.(3.22) should be able to expand at large 't Hooft coupling onto NS five branes wrapped on 5-spheres. The argument presented in [35] is that we can represent any partition of N by a Young diagram whose column lengths are the elements of the partition. In the D2-brane interpretation such diagram corresponds to a state with one D2-brane for each column, with the number of boxes in the column indicating the D0-brane charge carried by this D2-brane. Alternatively, in the NS five brane interpretation, the rows of the Young diagram correspond to the number of D0branes carried by each NS five brane. In this dual interpretation the number of NS five branes coincides with the dimension of the largest irreducible representation, and the D0-brane charge carried by the mth NS five brane is equal to the number of irreducible representations with a size greater than or equal to m. This partition of N is commonly referred as the dual partition. This is depicted in Figure 4. In our case the dual partition to the partition in eq. (3.21) describes p NS five branes, with the mth of them carrying D0-brane charge given by We can check that indeed N , as given by eq.(3.22), is reproduced from As we will discuss in Section 3.5, it is possible to construct BPS NS five branes wrapped on the 5-sphere of the non-Abelian T-dual geometry, located at r = 0 and at given positions in ρ, that carry D0-brane charge given by eq.(3.23).

Dielectric D2-branes
As we have discussed, Lin showed in [31] that in the presence of transverse RR 6-form and NS-NS 3-form fluxes, D0-branes polarise into concentric D2-branes with radii proportional to the D0-charge that they carry. This description is valid in the limit of large D0-charge. In this section we show that there are similar concentric D2-branes wrapped on the 2-sphere of the non-Abelian T-dual solution that are interpreted as D0-branes smeared on r and concentric 2-spheres of radii kπ. Let us consider a probe D2-brane wrapped on the 2-sphere of the non-Abelian T-dual solution, the radius of which depends on both ρ and r. In the presence of large gauge transformations this brane is dielectric and carries D0-brane charge. This can be shown from the couplings to the RR potentials in the WZ action, Thus, substituting in eq. (3.26) we see that the brane carries k units of D0-brane charge, In the particular case of the NATD solution this coupling reads (3.28) In order to find the equilibrium position of the brane we compute as well its DBI action: A D2 placed in this background has then a potential The equilibrium position of the brane, at which it carries zero energy, is then found minimising in ρ and r. We find that ρ = kπ and that r ∼ ∞. The radius of the brane is thus which corresponds to a metric where we have used that sinh r ∼ cosh r ∼ σ for large r. The asymptotic behaviour given by this equation corresponds to D0-branes smeared on r and an S 2 with radius kπ.

Dielectric NS five branes
The previous configuration in terms of concentric shells of smeared D0-branes is valid at weak coupling [35]. Alternatively, at strong coupling we should find configurations of NS five branes carrying D0-brane charge. In this section we study these configurations using the effective action describing NS five branes in Type IIA constructed in [47]. We refer the reader to this reference for more details on its construction. Consider a probe NS five brane wrapped on the 5-sphere of the non-Abelian T-dual solution 4 , This ansatz does not allow configurations of concentric NS five branes, since all branes have the same radius L. Instead, we will see that the NS five branes will be located at different positions in the ρ direction according to the D0-charge that they carry. In order to describe this brane it is enough to use the truncated DBI action in [47], where G 1 is the invariant field strength associated to the RR 1-form potential, through which we will introduce D0-brane charge in the configuration. Here c 0 is the worldvolume field that originates in the eleventh scalar of the M5-brane. In close analogy with the more usual BI field, we can induce D0-brane charge in the configuration takingċ 0 = 0. The relevant part of the WZ action reads in turn: where B 6 is such that Particularising to the non-Abelian T-dual background, where we have used that A 5 = −2 √ α ′ L 4 ρ 2 dΩ 5 after integrating F 6 in eq.(3.11). As already mentioned, in order to induce D0-brane charge in the NS five brane we switch onċ 0 = 0. Its conserved conjugate momentum will then give the D0-brane charge carried by the NS five brane. Indeed, substituting in the action (of a NS five brane in our conventions) we find, The corresponding Hamiltonian H = Pċ 0 − L reads, (3.41) Minimising in P we find that and This is in turn minimised for r = 0, where the NS five brane becomes in fact tensionless. The output of this analysis is that an NS five brane wrapped on the S 5 can be in equilibrium if it is located at the singularity of the background, r = 0, and at a position in the ρ-direction that satisfies that its D0-charge, P , is quantised, as imposed by eq. (3.42). This condition shows that NS five branes located at ρ = mπ carry D0-brane charge P = L 4 m 2 16α ′2 . In turn, the NS five branes defined by the dual partition of N should be located at ρ = π p(p + 1) − m(m − 1), such that the relation in eq.(3.23) is satisfied.
Below, we will use all the information collected in this section to propose a field theory dual to our non-Abelian T-dual solution of eq.(3.2) Indeed, using the expressions in eq. (2.12), we can derive the charge density that approximates our distribution of conducting plates taking the partition of N suggested by the quantised charges of the non-Abelian T-dual solution. This partition is given by eq.(3.21), as found out in Section 3.3. Doing this one finds that, (4.1) Plugging now this expression for Q(z) in the 'coarse-grained' potential given by eq. (2.11), and using that -see [41], where, as explained in [41], the η log 2L term can be added in order to regularise the potential, and not changing the expressionV of the background, we are able to reproduce the potential associated to the NATD solution, given by eq.(3.2). The result of this analysis is that the NATD solution arises as a 'coarse-grained' approximation of the Lin-Maldacena problem associated to it. This justifies the existence of the singularity at r = 0 (or σ = 1). However, an important aspect of the solution remains that obscures its interpretation as dual to a PWMM vacuum. As we have seen, the large σ asymptotics of the solution is that of D0-branes smeared on R 3 (η, Ω 2 ). Consistently with this, we have shown that it satisfies a smeared version of the fluctuated equations written by Lin in [31]. Our solution thus represents a smeared-deformation of the PWMM.
Following the ideas in [22,23], we would like to propose that the configuration obtained using non-Abelian T-duality-eq.(3.2)-represents a patch of a background with the right D0brane asymptotics of a faithful dual representation of the PWMM. Thus, in this 'completed' background the large σ and large η effects of smearing the D0 branes, as explained in Sections 3.1 and 3.2, are amended. The consequence of this will be to have a (0 + 1) dual field theory, a matrix model.
A simple way to achieve this is to take the charge distribution given by eq. (4.1) but with η living in a finite interval, η ∈ [0, L], as suggested by the regularised potential in eq.(4.2). Ending the density of SU (2) representations in this manner has the effect of recovering the asymptotic behaviour driven by V ∼ ησ 2 − 2 3 η 3 , and thus the large distance asymptotics of D0 branes.
The explicit potential generated by the charge distribution in eq. (4.1) with η ∈ [0, L] is, using eq.(2.11), One can check that this potential satisfies the Laplace equation (2.6), and that for finite L the asymptotics is that of the (ησ 2 − 2 3 η 3 ) term. Taking instead L very large we find, The first term in ∆V does not contribute toV nor to the background in eqs.(2.4)-(2.5). All subleading terms in the expansion satisfy the Laplace equation. Thus, taking L → ∞ we encounter the V NATD , which produces the solution in eq.(3.2), whose large distance behaviour is that of smeared D0 branes. This can be interpreted as the deformation of the dual PWMM with an irrelevant operator, which mirrors the ever-growing density of SU (2) representations, Q(z) = z in that system. This is an effect of focusing on the solution given by V NATD . On the other hand, the full completed solution obtained with the potential in eq.(4.3) replaces the ever-growing density of representations, by that generated with eqs.(2.11) and (4.1) up to a maximum value of η, η = L. This ameliorates the long distances behaviour, by using the background generated by the potential in eq.(4.3). We are working with a matrix model deformed by relevant operators. As anticipated in the papers [39], [40] and [41], this procedure does not resolve the small distance singularity. Indeed, what happens is that the singularity (originally at σ = 1 for the non-Abelian T-dual solution), moves to σ = 0 for the completed background obtained with the potential in eq.(4.3). In our case, this singularity could be screened by a black hole horizon adding temperature to the dual PWMM.
In summary, the non-Abelian T-dual solution should be thought of as a patch of a more general background. This completed solution, encoded by the potential in eq.(4.3), is dual to a vacuum of a (0 + 1) quantum field theory.

Connection with the AdS 7 × S 4 superstar
In the previous section we have found that the geometry in eq.(3.2) can be completed using the ideas of [39][40][41], in such a way that the asymptotics of the metric becomes that of D0-branes. Following this idea, the geometry of eqs.(3.2)-(3.7), encoded by the potential function V NATD of eq.(3.6), was completed by replacing the potential with that found in eq.(4.3).
This completion was motivated by the fact that the 'smeared D0-branes' asymptotics of the NATD solution was in conflict with its interpretation as dual to a vacuum of a (0 + 1) dimensional QFT such as the PWMM. Even if the completion in Section 4 gives rise to the right D0-brane asymptotics, a singularity, this time at σ = 0, still remains, which indicates that some effects of the QFT are not captured by the completed supergravity solution.
In this section, we will propose a different way of globally defining our non-Abelian Tdual geometry by 'growing dimensions' to the (0 + 1)-QFT. More specifically, we will show that we can connect our background with the Penrose limit of the superstar solution in [26], which is dual to a relevant deformation of the 6d CFT living in M5-branes.
We will start by studying in detail the Penrose limit of the AdS 7 × S 4 superstar solution constructed in [26], that describes giant gravitons backreacting on AdS 7 × S 4 . We will show that this geometry exactly agrees with our non-Abelian T-dual solution uplifted to eleven dimensions. A possible interpretation of non-Abelian T-duals as Penrose limits of AdS solutions was suggested in the seminal paper [8] by Sfetsos and Thompson. Our explicit realisation of this idea -the first in the literature-suggests that it may apply as well to other non-Abelian T-dual geometries.

Penrose limit of the half-BPS AdS 7 × S 4 superstar solution
As we have discussed before, by relaxing one of the boundary conditions of the Laplace equation we can still generate physically relevant solutions. An interesting example is the deformation by giant gravitons of AdS 7 × S 4 , known as the superstar solution [26]. Such a deformation breaks half of the supersymmetries of the vacuum solution and corresponds to a horizon free black hole whose singularity is interpreted as a condensate of giant gravitons. Interestingly, it was found that the origin of the singularity in the electrostatic description corresponds to the coarse grained limit of a distribution of conducting disks whose sizes are much bigger than the separation between them, such that far from the disks we can replace the discrete configuration by a conducting continuous surface of shape R(η) [41]. Such a description was shown to reproduce exactly the microscopic interpretation of the solution. Let us now proceed by describing the relationship between this solution and the non-Abelian T-dual solution studied in this work. For concreteness, we will only work with the elevendimensional metric, though the configuration is completed by a F 4 -field.
The metric of the eleven dimensional AdS 7 × S 4 superstar solution is described by [26], and dΩ 2 5 and dΩ 2 2 are the line elements of the five and two spheres respectively. To proceed consider the following scalings, After these scalings, the metric in eq. (5.1) reads, where,∆ = 1 +q r 4 sin 2 θ,H = 1 +q r 4 ,f = 1 +r 2H .
We study the behaviour of the metric in the neighbourhood ofr = 0, θ = 0. We take the limit L → ∞ and let,r Notice that we have rescaled the charge parameter in order to obtain a non trivial q-dependent solution in the limit. The leading order behaviour of the functions in eq. (5.5) in the limit L → ∞ is,∆ The metric becomes, (5.8) Redefining r 2 5 = 4(σ 2 − 1), r 2 = 2η and setting q = 4, the metric in eq. (5.8) can be written as, In summary, by focussing (or Taylor expanding) the solution in eq. (5.4) near the singularity atr = 0 along θ = 0, with the appropriate scalings in eq. (5.6), we find the eleven dimensional uplift of the non-Abelian T-dual solution in eq. (3.7). Since Taylor expanding a solution does not (in general) produce a new solution of the equations of motion it is worth to investigate the nature of the above limit. In the following we will show that the procedure discussed above corresponds to a Penrose limit. Consider a null geodesic γ parametrised by the coordinates (t, x 11 ) that sits at the point θ = 0,r = 0 of the metric in eq.(5.4). We approach this point by takingr = r 0 + r 5 L and θ = θ 0 + r 2 2L . We will be interested in the case of r 0 and θ 0 vanishing, and L → ∞. Before formally taking the limit of large L, we consider the induced metric for this geodesic. It reads, -24 -Notice that in the limit r 0 → 0, the induced metric in eq. (5.10) and the geodesic solution in eq.(5.12) are ill defined as H 0 → ∞. This is expected since we are considering a geodesic in the neighbourhood of a singularity. However, there is a nice conspiracy in the limit L → ∞ such that the induced metric in eq. (5.10) is well defined 5 . Indeed by taking the limit of large L and choosing null coordinates such that 13) in the limit in which r 0 → 0 and θ 0 → 0 with L → ∞ and also scalingq = q L 2 -the content of eq.(5.6), the metric in eq. (5.4) takes the form of a pp-wave in Rosen coordinates. Namely where∆,f are defined in eq.(5.7) We therefore see that the non-Abelian T-dual background corresponds to the Penrose limit of the half-BPS AdS 7 × S 4 superstar. Notice that the nontrivial rescaling in the charge parameter in eq. (5.6), is necessary in order to keep the effect of the parameterq that controls a relevant deformation on the (0, 2) CFT. Had we rescaled theq parameter in a different way such that the limit L → ∞ remains finite we would have obtained the maximally supersymmetric pp-wave in eleven dimensions. Indeed, settingq = 0 in the metric of eq. (5.14), we obtain the maximally supersymmetric pp-wave, as expected.
To make our argument stronger, notice that the null geodesic chosen in eq. (5.6) is consistent (up-to typos) with the analysis performed in [48], where the Penrose limit of the superstar solution in eq.(5.1) was studied in the context of the LLM formalism. In this approach the Penrose limit is taken by blowing up the neighbourhood around the edge of a droplet describing the superstar solution. The upshot of this procedure is the potential function in eq. (3.6) which corresponds to the non-Abelian T-dual solution.
In summary, the Penrose limit of the AdS 7 × S 4 superstar coincides with the uplift to eleven dimensional supergravity of the geometry in eq.(3.7). Using Holography, this tells us that the CFT dual to our background in eq.(3.2) could be thought of as a particular sector of the (0, 2) SCFT with relevant deformations.

Connection with the Sfetsos-Thompson solution
In this section we show that our non-Abelian T-dual solution (obtained upon dualisation on an S 3 ⊂ AdS 5 ) can also be obtained with a Wick rotation on the Sfetsos-Thompson solution [8] (obtained upon dualisation on an S 3 ⊂ S 5 ). This observation, combined with the previous interpretation of our non-Abelian T-dual solution as a Penrose limit, suggests that it should be possible to interpret the Sfetsos-Thompson solution as the result of taking a Penrose limit on a globally well-defined N = 2 geometry.
Indeed, the Sfetsos-Thompson solution is an example of an AdS 5 background belonging to the Gaiotto-Maldacena class of geometries [49]. Inspired by this, an interpretation to its dual CFT could be given in terms of an infinitely long linear quiver [22]. Starting with the AdS 5 × S 5 solution, ds 2 =L 2 ds 2 AdS 5 + L 2 4dα 2 + 4 sin 2 αdβ 2 + cos 2 α(ω 2 1 + ω 2 2 + ω 2 3 ) , where 4L 4 = πg s N α ′2 and ω i are the SU (2) left-invariant Maurer-Cartan forms, the Sfetsos-Thompson solution was obtained acting with non-Abelian T-duality on one of the two SU (2) isometries of the internal space. The resulting background is given by [8], A potential function satisfying the Laplace equation (2.6) can be written, in this case with a charge density given by λ(η) = σ∂ σ V (σ, η)| σ=0 . After changing variables as, ρ = 2L 2 α ′ η, σ = sin α, (6.3) the potential is Let us now Wick rotate the solution in eq. (6.2), where ds 2 AdS 5 = − cosh 2r dt 2 + dr 2 + sinh 2r dΩ 2 3 , (6.5) in global coordinates, by performing the following analytical continuation in the coordinates: α → −ir − π 2 , t → φ andr → iθ. The effect of the above transformations is to replace , cos α → −i sinh r, ds 2 AdS 5 → −dΩ 2 5 , (6.6) in eq. (6.2), where dΩ 2 5 is given by, dΩ 2 5 = cos 2 θdφ 2 + dθ 2 + sin 2 θdΩ 2 3 . (6.7) In order to get the metric with the correct signature we further impose L → 2iL. In this way one arrives at the non-Abelian T-dual solution in eq. (3.7). Notice that the relation between the solutions in eqs. (3.7) and (6.2) is an explicit ten dimensional realisation of the analytical continuation procedure discussed in [33]. We would like to stress that even if both solutions in eqs. (3.7) and (6.2) solve the same Laplace equation (2.6), they differ in the boundary conditions they satisfy. Notice as well that it is necessary to change the sign of the potential, V ST → −V N AT D , to ensure that the metric gets the correct signature.

Conclusions and future directions
For the benefit of the reader, let us start with a brief summary of this work.
In this paper, we constructed and studied a new background, solution to the equations of motion of Type IIA supergravity. It was obtained by the application of non-Abelian Tduality on an SU (2) inside the SO(2, 4) isometry of AdS 5 × S 5 . In Section 3, we worked-out the asymptotics, conserved charges and interesting probes of the space-time. A proposal for a quantum field theory dual in terms of a given BMN vacuum was then given in Section 4. The precise way in which we were able to define a consistent field theory dual associated to our solution provided a very explicit realisation of the idea that solutions generated by non-Abelian T-duality correspond to a patch of a more generic background. The quantum field theory proposed in Section 4 to inform the completion of our Type IIA background is (0 + 1)-dimensional, a matrix model obtained by a relevant deformation of the BMN Plane Wave Matrix Model. As we emphasised the background is still singular at short distances (σ = 0), being this a feature of the procedure of [39], that we adopted. Meanwhile, in Section 5 a different completion of the non-Abelian T-dual background and its associated field theory, was proposed in terms of a 6d (0,2) SCFT, also deformed by relevant operators. The material in Section 5 explicitly realises (by the first time) a conjecture that appeared in various previous instances in the literature. Namely, our new background is the Penrose limit of a well-known solution, the Superstar in AdS 7 × S 4 . Finally, in Section 6 we connected our Type IIA background with that constructed in [8].
The relevance of the material we have presented immediately suggests numerous new avenues of research. Let us mention some of them here.
It would be interesting to extend the realisation of the non-Abelian T-dual solution as the Penrose limit of the AdS 7 ×S 4 solution to other examples. It would also be nice to study more exhaustively the relation between our non-Abelian T-dual solution and the Sfetsos-Thompson solution presented in Section 6. This would amount to a way of connecting Gaiotto-Maldacena backgrounds (dual to Gaiotto field theories) with Lin-Maldacena solutions (dual to SU (2|4) quantum field theories).
It is clear that-by construction-the completion proposed in Section 4 to get a welldefined (0+1) dual QFT is just one among many possible ones. The only reason we completed in the way presented here is to ease the calculations. But this in turn, suggests that there exist a host of different quantum field theories and associated string duals that share a common sector, represented by the solution obtained via non-Abelian T-duality. It would be interesting to study this universality in more detail, which points to a generalised form of Eguchi-Kawai mechanism of reduction [50].
It would also be of interest to combine our solution and its completion with the blackhole geometry discussed in Section 3. The cloaking of the singularity makes it a perfect arena to study finite temperature effects. A computer simulation of the matrix model can be compared with the predictions calculated from our Type IIA background supplemented by the black hole. Indeed, following the lead of the papers [51], a simulation for the PWMM at finite temperature should be feasible. Then, correlators calculated with the string background could be checked against lattice simulations.
There are also a number of points of more technical nature that it would be nice to clarify. One of them is the study of the multipole expansion that follows from sub-leading terms to eq.(2.8), and the similar 1/L-sub-leading terms of eq.(4.4). The study of the different terms in the expansion, might contain interesting structure. It would also be of interest to study some sectors of the (0,2) SCFT in six dimensions in terms of the PWMM. Similarly, the fact that our completion in Section 5 involves M5-branes points to a feature that threads many different solutions obtained via non-Abelian T-duality. We hope to report on these and other interesting questions soon.