Comments on Squashed-sphere Partition Functions

We study the partition function of odd-dimensional conformal field theories placed on spheres with a squashed metric. We establish that the round sphere provides a local extremum for the free energy which, in general, is not a global extremum. In addition, we show that the leading quadratic correction to the free energy around this extremum is proportional to the coefficient, $C_T$, determining the two-point function of the energy-momentum tensor in the CFT. For three-dimensional CFTs, we compute explicitly this proportionality constant for a class of squashing deformations which preserve an $SU(2)\times U(1)$ isometry group on the sphere. In addition, we evaluate the free energy as a function of the squashing parameter for theories of free bosons, free fermions, as well as CFTs holographically dual to Einstein gravity with a negative cosmological constant. We observe that, after suitable normalization, the dependence of the free energy on the squashing parameter for all these theories is nearly universal for a large region of parameter space and is well approximated by a simple quadratic function arising from holography. We generalize our results to five-dimensional CFTs and, in this context, we also study theories holographically dual to six-dimensional Gauss-Bonnet gravity.


Introduction and summary of results
Coupling a quantum field theory to background fields is a powerful tool which can be utilized to gain insights into the dynamics of the theory. This idea has been used often in studying quantum field theories, and is typically most powerful when the QFT at hand enjoys extra symmetries which constrain its dynamics. A well-known example that has gathered a lot of attention recently is coupling supersymmetric theories in various dimensions to background metrics and background gauge fields for the global symmetries of the theory. One can then often rely on the powerful tool of supersymmetric localization to extract many exact results for large classes of QFTs -see [1] for a review. In the absence of supersymmetry, one loses much of the technical control but, as we show below, if the theory at hand possesses conformal symmetry, some progress can still be made.
Our interest here is in studying Euclidean conformal field theories in odd dimensions, and couple them to a curved background metric. Every CFT has a conserved energy momentum tensor which naturally couples to the background metric and, thus, the type of deformations we study have a universal character and exist for any CFT. We choose to study theories in odd dimensions since these are not plagued by conformal anomalies. In addition, there is a natural class of deformations of the Einstein metric on odd-dimensional spheres which allow for many explicit computations.
The main focus of our work is to understand how the partition function or, equivalently, the free energy of general CFTs depends on the background metric. To achieve this, we can place the CFT on a conformally flat background metricḡ µν , a prime example being the round d-dimensional sphere S d , and switch on a "squashing" deformation of this metric controlled by the symmetric tensor h µν where is a real parameter. We can then study the free energy 1 of the CFT as a function of . In Section 2, we show that, for any odd-dimensional CFT, this function has the following series expansion for small where andG(ḡ µν , h µν , d) is a function which depends (see (2.17) below) on the corresponding conformally flat metricḡ µν , the metric deformation h µν and the spacetime dimension, d. 2 The real number C T is theory dependent and positive for unitary CFTs. It is often referred to as the "central charge" and determines the two-point function of the energy-momentum tensor. In particular, for CFTs in Euclidean flat space R d with coordinates X a one has [4] 3 where S d ≡ 2π d/2 /Γ[d/2] is the area of the unit (d − 1)-sphere and the tensorial object I ab,cd (X) reads I ab,cd (X) ≡ 1 2 (I ac (X)I bd (X) + I ad (X)I bc (X)) − 1 d δ ab δ cd , where I ac (X) ≡ δ ac − 2 X a X c |X| 2 . (1.4) Computing the explicit form of F (0) and the higher order terms in (1.1) is hard for general squashing deformations. To make further progress, we focus on squashed spheres S d , which are characterized by simple explicit metrics with a large isometry group. The squashed spheres we study are singled out by the fact that, just like the usual odd-dimensional round spheres S d , they can be described as "Hopf fibrations" over the complex projective spaces CP k (we use the notation k ≡ (d − 1)/2), i.e., The metric we put on these squashed spheres is given by (1. 6) where ds 2 CP k is the Einstein metric on the complex projective space CP k normalized in a way such that R ab = g ab , ψ is a periodic coordinate parametrizing the S 1 and J = dA CP k , where J is the 1 In this work we focus on the real part of the free energy of odd-dimensional CFTs and thus avoid the subtleties associated with contact terms that can affect its imaginary part, see [2,3] 2 A similar formula has been derived for entanglement entropies of deformed spherical and planar regions. In Section 5.3 we provide some comments on the connection between these two results. 3 Note that our convention for C T is the same as in [5], and differs from that in e.g., [4,6] by a factor of 1/S 2 d = Γ[d/2] 2 /(4π d ). For example, in our convention, C T = d/(d − 1) for a free real scalar field.
Kähler form on CP d− 1 2 -see Appendix B for more details. Note that we find it convenient to parametrize the squashing deformation in (1.1) as = − α 1+α with α ≥ −1. The metric on the usual round sphere is restored for α = = 0. For non-vanishing α, the metric in (1.6) preserves an SU ( d+1 2 ) × U (1) subgroup of the SO(d + 1) isometry group of the round d-dimensional sphere. The free energy of an odd-dimensional CFT on the round sphere is a useful quantity since it is conjectured to be monotonic along unitary RG flows. More precisely the quantity conjectured to decrease along RG flows for odd-dimensional QFTs is [7] F = (−1) (d−1)/2 log |Z S d | . (1.7) This F -theorem was proven for three-dimensional CFTs in [8] using entanglement entropy techniques (see [9] for a review and further references), but remains a conjecture for higher-dimensional theories. The squashing deformations of the round sphere metric that we study here, break conformal invariance in the CFT and thus induce an RG flow. This flow is however outside the scope of the established F -theorem, since the induced CFT perturbation is driven by the energy-momentum tensor, which is an operator of spin 2. Nevertheless, it is still natural to expect that the free energy should be a local maximum in the space of deformations controlled by the parameter in (1.1). This expectation indeed bears out. For CFTs in three and five dimensions, and the particular metric deformation in (1.6), we are able to explicitly compute the quadratic correction toF in the small expansion, and find 4F For unitary CFTs, the central charge C T is positive and thus the free energy on the round spherẽ F(0) is a local maximum in the space of squashing deformations given by (1.6).
In addition to the general results discussed above, we are able to compute the functionF( ) for general values of the squashing parameter in various theories in three and five dimensions. In three-dimensional CFTs, we can compute this free energy for a conformally coupled scalar and a free Dirac fermion by utilizing the explicitly known eigenvalues of the Laplacian and the Dirac operator for the metric in (1.6). The calculation then boils down to an evaluation of a Gaussian path integral which we are able to carry out numerically using an approach similar to the one employed in [10][11][12]. These explicit calculations are in harmony with the general results in (1.1) and (1.8). For threedimensional CFTs with weakly coupled gravity duals, we have access to analytic results forF( ). The CFT at hand placed on the background with metric (1.6) is described holographically by the well-known AdS-Taub-NUT solution [13,14]. Via the standard holographic dictionary, the partition function of the CFT is captured by the properly regularized on-shell gravitational action for such solution. The result is quite simple and somewhat surprising. The cubic and higher order terms in in (1.1) vanish, and C T completely controls the simple quadratic dependence of the squashed sphere partition function on through the expressions in (1.1), (1.2), and (1.8). We should stress that the free energy for CFTs in this squashed sphere background has been studied before in [15], where the free energy of the free and interacting O(N ) model was discussed in the context of higher-spin holography -see also [10][11][12] for a more recent discussion. It is also worth pointing out that the partition function of three-dimensional N = 2 SCFTs on the squashed sphere with the metric (1.6) was studied in a number of papers using the tools of supersymmetric localization -see for example [16,17]. In the supersymmetric context, in addition to the background metric, one has to turn on a non-trivial background gauge field for the U (1) R-symmetry current present in every such CFT. Due to this reason, the results in [18], although similar in spirit, are different from ours since we do not use supersymmetry and the only background field in our setup is the metric.
To confirm the general structure of the free energy, and the result in (1.8), we also evaluate explicitly the functionF( ) in five dimensions for two particular classes of CFTs. For a fivedimensional conformally coupled scalar, we can again explore the high-symmetry of the metric in (1.6) to evaluate explicitly the eigenvalues of the Laplacian and subsequently study numerically the resulting path integral. We find excellent agreement with the general formulas in (1.1), (1.2), and (1.8). In addition to this free CFT, we also employ a bottom-up holographic approach and use Gauss-Bonnet gravity with a negative cosmological constant as a dual description to a large class of putative five-dimensional CFTs. This class of gravitational theories admit AdS-Taub-NUT type solutions which have the squashed metric in (1.6) as an asymptotic boundary. Using holographic technology, we are able to compute explicitly the squashed sphere free energy, and find that it is a cubic function of , given in (4.20) below, with coefficients determined by the central charge C T and one of the coefficients, usually dubbed t 2 , determining the three-point function of the energymomentum tensor.
The outline of this paper is as follows. We begin in the next Section with a general discussion of the partition function of odd-dimensional CFTs on deformed conformally flat manifolds. In Section 3, we study the partition function of CFTs on a squashed three sphere and compute this function explicitly for free boson and fermion theories as well as for CFTs with a weakly coupled Einstein holographic dual. In Section 4, we extend this analysis to the case of five-dimensional CFTs placed on a squashed five sphere and compute the partition function explicitly for a free boson as well as for CFTs dual to six-dimensional Gauss-Bonnet gravity. We conclude in Section 5 with some comments on supersymmetric theories on the squashed sphere as well as on an extension of our results to more general deformations of the round sphere. We also discuss the connection between our study and recent results on the shape-dependence of entanglement entropy in CFTs. In the four appendices, we summarize the structure of three-point correlation functions of the energymomentum tensor in CFTs, present details on the holographic evaluation of the partition function, the spectrum of the scalar Laplacian on squashed spheres, as well as a summary of the numerical approach to calculating the partition function of free theories.

Partition functions on deformed manifolds
Let us consider the Euclidean partition function of a general CFT on a d-dimensional 5 manifold M with metric g µν where ϕ stands schematically for the set of dynamical fields in the theory. We wish to understand how the free energy F ≡ − log Z changes under small deformations of the metric. Let us parametrize 5 We will soon restrict d to be an odd integer.
such deformations as whereḡ µν is some arbitrary reference metric on M, is a real constant, and h µν encodes the geometry of the perturbation. If we assume that | | 1, we can expand the action in a power series in as follows where, as usual with g andḡ, we denote the determinant of the respective metrics, and we have defined h µν ≡ −ḡ µαḡνβ h αβ , h ≡ḡ µν h µν . 6 We have also used the definition of the stress-energy tensor Hence, the partition function on the deformed background reads Now, let us further assume that the free energy is an analytic function of , i.e., that we can expand it as . A straightforward calculation of the leading correction to F(0) yields Whenever the one point function of the energy momentum tensor in the CFT on the undeformed space vanishes one has F (0) = 0, and thus the background associated toḡ µν is a local extremum of the free energy. From now on we will assume that (2.8) holds. This is the case for odd-dimensional CFTs placed on a background metricḡ µν on M which is conformally flat. An example of this situation which will be of particular interest in our work is the round sphere in odd dimensions. The fact that F (0) = 0 for small perturbations around the round sphere is in the spirit of the F -theorem in three dimensions [7,19]. Note, however, that the two statements are not identical since we are considering deformations of the CFT with the energy momentum tensor, which is not a scalar operator. Assuming that (2.8) holds, one finds the following expression for the coefficient of the second order correction to the free energy Analogously, the coefficient of the third order correction reads If the CFT at hand lives in an odd number of dimensions, and the metricḡ µν on M is conformally flat, one can show that the terms underlined with a bracket in (2.9) and (2.10) vanish. This is due to the fact that one point functions of local operators in CFTs should vanish in the absence of conformal anomalies. Since we are interested in deformations of odd-dimensional round spheres, from now on we will ignore the underlined terms in (2.9) and (2.10). In that case, the leading correction to the free energy is quadratic in with coefficient given by Observe that ifḡ µν is the flat metric on R d , (2.8) is automatically satisfied and C T controls the leading correction to F(0), as it is obvious from (1.3). Now let us take M to be a (locally) conformally flat manifold, i.e., let us assume it to be related to R d through a conformal mapping Then, the metrics of both spaces are related through where f * is the pullback and Ω 2 (x) the corresponding conformal factor. In odd-dimensional theories, the stress tensors for the theories on M and R d are related through 7 and where we denoted the coordinates in M and R d by {x µ } and {X a } respectively. For the correlator in the integrand of (2.11), one finds This expression allows us to write the two-point function on M in terms of the two-point function on R d , whose explicit expression is given by (1.3).
To recap, we find that on any odd-dimensional manifold M for which the map f in (2.12) exists, the free energy behaves, under a small deformation of the metric, as where (2.17) and where in the last term it is understood that we need to write all contributions as functions of {x} and {y} using the map f defined in (2.12).
In the rest of the paper, we will focus on three-and five-dimensional squashed spheres, which provide a class of deformed conformally flat manifolds. For these examples we will be able to perform the integrals in (2.17) analytically and thus obtain a simple explicit expression for F (0). 8 A few comments about the third order correction F (0) are in order. The terms on the first and third lines of (2.10) which are not underlined depend on the particular CFT of interest only through the constants C T , t 2 , and t 4 appearing in the three-point function of the energy momentum tensor -see appendix A for details. Thus, these contributions are in principle computable for general CFTs. However, the one on the second line depends on the details of the CFT at hand and, in particular, on OPE coefficients of local operators with the energy momentum tensor. This leads to the conclusion that it seems hard to obtain an explicit expression for F (0) valid for general CFTs. In the examples discussed below, we will be able to estimate F (0) for squashed spheres in some specific CFTs, including free theories as well as theories with a weakly curved holographic dual.
We would like to emphasize that in our discussion we have treated the free energy as a real quantity. In Euclidean signature this is imprecise and in general one should worry about possible contact terms with imaginary coefficients in the energy-momentum tensor correlators. 9 These type of contact terms are discussed in some detail in [2,3] for three-dimensional CFTs on curved manifolds. In this context the relevant finite contact term for our discussion has an imaginary coefficient and is given by where ω µ is the spin connection on M. This gravitational Chern-Simons term can be added to the Lagrangian with an arbitrary integer coefficient and will in general affect the imaginary part of the free energy of the CFT at hand. One can also imagine that a similar gravitational Chern-Simons term built out of the spin connection can be added also to CFTs in higher odd dimensions. We expect that these counterterms will in general have an imaginary coefficient in Euclidean signature and thus they will not affect the real part of the squashed sphere free energy which is the main focus of our work.

Three-dimensional squashed spheres
In this Section, we study the free energy of three-dimensional conformal field theories on a squashed sphere as a function of the squashing parameter. In this case, the metric in (1.6) becomes 10 ds 2 where ds 2 CP 1 is the usual round metric on S 2 , 2) and the coordinate range is 0 ≤ θ ≤ π, 0 ≤ φ ≤ 2π, 0 ≤ψ ≤ 4π. Note that in this parametrization we have the following relation for the small parameter defined in (2.2) At first order in α (which is also first order in ), we can write (3.1) as in (2.2), whereḡ µν is the metric of the round sphere  Clearly, this is a very symmetric squashing deformation of the metric on the round sphere, which will allow for the explicit results below. 10 Note that in the introduction and in appendix B we write ds 2 To evaluate the integral in (2.17), we need the explicit form of transformations in (2.13) and (2.14). To find this, we employ a conformal mapping (f −1 ) from R 3 to S 3 as follows In these coordinates, the metric of the round S 3 (3.4) is related to the flat metric on R 3 through With these ingredients, we are ready to evaluate F S 3 α (0) using (2.17). A somewhat tedious, but otherwise straightforward, calculation yields Notice that for unitary CFTs, C T is a positive number, which in turn implies that the partition function on the round S 3 is a local maximum in the space of squashing deformations considered here. Similar observations were made in the case of CFTs with an Einstein holographic dual, a free scalar theory and the interacting O(N ) model in three dimensions in [10,15]. Our result shows that this holds for general CFTs and, furthermore, that the leading correction to the round sphere result, which is always quadratic in the deformation, is controlled by the central charge C T (3.8).
From these results, we can conclude that for general three-dimensional CFTs, the free energy on a squashed three-sphere with the metric (3.1) is given by for small values of the squashing parameter α (or ). In the following subsections, we study F S 3 α for general values of α, i.e., beyond the quadratic approximation, for free CFTs as well as for theories with weakly coupled gravitational duals. The results we find, in the regime of small α, are in harmony with the general result in (3.9).

Conformally-coupled scalar
Let us consider a conformally-coupled scalar field φ. In order to compute the free energy F sc, S 3 α , we need to evaluate the partition function on the squashed-sphere background (3.1). 11 After a Gaussian integration, we find where we have introduced the UV cutoff Λ. In order to evaluate (3.11) on the squashed-sphere (3.1), we need the eigenvalues of the conformal Laplacian on this manifold. These were used in a similar context in [10,22], and are given by 12 with n = 1, 2, . . . and q = 0, 1, . . . , n − 1.
The partition function Z sc has UV divergences and needs to be regularized. While it is possible to use zeta-function regularization as in [15,22], this leads to somewhat cumbersome calculations that are non-trivial to extend to other CFTs, like the free fermion considered in the following subsection, higher-dimensions, or to spheres with more complicated deformations. Therefore, we will use the numerical regularization methods employed in [10][11][12]. The essential idea behind these is to use a "soft regulator", which makes it possible to split the free energy into UV-and IR-dependent parts. The IR-part consists of a convergent sum which can be straightforwardly calculated numerically, while the UV one contains all the divergences. The only counter-terms that can be added to the action to control these divergences are where Λ is the cutoff and A and B are real coefficients which depend on the details of the theory at hand. The regularization procedure amounts to performing a numerical fit of the on-shell action plus the counterterms in (3.13). From this, we can determine the values of A and B which allow us to remove the unwanted divergences. More details about the numerical procedure we used can be found in appendix D (see also [12]). After subtracting the divergences, we obtain the regularized free energy as a function of the squashing parameter α. We can perform this calculation for arbitrary values of α and our results for the free energy F sc, S 3 α are depicted 13 in Figure 1. Two features of our numerical results are immediately noticeable. First, the round sphere at α = 0 leads to a local maximum of the free energy. For the large range of values of α that we studied, this is also a global maximum and it is natural to conjecture that this is true for all allowed values in the range α > −1. Second, for large positive values of α the free energy becomes a linear function. 14 This linear behavior for large α was also observed and discussed in [15,23]. In particular, in [23], it was argued that in the limit α 1 the free energy is given by A fit to our numerical curve yields instead 12 For more details, see appendix C. 13 As explained in Section 3.4 we find it convenient to subtract the free energy of the round sphere and normalize the result by dividing by C T . 14 We are grateful to Louise Anderson and Paul McFadden for a useful correspondence on this topic.
It is clear that the leading O(α) term in our numerical results agrees very well with the result in [23]. This is not the case, however, of the subleading contributions. It would be interesting to identify the source of these discrepancies, which do not seem attributable to numerical errors. 15 Figure 1: The free energy for a conformally coupled scalar on the squashed S 3 as a function of the squashing parameter α.
We provide further comments on the squashed sphere free energy for the free scalar and how it compares to the free energy of free fermions and holographic theories in Section 3.4. However, we can readily perform some consistency checks of our numerical procedure. First, we can evaluate the free energy for α = 0. This should correspond to the well-known result for the round three-sphere free energy, given by [7] F sc, From our numerics we find F sc, S 3 0.063807, in very good agreement with (3.16). According to our general result (2.16), the first derivative of the free energy with respect to the squashing parameter evaluated at α = 0 should vanish. Numerically, we find F sc, S 3 α (0) 2.25 · 10 −11 , which is in clear agreement with our expectations. In addition, according to (3.9), we expect to find −48/π 2 F sc, S 3 α (0) = C T , sc = 3/2. Numerically, we obtain −48/π 2 F sc, S 3 α (0) 1.5000004, which is clearly in very good agreement with the general result (3.9).

Free fermion
Let us now derive similar results for a free Dirac fermion. The partition function in this case is where / D is the Dirac operator on the corresponding curved background. The free energy can be found by solving the Gaussian integral where, once again, we introduce an energy cutoff Λ. We are interested in evaluating the determinant in (3.18) for the squashed sphere metric in (3.1). To this end we need the eigenvalues of the Dirac operator on this background. These have been found in [22,24,25], and read where n and q are integers. For the positive branch, denoted by +, n goes from 1 to ∞ and q from 0 to n, while for the negative branch, denoted by −, n goes from 2 to ∞ and q from 1 to n − 1. In evaluating the determinant in (3.18), we have to sum over both branches of the eigenvalues.
To compute the free energy as a function of α, we apply the same numerical procedure used in Section 3.1 (and described in Appendix D) for the free scalar. Our numerical results are presented in Figure 2. We again observe that the round sphere partition function at α = 0 is a local maximum. However, it is clearly not a global maximum since, interestingly, for large positive values of α we find that the free energy increases linearly. Including subleading terms, we find Just like for the scalar, an analytic expression was found in [23] for F f, S 3 α in this regime. This reads which, in this case, agrees very well with our numerical result 16 -see also Figure 2. In contrast to the scalar case, the fermion free energy reaches a (local) minimum at an intermediate value of α, given by α min 4.6228 .
To gain further confidence in our numerical results, we can perform some consistency checks. For α = 0, the partition function for a free Dirac fermion was evaluated analytically in [7] and reads With our numerical method for α = 0, we find F f, S 3 0.21895949, which is in excellent agreement with (3.23). The first and second derivatives of the free energy with respect to α, evaluated at α = 0 are given by F f, S 3 α (0) = 9.01 · 10 −7 and −48/π 2 F f, S 3 α (0) = 2.99999, respectively. After taking into account that for the free Dirac fermion we have C T , f = 3, these values are in very good agreement with our general results (2.16) and (3.9).
The free energy for a free Dirac fermion on the squashed S 3 as a function of the squashing parameter α. The round sphere is not a global maximum, as the free energy grows unbounded after reaching a local minimum for α min 4.6228. For large enough α, the free energy is given by Equation (3.21) [23], which we also plot (gray dashed line).

Holography
Let us now study the squashed sphere free energy of three-dimensional CFTs with an Einstein gravity holographic dual using the standard tools of AdS/CFT. In order to do so, let us consider an Euclidean bulk theory of the form where G is the Newton constant and the parameter sets the scale of the Anti-de Sitter (AdS) vacuum solution of the above action. The Euclidean AdS 4 solution of (3.24) is dual to the vacuum state of some three-dimensional CFT. 17 The gravitational dual of this CFT placed on the squashed sphere with metric (3.1) should be a gravitational background which is asymptotically locally (Euclidean) AdS 4 with a boundary metric given by (3.1). There is indeed a class of analytic solutions to the equations of motion derived from (3.24) which posses the required asymptotic structure. These are the so-called AdS-Taub-NUT and AdS-Taub-Bolt metrics which were studied in a holographic context in [13,14] (see also [12] for a recent account and more complete list of references). Thus, we can apply standard holographic tools to these analytic solutions to extract information about the dual CFT theory on S 3 α .
The AdS-Taub-NUT solution is the one which is thermodynamically stable for values of α less than α HP = 6 + 2 √ 10 12.3246, and thus dominates the gravitational partition function in that range. At α HP , there is a Hawking-Page phase transition and the thermodynamically favored solution is AdS-Taub-Bolt, which has a non-trivial two-cycle [27,28]. In addition, one should note that the AdS-Taub-Bolt solutions do not exist for α < α crit ≡ 5 + 3 √ 3 ≈ 10.1962 (see [12] for a more detailed discussion). Our interest here is on relatively small values of α, since we would like to study deformations that smoothly connect to the AdS 4 vacuum solution. Therefore, from now on we will focus on the AdS-Taub-NUT solution, which is reviewed in Appendix B. As we show there, the free energy of this background as a function of α can be computed via holographic renormalization and reads [27] where the label "E" stands for Einstein gravity.
To make contact with the field theory around (3.9), we need to express this result in terms of the "central charge" C T . For CFTs dual to Einstein gravity with the action (3.24), one has (see e.g., [6] and Footnote 3) Using this and (3.25), it is straightforward to show that the general field theory results (2.16) and (3.9) are obeyed. This is a non-trivial quantitative check of the holographic interpretation of the AdS-Taub-NUT solution. The free energy in (3.25) takes a more suggestive form when rewritten in terms of the parameter Notice that this is an exact result valid for finite 18 values of α (or ). In particular this result implies that for three-dimensional CFTs with holographic duals captured by the action in (3.24), the quadratic approximation in (2.16) and (3.9) is exact, and all higher-order corrections in powers of vanish. Notice also that α = 0 is a global maximum for the gravitational free energy (3.27), i.e., the free energy is extremized when the boundary is the round S 3 .

Comparison
It is time that we pause and compare the results for the free energies computed so far. We find it illuminating to study the dependence of the free energy on the squashing parameter by subtracting the value of the free energy on the round sphere and dividing by C T , i.e., by considering the quantity Indeed, it is natural to normalize the free energy by the central charge C T , given that, as shown in Section 2, this quantity controls the leading term in the expansion around α = 0. Moreover, it is clear from (3.9) that the round sphere value is a local maximum of F S 3 α . In Figure  3, we plot this normalized free energy as a function of α for the three theories considered in the previous subsections. In addition to that, we present the sum of the free energy of one free boson and one free Dirac fermion. It is unclear to us why the result for this combination is so similar to the one for general CFTs with a holographic dual. The holographic result for the free energy in (3.27) is a simple parabola when parametrized in terms of the squashing parameter . To this end in Figure 4 we plot the free energies for the theories discussed above as a function of . 3] for holographic theories dual to Einstein gravity (green solid line), a free scalar (blue dots), a free fermion (red dots) and the combination of both (orange solid line).
Some further comments are in order. The quadratic holographic result in (3.27) is a very good approximation for the free boson and fermion in the range ∈ (−0.5, 0.5). This is better than the expected resemblance controlled by the leading order quadratic result. To quantify this better, we can study the subleading cubic correction to F S 3 where q ≡ F (0) 6C T , as can be seen from (2.16). For the free scalar and free fermion, we can find the coefficient q from our numerical results. A conservative estimate leads to |q| < 10 −3 . Clearly, we have q = 0 for theories with holographic duals. This suggests that the cubic correction in (3.28) is either vanishing or has a very small coefficient for general CFTs. Notice that it seems hard to calculate the number q for general CFTs, since this coefficient is determined by the integral in equation (2.10). The various correlation functions (in particular the two-point function terms) appearing in the integrand of (2.10) are theory-specific and, in general, non-vanishing. Thus, the only way to have q = 0 is to either have a miraculous cancelation between the five terms in (2.10), or to obtain a zero value of the integrand due to the contractions with the very special form of the perturbation h µν in (3.5). It would be very interesting to establish whether q vanishes for general theories, or is simply some small non-zero number.
The fact that C T provides the "right" normalization for the free energy -in the sense that the leading contribution is universally controlled by this quantity, and that curves corresponding to very different (presumably general) theories collapse to remarkably close curves -is reminiscent of the behavior found in [29,30] for the dependence of the universal contribution to the entanglement entropy of three-dimensional CFTs from corner regions, customarily denoted a(θ), where θ is the corner opening angle. In this setup, the available results for a(θ)/C T , corresponding to free fields, holographic theories and Wilson-Fisher fixed points of the O(N ) model (for N = 1, 2, 3), were all shown to give rise to extremely similar curves for general values of θ. Analogously to our findings for the free energy of the nearly-round squashed sphere, such agreement was observed to be exact in the almost smooth corner limit. 19 We discuss the connection between these two, seemingly unrelated, quantities further in Section 5.3.
Another notable feature of the numerical results in Figure 3 and Figure 4 is that the free fermion and free boson curves serve as an "envelope" for the holographic result and the combination of both free fields. In fact, one can check that for a free CFT consisting of any number of free bosons, N b , and fermions, N f , the normalized squashed sphere free energy lies between the red and the blue curves in Figures 3 and 4. We are not able to provide a rigorous explanation of this curious fact but it is natural to offer some speculations. The higher-order terms in the expansion of F S 3 α around α = 0 (or equivalently = 0) are controlled by some combination of integrated n-point functions of the energy-momentum tensor. Could it be that the free boson and fermion provide an upper an lower bound for such integrated correlators? A hint that this might be true is provided by the three-point function of the energy-momentum tensor. 20 This correlator in general dimension d has three independent tensor structures and is entirely fixed by three independent constants, one of which is C T and the other two are often dubbed t 2 and t 4 . For kinematical reasons, in d = 3 one finds that t 2 = 0. The constant t 4 is bounded above by the value for the free boson, t sc 4 = 4, and bounded below by the free fermion value, t f 4 = −4, see [4,6,34,35]. 21 For CFTs with Einstein holographic duals, one finds t 4 = 0 [6,35]. Of course, the three-point function coefficients C T and t 4 should not determine the full behavior of the function F S 3 α , since higher point functions of T µν are not simply determined by constants, but are complicated (and in general not known) functions of the conformal invariant cross ratios. Nevertheless, it is tempting to speculate that the free boson and fermion curves in Figures 3 and 4 provide some analogue of the Hofman-Maldacena bounds [6,35], which would be valid for odd-dimensional CFTs with a broken conformal invariance induced by the squashing deformation controlled by α.
It is clear from the above discussion that it would be most interesting to have access to analytic expressions for F S 3 α for some CFTs (including free theories) beyond the holographic result in (3.27). While we shall not pursue this goal in the current work, let us mention in passing that a suggestive (approximate) relation between the free energies of the free boson and fermion involving the holographic result can be found as follows. Let us define , which is simply obtained by subtracting the holographic result in (3.27) from the normalized free energy. Then, within the limits of our numerical calculations, the relation T sc, α −π 2 T f, α /6 holds with a precision of ∼ 99% or better for all positive values of α up to α ∼ 8. 22
(4.1) 20 See appendix A for a review on the form of such three-point functions. 21 Notice that the extensive quantities are C T and C T t 4 (as well as C T t 2 when d ≥ 4). 22 Note that the factor ζ(2) = π 2 /6 ≈ 1.645 has been obtained by inspection and is subject to the uncertainty associated with our numerical procedure. A more precise determination of the proportionality constant might lead to a less charismatic value such as 1.646 or 1.644. Note also that, while it is remarkable that such a relation holds for not so small values of α, our numerics suggests that it stops being a good approximation for α ≥ 9. 23 Observe that, as opposed to the three-dimensional case studied in Section 3, here we use the same coordinates as in appendix B.
In this coordinates, which satisfy 0 ≤ ξ 1,2 ≤ π/2, −2π ≤ ψ 1 ≤ 0, 0 ≤ ψ ≤ 2π and 0 ≤ ψ 2 ≤ 3π, the metric of the round S 5 is (see Appendix B for details) ds 2 S 5 = dξ 2 2 + dψ 2 + sin 2 ξ 2 dξ 2 1 + dψ 2 1 /2 − dψ 1 cos(2ξ 1 ) (dψ 1 /2 + dψ 2 + dψ) (4.2) +dψ 2 (dψ 2 + 2dψ) + dψ 1 (dψ 2 + dψ)] , and is related to the usual flat metric on R 5 through The only non-vanishing component of the inverse metric perturbation h µν reads Our main interest is in evaluating the free energy of CFTs on this squashed sphere background. As discussed around (1.7), it is convenient in the context of the F -theorem to discuss the quantityF S 5 α instead, which for d = 5 differs by a sign from the usual definition of the free energy. In complete analogy with the d = 3 case studied in Section 3, we can now evaluateF S 5 α (0) using the general integral in (2.17). The calculation is tedious, but the final result takes the simple form This implies that, for general CFTs, the free energy on a squashed five-sphere is given bỹ for small values of the squashing parameter. We now proceed beyond this leading order approximation and evaluate the free energy for general values of α in some specific CFTs.

Conformally-coupled scalar
The partition function for a conformally coupled scalar field in a general five-dimensional background is given by As in the three-dimensional case, this leads to the following expression for the free energy: To evaluate this determinant, we need the eigenvalues of the conformal Laplacian operator as well as their multiplicities. In Appendix C, these are worked out for a general d-dimensional squashed sphere for odd values of d. The result for the eigenvalues, λ n,q , and the multiplicities, m n,q , in d = 5 reads λ n,q = (n − 1)(n + 3) + α(n − 1 − 2q) 2 + 3(5 + 6α) 4(1 + α) , (4.9) m n,q = (n + 1) (q + 1)(n − q) 2 , (4.10) where the integers n and q obey n ≥ 1 and 0 ≤ q ≤ n − 1. Plugging these eigenvalues into the free energy (4.8), gives rise to a diverging sum, on which we apply the same numerical regularization procedure as in the three-dimensional case. In particular, we use a soft regulator to separate the divergent and convergent parts. In five dimensions, the divergent terms are Subtracting these counterterms with appropriate coefficients, we obtain the regularized free energy as a function of the squashing parameter α. In Figure 5, we plot our results for the free energy for a broad range of values of α. Some additional comments on this function are provided in Section 4.3. As a consistency check of our numerical results, we can once again compare with the result for the free energy on a round S 5 , i.e at α = 0 [7] F sc, S 5 = 1 which is in excellent agreement with our numerical result. Similarly, we find thatF sc, S 5 α (0) ∼ 10 −10 , which is in harmony with the general result in (4.6). The five-dimensional free scalar has C T = 5/4, which is again in excellent agreement with (4.6) and our numerical estimate − 320 3π 2F sc, S 5 α (0) = 1.249995 ≈ C T .

Gauss-Bonnet holography
Another class of CFTs for which the squashed sphere partition function can be explicitly computed is provided by conformal theories with weakly coupled holographic duals. We will again adopt a bottom-up approach and will consider five-dimensional CFTs dual to six-dimensional Gauss-Bonnet gravity in the bulk. The Euclidean action of this theory is given by where X 4 = R µνρσ R µνρσ − 4R µν R µν + R 2 is the dimensionally-extended Euler density of fourmanifolds. The Gauss-Bonnet term is topological in four bulk dimensions, i.e., for d = 3, and does not modify the gravitational equations of motion for d ≤ 3. It does, however, contribute to the equations of motion for d > 3. Interestingly, those equations are still second order, 24 which in turn allows for the construction of many non-trivial analytic solutions of this theory in various dimensions. In particular, AdS-Taub-NUT and Bolt solutions to (4.13) have been constructed not only for d = 5, but for general odd values of d in [37] (when λ GB = 0) and in [38] (for general values of λ GB ). The boundary metric of the class of NUT solutions relevant for our purposes is precisely that of a squashed-sphere S 5 α (1.6) with (see Appendix B for details) where we have defined and the Gauss-Bonnet coupling is constrained to obey λ GB ≤ 1/4. 25 The free energy of these solutions was computed in [39,40] and can be written as 26 (4.16) The central charge C T of the dual CFT can be computed using standard holographic tools, and reads [6] (4.17) Using (4.16) and (4.17), it is easy to show that the general relation in (4.6) is obeyed. In the limit λ GB → 0, we have f ∞ → 1, and we obtain the Einstein gravity result Observe that for λ GB = 0, we have two dimensionless quantities, namely 4 /G and λ GB , which appear in the expression for the free energy (4.16). In particular, this means that in general we cannot write F GB, S 5 α in terms of C T alone, like in the case of Einstein gravity. To understand better the result in (4.16) in terms of field theory quantities, we can use some of the holographic results in [6]. In particular, it is known that the three-point function of the energy-momentum tensor of CFTs dual to Gauss-Bonnet gravity has a special structure which leads to the vanishing of one of the three independent constant coefficients, namely, t 4 = 0. In addition to C T in (4.17), the other non-vanishing coefficient which fully determines the T µν three-point function is With this at hand, it is possible to express (4.16) in terms of the squashing parameter and the coefficients C T and t 2 as Notice that the leading order quadratic term in (4.6) is supplemented only by a cubic term in .
In the context of holography, λ GB cannot take arbitrary values. Indeed, in general, it is possible to find lower and upper bounds on the couplings of higher-derivative gravities by using certain physically sensible criteria on the corresponding dual theories -see e.g., [6,35,[41][42][43][44][45][46]. The most stringent known bounds in the case of holographic theories dual to (d + 1)-dimensional Gauss-Bonnet gravity come from requiring that the energy flux is positive everywhere for certain scattering processes or, alternatively, from imposing that causality is respected in the dual CFT [6]. In the case of interest for us, i.e., for d = 5, the corresponding constraints read (4.21) In Figure 5 we have plotted (F GB, S 5 α −F GB, S 5 )/C T for all the values of λ GB allowed by (4.21).

Comparison
After computing the squashed sphere free energy for a conformally coupled scalar field and a oneparameter family of holographic theories, it is instructive to compare the results. As in Section 3.4, it is useful to perform this comparison by plotting the quantity (F S 5 α −F S 5 )/C T as a function of the squashing parameter. We present such plots in Figure 5 and Figure 6.
Similarly to the three-dimensional setup, the round sphere free energyF S 5 is a local maximum of the functionF S 5 α for all theories discussed above. This is expected, and follows from the general result in (4.6). It is clear however from Figure 6, and from the analytic holographic result in (4.20), that the round sphere free energy is not a global maximum. Another interesting feature of the curves in Figures 5 and 6 is that the analytic holographic results are very close to the numerical curve for the free boson in a surprisingly large region of values for the squashing parameter. In this case, it looks like the leading quadratic term produces an approximation which is considerably worse than its three-dimensional analog. It seems, however, that the Gauss-Bonnet result (4.20), which contains a cubic term in , does a much better job. Figure 5: A plot of (F S 5 α −F S 5 )/C T as a function of α for a free scalar (blue dots), and for holographic theories dual to Gauss-Bonnet gravity (green shaded region). The holographic results include the Einstein theory (green) and the limiting cases (4.21) allowed by causality in the dual CFT (purple and orange).
Some additional speculations and conjectures naturally emerge from our results in Figures 5 and  6, and some of the results in [6], valid at the conformal fixed point at α = 0. In particular it was shown in [6] that the following inequalities are obeyed by the coefficients t 2 and t 4 , which control the three-point function of the energy-momentum tensor in unitary CFTs For CFTs dual to the Gauss-Bonnet theory defined by the action in (4.13), one finds that t 4 = 0, and t 2 is given by (4.19). This, combined with the bounds in (4.22), leads to the bound on λ GB in (4.21). This bound on λ GB at the conformal fixed point translates into the shaded green region on the plots in Figures 5 and 6. This fact naturally suggests that the unitarity bounds on t 2 and t 4 in (4.22) may extend away from the α = 0 conformal point and ultimately result in a set of "enveloping" curves which constrain the function (F S 5 α −F S 5 )/C T for general CFTs on a squashed sphere. Notice that the fact that the numerical results for the free boson lie outside the shaded region in Figures 5 and 6 is not incompatible with this conjecture, since we have t 4 = 12 for the free boson, which is greater than the t 4 = 0 value for Gauss-Bonnet. In fact, if our conjecture is true, Figure 6: The same quantity as in Figure 5 as a function of . We have also included the leading universal quadratic term (dashed pink line).

Discussion
In this section, we offer a short discussion on several topics that are related to our main results presented above. In particular, we discuss the difference between the supersymmetric free energy evaluated by supersymmetric localization and the non-supersymmetric free energy studied in this work. We also study the effect on the free energy of a generalization of the squashed sphere metric in (3.1) which preserves less isometry and depends on one more squashing parameter. Finally we comment on the resemblance between our results and recent studies of entanglement entropy for deformed entangling regions. We conclude by pointing out several interesting avenues for further study.

Three-dimensional N = 2 SCFTs on the squashed S 3
As mentioned in Section 1, supersymmetric localization is a powerful tool to extract results for supersymmetric QFT observables on curved manifolds. In particular, it is possible to compute the path integral of three-dimensional N = 2 SCFTs on the squashed S 3 with the metric (3.1) with this method -see [16,17]. In addition to the non-trivial background metric, in this setup it is required to turn on a background field for the gauge field that couples to the superconformal U (1) R-current present in every such theory. This modifies the path integral and makes the resulting supersymmetric free energy, which we will denote by F, not equivalent to the intrinsically nonsupersymmetric free energy F discussed throughout this work. To illustrate this difference, we focus on the simple example of a theory of a single N = 2 chiral multiplet, X. The bosonic content of this multiplet is simply given by a complex scalar and corresponding complex fermion. The non-supersymmetric free energy F( ) for a free N = 2 chiral multiplet is thus easily obtained by combining our numerical results in Section 3.1 and 3.2, and is plotted in Figure 7. On the other hand, the supersymmetric free energy of a CFT involving only a chiral multiplet on the squashed sphere with metric (3.1) is computed by localization and found to depend only on the R-charge of the chiral multiplet, ∆ -see [17] and [18,47]. The result is compactly written as the following integral where, following [47], we have used a new definition of the squashing parameter given by the relations 27 The integral in (5.1) can be evaluated explicitly for the free chiral multiplet with ∆ = 1/2. Using this value for ∆, we have evaluated numerically the integral in (5.1) as a function of b (and thus ) and the results are plotted in Figure 7. We have used our usual normalization and the fact that for a free chiral multiplet one has C T = 6. It is clear from this plot that the non-supersymmetric and supersymmetric free energies have quite different behavior as a function of the squashing parameter . 28 Another well-known CFT for which the integral in (5.1) can be evaluated explicitly is the theory of a chiral multiplet with a holomorphic superpotential W = X 3 . This is a supersymmetric relevant deformation that leads to an interacting fixed point in the IR, known as the critical Wess-Zumino CFT -see, e.g., [5]. In this interacting theory, the chiral operator is protected and has ∆ = 2/3. We note in passing that C T for the critical Wess-Zumino model can also be evaluated by localization (or via the conformal bootstrap) and one has C T = 32 81 16 − 9 √ 3/π ≈ 4.36071 [5,47,48]. We have evaluated numerically the dependence of the free energy F on for this interacting theory, and we find that (using our normalization with C T ) this function is almost indistinguishable (within the chosen resolution) from the supersymmetric curve for the free chiral multiplet in Figure 7. We emphasize that the two functions are not the same, as is clear from (5.1), they just take very similar numerical values for a large range of the parameter .
Finally, let us mention that the supersymmetric squashed sphere partition function can also be computed for three-dimensional N = 2 SCFTs with a holographic dual -see [17,49]. Using the same conventions and notation as in Section 3.3 one finds This can be easily converted to a function of the squashing parameter used throughout this work This simple linear function of is clearly very different from the quadratic holographic result in (3.27) derived for non-supersymmetric CFTs with a gravity dual. The difference should be attributed to the fact that in the supersymmetric setup one has to turn on a particular fixed combination of background fields in the QFT which then induce certain couplings in the bulk four-dimensional supergravity that modify the supergravity on-shell action.
Observe that (5.5) reduces to the metric of S 3 α in (3.1) when β = 0. The metric in (5.5) preserves only an SU (2) subgroup of the SO(4) isometry group of the round S 3 .
The eigenvalues of the scalar Laplacian for the metric in (5.5) can be computed numerically -see [12,50]. Therefore, the free energy for a conformally coupled free scalar field can be found numerically using a method similar to the one employed in Section 3.1 for the single squashed sphere with β = 0. The details of this analysis are presented in [12] for a large range of the parameters α and β. Here we would like to point out some salient features of these results. First, following an argument completely analogous to the one in Section 2, one can show that the round sphere free energy is a local extremum of the function F(α, β), i.e., ∂ α F| α=β=0 = ∂ β F| α=β=0 = 0. The numerical results in [12] are in harmony with this general statement to a very good accuracy. In addition to that, we can extract the following leading order behavior of the free energy for small values of α and β. The numerical results are in a very good agreement with the following formula (we use C T = 3/2 for the free scalar) Clearly this is very similar to the result in (3.9) and it is natural to conjecture that (5.7) is valid for general CFTs. This conjecture is supported also by the holographic results in [12]. In particular, a two-parameter family of AdS-Taub-NUT solutions of the Einstein theory with a negative cosmological constant were numerically constructed in [12]. These solutions have the metric (5.5) as an asymptotic boundary. Using holographic renormalization, one can extract the finite part of the on-shell gravitational action and, as in Section 3.3, this is interpreted as the free energy of the dual CFT. These numerical holographic results are again in good agreement with (5.7). Finally, we would like to emphasize that if (5.7) holds for general CFTs, it implies that, for unitary theories, the round sphere free energy is a local maximum of F(α, β) in the space of squashing deformations given by (5.5).

Connection to entanglement entropy
In Section 3.4, we pointed out the similarity between our results for the dependence of F S 3 α on the squashing parameter, and those corresponding to the universal contribution to the entanglement entropy of corner regions in three-dimensional CFTs. In both situations, C T controls the leading (quadratic) correction to the most symmetric configuration (round sphere and absence of corner respectively), and in both, the curves corresponding to very different theories collapse to an almostuniversal curve when normalized by this quantity. The connection to entanglement entropy results can, in fact, be made more precise, as we explain now.
On one hand, as shown in [51,52], the universal contribution to the entanglement entropy of a region bounded by a (d − 2)-sphere is equivalent to the free energy of the corresponding theory on S d for any odd-dimensional CFT, namely On the other hand, a lot of effort has been invested recently in understanding how the entanglement entropy of planar or spherical entangling regions changes under small geometric deformations -see [29][30][31][32][33][53][54][55][56][57][58][59][60][61][62][63][64] and references therein. The result of main interest for us (proven in full generality in [33]) states that, given a general CFT in arbitrary (odd or even) dimensions, if we consider a deformed planar or spherical entangling surface, the leading correction to the entanglement entropy universal term is quadratic in the deformation (schematically) and proportional to the central charge, C T , controlling the stress-tensor two-point function, i.e., where the proportionality constant depends on the explicit details of the geometric deformation. The similarity between (5.9), (5.10) and (2.16), (2.17) is manifest, and particularly striking in the case of a deformed spherical entangling surface S d−2 and a deformed spherical background S d , in which case (5.8) implies that the leading contributions in (2.16) and (5.9) exactly coincide. Note, however, that the identity (5.8) is in principle only valid for round spheres S d−2 and S d . The fact that the leading correction to the free energy, under geometric perturbations, in odd-dimensional spheres is also controlled by C T , suggests that the relation (5.8) might somehow be extended to a relation between the free energy of CFTs in deformed spheres S d , and the entanglement entropy of deformed S d−2 entangling surfaces. This looks like an intriguing possibility worth further study. Note that a d-dimensional sphere can be deformed in more ways than a (d − 2)-dimensional one, which means that any putative relation generalizing (5.8) along this lines should no longer be a bijection.
Observe also that (5.9) and (5.10) capture, in particular, the case of an almost-smooth corner entangling region, for which the deformation parameter is naturally identified with the corner opening angle, = θ − π. 29 With this identification, the analogy established in Section 3.4 between the corner function as a function of the opening angle and the squashed-sphere as a function of the deformation parameter becomes clear.

Outlook
It is clear that in this work we have only scratched the surface in understanding how the partition function of odd-dimensional CFTs placed on deformed conformally flat manifolds depends on the deformations parameters. There are some clear avenues for further work that will shed more light on this issue.
In Sections 3.1, 3.2 and 4.1, we computed the squashed sphere free energy of free theories numerically. It will be desirable to find an analytic method to compute these partition functions. The zeta-function approach employed in [15] might be a good line of attack. It will also be interesting to extend the five-dimensional results of Section 4.1 to free Dirac fermions. To this end one will need to find the eigenvalues of the Dirac operator on the squashed sphere in (1.6) generalizing the results of [22,24]. It will of course be most welcome if we have tools to access this squashed sphere free energy for interacting non-supersymmetric CFTs. These explicit calculations will allow us to gain further insight into squashed sphere partition functions and will also serve as a test ground for the various conjectures and speculations we made in this work. In particular, it would be interesting to understand whether the speculations about a possible extension of the Hofman-Maldacena type bounds for CFTs on squashed sphere discussed in Section 3.4 have any merit. Perhaps a more modest goal will be to find a way to compute analytically the third order correction to the free energy in the small expansion, F (0), in (2.10). The numerical results for free theories in three dimensions suggest that this coefficient is a very small number and it is desirable to understand why this is the case. Finally, it would be very interesting to make the tantalizing connections with entanglement entropy for non-spherical regions that we outlined in Section 5.3 more precise, and to understand if there is a generalized version of the relation (5.8) between partition functions and entanglement entropies in odd dimensions for deformed spheres.
The tensor t a b ,c d ,ef (W ) is given by the unwieldy expression where (A, B, C, D, E) are constants. Only three of these constants are independent. This is made manifest by the following relations There is an alternative parametrization of these three independent constants which is convenient for the analysis of unitarity bounds in [35] where S d−1 was defined below (1.3). The important upshot from these explicit formulae is that the 3-point function of the energymomentum tensor of CFTs in general dimension d is completely determined by conformal symmetry up to three real constants -C T , t 2 , and t 4 . The constant C T is the same as the one appearing in the 2-point function of the energy-momentum tensor (1.3). For low values of the dimension d, there are kinematical accidents that further reduce the number of independent constants. For example, in d = 3 one finds that t 2 = 0, which in turn implies that e.g., A can be written in terms of B and C.
Finally let us point out that the 3-point function of the energy-momentum tensor on a conformally flat manifold can be expressed in terms of the correlator in flat space using the expression where Ω(x) and M ab µν are defined in (2.13) and (2.14).

B AdS-Taub-NUT solutions in Gauss-Bonnet gravity
An important part in our analysis is played by the so-called Euclidean AdS-Taub-NUT solution of four-dimensional general relativity with a negative cosmological and its generalizations to higherdimensions -see [13,27] and references therein. The metric on this space can be thought of as an S 1 fibration over a (d − 1)-dimensional base space, B, equipped with an Einstein-Kähler metric g B . The (d + 1)-dimensional metric can be written in general as where n is the so-called NUT charge, τ is the coordinate of the fibre S 1 and A is related to the Kähler form J on B through J = dA. When the dimension of the set of fixed points of the S 1 isometry is less than (d − 1), the solution describes a NUT, while when it is equal to (d − 1) it describes a Bolt. 30 The Bolt solutions are not smoothly connected to the AdS d+1 metric and will not be considered in this work. The explicit form of the function F (r) depends on the dimension d and on the particular gravitational theory of interest.
Here we focus on metrics for which the base space is the complex projective plane, i.e., B = CP d−1 2 . We can write the metric on that space in the following recursive form (we use the notation k ≡ d−1 2 ) where A CP k−1 is related to the Kähler form of CP k−1 as J CP k−1 = dA CP k−1 and ξ k and ψ k are the new coordinates of CP k with respect to CP k−1 . The normalization of (B.2) is chosen such that the Ricci tensor obeys R ab = g ab . The 1-form A CP k is given by It is clear from these formulae that it is possible to reconstruct the Kähler form and the metric of CP k from the analogous information for CP k−1 in a recursive manner starting from the k = 1 case corresponding to the usual round metric on S 2 ds 2 CP 1 = 4 dξ 2 1 + sin 2 ξ 1 cos 2 ξ 1 dψ 2 1 , A CP 1 = 4 sin 2 ξ 1 dψ 1 . (B.4) While NUT and Bolt solutions in various dimensions have been constructed in the literature for many classes of theories, our interest here is in even-dimensional Gauss-Bonnet gravities, whose Euclidean action is given by , where G is the Newton constant, sets the scale of the cosmological constant, λ GB is the Gauss-Bonnet coupling and X 4 is the Euler density defined below (4.13). The not very illuminating form of F (r) for the NUT solutions with a CP d−1 2 base space for general even-dimensional Gauss-Bonnet gravity can be found in [38] (see also references therein). For our purposes, it suffices to point out that the asymptotic metric of such solutions can be recast in the form where the boundary metric reads and where we made the coordinate redefinitions τ → ψ · n(d + 1) and r → r/ √ d + 1 and defined

B.1 Holographic free energy
To make the variational problem associated with the functional (B.5) well-defined one has to add the following boundary term [66][67][68] where h is the determinant of the induced metric h ij on the boundary, 31 K ij is the second fundamental form associated to the vector normal to ∂M , K its trace, G ij the Einstein tensor associated to h ij and J the trace of the boundary tensor This is defined as h ij = t i µ t j ν g µν , where t i µ = ∂y i /∂x µ and the y i are coordinates on ∂M. Similarly, K ij = t i µ t j ν K µν , where K µν = ∇ µ n ν , with n ν the unit vector orthogonal to the boundary.
For holographic applications, one would like to have solutions to this Gauss-Bonnet gravity theory with a finite on-shell action. To achieve this, we can rely on the standard techniques of holographic renormalization [69] and supplement the original action (B.5) with the following counterterms [27,28,70,71]

C Conformal Laplace operator on squashed spheres
In this appendix, we construct the eigenvalues of the conformal Laplacian on a general-dimension squashed sphere with the metric (1.6). Given an arbitrary Riemannian manifold, the conformal Laplace operator is defined as where R and ∇ 2 are the Ricci scalar and the usual Laplace operator respectively. The eigenvalues of D and their corresponding degeneracies on the round unit-radius sphere S d are respectively given by (see for example [7]) where ≥ 0. Now, the Laplacian on the squashed-sphere S d α can be written in terms of the round-sphere Laplacian as Using the expression for the Ricci scalar associated to (1.6), it is possible to show (see below) that the eigenvalues of D on S d α read where ≥ 0 and 0 ≤ q ≤ . Their degeneracy is in turn As a non-trivial consistency check of this formula, it is easy to show that q=0 m d, ,q = m d, , as expected. Similarly, one also finds that λ d, ,q = λ d, for α = 0. The expressions in (C.5) and (C.6) can be deduced by using the fact that the squashed sphere we are considering preserves the isometries of CP d−1 2 , which in turn is a coset space .
( (C.11) The general rule for decomposing the ( , 0, . . . , 0) representation of SO(d+1) is the following. First, one has to remove one of the Dynkin labels, since the rank of SU ( d+1 2 ) is one lower than the rank of SO(d + 1). Then one writes all possible decompositions of the non-negative integer into a pair of non-negative integers p 1 and p 2 such that = p 1 + p 2 . Finally, one writes all possible Dynkin labels with p 1 as the first entry and p 2 as the last with zeroes in between. One can then identify the integer q in equations (C.5) and (C.6) with p 1 .

D Numerical regularization
In this appendix, we give more details on the numerical techniques utilized in the free field theory calculations throughout the paper. As already mentioned, our regularization scheme is based on the same technique as the one used in [10] (see also [12]).
The free energy of the free scalar and fermion in arbitrary dimensions can be written as where D is the conformal Laplace operator or the Dirac operator respectively, f = 1 for fermions, f = 2 for scalar fields, and Λ is an energy cutoff. The free energy will, in general, be a UV divergent quantity which needs to be regularized. To this end, we use a heat-kernel type regulator [10,73] where we have denoted the eigenvalues of D by λ i . This expression yields the determinant for modes whose energies are less than a "soft" cutoff Λ, while cutting off the sum exponentially above this value. In particular, for λ Λ f , one finds The integral can now be split into two pieces, one with low energy modes (IR) and another with high energy modes (UV) (D.6) Here, Γ(a, z) is the incomplete Euler Gamma function, δ is an arbitrary positive real number that we can change to get a better convergence, and m i is the multiplicity of the eigenvalue λ i . The sum over the IR modes converges for large λ i , and can therefore be done numerically if the maximum number of eigenvalues is chosen correctly. The divergences are all contained in det U V . These have to be controlled and subtracted. For the case of the free scalar field in three dimensions, one can apply the Euler-MacLaurin formula  to the infinite sums to obtain these analytically [10]. In more general cases, the Euler-MacLaurin formula does not give an estimate of the infinite sum. However, for a given value of t, this sum can be obtained numerically, as it converges fast enough if we take into account enough eigenvalues. The maximum number of eigenvalues that we will sum over will be called n max . The next step in the algorithm consists in evaluating this converging sum for different values of t, from a starting value t init , to the final value δ, with step size ∆t. With this sum evaluated for a large range of t's, it is possible to do a numerical fit to the "data points". This gives the dependence of the sum on the energy cutoff. In order to know which function we need to use for the fit, it is necessary to know the behavior of the integrand as a function of t. Since t ∼ Λ −2 , we know that it should have the same behavior as the on-shell action in terms of the cutoff powers -see e.g., (3.13). Therefore, we perform a fit with the function fit(t) = A 1 t (d+1)/2 + A 2 t d/2 + A 3 t (d−1)/2 + · · · + A d+1 t 1/2 . (D.7) We expect all the A i with i odd to vanish -i.e., to be very small numbers in our numerical calculations, so that the leading term of the fit contains the volume divergence, the subleading one the first curvature divergence etc. In each case, the coefficients that should vanish are used as a precision check of our numerical work. If these are really small, the fit can be trusted. To give a flavor of our results for the coefficients in (D.7) in Table 1, we present explicitly the values for A i for different values of α in d = 3 for the free fermion field. Notice that for the free fermion the eigenvalues in (3.19) are labelled by the integers, and thus the sum over i splits into 2 sums, one over q and another over n in (3.19). The multiplicity in this case is equal to n. To obtain the values in Table 1, we take for det U V , n max = 3101, and we start our fit from t init = 10 −6 with step size ∆t = 10 −6 . For the parameter δ we choose δ = 10 −3 . After the fit function in (D.7) is obtained in this way, we have good control over the divergences in det U V which we then subtract and, after that, evaluate the integral in (D.6). The result is then used in (D.5) and (D.1) to obtain the finite regularized value of the free energy.