Perturbative solutions of N=1∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}={1}^{\ast } $$\end{document} holography on S4

We apply the recently proposed perturbative technique to solve the supergravity BPS equations of N=1∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}={1}^{\ast } $$\end{document} theories put on S4. In particular, we have calculated the coefficients of the leading quartic terms exactly, in the expression of the universal part for the holographic free energy as a function of the mass parameters. We also report on the coefficients of higher order terms up to 10th order, which are computed numerically.


Introduction
According to AdS/CFT correspondence [1], we have a large number of duality pairs where we have on one side strongly coupled gauge field theories and on the other side supergravity solutions with an anti-de Sitter factor in string/M-theory. Thanks to the breakthrough made in [2], using the supersymmetric localization technique, at least for theories with enough number of supersymmetries and appropriate coupling to background geometry, one can reduce a certain class of path integrals to a finite-dimensional matrix integral. When the field theory does have a holographic dual, the large-N limit of the matrix integral should match the counterpart quantities, e.g. Wilson loops and free energy, in supergravity. This program has been successfully applied to many examples in various dimensions, see for instance [3] for more details and references.
In this paper we are interested in non-conformal deformations of AdS/CFT. Within the lower-dimensional supergravity theories, which are related to 10/11 dimensional supergravity through a consistent truncation, we have scalar fields which are dual to relevant and marginal operators in the dual theory. In particular, those scalars which are dual to mass terms can be identified and the associated BPS flow equations can be written down. These BPS equations are first-order nonlinear differential equations, and exact solutions are usually not available. When one evaluates the holographically renormalized action for regular solutions, it is expected to match the large-N limit of mass-deformed free energy for the dual field theory. This task is sometimes called "precision holography" for mass-deformed conformal field theories. It started in [4], where exact supergravity BPS solutions to massdeformed ABJM model [5] were constructed and its renormalized supergravity action was shown to agree with the large-N limit of localization computation with supersymmetric mass terms. Then the authors of [6] tackled a similar problem of comparing the so-called N = 2 * super Yang-Mills theory in D = 4 (which can be obtained by giving mass to JHEP07(2019)169 N = 4 super Yang-Mills theory in a way compatible with N = 2 supersymmetry), and their supergravity dual solutions. By analyzing the numerically constructed supergravity solutions, it was argued in [6] that the evaluated action matches the prediction of localization computation. See also [7] for the uplift of the mass-deformed supergravity solutions from D = 5 to D = 10 and its application.
In the construction of supergravity solutions and their holographic analysis, the crucial information we need is a relation between the integration constants in the UV expansion, constrained by regularity at IR. In earlier works like [6,8,9], such relations were inferred from numerical solutions. It is certainly not satisfactory, especially when analytic expressions are available from the localization computations. In a recent work [10], one of the present authors proposed a perturbative method using which one can analytically extract the aforementioned relation in a series expansion form. In [10], three non-conformal holography problems were addressed: mass-deformed ABJM, the mass deformation of Brandhuber-Oz superconformal field theory in D = 5 [11], and also the N = 2 * deformation of N = 4, D = 4 super Yang-Mills. For the first two examples, one can either find exact supergravity solutions [4], or at least the series form of the relation between integration constants can be summed [10]. On the other hand, for the case of mass-deformed AdS 5 , the linearized equations involve log and polylogarithm terms which render explicit integrations difficult, and [10] resorted to approximation using series expansion up to high but finite orders. It is the purpose of this paper to improve and extend the study of massdeformed N = 4, D = 4 super Yang-Mills in the gravity dual side. In particular, we tackle the problem of N = 1 mass deformations, which is also called N = 1 * models. The dual gravity action was already constructed as a consistent truncation of D = 5 maximal SO(6)gauged supergravity by Bobev et al. in [8]. Combining both brute-force integration and also Padé approximants, we report a series expansion form of gravity side free energy, containing two more coefficients than the numerical results of [8]. We also push the integration analytically to 3rd order, while in [10] we used approximation already at 3rd order.
The plan of this report is as follows. In section 2 we present the five-dimensional Einstein-scalar action we are interested in, and also its associated BPS equations. In section 3, we describe the general feature of our perturbative approach when applied to holographic N = 1 * system. We also report on the analytic results up to 3rd-order perturbation, for general three-mass case. In section 4 we consider three special cases of single mass models, and report on the expression of holographic free energy obtained from perturbation up to 9th order. We conclude with discussions in section 5.

Summary of previous works
The aim of this paper is to re-visit the analysis of BPS equations presented in [8]. The action derived there is for a five-dimensional Einstein gravity which additionally contains ten (real) 1 scalar fields. In terms of AdS/CFT, those scalars are dual to various operators JHEP07(2019)169 appearing in the action of N = 4 super Yang-Mills theory when we put the theory on S 4 and turn on mass terms of three chiral multiplets and keep N = 1 supersymmetry. In this article we will closely follow and use the formulae in [8], with minor notational changes. 2 Thanks to supersymmetry, the Lagrangian density can be written succinctly in the following manner. and W is the same as W, except for the replacement of z a byzā.
Note that here η 1 , η 2 are two real scalar fields, while z a ,z b (a,b = 1, 2, 3, 4) together originally constitute four complex scalars. Since we consider supergravity in Euclidean signature, z a andzā are not mutually conjugate any more and we treat them as independent real scalars. We choose conformal gauge for the metric convention, (

2.4)
Then the AdS vacuum solution has unit radius, with the following scalars and warp factor The superpotential W and the Kähler potential K carry the information of supersymmetry transformations and eventually determine the BPS equations which provide firstorder differential relations for the scalar fields and warp factor. In the conformal gauge eq. (2.4), the BPS equations can be written as follows. 3 These equations are equivalent to eq. (4.14) in [8]. Here they are re-arranged to better suit our perturbative prescription.

JHEP07(2019)169
One may try to solve these equations near r = 1 (UV). Then it turns out that generally the solutions contain eight integration constants. For more details readers are referred to eq. (4.19) in [8], and also (B.6) of this paper. In terms of ρ = 2 tanh −1 r, the UV expansion (ρ → ∞) contains As explained in [8], µ 1 , µ 2 , µ 3 are interpreted as sources for the mass of three chiral multiplets in the field theory, v 1 , v 2 , v 3 are the expectation values of mass term operators, w is dual to gaugino expectation value, and s is to be identified with the Yang-Mills coupling constant. We also note that, to be precise the holographic dictionary identifies µ i = ±im i a [6,8], where m i is the mass of chiral multiplets and a is the radius of S 4 where the gauge field theory is put on. The crucial information we need is how v i are determined as functions of µ i , once we demand IR regularity (r = 0). More specifically, according to the holographic computation in [8], the holographic free energy F satisfies It turns out that the localization side computation for F contains a scheme-dependent factor, which can be removed when we take third or higher derivative. Integration of the above expression with F = F = F = 0 is thus called the universal part. In this paper from now on we will always mean the universal part, when we refer to F .

Perturbation for general solutions
In our perturbative approach we take the Euclidean AdS (i.e. hyperbolic space) as a zerothorder reference solution. At first, we turn on scalars z a ,zb, while A, η are still at vacuum configuration. Then at the second order, through gravitational and inter-scalar interactions we begin to have non-vacuum values for warp factor A and scalar field η. It turns out that via appropriate choice of the perturbative parameter , without losing generality we may set

JHEP07(2019)169
We substitute this perturbative expansion into the BPS equations and choose the upper sign for concreteness. Demanding that the equations should be satisfied for arbitrary , and doing some algebraic manipulation, one obtains the following form of equations.
The right-hand-side expressions here contain some rational functions of r, and also func- We note that these equations are always homogeneous: if we assign weight k to functions z and the in-homogeneous parts are known functions determined from lower orders of . Solving the first two coupled differential equations is thus in principle straightforward. Then this result can be substituted into H (2n) , determining η (2n) algebraically. A (2n) is determined through integration of Σ (2n) , which does not contain A (2n) and the rest of perturbative functions are determined already from previous steps. This way, we can determine the solution iteratively to any higher orders, in principle. When we perform the integration, a guiding principle is that the scalar fields should vanish at r = 1 (in the UV) because we want the solution should be asymptotically AdS. We also demand all functions are regular (i.e. non-divergent) at r = 0 (in the IR).
Let us now illustrate how this scenario leads to explicit solutions. At leading non-trivial order, n = 1, and we obtain Ξ (1) 3 .

(3.3)
For Ξ (1) the expression is essentially the same, except for exchange of z (1) vs.z (1) etc. We see that at this order we just have a coupled system of linear differential equations. General solutions should contain eight integration constants, and one can write them as follows z (1) 1 = c 1 u 1 +c 5 u 3 + c 5 +c 1 u 5 /r,z
The corrections to warp factor start to appear at O( 2 ). Substituting our O( ) solutions, we obtain We need to integrate it with boundary condition A(r = 1) = 0. The result is rather messy, containing polylogarithms. Explicit result can be found in the appendix. Then the real scalars η 1 , η 2 are determined algebraically, and η Again their explicit forms are relegated to the appendix. At O( 3 ), we substitute our O( 2 ) solutions which play a role of in-homogeneous terms which we need to integrate. Although it is not impossible, the results are quite messy and we do not try to present the result in full detail in this paper. However their precise UV asymptotic behavior can be more easily studied, as the derivative of O( 3 ) solution is JHEP07(2019)169 determined by O( 2 ) data through BPS equations. We have managed to determine v i as functions of µ i . The result is as follows, We can also calculate the gaugino condensate holographically, using our solutions. According to the analysis of [8], it is given by w(µ) which appears in the last line of UV expansion in (2.7). Our O( 3 ) result gives (3.11) which is in agreement with the conjecture made in [8], based on numerical solutions.

Further analysis of single mass models
Now we specialize to three special sub-sectors of the general N = 1 * models, following [8].
On the gauge theory side, we consider first the N = 2 * theory where we make a hypermultiplet massive in super Yang-Mills theory. Then we also consider two special cases of N = 1 * deformations.

N = 2 * model
Compared to the undeformed N = 4 super Yang-Mills, here we give the same non-zero mass to two chiral multiplets (i.e. a hypermultiplet). The relevant supergravity truncation was constructed earlier in [6]. If we start with the BPS system of the general N = 1 * cases, we set z 2 = z 4 = 0,z 2 =z 4 = 0, z := z 1 = z 3 ,z :=z 1 =z 3 , η 1 := η, η 2 = 1. Then the action simplifies to We recall the expansion of the solutions near UV, done in the Fefferman-Graham coordinate ρ = 2 tanh −1 r. In particular, the parameters v, µ are defined in terms of the ρ → ∞ asymptotic behavior of scalar fields. which agrees exactly with the large-N limit of localization computation [12]. This problem was re-visited in [10] using the perturbative technique we employ in this paper. We note that our refined result up to O( 3 ) in eq. (3.10) is consistent with the above formula, when e.g. we set µ 3 = 0, µ : Approximation using a truncated series expansion at IR (r = 0), to solve the BPS equations, was reported already in [10]. For the results reported in this paper, we have used an improved method: the BPS equations are solved by series expansion up to certain order, and then the remaining parts are replaced by Padé approximation for substitution to higher -order equations. We have performed the computation up to O( 9 ), and find that the result is v(µ) = −2µ + 1.00017µ 3 + 0.500022µ 5 + 0.333344µ 7 + 0.250378µ 9 + · · · . (4.5) It is confirmed that our perturbative solution agrees nicely with eq. (4.4).
In order to illustrate the reliability of our method, in figure 1 we show how the values of the expansion coefficients of µ 3 , µ 5 in eq. (4.5), extracted from the limiting behavior of z 1 , converge as we increase the order of truncated series solution to the BPS equations. We obtain similar results for the single mass case with N = 1 mass deformation.

Equal mass N = 1 * model
In this case one gives the same non-zero mass to all three chiral multiplets in N = 4 super Yang-Mills theory. On the supergravity side it is implemented by z 2 = z 3 = −z 4 , z 2 =z 3 = −z 4 , η 1 = 1, η 2 = 1. Then the action reduces to (4.12) One finds that the UV expansion of the BPS equations leads to the following results, to all expansion coefficients at higher orders in . We believe it is due to essentially the same difficulty that the authors of [8] could extract the coefficients only upto µ 4 in eq. (5.12) for this case.
In order to isolate fields with different UV asymptotics, we solve the BPS equations for re-defined functions y 1 ≡ z 1 + z 2 , y 2 ≡ z 1 − 3z 2 (and analogously for barred ones). Then the function v(µ) can be extracted from y 1 andȳ 1 . Due to the problem of zero modes, it turns out that the UV behavior and equivalently the Taylor coefficients of v(µ) are best extracted for the truncated series expansion in r with intermediate lengths. We find that the UV limits are stable for the coefficients of µ 3 when we truncate to the range of 300−400 orders in r. For higher orders of perturbation, the limits are best taken within the range of truncation to 150−250 (µ 5 ), and 110−150 (µ 7 ). See figure 2. We report, with a truncation upto 150 orders in r and Padé approximation, (4.14) Then the free energy is F/N 2 = −0.0433755µ 4 + 0.000693235µ 6 − 0.0000304161µ 8 + · · · . (4.15) Using our results we can also compute the gaugino condensate, w(µ). Note that this UV parameter vanishes for both N = 2 * and N = 1 * single mass deformation. Based on numerical solutions, it was conjectured w(µ) = 2µ 3 , precisely. Our analysis also supports this conjecture, and for instance we have witnessed that the coefficient of µ 7 term in w(µ) can be suppressed to less than 10 −4 .

Discussions
In this paper we have re-visited the problem of studying N = 1 * theories put on S 4 , using the perturbative prescription recently advocated in [10]. Presently our goal is to calculate the holographic free energy F . In a series expansion form, F is an even function of mass parameters µ i , starting with µ 4 . We have calculated terms up to µ 10 for single mass model, and µ 8 for equal mass model. Unfortunately however, with our current technology of supersymmetric localization, it is not feasible to calculate the counterparts in gauge field theory since S 4 localization needs N = 2 supersymmetry. We hope, one day, our predictions can be confirmed by bona fide field theory calculations.