Exact holographic RG flows and the $A_{1}\times A_{1}$ Toda chain

We construct analytic solutions of Einstein gravity coupled to a dilaton field with a potential given by a sum of two exponentials, by rewriting the equations of motion in terms of an integrable Toda chain. These solutions can be interpreted as domain walls interpolating between different asymptotics, and as such they can have interesting applications in holography. In some cases, we can construct a solution which interpolates between an AdS fixed point in the UV limit and a hyperscaling violating boundary in the IR region. We also find analytic black brane solutions at finite temperature. We discuss the properties of the solutions and the interpretation in terms of RG flow.


Introduction
The notion of the renormalization group has dominated our thinking about quantum field theories and statistical systems since its elaboration by K. Wilson. It has been a paradigmchanging idea, providing a unified and systematic picture to understand the dynamics of systems with many degrees of freedoms. Quantum field theories are now not seen as isolated items, but as classes connected by the RG flow. The method is very powerful, however in practice one can usually only determine the structure of the flow in a perturbative expansion around a fixed point that is weakly coupled. Luckily there are many interesting theories for which the program can be implemented, including of course Yang-Mills and QCD that are asymptotically free, the Wilson-Fisher fixed point of λφ 4 in dimension 4 − , etc., but strongly-coupled examples are few.
The holographic correspondence gives a description of a class of strongly-coupled field theories in terms of a dual weakly-coupled gravitational description. The RG flow in this description is geometrized, and corresponds to a gravitational solution with particular asymptotic properties; the holographic direction corresponds to the energy scale, and so the Hamiltonian evolution in this direction can be put in correspondence with the evolution of the system under the change of the RG scale [1]. These solutions are often described as "domain wall" solutions, as they interpolate between two asymptotic regions, each of them being a solution on its own and corresponding to a given theory [2]- [4].
This aspect of the holographic correspondence has been intensely explored, but some issues have so far eluded a complete resolution; for instance, the precise relation between the holographic and the Wilsonian scheme [5,6], and the fact that Einstein equations are second order while RG equations are first-order; even though the equations can be cast in first-order form using the Hamilton-Jacobi formalism, there seems still to be a mismatch in that the couplings of the field theory are promoted to fluctuating fields in the gravity description; recently an attempt to solve this problem has been made [7] with the notion of a "quantum RG flow", that arises from the classical one after integrating out the doubletrace operators. This flow may be more complex than the one typically arising in QFT, and these possibilities have only been considered in the last few years [8].
The matching of RG equations with the gravity equations has been precisely formulated only for the solutions that have the Poincaré invariance at the boundary, and correspond to the vacuum of the dual field theory. In these cases, one can show that the equations of motion can be expressed in terms of a superpotential, which is determined by the scalar potential (up to the choice of some integration constants), and is in one-to-one correspondence with the beta function of the dual theory. But it is not yet known how to extend this procedure to non-Poincaré invariant states, in particular finite temperature/density.
In this paper we will not be concerned with these conceptual issues. Rather, in the spirit of "bottom-up" holography, we remark that there are some relatively simple models of gravity coupled to a scalar field, that allow for interesting examples of RG flows that can be found analytically. This is due to the fact that the Einstein equations, with the Ansatz that corresponds to domain-walls solutions, reduce to dynamical equations that are completely integrable (they can be reduced to the equations of a Toda chain). It is -2 -

JHEP05(2019)117
remarkable that such analytic domain-wall solutions can be found not only for the vacuum but also at finite temperature. The Toda equations arise when the potential is a sum of exponentials, when the coefficients in the exponent satisfy certain relations. The case of a single exponential term had been solved in [9]. The Chamblin-Reall black brane solution has found applications to the study of the dynamics of a non-conformal strongly coupled plasma in the hydro regime and beyond [10,11]. However the single exponential case has some limitations: the potential is monotonous, so it does not have a minimum but only run-away solutions; the RG flow in this case does not start from a UV fixed point, and in fact the dual theory is not well-behaved in the UV. Moreover, the single exponential has a definite sign. This precludes the possibility of studying cases in which the potential changes sign along the solution. The general analysis of RG flow with a single scalar, performed in [8], only deals with potentials that are definite negative; the general case of a potential with zeros is less explored.
In this paper we consider the next simplest case, with a sum of two exponentials. Integrability leaves one of the two exponents unfixed, so we have one free parameter, like in the case of a single exponential. This case already allows to overcome the limitations just mentioned. We will leave the consideration of more complicated potentials for future investigations. Qualitatively, the properties of the potential is that it has a minimum with V min < 0, it goes from being negative and vanishing at φ → −∞ to positive and diverging at φ → +∞.
It turns out that for our choice of the potential, the integrable Toda chain is associated to the Lie algebra is A 1 × A 1 [12,13], so it is a particularly simple case and the solutions can be given explicitly in terms of elementary functions.
The solutions we find depend on a certain number of parameters. One parameter distinguishes between vacuum and nonvacuum solutions. The vacuum solutions are Poincaré invariant. Turning on the parameter gives a deformation to non-vacuum solutions, in which Poincaré symmetry is broken, and horizons can be formed. One more constant is the "energy" of the solutions (more exactly the energy in the associated Toda chain description) and determines the type of solutions; there are four general classes of solutions -depending on the type of functions that appear, we will call them the sinh-class, sin-class, linear class and cosh-class of solutions (see eqs. (2.33)). We will consider mainly the sinh-class. Two other parameters determine the position of singularities, whose presence requires to split the domain on the radial variable in three regions; correspondingly we have "left", "middle" and "right" solutions. In the left and middle solutions the dilaton interpolates between +∞ and −∞, whereas the right solution is bouncing: the dilaton starts at −∞, goes to a maximal value and then goes back to −∞. It should be noted that the singularities are not just coordinate singularities; the scalar curvature usually diverges at the end points of the branches, except for the right end of the middle branch.
Coming to the domain wall coordinates we explore the solutions in the holographic framework. We define the energy scale A= e A , where A is the scale factor of the domain wall metric and the running coupling as λ = e φ through the dilaton φ. It is interesting that, in spite of the fact that the dilaton has a similar behaviour in the left and middle solutions, the scale factor A has rather different behaviour on these solutions, namely, on the left -3 -

JHEP05(2019)117
solution when the dilaton varies from −∞ to a special value φ s , the energy scale increases, but after the dilaton passes the special value φ s , the scale factor starts to decrease down. This non-monotonic behavior of A precludes the possibility of a holographic interpretation of this branch of solutions. For the middle solution, when the dilaton increases from −∞ to ∞, the scale factor decreases from large positive values to large negative values. This behaviour corresponds to the running coupling in an asymptotically UV free theory, so it is of interest from the point of view of possible applications in QCD. In the right solution the scale factor is increasing from minus infinity to plus infinity, but the dilaton bounces back after reaching a maximal value φ max , as already mentioned. These behaviours are illustrated in figure 8.
There is also a special case when the points of singularities coincide. In this case the singularities are cancelled and the solution has a smooth AdS boundary in the UV limit and a hyperscaling violating boundary in the IR region. The dilaton supporting this geometry runs from a constant value in the UV to −∞ in the IR.
As we mentioned before there are other solutions, in particular this one from the linear class describes the opposite flow, from hyperscaling-violating in the UV to AdS in the IR.
When we turn to non-vacuum solutions, the deformation can give solutions with a horizon or without. For the cases having horizons, the parameter characterizing the deformation is related to the temperature. It happens that the black brane solutions can be constructed only from those vacuum solutions which are defined for u → ±∞. For these solutions the dilaton potential evaluated on-shell is bounded from above and therefore the solutions obey Gubser's criterion [22]. In particular, for the solutions defined for u → +∞ we change the scaling properties of these solutions in the IR regime varying the temperature, meanwhile the UV behavior does not change.
In accordance with Gubser's criterion, we find finite temperature solutions only for the "regular" vacuum flow. In particular, we do not have any black brane solution with AdS UV asymptotics, as this is a singular flow. The finite temperature generalization of the vacuum solution with coinciding singularities yields to be just a AdS-Schwarzschild black brane, i.e. it does not flow anymore.
The interpretation of the solutions in terms of RG flows is clarified by using the first order variables X, Y [14] that have a direct relation to the superpotential and the beta function, when a dual theory interpretation is possible. We discuss how our solutions describe flows between different attractor points in the phase space, i.e. the (φ, X) plane, but in order to recover all possible flows we also have to include other classes of solutions, namely the linear and cosh-solutions. This is illustrated in figures 9 and 10. In figure 16 we show possible flows at non-zero temperature.
The paper is organized as follows. In section 2 we describe the holographic gravity model, the ansatz for the metric and the dilaton that leads to a special mechanical model, that can be explicitly integrated. Here (section 2.5) we also introduce some general relations for the holographic RG flow. In section 3 we describe in details the vacuum solutions and give their interpretation in terms of the holographic RG flow. In section 4 we study nonvacuum solutions, derive a black brane and study the black brane solution as a holographic RG flow at finite temperature. In the appendix, we collect information about curvature Figure 1. The behaviour of the potential V (φ) for C 1 < 0, C 2 > 0.
invariants for our background in section A, some formula about the dilaton field in section B and details about the superpotential for vacuum case in section C.

The setup
We consider a holographic model with gravity coupled to a dilaton field, and the dilaton potential is taken as a sum of two exponential functions. As explained in the introduction, the choice of the dilaton potential is motivated by studies of models with only one exponential function in the potential [10,14].

The holographic gravity model
The holographic model is governed by an action of the form where C i , k i with i = 1, 2 are constants. We choose the constants C 1 and C 2 as follows In this case the potential has a minimum and regions of positive and negative sign, as shown in figure 1. The model (2.1)-(2.2) is defined on a 5-dim manifold M equipped with the metric where A = A(u), B = B(u) and C = C(u) are some smooth functions. We make an ansatz in which these functions, as well as the dilaton, depend only on the u-coordinate. The equations of motion which follow from the action (2.1) read

JHEP05(2019)117
The equation for the scalar field is Taking into account the relations (A.2)-(A.6), (B.1), the Lagrangian can be reduced to the following form, up to total derivative terms: where we denote˙≡ d du . For convenience we redefine the functions as and the so-called lapse function In what follows we use the harmonic gauge, i.e.
The explicit form of E.O.M. is presented in appendix A.1 in the harmonic gauge C = A + 3B.

Mechanical model
The Lagrangian (2.7) with the help of (2.8)-(2.11) can be taken to the following form (see [21] and refs. therein) where the minisupermetric G ij on the target space M reads (2.14) The corresponding energy constraint is given by However, the system (2.12)-(2.15) is still difficult to work with, because of the degrees of freedom are coupled in the potential. The idea of the next section is to separate the equations of motion rotating the system [18,19].

Integration of the mechanical model
To perform the following calculations it is useful to present the mechanical Lagrangian (2.12) in the form where V and W are some vectors on the target space (we suppose that the original basis is (e 1 , e 2 , e 3 )), the brackets denote a scalar product on the target space M with the metric (2.14).
Then the components of V and W are defined by The method we will use applies generally to systems with interactions of the form (2.16), with arbitrarily many exponentials [18][19][20]. It turns out that the system is integrable if the vectors V and W can be identified with the root vectors of a Lie algebra (for more details see [12,13,18,19]). Then the system becomes a Toda chain. In our case, we can see that the scalar products of the vectors are: These correspond to the roots of a Lie algebra only if V and W are orthogonal, in which case we have the A 1 × A 1 Toda model. The integrability requirement puts a restriction on the dilaton couplings: (2.20) We also suppose that V is a timelike vector and W is a spacelike one (notice that the symmetry between k 1 and k 2 is broken by the choice of the signs of the coefficients in the potential). We rename and we have the condition We can find a basis of orthonormal vectors with respect to the metric G ij :

JHEP05(2019)117
The new basis is related to the old one by a Lorentz transformation: 25) and the coordinates transform as The Lagrangian and the energy constraint (2.12)-(2.15) in the new basis take the form From (2.27)-(2.28) we see that mechanical variables are decoupled, so the equations of motion following from the Lagrangian (2.27) arë where s = 1, 2 and we introduce a notation for the scalar product The eqs. (2.29) with s = 1, 2 are two decoupled Liouville equations. The solution can be given explicitly: with the functions where u 0s , E s are constants of integration, s = 1, 2. The E s are the conserved energies of the two decoupled Liouville modes. In our case, with C 1 < 0 and C 2 > 0, we have η ss C s > 0, so the relevant solutions are given by the first three lines. We will analyze in detail the sinh-solutions, which are the most interesting from the RG point of view. The solutions to eq. (2.30) is with constants of integration p 3 , q 3 . Having solved the system, we can go back to the original variables applying the transformations (2.26) with components and obtain with the parameters α i , β i defined using (2.25) The parameters α i satisfy the following conditions The constants E 1 , E 2 and α i are related by the constraint (2.28) which reads with E 1 < 0, E 2 > 0 and the value of E 1 is bounded as |E 1 | = E 2 + 2(α 1 ) 2 3 . Then the functions F 1 and F 2 in (2.43)-(2.46) are given by where 0 < k < 4/3, C 1 and C 2 are given by (2.2)-(2.3). The generic form of the metric that solves EOM following from (2.1) reads where y = (y 1 , y 2 , y 3 ).
The dilaton potential evaluated on the solution becomes (2.51) We note that we have two more solutions for our choice of the potential (2.2) with C 1 < 0 and C 2 > 0 governed by F 1 and F 2 from the second and third branches of (2.33). The solutions differ by the range of integrating constants, i.e. for E 1 > 0 and E 2 < 0 (opposite signes to (2.48)-(2.49)), we have ) , (2.52) while taking E 1 = E 2 = 0 one has

Solutions as RG flows
Let's briefly discuss the gravity solutions as holographic RG flows. It is useful to come to so-called domain wall coordinates. In this coordinates a general form of the non-vacuum solutions is which covers the vacuum case with f (w) = 1. Both temperature and vacuum solutions are characterized by the scale factor e A(w) , that measures the field theory energy scale, the blackening function f (w) and by a scalar field profile φ(w) which is interpreted as the running coupling. If we define we get the system of first order differential equations For the vacuum case we get

The metric and the dilaton for vacuum exact solutions
The vacuum solutions are those that preserve the Poincaré invariance at the boundary. As can be seen from (2.50), this requires that α 1 = α 2 = 0, so in this case the two Liouville energies have to be: Owing to (3.1) the metric (2.50) takes the form and the dilaton is given by In the vacuum case, due to (3.1) the following relation holds: From the form of the functions F 1 and F 2 given by (2.48)-(2.49) we can see that the solution will have coordinate singularities at the points u 01 , u 02 . We have to consider three coordinate charts. Let us take u 02 < u 01 . Then the charts are

3-branch solution
The degenerate case with u 01 = u 02 = u 0 requires only two charts: In figure 2 we plot the solutions for the dilaton φ (3.3) for each choice of the coordinates (3.5)-(3.7).
From figure 2 we see that for u → u 01 ± the dilaton tends to −∞, while for u → u 02 ± the dilaton asymptotes to +∞. It can be verified by direct calculations.

Asymptotics of the metric and the dilaton
We present here the asymptotics of the dilaton (3.3) and the metric (3.2) with (3.5)-(3.7) near the boundaries.
• The left solution with u < u 02 at u → −∞: where we use a new coordinate z ∼ 4+3k 3µ 1 e at u → u 02 − : with the radial coordinate defined by From the relation for the scalar curvature (A.21) with u → u 02 − it can be seen that the solution has a non-removable singularity at the point u 02 .
• The middle solution with u 02 < u < u 01 at u → u 02 + the asymptotics of the dilaton and the metric are the same as in (3.12)-(3.13). As in the case of the left solution, one can see from the scalar curvature (A.21) that the middle solution has also a singularity at the point u 02 ; -13 -

JHEP05(2019)117
at u → u 01 − : where the conformal coordinate is Notice that as u → u 01 , the conformal coordinate z → 0 if 0 < k < 2/3 and z → ∞ if 2/3 < k < 4/3. These asymptotics are the same as for the solutions in a single exponential potential [9]. It is worth nothing the scalar curvature of the middle solution (A.23) has a regular behaviour at u 01 .
• The right solution with u > u 01 at u → u 01 + the asymptotics are as in (3.15)-(3.16).
at u → +∞: The asymptotics are the same as in (3.10)-(3.11). Even though the dilaton goes to −∞, these are different than the asymptotics for a single exponential. The scalar curvature (A.25) goes to +∞ with u → +∞. Now let us turn the discussion to the dilaton potential, that can be written on solutions as (3.20) In figure 3 we plot the dilaton potential on the solutions for φ. From figure 3C one can see that for the solutions defined for u > u 01 there is a turning point V s = V (φ) with φ = φ s , where the potential stops and goes back to zero.
One can find the stop points of the potentials analytically. The stop point u s can be calculated as follows. The dilaton should obey φ = 0, (3.21) this gives rise to Figure 3. The dilaton potential plotted on solutions for φ. A) The potential in the left-solution varies from V (φ 1 ) = 0 at φ 1 = −∞, reaches its minimum and then goes to V (φ 2 ) = +∞ at φ 2 (u 02 − ) = +∞ (thick lines). B) The potential on the middle solution varies from V (φ 2 ) = +∞ at φ 2 = φ(u 02 + ) = +∞ then reaches the minimal value V min < 0 and goes to V (φ 1 ) = 0 Note, that the value of the potential at the turning point V s = V (φ s ) (where φ s = φ(u s )) doesn't coincide with the extremal point of the potential. The extremal point of V (φ) can also be computed analytically. Using and taking into account that C 1 C 2 < 0, we can find We show the behaviour of the potential as a function of φ and the dilaton as a function of u on the same plot in figure 4. In this picture we draw the dependences for all solutions.

Special case u 01 = u 02 , solutions with AdS boundary
Let us see the features of the solution given by (3.2)-(3.3) with (3.9) and u 01 = u 02 = u 0 . In figure 5 the behavior of the dilaton (3.3) with u 01 = u 02 = u 0 is shown. From this picture one can see that the dilaton tends to −∞ as u → ±∞.
As for the previous solutions we present the asymptotics in the conformal coordinates.
• In the limit with u → ±∞    where the conformal radial coordinate is given by z ∼ ∓ 4+3k 3µ 1 e ∓ 3µ 1 u 4+3k . So we come to the same asymptotics as for the left (3.10) and right (3.18) solutions with u → ±∞ for the 3-branch case. The scalar curvature has also a common behaviour with the right and left solutions, namely it goes to infinity (A.29) with u → +∞.
• For u → u 0 one gets has the following form of the metric (3.2) where the conformal "radial" coordinate is defined as z = 4(u − u 0 ) 1/4 and z → 0 as u → u 0 . In (3.28) one can easily recognize the 5d AdS metric supported by the constant dilaton which coincides with the minimum of the potential (3.25). As expected the scalar curvature of this solution with u → u 0 has a constant value (A.28).
We note that the figure 5 is an agreement with the calculations of the asymptotics for the dilaton. The potential of the dilaton as a function of φ with u 01 = u 02 = u 0 is presented in figure 5B. From this plot we observe the existence of the turning point of the two-branch solution. The equation for the stop point (3.23) in this case becomes and has a solution u s = u 0 . (3.31) As already observed, in this case we find that the stop point coincides with the extremal point of the potential, φ s = φ c = φ(u 0 ).

Details of RG flow for vacuum solutions
The domain wall form (2.57) of the vacuum solutions looks To come to (3.32) we use the change of variables for We can represent the scale factor of the domain wall (3.32) as follows so the energy scale A is (3.35) The running coupling is defined through the dilaton (3.3) reads (3.36) The function X (2.59) is represented as follows with λ given by (3.36), A -by (3.34). It is useful to explore the evolution of the rescaled β-function X (3.37) as a function of log of the running coupling (3.36), which is the dilaton φ. In figure 6 we present the parametric plots of the behaviour of X (3.37) using solutions for φ (2.46) and A (3.34) with the dependence on the parameter u (3.5)-(3.7). In these plots we fix the shape of the potential putting C 1 = −C 2 = −2 and k = 1 while we vary the constants |E 1 | = |E 2 | (labeled as E in figure 6). From picture 6A and 6C we see that the holographic β-functions at zero temperature constructed on the left u < u 02 and right solutions u > u 01 can take both negative and positive values. As for the case 6B, which corresponds to the middle solution, X is always negative. It is also of our interest to see the function X(φ) plotted on the solutions with u 01 = u 02 = u 0 , particularly for u 0 = 0. In figure 7 we show the function X(φ) on the dilaton solution (3.3) with u 01 = u 02 . We see that the behaviour of the function X(φ) is the same for all values of E 2 . From figures 6-7 we see that X are regular except those plotted on the left solutions on figure 6A at some points φ s where X(φ s ) takes an infinite value. From . The X(φ) function for the dilaton solution with u 01 = u 02 = 0, the potential fixed as Eq. (3.38) defines the singular point u s , which is related to φ s . The singular point φ s coincides with the point where A = 0, and its position depends on E, k, u 01 and u 02 , see figure 8A. From this picture we see that for the left solution the dilaton varies from −∞ to a special value φ s , the scale factor is non-monotonic function increases, but after the dilaton passes this special value φ s , the scale factor starts to decrease. From figure 8B we observe that the scale factor A has a monotonic behaviour, decreasing with respect to dilaton running from −∞ to +∞. Figure 8C it is demonstrated by orange curves that the scale factor of the right solution decreases from +∞ to −∞ for all values of the dilaton, which runs from −∞ to some constant value and then goes back to −∞. We also present that the behavior of the scale factor A on φ for the solution with u 01 = u 02 = 0 by the brown curve in figure 8C. We see that the scale factor for this solution starts to decrease from +∞ with some constant value of the dilaton, then, passing some value of A, both the scale factor and the dilaton tend to −∞.
In the section 2.5 it was already said that the system (2.61)-(2.62) on the vacuum solutions reduces to eq. (2.63), which is quite simple and one can treat it in the general form. We remind that the dilaton potential is given by (2.2) and in the cases presented at the plots below with k = k 1 and k 2 = 16 9k 1 (a constraint from the solution) and C 1 = −2, C 2 = 2. The values of the dilaton coupling constant have the following restriction 0 < k < 4/3.
Taking into account (3.39) we have potential (3.39). We also impose figures 6 and 7 on figure 9A and see that they can partially cover the plot. In figure 9A the red lines correspond to the left solution, see figure 6A, the green lines are those from the middle solution, see figure 6B, the orange ones correspond to the right solution, see figure 6C. In figure 9A we observe that the function X corresponding to the right solution (the orange curves) interpolates between X c2 and 1. The X for the middle solution (the green curves) interpolates between X c1 and X c2 in figure 9. The dark red curves start at 1 and go to +∞ as φ → φ s − 0, some of these lines have local minimum and are located very close to the brown curve corresponding to the vacuum flow for the two branch solution u 01 = u 02 . The flow with local maximum in figure 9A corresponds to our solutions with small difference between u 01 and u 02 . The darker red lines in the right bottom part of 9A start at X c1 and go to −∞ when φ → φ s +0.  figure 9B we present the stream of (2.63) by blue lines, X(φ) on φ plotted using (2.52)-(2.53) by grey lines and X(φ) related to (2.54) by the dark red dashed lines, correspondingly. We see that these lines partially fit the solutions to equation (2.63). We show separately the behaviour of X(φ) with (2.54) in figure 9C. From figures 9A and 9B we see that X(φ) corresponding to the linear solutions can be considered as a boundary between X(φ) corresponding to "sin"solutions and X(φ) goverened by the left sinh-solutions (i.e. sinh-solutions with u < u 02 ).
However, we still observe in figures 9B and 9C that some regions of the figure are not covered. If we look at the equation (2.63) it is evidently that it is invariant under the -20 - In figure 9D we draw X(φ) as a solution to (2.63) (blue lines) and with the help of the analytical solution related to (2.55)-(2.56) (the curves of other colors), that perfectly complete the necessary part of the stream.
Let's look at figures 9A and 9B more closely. On these pictures we see a special point P , which is a stationary point of our potential V (φ) = 0. We zoom the scale near the point P in figure 10 to see the behaviour of the vacuum X(φ) functions near P in details. As in the previous figures we show the stream of eq. (2.63) by blue curves and by the other color curves the functions X(φ) plotted using analytic solutions for A and φ.
It worth to be noted that in figures 9 and 10 the arrows always show the direction of decreasing scale factor A, see figure 8, that corresponds to the flow from UV fixed point to IR.

The running coupling λ = e φ on the energy scale
It is important to know the behaviour of the running coupling λ = e φ on the energy scale A= exp A. One can trace this using the analytic solutions for φ and A. For the vacuum case both of them depend on the constant of integration E 1 , E 2 , u 01 , u 02 and on the shape of the potential, which is defined by C 1 , C 2 and k.
Let us see how the energy scale A depends on u. In in figure 11A we observe that A is non-monotonic function for the left solutions. First it increases, but near u 02 it decreases that can be clearly seen in the zoomed picture figure 11B. This non-monotonic behaviour is also read from figure 8 for the behaviour of the energy A on the dilaton for the left branch.
In figure 12 the dependence of the coupling constant λ = exp φ on the scale A is presented for the zero-temperature solutions with (3.5)-(3.7) and figure 12C with (3.9). We recall that on the left and middle solutions the dilaton interpolates between ±∞ and for the right solution it starts at −∞ goes to the maximal value φ max and then goes back to −∞ (we have observed this behaviour already in the previous sections). In spite of that the dilaton has the similar behaviour on the left and middle solutions, the scale factor has rather different behaviour on these solutions, namely, on the middle solution the scale factor monotonically decreases from large positive values to zero as the dilaton runs from +∞ to −∞. The right solution is a bouncing solution with the decreasing scale factor. These behaviours of the scale factors and the dilaton are reflected in the plots for the dependence of the running coupling on the energy scale.
We can summarize the results for the running coupling in the following form • as expected for solutions with u < u 02 , see figure 12A, where the dilaton tends to −∞ in the IR region and we have IR-free theory, while in the UV region the effective coupling λ → +∞.
• In figure 12B we see that the dependence of λ on A plotted on the middle solutions with u 02 < u < u 01 mimics the QCD behavior.
• The running coupling plotted as a function of the energy scale at figure 12C for the solutions with (3.7) shows that it is UV free as well as IR free theory, i.e. the running coupling has the form of the hill.

The metric and the dilaton
The metric and the dilaton solutions to the model (2.1) in the non-vacuum case are where the functions F 1 and F 2 are given by (2.48)-(2.49) as before. Just as in the vacuum case, we need to separate the solutions in the branches (3.5)-(3.9) with respect to values of u 01 and u 02 . We note that the factors of α 1 in the metric (4.1) break the Poincaré symmetry. However, below we will see that the presence of this parameter allows us to define a horizon and to construct black branes. We also recall that E 1 and E 2 must obey the constraint The condition (4.3) allows to tune parameters thus we have two additional regimes respectively. This leads to new dynamics of the dilaton (4.2). Particularly, the dilaton can be constant, that is inapplicable for the vacuum case where it always holds µ 2 = 4 3k µ 1 .  We illustrate the behaviour of the dilaton solution for the non-vacuum case in figure 13 for branches (3.5)-(3.7) with u 01 = u 02 . As it was expected the condition (4.3) changes the behaviour of the dilaton. In figure 13B we plot the dilaton keeping the same shape of the potential and the value of α 1 as for figure 13A, but the values of the parameter E 2 are changed with respect to (4.5). In figure 13C we again save the form of the potential, use the same value of E 2 as in 13A, take bigger α 1 and show that the opposite sign of this parameter doesn't change the asymptotics of the dilaton. From figure 13C one observes that for u → ±∞ φ can tend to +∞, while for the vacuum case φ → −∞ as u → ±∞. In figure 13C we also see the behaviour of the dilaton doesn't depend on the sign of the parameter α.
In figure 14 we plot the dependences of the dilaton solution with u 01 = u 02 = u 0 on u (3.9). As for the 3-branch solutions from figure 14B we see that the dilaton can be constant with the appropriate choice of parameters agreed with (4.3) and can change its asymptotics from −∞ to +∞ at u → ±∞, see figure 14C.
In figure 15 we draw the potential V as a function of φ and the dependence u on φ. The functions u(φ) are different on the different branches. The different values of u corresponding to the same φ are indicated by points at the vertical lines. We see that  the function u(φ) for |α 1 | < α 1 cr is double-valued on the right branch, and for |α 1 | > α 1 cr it is double-valued on the left branch, and for |α 1 | = α 1 cr the both functions one-to-one functions for all branches, and φ < φ 0 for the right branch and φ > φ 0 at the left branch.
The dilaton potential which is plotted in figure 15 on the right solutions for the dilaton is bounded above like for the vacuum right solutions.

The black brane solutions
Now we are going to find a black brane representation of the solution (4.1)-(4.2) defined for (3.7) with u 01 = u 02 . The metric (4.1) can be rewritten in the following form with the constant C given by the function X (u) written as  where we took into account the relations for µ 1 and µ 2 (2.48)-(2.49). We see that for 0 < k < 4/3 one has κ < 0.
To have the black brane solutions we need to remove a conical singularity in the metric (4.6) (since with α 1 < 0 we get zero in front of dt 2 ). Taking e 4 3 α 1 u = ρ, t = iτ, (4.10) so the metric reads where X (u) → 1 for ρ → 0 as u → ∞. Therefore, there is no conic singularity if the following constraint is satisfied We also fix the periodicity 4 3C 3/2 α 1 β = 2π. (4.13) Plugging (4.9) in (4.12) we come to the condition to the parameters , E 2 = 6k 2 (α 1 ) 2 16 − 9k 2 , (4.14) that corresponds to µ 1 = µ 2 = µ with Therefore, we get the black brane, if (4.12) is satisfied, and the temperature is Here C is taken with the constraint (4.15), see (4.19) below. Under the condition (4.12) the black brane metric has the form where C and X are given by  with the horizon located at u = +∞ and the near-horizon expansion of X (u) is We note that the boundary is at u 01 and we fixed u 01 = 0 to have f = 1 at this boundary. One can check that null geodesics we have the correct behaviour. Null geodesics imply ds 2 = 0, (4.21) i.e. for the light moving in the radial direction This calculation confirms that we have the horizon at u = +∞. The scalar curvature and the Kretschmann scalar near horizon u → +∞ are We note that with respect to the constraint to absence of the conic singularity (4.11) both the scalar curvature (4.24) and Kretschmann scalar (4.25) tend to zero with u → +∞.

The Gubser bound
The dilaton supporting the geometry (4.17) reads and takes the constant value near horizon Now one can check if the Gubser's bound [22] for asymptotically non-AdS solutions holds where V UV is the value of V (φ) at ultraviolet fixed point. Since V UV = 0 with the dilaton φ → −∞ at the UV point the constraint (4.30) comes to be (4.28).
In the UV limit, i.e. near u 01 = 0, the solutions turns to have the asymptotics as the Chamblin-Reall solution governed by the single exponential potential We note that we can construct a black brane background for the left solutions (3.5) with u 01 = u 02 assuming that the horizon is located at u = −∞, and the parameter α 1 is positive.

Special case u 01 = u 02 , AdS black brane
Now we turn to the special case of the non-vacuum solutions with u 01 = u 02 = 0 with u > 0. The construction of a black brane metric is the same as presented before for solutions with u 01 = u 02 with the horizon located at u = +∞. The metric (4.17) has the form where we took into account µ = − 4 3 α 1 and the constant C reads Due to the constraint µ 1 = µ 2 and u 01 = u 02 the dilaton (4.2) becomes constant The curvature of the metric (4.33) is negative and reads

JHEP05(2019)117
Doing the change of coordinates one gets the usual form for the 5d AdS black brane For the dilaton potential we have the saturation of the Gubser's bound (4.30) that is in agreement with the suggestion from [22], since the solution (4.33)-(4.35) is antide Sitter black brane.
One can summarize our studies on non-vacuum solutions as follows. For u → ±∞ the scalar curvature and the dilaton for the right and left solutions can be constant for arbitrary value of the temperature which at the same time defines the constant E 2 (4.5). However, the left solutions have a special point u 02 at which the scalar curvature has a non-removable singularity, while the scalar curvature of the right solutions is regular at its special point u 01 . Finally, from (4.29) and (4.40) we see that the dilaton potential calculated on-shell is bounded only for the right solutions and special solutions with u 01 = u 02 . Therefore, these solutions can satisfy Gubser's criterion [22].  The finite temperature case is described by the additional variable Y defined through the blackening function (2.60). Here we deal with the system (2.61)-(2.62), which seems to be rather complicated comparing to the zero-T case and one has to apply analytic solutions for φ and A with the coordinate u > u 01 to show the RG flow.

JHEP05(2019)117
To see the behaviour of the RG flow at finite temperature it is useful to plot X and Y as functions of φ for the analytical solution. We remind the horizon of the black brane at u → +∞ in this case is defined for α 1 < 0. In figure 16 we draw X and Y as functions of φ for different values of the negative α 1 . From figure 16A we see that the behaviour of X is changed by α 1 , so X becomes to be negative in the finite-T case. As for the Y function one can see from figure 16B that it is positive. In figure 16A we observe that the X function has stop points at some values of φ. These values of φ correspond to the asymptotics of the dilaton at the horizon φ h (4.27). One can show that X given by Then the value of X at the horizon is We see that as µ → 0. As for the Y -function (2.60), it takes infinite values at the horizon that can be read from figure 16B.
In figure 17 we present the dependences of X and Y on the energy scale A. In figure 17A we again observe stop points of X at some A (4.47) with u → +∞.  Now let us look what happens with the behaviour of the running coupling on the energy scale at finite temperature. As for the vacuum case we have the parametric dependence on E 1 , E 2 , α 1 , the position of poles u 01 and u 02 . We note that E 1 , E 2 and k are not independent. We have also seen before the temperature is related to the parameter α 1 (4.12). As in the previous section, section 3.3.2, to have an insight to possible behavior of the running coupling on the constructed solutions, we start with presenting the behavior of the energy scale A as a function of the coordinate u with (3.5)-(3.7) and then incorporated found in section 4, see figure 13, the behavior of the dilaton as a function of u.
In figure 18 we see that A tends to some constant value at u → ∞ that can be supported by (4.47).
In figure 19 we present the parametric dependence of the running coupling λ as a function of the energy scale A for the black brane solutions with u 01 = 0, u 02 = −1 and  different values of the parameter α 1 . We see that the IR dynamics is changed by α 1 . This comes from the fact the dilaton can change the asymptotics from −∞ to some constant value with (4.15) as has been explained in section 4. So we observe the increasing coupling constant in the IR regions. For A→ +∞ the running coupling λ goes to 0 for all chosen parameters providing the UV freedom. Or in others words, we can mimic the QCD RG flow for negative α 1 .

Free energy
The free energy corresponding to the black brane solution is given by the renormalized on-shell action. This can be computed directly, but the computation is simplified by the -33 -

JHEP05(2019)117
following observations. We use the domain wall coordinates here. First, the trace of the Einstein equations gives so that the bulk Lagrangian on-shell is (4.52) The finite temperature potential equations imply that and from this relation one can see that the bulk term is a total derivative The normal vector to the cutoff surface w = is n w = √ f , n i = 0. The extrinsic curvature reads The on-shell Einstein action with regularization while the Gibbons-Hawking term with regularization reads We need to evaluate it on the regular black brane solutions given by eqs. (4.17)-(4.19), and (4.26). Taking into account the change of coordinates dw = e 4A f du, the action becomes The expansion of the scale factor A near u ∼ 0 reads with log(1 − e 2µu 02 ), (4.61) Plugging (4.60) in (4.59) we obtain the regularised on-shell action The regularized on-shell is divergent and one needs to add the counterterms before removing the cut-off. The counterterms for the general dilaton-gravity system have been derived in [23]. For our homogeneous solutions the only relevant term is where U is any function that satisfies the equation of the zero-temperature superpotential. We use then the superpotential appropriate for the regular solution, i.e. (see appendix C) The asymptotics of the dilaton is given by (4.68) Using the asymptotics (4.60)-(4.62) and (4.66)-(4.68) we find The renormalized action is then The free energy can be computed through the renormalized on-shell action (4.70); the difference between the free energy of the black brane solution and the free energy of the vacuum, obtained at µ = 0, is On the other hand one can calculate the free energy using black brane thermodynamics, that involves the following relation for the free energy, the entropy density and the temperature of the black brane dF = −sdT.  Performing integration in (4.78) with T given by (4.80) we obtain for the free energy in agreement with the result from the renormalized action.

JHEP05(2019)117
The temperature and the free energy exhibit a qualitatively different behavior depending on the value of k. In the case 0 < k < 2 3 , both are monotonic functions, so the black brane solution exists for any temperature, and there is no phase transition, so the black brane is thermodynamically favored over the thermal gas solution.
In the other case, for 2 3 < k < 4 3 , the temperature increases as function of µ up to a maximum value and then decreases back to zero. The maximum is attained at The free energy is initally negative, but changes sign at We notice that considering u 02 → 0, we get Λ → 0 and the free energy comes to the dependence on the scale disappears, consistently with the fact that the solution becomes the AdS black brane, with the constant dilaton, and there is no phase transition.
In figure 20 we present the behaviour of the free energy on the temperature for different values of k. As we discussed above for 0 < k < 2 3 the free energy is monotonic, while for 2 3 < k < 4 3 the free energy decreases up to T max and then starts to increase, while the temperature is decreasing. This phase diagram cannot be complete, as the free energy would be discontinuous at T max . We do not know whether this means that the dual theory does not make sense for this range of k, or there is some other solution that connects to the ones that we know and restores the continuity. It is worth to notice that a similar behaviour of the free energy was observed in [24] for black hole solutions which are the finite temperature generalizations of the bouncing vacuum solutions.

Conclusion and discussion
We have presented some analytic solutions of Einstein equations coupled to a dilaton field, at zero and non-zero temperatures, that correspond to holographic RG flows between different fixed points. The non-trivial form of the potential allows for a rich variety of different behaviours for the dilaton and the scale factor, corresponding to a coupling that can run to zero or to infinity in the UV and in the IR. We mainly considered solutions which are governed by sinh-functions. However we have two more classes of solutions for our choice of the potential, namely the sin-class and linear class. All solutions depend on two parameters, u 01 and u 02 , splitting the solutions into three branches: "left", "middle" and "right". Moreover, a special solution appears when u 01 = u 02 . In the zero-temperature case the dilaton flows from +∞ to −∞ both for the left and middle solutions. However, the scale factor of the left solution is non-monotonic function while the scale factor of the middle solution monotonically decreases. Correspondingly, the left vacuum solution does not have a holographic interpretation and the middle one is supposed to describe a holographic RG flow that mimics QCD behaviour. The right solution at T = 0 has a bouncing dilaton and a scale factor, that monotonically decreases, i.e. it has both IR and UV free limits. We note that the left, middle and right solutions interpolate between the hyperscaling violating boundaries. A special family of solutions with u 01 = u 02 interpolates between an AdS boundary in the UV limit given by (3.28) and the hyperscaling violating boundary in the IR limit. Another solution with an AdS boundary but in the IR limit (and the hyperscaling violating boundary in the UV limit) is the one from the linear class. For the right vacuum solution we could find its black brane analogue, which again interpolates between the hyperscaling violating boundaries with a constant dilaton at the horizon. The special solution with coinciding point turns to be AdS-Schwarzschild black hole with the corresponding constant dilaton. We also showed that the finite temperature solutions can behave very differently from the corresponding zero-temperature ones. The next step in our program is to understand better these solutions from the point of view of the putative dual field theory. From the phenomenological perspective, we would like to be able to apply our model to real QCD (or Yang-Mills, since we do not introduce quarks). As mentioned in the introduction, the model with a single exponential has undesirable features in this respect. It corresponds to a dual theory with a beta function that is linear in the coupling. The model studied in this paper is an improvement in this direction, since some of the flows describe a theory that is asymptotically free in the UV, with a scale Λ, determined by the boundary conditions of the fields, that can be thought of as Λ QCD , and the coupling grows large in the IR.
It would be interesting to determine the spectrum of fluctuations around these solutions, in particular the quasi-normal modes of the black brane solutions, and also to extend our ansatz in order to look for analytic time dependent solutions that describe out-of-equilibrium dynamics.
It would also be worthwile to understand if the method we used, that reduces the Einstein equations to an integrable Toda chain, could be used for other case. In the JHEP05(2019)117 model of two-exponential potentials that we used, there is a constraint (2.21)-(2.22) on the exponents of the potential; this restricts considerably the space of parameters that we can explore. We plan to consider extensions of the present work to different forms of the potential. For instance, [25] considers a potential with three exponential functions, that can appear after a reduction in the bosonic sector of gauged supergravity (a study of black holes for a particular case of such a potential was done in [26]). Having more terms in the potential allows us to have a richer landscape of critical points, including also de Sitter solutions in addition to AdS and hyperscaling geometries. In general, we would like to have a classification of all cases that can be solved with the techniques used here. One way to approach this problem could be to reverse-engineer the dilaton potential and understand if it can come from some dimensional reduction.
It would also be interesting to consider the coupling of a Maxwell field in the bulk and find charged black hole solutions. tions. AG is also very grateful for the warm hospitality DESY Theory Group, APCTP (Pohang) and Sogang University (Seoul). I.A. is supported by the "BASIS" Science Foundation (grant No. 18-1-1-80-4).

A The curvature invariants of the background
The ansatz for the metric is given by The non-zero components of the Ricci tensor are The scalar curvature reads R = −2e −2C Ȧ 2 + 3ȦḂ + 6(Ḃ) 2 −ȦĊ − 3ḂĊ +Ä + 3B . where the functions A 1 = B 1 and α 1 is a constant, which is equal to zero for the vacuum case. Now the scalar curvature (A.5) reads that covers the vacuum case with α 1 = 0 The Kretschmann scalar is defined as follows Pluggin the metric (A.1) we have One can write down the Kretschmann scalar of the non-vacuum solution (α 1 = 0) Taking into account that C = 4A 1 and α 1 = 0 for the vacuum solution one obtains
• For the right solution one has u → u 01 + The scalar curvature matches with that one (A.23)-(A.24) for the middle solution near the point u 01 .
(A.27) Here one has to study the behaviour of the scalar curvature in the limits of small u and u → +∞.
For u → u 0 one has that is an agreement with (A.25). The corresponding expression in the scalar coordinates matches with (A.26).