Phase structures emerging from holography with Einstein gravity - dilaton models at finite temperature

Asymptotic AdS Riemann space-times in five dimensions with a black brane (horizon) sourced by a fully back-reacted scalar field (dilaton) offer - via the holographic dictionary - various options for the thermodynamics of the flat four-dimensional boundary theory, uncovering Hawking-Page, first-order and second-order phase transitions up to a cross-over or featureless behavior. The relation of these phase structures to the dilaton potential is clarified and illustrating examples are presented. Having in mind applications to QCD we study vector mesons in the probe limit with the goal to figure out conditions for forming Regge type series of radial excitations and address the issue of meson melting.


Introduction
The advent of the insight into AdS/CFT correspondence [1][2][3] offered the option of having an alternative access to strongly coupled systems, e.g. to various facets of QCD in the nonpertubative regime, for instance. Phenomenologically interesting problems, e.g. the hadron spectrum or properties of the quark-gluon plasma, become treatable within a framework called holographic approaches. Most desirably would be to have a holographic QCD dual at our disposal, from which statements on QCD-related quantities can be derived in a unique manner. However, such a dual is presently not available [4][5][6][7]. Therefore, in practice, field-theory quantities in a five dimensional asymptotic Anti-de Sitter (AdS) are often related to observables (or expectation values of operators) in four dimensional Minkowski space-time in the spirit of the field-operator duality ( [8], see also §5.3 in [9], §10.3 in [10], for instance). Top-down approaches attempt to use input from string theory constructions or elements thereof. These are to be contrasted with bottom-up approaches which aim at starting with an appropriate model on the field-theory side to mimic certain selected features of the boundary theory side, with the latter being connected to the quantum field theory in Minkowski space, while the former includes the dynamics in the bulk. Besides early emphasis on accessing principal features of strongly coupled systems with many extensions to higher or lower dimensions than mentioned above, one can also take the attitude of adjusting sufficiently simple and thus transparent bottom-up models to a certain input and then employ them for predictions. Of course, the predictive power becomes a relevant issue here. Moreover, the foundations of the AdS/CFT correspondence, namely a very large number N c of gauge degrees of freedom and a very large 't Hooft coupling, are often argued only to hold under special conditions, too. For instance, w.r.t. QCD one knows [11] that certain thermodynamic observables of the Yang-Mills gauge theory obey the proper scaling with N c and hopes that the physics case of N c = 3 is adequately captured. Holographic modeling of QCD related problems became popular due to some particularly striking findings. Among them are the Regge type spectrum of hadronic and glue ball states, e.g. within the soft-wall model [7], the famous ratio of shear viscosity η to entropy density s, η/s = 1/4π [12] , and the phase diagram with a critical point [13,14], to mention a few ones. Besides gravity, the dilaton plays an important role as a breaker of the conformal symmetry since it introduces an energy scale. Obviously, the holographic models with gravity coupled to and sourced by a dilaton field -including the negative cosmological constant to ensure the asymptotic AdS geometry -represent some minimalistic set-up. To be specific, we restrict ourselves here to Einstein gravity. What remains is fixing the dilaton self-interaction. This may refer to roots in string theory, as recently put forward, e.g. in [15,16] in a topdown approach, or to shape the dilaton self-interaction -encoded in the dilaton potentialby reproducing a certain set of wanted results within the dilaton engineering to reproduce Lattice QCD thermodynamics results of [17,18]. The resulting set-up is called Einstein gravity -dilaton model. It continues numerous previous studies in cosmology, most notably in inflationary scenarios. Analogously, in holography such gravity-dilaton models enjoy some popularity due to their conceptual simplicity. There is an overwhelming number of studies, e.g. [15,, based on that type of model: [19,20] represent an in-depth analysis and review of the model in detail and [21] touches the issue of consistency. References [15, focus on thermodynamics, where [22] can be considered to be the prototype of modeling thermodynamics with a dual black hole. The authors of [23,24] give a fundamental bound of the speed of sound and [25] derived a relation between the speed of sound and the single heavy quark free energy. Discussions about different phase structures can be found in [26][27][28][29], and for quantitative comparisons and parameter fits to results on Lattice QCD thermodynamics in case of vanishing and finite baryo-chemical potential we refer to [30][31][32][33][34]. There are investigations about the temperature dependence and the behavior during phase transitions of related quantities, e.g. [35] calculates string tension at finite temperature, [36] chooses an approach based on the beta function, [37][38][39] deal with electric and magnetic quantities, [15,40] calculate the Debye screening mass; transport coefficients and bulk viscosities are the topics of [41,42] and [43,44], respectively. A holographic approach to the broad field of hydrodynamics and thermalization is given e.g. by [45][46][47] within the gravity -dilaton model class.
The soft-wall model [7] developed into a role model for computing particle spectra holographically. While in the original model the metric background is fixed by ad-hoc ansätze, the main idea of its generalizations [48][49][50][51][52][53][54][55][56][57][58][59][60] is to obtain the metric background as a solution of the Einstein equations. The particles are considered then as test particles. [48][49][50][51][52] give particular attention to Regge type spectra, [54] investigates the case of finite temperature and [55,56] focus on chiral symmetry breaking. Due to the various applications of Einstein gravity -dilaton models (see as well [61] for a braneless approach, [62] for fluctuating branes, [63] for a real-time formulation, [64] for cosmological discussions, [65] for scalar condensates or [66] for a generalization to higher dimensions) this list does not purport to be complete. Finite temperature effects are generated by plugging a black hole in the originally AdS and deform it accordingly. Thus, a Hawking temperature and a Bekenstein-Hawking entropy density link to thermodynamics. In the present paper, we also stay within such a framework: holographic gravity-dilaton model with the goal to elaborate the emerging thermodynamics w.r.t. conditions for the dilaton potential to catch certain phase structures with relevance to QCD. 1 The Columbia plot (cf. [68] for an updated version) provides several options for 2+1 flavor QCD: in dependence on the quark masses, first-and secondorder phase transitions may show up as well as a cross-over and some others. We try to answer the question which properties the dilaton potential must have to enable these phase properties related to deconfinement and chiral restoration in QCD. On top of thermodynamic aspects we consider holographic vector mesons in the probe limit. That is, the gravity and dilaton background resulting from the field equations and equation of motion governs -besides the thermodynamic features -the behavior of vector mesons. We thus extend our previous studies [69] and investigate to which extent the disappearance of vector mesons as a possible indicator of deconfinement occurs at the QCD cross-over temperature. This is important for a holographic realization of the thermo-statistical interpretation of hadron multiplicities in ultra-relativistic heavy-ion collisions, e.g. at LHC [70]. The ultimate goal of such investigations, which however is beyond the scope of the present paper, is an extension to non-zero chemical potential, e.g. to address issues of the QCD phase diagram and the chemical freeze-out curve therein. Our paper is organized as follows. In Section 2, we recall the holographic settings, that is the gravity-dilaton model, its field equations and equation of motion as well as the emerging thermodynamics and the access to phase structure. Also, the holographic description of vector mesons is recalled. Section 3 deals with thermodynamic scenarios by presenting a series of selected examples of transition types characterized by entropy density, sound velocity and pressure. This is supplemented by showing the Schrödinger equivalent potential which governs the existence or non-existence of vector mesons. After some general remarks on shaping the dilaton potential, we try to elucidate the conditions on the dilaton potential to enforce a first-order phase transition or a cross-over. In Section 4, we explain a relation between the Schrödinger equivalent potential at zero temperature and the thermodynamic features. Both ones are linked by the field equations; details can be found in Appendices A and B. In the second part of Section 4, we reverse our view: instead of starting with a dilaton potential we model a certain shape of the Schrödinger equivalent potential which allows for a certain wanted hadron spectrum (ideally of a Regge type) and derive -again via field equations -the resulting dilaton potential. 2 We conclude that part by considering the vector meson melting upon temperature increase. The summary and a discussion of possible extensions towards the goal of a consistent scenario of QCD thermodynamics with the chemical freeze-out model in the LHC energy regime can be found in Section 5.

Thermodynamics from Einstein gravity -dilaton model
We consider the action over the five-dimensional Riemann space-time with special ansatz of the metric given by the line element squared as where A denotes the warp factor with A(z) → −2 ln(z/L) as z → 0 to ensure an asymptotic Anti-de Sitter (AdS) and f denotes the blackness function with f (z) = 1 − O(z 4 ) which encodes the temperature via with the horizon position z H and the simple zero f (z H ) = 0. The vacuum case, T = 0, is equivalent to f = 1. The dilaton Φ in the action (2.1) is a dimensionless real-valued scalar bulk field. Its potential V (Φ) has the asymptotic small-Φ form where the first term refers to the negative cosmological constant and the second one has to obey the Breitenlohner-Freedman (BF) bound −4 ≤ m 2 Φ L 2 ≤ 0 [71,72]; L sets a scale, as κ in (2.1), to make the action dimensionless in natural units. From (2.1) and (2.2) the field equations follow and the equation of motion which is redundant since it follows from (2.6) with (2.4, 2.5). A prime means derivative w.r.t. the bulk coordinate z. Equations (2.4-2.6) can be solved for a given V (Φ) with the above side conditions. In such a way, the dilaton potential V (Φ) determines the temperature via (2.3) and the speed of sound squared via } stands for the entropy density. The quantities T , c 2 s and s refer to the boundary theory according to the holographic dictionary. The pressure is calculated via p = dT s(T ) with the side condition p(T = 0) = 0.

Vector mesons in probe limit
Additionally, we study the behavior of vector mesons in the probe limit, i.e. they are not back-reacted. To capture their features, we use the standard action [7] where F 2 is the squared field strength tensor of a U (1) vector field. The equation of motion follows, after some manipulations [73], as one-dimensional Schrödinger type equation [69] with the Schrödinger equivalent potential For T = 0, the tortoise coordinate ξ becomes z and only U 0 (z) is relevant. m n denotes the masses of normalizable modes as solutions of (2.10) with n = 0, 1, 2, · · · as the quantum number of radial excitations. In general, A and Φ depend on both, z and z H . To distinguish the vacuum case (T = 0: we add the index 0 to A and Φ, as above for U 0 . By employing the field equations, we find the following relation for U 0 (z): In Section 3.1 we study the vector meson spectrum over the background determined by solutions of (2.4-2.7) within U 0 from (2.13), while in Section 4.2 an ansatz for U 0 (z), which facilitates a certain mass spectrum, is used as an input for (2.4-2.7) to figure out the related thermodynamics and phase structure.  3 Thermodynamic scenarios

Selected examples
To illustrate the systematics of the thermodynamics related to the dilaton potential V we choose the three-parameter ansatz [22,74] − L 2 V (Φ) = 12 cosh(γΦ) + aΦ 2 + bΦ 4 (3.1) because we can go through several thermodynamic scenarios (see Fig. 1) by changing the values of parameters γ, a and b, which are related by −L 2 m 2 Φ = 12γ 2 + 2a to the dilaton mass parameter m 2 Φ . As already remarked in [22], quite featureless dilaton potentials V (Φ) can lead to fairly different thermodynamic features. Since the field equations (2.4-2.6) can be rearranged to display only a sensitivity to ∂ Φ V /V as a function of Φ, we plot this key quantity in the left column of Fig. 1. Example a) does not exhibit any features: T (z H ) is monotonously decreasing, c 2 s (T ) is increasing, and U 0 has no minimum, meaning that normalisable modes as vector mesons do not exist at all. Examples b) -d) exhibit ∂ Φ V /V with a pronounced maximum which becomes gradually higher. In case b), T (z H ) is monotonously dropping, albeit with a shallow near-flat section; it causes a pronounced minimum of the sound velocity squared; U 0 does not allow for any normalizable states due to the lacking minimum. That case is classified as cross-over. Case c) features a lifted maximum of ∂ Φ V /V , resulting in a flat section of T (z H ), and the sound velocity drops at a certain temperature to zero, thus representing an example of a second-order phase transition. The Schrödinger potential U 0 (z) has here a very shallow minimum, i.e. vector mesons as normalizable modes cannot show up. Lifting the maximum of ∂ Φ V /V further (case d)), T (z H ) shows a local minimum that is connected to an inflection point, implying metastable states, unstable states and spinodales as well. That becomes most evident by the sound velocity squared c 2 s (T ) and the pressure p(T ) (see Fig. 2): metastable states are depicted by thin solid curve sections and the unstable ones by dotted sections. Clearly, states with c 2 s < 0 cannot be realized in nature. All these features classify a first-order phase transition, for which U 0 (z) exhibits a local minimum; it is -for the given parameters -still too shallow to accommodate vector mesons. Only if ∂ Φ V /V exceeds 2/3 (case e)), the global minimum of U 0 (z) allows for meson states, depicted by horizontal lines. At the same time, T (z H ) has a global minimum pointing to a Hawking-Page (HP) phase transition. Let be T min = min T at z min H . Then, the branch for z H < z min H is stable (for p > 0) and metastable (for p < 0), while the branch for z H > z min H is unstable, and its free energy is above the thermal gas solution (see [26] for the related construction) which applies for 0 < T < T min < T c < ∞, where T c (slightly above T min ) is the first-order phase transition temperature. (While T min follows as minimizer of (2.3) over z H , T c is obtained by the loop construction sketched in the rightmost graphs of Fig. 2.) The velocity of sound drops to zero at T c . Figures 3 and 4 exhibit that region in the γ-a-b parameter space, where a first-order or HP phase transition occurs; and in Fig. 4, curves on which a second-order phase transition happens are depicted too. Thanks to the three-parameter potential (3.1), the visualization of the parameter space structure is quite straightforward, while multi-parameter ansätze may lead to intricate structures. In the present case, the HP transition happens for γ ≥ 0.8 for all BF permitted values of a and b as well, as exhibited in the left panel of Fig. 4. For γ < 0.8, the blue areas depict the first-order phase transition, which are limited by the second-order transition (red curves). Further left (gray parts of the panels with constant γ or a or b cross sections), a cross-over occurs, which turns smoothly into a featureless behavior for smaller values of b.

Shaping the dilaton potential
We refrain here from a quantitative comparison with QCD thermodynamics and refer the interested reader to [15,22,26,32,33], for example, where the thermodynamics of Yang-Mills and 2+1 flavor QCD with physical quark masses is considered. The essence is to employ multi-parameter ansätze of the dilaton potential with the aim to reproduce the Lattice QCD thermodynamics results as good as possible. 3 Information of QCD is thus imported and mapped in a cumulative manner on V (Φ) without explicit reference to quarks and gluons and their masses, colors, flavors, couplings etc. 4 Having successfully accomplished the shaping of V (Φ), one can proceed to derive further quantities, such as viscosities [43,44], diffusion constants [41,42] etc. as predictions. By extending the model (2.1) by further fields, e.g. a Maxwell type U (1) gauge field [13,33,81], one may address non-zero chemical potential effects to access the phase diagram and issues of a critical point. Here, susceptibilities serve as crucial further information to be imported from QCD. However, our present goal here is to answer the question whether one can describe -within the modeling (2.1, 2.9) -at the same time QCD thermodynamic features, i.e. a cross-over at about 155 MeV [17,18], and a proper in-medium behavior of hadrons, i.e. vector mesons in probe limit as representatives thereof. By a proper "in-medium behavior" we mean that  . Plot of parameter ranges of the dilaton potential (3.1) for first-order phase transition (blue areas) or HP transition (yellow areas) for a series of selected constant values of γ = 0, · · · , 1 (left), a = −6, · · · , 2 (middle) and b = 0, · · · , 0.4 (right). The red curves mark the second-order phase transition ranges; beyond, the transition turns into a cross-over, followed by a featureless behavior (dark gray areas). The BF bound restricts the values of a to the strip −6γ 2 < a < 2 − 6γ 2 (depicted by the two black solid curves in left and middle panels). The panels continue to larger values of b and γ without changes.
at chemical freeze-out temperature of about 155 MeV [70], hadrons do exist with hardly medium-modified properties. Otherwise, the famous thermo-statistical interpretation of hadron multiplicities in ultra-relativistic heavy-ion collisions [70] would be invalidated.
The answer to the posed problem seems to be negative. Hints come, for example, from [73,[82][83][84][85], where the melting (disappearance) of hadrons was found to happen at temperatures significantly below 155 MeV. While [69] offers an avenue to remedy such an insanity, the given framework of (2.1, 2.9) seems to be too restricted and calls for extensions. Leaving the latter ones for separate work, we try to find a loophole to join the cross-over thermodynamics and suitable vector meson states. Prior to that, however, we attempt to clarify the systematics of thermodynamic features in the spirit of the last column of Tab. 1 in relation to the dilaton potential.

Beyond the adiabatic approximation
The authors of [22] derived the relation (henceforth called Gubser's adiabatic criterion) where '≈' indicates the validity in adiabatic approximation. The formula implies that for ∂ Φ V /V > 2/3 the sound velocity becomes imaginary, thus pointing to a first-order phase transition, either as a standard constructionà la example d) or the HP transitionà la example e) in To formalize these findings we derive in Appendix A the relation Given the facts that (i) T (z H → 0) → 1 πz H [23,24], (ii) the monotonous behavior of Φ(z, z H ) as a function of z with Φ > 0, and (iii) the above quoted asymptotic behavior of V (Φ) at small Φ (implying ∂ Φ V /V = m 2 Φ L 2 Φ/12), one recognizes from the first line of (3.3) that the slope dT /dz H , which is negative at small z H , can turn into a positive one, once ∂ Φ V /V > 2/3 is reached, indicating a local or global minimum of T (z H ). Understanding "adiabatic approximation" as a situation where ∂A/∂z H , ∂A /∂z H and ∂Φ/∂z H are small, one thus recovers Gubser's adiabatic criterion. Otherwise, the second line of (3.3) provides corrections. In fact, in example d), ∂ Φ V /V stays below the 2/3 line but facilitates a first-order phase transition. (The relation of T (z H ) to transitions is discussed in [19,20]: in essence, a minimum of T (z H ) points to a first-order phase transition, since s(z H ) is a monotonously increasing function.) The corrections give eventually a border line (called T min curve) for each type of dilaton potential which is determined by calculating the minimal value of ∂ Φ V /V (as a function of Φ and depending on all parameters denoted shortly by p ) such that T forms a minimum. Systematic numerical analyses with the dilaton potential (3.1) show that this line is shifted down if ∂ Φ V /V as a function of Φ becomes steeper when varying the parameters p in V (Φ; p). The right part of Figure 5 shows an example of such a line and the dependence of the difference between Gubser's criterion and the T min curve as a function of Φ. This is further visualized in Fig. 6 for two special parameters, b = 0 and 2, in the projections on the γ-a plane: the offset of the regions determined by (3.2) and the true onset of a first-order phase transition increases with b; in addition, the region where the Schrödinger equivalent potential U 0 displays a minimum is shown by the red curves.

Adiabatic approximation to Schrödinger potential
The relationship between the situation of T = 0 and features at T > 0 has been stressed in [19,20]. Here, we envisage a relation of U 0 (z) and T (z H ) in adiabatic approximation. In Appendix B we derive the relation 9U 0 ∼ = 17 12 It is valid if, in a decomposition A(z, z H ) = A 0 (z) + a(z, z H ) and Φ(z, z H ) = Φ 0 (z) + ϕ(z, z H ), the terms a(z, z H ) and ϕ(z, z H ) are sub-leading and can be neglected. The Chamblin-Reall solution [86] with V (Φ) ∝ exp{γΦ} is an example, where a = ϕ = 0 can

Requiring a minimum of U 0 (z)
The above examples demonstrate that for many parameter choices of the dilaton potential V (Φ) the Schrödinger potential U 0 (z) does not exhibit a minimum and thus does not allow for modes which can be interpreted as vector mesons. Instead of deriving the background (warp factor and dilaton profile at T = 0) from given dilaton potential, we start now with an ansatz for U 0 (z) such as to have a minimum. Assuming the latter one is sufficiently deep, normalizable modes would then be expected. Our ansatz for demonstrative purposes is where the first term comes from the asymptotic warp factor at z → 0; the second term facilitates the required minimum at z min /L = (3/2p) withẑ ≡ z/L and c 1 = 1 (due to AdS behavior at boundary z → 0) and ) (due to the assumption that (4.2) is globally valid). The field equations (2.4-2.6) must be solved numerically to get A(z), Φ(z) and V (Φ). The same V (Φ), which is supposed to be independent of T , is then used to derive U T (z, z H ) and T (z H ). Figure 7 exhibits such solutions for p = 0.5, 1 and 2, where the latter value reproduces the soft-wall model [7] with a strictly linear Regge type spectrum L 2 m 2 n = 4(n + 1) for n = 0, 1, 2 · · · . The left panel is for U 0 (z) according to (4.2), while the middle panel shows T (z H ); the right panel displays ∂ Φ V /V as a function of Φ. There is a striking similarity of the curves T (z H ) and U 0 (z) (isotopy referring to the monotony behavior) which we interpret as follows: a minimum of U 0 (z) is related to a minimum of T (z H ), i.e. a first-order phase transition -here a HP transition, since it is a global minimum. The behavior of ∂ Φ V /V as a function of Φ is in agreement with our assessments in Sections 3.1 and 3.3. To get an idea on the related scale, the left-hand plot of Fig. 8 shows the first three states, that is L 2 m 2 n as a function of p for n = 0 (ground state) and n = 1, 2 (first two radial excitations). For comparison, the middle plot of Fig. 8 displays LT min also as a function of p. Both figures can be combined to T min /m n as a function of p (right-hand plot of Fig. 8). Having in mind applications to QCD and identifying T min with 155 MeV (see above) and m 0 with the ρ meson ground state mass of 770 MeV, one arrives at T min /m 0 ≈ 0.20, i.e. a value not too far from the range of values shown in the right-hand plot of Fig. 8. However, 2+1 flavor QCD with physical quark masses does not provide a first-order phase transition. Insofar, the present set-up is more appropriate for 2+1 flavor QCD in the chiral limit, which in fact enjoys a first-order phase transition [68], but detailed information on  thermodynamic quantities as well as the vector meson spectrum is lacking (cf. [87] for a search for the delineation curve in the Columbia plot where cross-over and first-order phase transitions touch each other; very preliminary first estimates [88] point to a ratio of T c /m 0 in the same order of magnitude as for the case of physical quark masses). The temperature dependence of U T (not shown) is such to cause the instantaneous disappearance (melting) of vector meson states at T c ≈ T min (see Fig. 9). Denoting the disappearance temperature by T dis and choosing m 0 = m ρ as scale, Fig. 9 translates into T dis = m ρ LT dis Lm 0 , meaning T dis = 116 MeV for LT dis = 0.27 (0.3) and T m 0 = 1.7 (2). That value of T dis does not apply to 2+1 flavour QCD with physical quark masses, since, as stressed above, the present scenario is more suitable for the chiral limit where a first-order phase transition occurs.

Summary and discussion
We focus here on the QCD relevant first-order, second-order and cross-over transitions. Obviously, more complicated structures are possible, e.g. a sequence or nested first-order transitions for functions T (z H ) with multiple local minima (see [26] for the case of a double transition). These require further shaping of the dilaton potential. This can be easily done by combining the elements of our systematics presented in this paper: an extreme point of ∂ Φ V /V generates an inflection point of T (z H ) which points to a cross-over or a second-order phase transition (if it is a horizontal turning point) or to a first-order phase transition if the temperature exhibits additional extreme points which can be controlled by the altitude of ∂ Φ V /V . We did not touch such issues as good and bad curvature singularities [89], adding further (e.g. charged scalar) fields which can bridge to order parameters and/or condensation [90], larger classes of dilaton potentials (e.g. Liouville potentials or linear combination thereof [91]) and fluctuations. Unfortunately, the Einstein gravity -dilaton model seems to be not flexible enough to allow simultaneously for a cross-over and meson states in the probe limit because the existence of the latter ones requires a minimum in T (z H ). The obvious idea to construct a dilaton potential such that T (z H ) has a minimum at a horizon z * H with T (z * H ) being small (to ensure the existence of the mesons) and a cross-over at the QCD critical temperature of about 155 MeV does not work very well, since all mesons states disappear already at the minimum temperature T (z * H ). A first step further on the road to a fully consistent approach could be the consideration of a U (1) Maxwell type gauge field. Such a field has been used to address the question of the behavior of vector mesons in relation to the thermodynamics: ρ mesons are described through the U (1) gauge field (see (2.9)) and putting together the actions (2.1) and (2.9) would yield a model with full back reaction from the mesons to the gravity background. Since QCD thermodynamics is not driven by vector mesons alone, another step is adding flavor by including the pseudo-scalar and scalar sectors via the bulk fundamental fields and its vacuum expectation values. Some works point directly in this direction: the authors of [92] introduce a second scalar field (glue ball field) and solve the field equations for the case T = 0; many other investigate the behavior of hadron species in a given background without back reaction (see e.g. [73,82,[93][94][95][96][97][98][99]). [100,101] give a study of phase transitions in relation to a flavor containing model with given metric background. Bringing the characteristic features of the mentioned works together would be an improvement. The Holy Grail would be a model with parameters steering quark masses and condensates for different flavors separately with proper hadron spectra in all (scalar, pseudo-scalar, vector, axial-vector, tensor and axial-tensor) sectors. All mentioned extensions point to the leading question, which framework is needed to have a QCD consistent thermodynamics and proper in-medium modifications of the hadron species. However, increasing the variety of a model means increasing its complexity and requires much follow-up work. By varying a second parameter (here p 3 = −0.01 · · · 0.01), the T min curve becomes a strip (red area), where a negative (positive) value of p 3 belongs to the lower (upper) part, since increasing p 3 means increasing the slope of ∂ Φ V /V and therefore increasing also the difference between Gubser's adiabatic criterion and the T min curve. If ∂ Φ V /V as a function of Φ with given parameter set p has a section with the corresponding T min area (here only displayed of a p 3 interval) then a first-order or HP phase transition is facilitated. Right: Border lines of parameter regions in −L 2 m 2 Φ = p 2 /24 vs. p 4 space, where either T (z H ) has a minimum (right to blue curve), or U 0 (z) has a minimum (right to red curve), or ∂ Φ V /V is greater than 2/3 (right to green curve) for the dilaton potential (C.1) with p 3 = 0 and p i≥5 = 0; if m Φ and therefore the slope of ∂ Φ V /V increase, the difference between Gubser's adiabatic criterion and the T min curve becomes larger.
B Derivation of (4.1) We start with (2.13). To relate U 0 with T , we evaluate (A.1, A.2, A.5, A.8) in leading order, i.e. neglecting a and ϕ: These four equations allow for eliminating all A's and Φ's in (2.13). Finally we arrive at (4.1).

C Towards a systematic dilaton potential expansion analysis
Another useful form of the dilaton potential is since ∂ Φ V /V runs over all BF permitted and AdS conform polynomials if p = (p 2 , p 3 , · · · ) runs over all vectors. The case p i≥3 = 0 characterizes the leading order (straight lines) in the spirit of an expansion of ∂ Φ V /V in powers of Φ. The left panel in Fig. 10 shows the influence of the varying p 3 -depending term (red strip) on the T min curve (black curve) which is generated with running p 2 analogously to Fig. 5. The parameter regions in the −L 2 m 2 Φ vs. p 4 plane, where T (z H ) (blue curve) or U 0 (z) (red curve) exhibits a minimum, as well as the area, where ∂ Φ V /V is greater than 2/3 (green curve), are also shown (see right panel and compare with Fig. 6). In such a manner one can study, piece by piece, the impact of the individual terms in (C.1) on the issue of phase structure and capabilities to permit vector meson modes in the probe limit.