This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

The quest for a conifold conformal order

Alex Buchel
Department of Physics and Astronomy
University of Western Ontario
London, Ontario N6A 5B7, Canada
Perimeter Institute for Theoretical Physics
Waterloo, Ontario N2J 2W9, Canada
(April 29, 2022)

The quest for a conifold conformal order
Alex Buchel Department of Physics and Astronomy University of Western Ontario London, Ontario N6A 5B7, Canada Perimeter Institute for Theoretical Physics Waterloo, Ontario N2J 2W9, Canada


Abstract

The holographic duality between cascading gauge theory and type IIB supergravity on warped deformed conifold with fluxes reveals exotic thermal phases with nonzero expectation values of certain operators, persistent to high temperatures. These phases, in the limit of vanishing the strong coupling scale of the cascading gauge theory, would realize thermal ordered conformal phases in 3,1{\mathbb{R}}^{3,1} relativistic QFT. We find that the dual Klebanov-Strassler/Klebanov-Tseytlin black branes in this limit are outside the regime of the supergravity approximation, rendering the construction of such conformal ordered states unreliable. While we have been able to construct conformal order in phenomenologically deformed effective theory of type IIB supergravity reduced on warped deformed conifold with fluxes, the removal of the deformation parameter causes the destruction of the thermal conformal ordered phases. Once again, we find that the holographic models with the conformal ordered phases are in the String Theory swampland.


April 29, 2022

1 Introduction

Conformal order stands for exotic thermal phases of conformal field theories (CFTs), characterized with nonzero one-point correlation function(s) of certain operator(s) [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. For a CFTd+1 in Minkowski space-time d,1{\mathbb{R}}^{d,1} the existence of the ordered phases implies that there are at least two distinct thermal phases:

Td+1=𝒞×{1,TΔi𝒪Δi=0,κ,TΔi𝒪Δi=γi0,\frac{{\cal F}}{T^{d+1}}=-{\cal C}\ \times\ \begin{cases}1,\qquad T^{-\Delta_{i}}\langle{\cal O}_{\Delta_{i}}\rangle=0\,,\\ \kappa,\qquad T^{-\Delta_{i}}\langle{\cal O}_{\Delta_{i}}\rangle=\gamma_{i}\neq 0\,,\ \end{cases} (1.1)

where {\cal F} is the free energy density, TT is the temperature, 𝒞{\cal C} is a positive constant proportional to the central charge of the theory, and {𝒪Δi}\{{\cal O}_{\Delta_{i}}\} is the set of the order parameters with the conformal dimension spectrum {Δi}\{\Delta_{i}\}. The parameters κ\kappa and {γi}\{\gamma_{i}\} characterizing the thermodynamics of the ordered phase are necessarily constants. Conformal order can realize spontaneous breaking of discrete [2, 1, 3] of continuous [10] global symmetries; but it does not have to be the case: in the model we discussed in [9], the conformal order parameter is not associated with spontaneous breaking of any global symmetry111This is strictly true in the context of the effective Kaluza-Klein reduced holographic models. The nonvanishing scalar expectation value would signal the spontaneous breaking of the global symmetry of the compactification manifold upon uplift to 10d supergravity. Unfortunately, no model with a thermal conformal order has been constructed in top-down holography yet..

From (1.1), note that when κ>1\kappa>1 (κ<1\kappa<1), the symmetry broken phase dominates (is subdominant) both in the canonical and the microcanonical ensembles. Irrespectively of the value, provided κ>0\kappa>0, the symmetry broken phase is thermodynamically stable. It is difficult to compute directly in a CFT the values {κ,γi}\{\kappa,\gamma_{i}\}, thus establishing the presence and the (in)stability of the ordered phase. Rather, the authors of [1, 6, 7, 8] established the instability of the disordered thermal phases in discussed CFTs. The condensation of the identified unstable mode then leads to 𝒪Δi0\langle{\cal O}_{\Delta_{i}}\rangle\neq 0 for the new equilibrium thermal state — the conformal order.

Conformal order states are very interesting in the context of holography [11, 12], as they imply the existence of the dual black branes in a Poincare patch of asymptotically AdSd+2AdS_{d+2} bulk geometry that violate the no-hair theorem. In [9] we proved a theorem that the disordered conformal thermal states are always stable in dual holographic models of Einstein gravity with multiple scalars. Thus, the mechanism for the conformal order presented in [1] is not viable in these holographic models.

The first holographic model of the conformal order was discovered (though not appreciated in this context prior to the QFT construction [1]) purely by accident in [2]. The general framework for constructing holographic conformal order, the phenomenologically deformed effective theory, was presented in [5] — we now review those arguments as they will be utilized in this paper222All the known constructions of the holographic conformal order can be understood within this framework.. Consider a top-down holographic model, dual to333Extensions to AdSd+2/CFTd+1 models with multiple order parameters 𝒪Δi{\cal O}_{\Delta_{i}} is trivial — such models will be studied in section 3. a CFT4 with a single operator 𝒪Δ{\cal O}_{\Delta} of a dimension Δ\Delta. The five-dimensional gravitational bulk effective action takes the form

S5=c2π2L35vol5{R12(ϕ)2𝒫[ϕ]},S_{5}=\frac{c}{2\pi^{2}L^{3}}\int_{{\cal M}_{5}}{\rm vol}_{{\cal M}_{5}}\biggl{\{}R-\frac{1}{2}(\nabla\phi)^{2}-{\cal P}[\phi]\biggr{\}}\,, (1.2)

where cc is the central charge of the boundary CFT, ϕ\phi is the gravitational bulk scalar dual to the order parameter 𝒪Δ{\cal O}_{\Delta}, and the scalar potential 𝒫[ϕ]{\cal P}[\phi] is

𝒫[ϕ]=12L2+Δ(Δ4)2L2ϕ2+𝒪(ϕ3).{\cal P}[\phi]=-\frac{12}{L^{2}}+\frac{\Delta(\Delta-4)}{2L^{2}}\phi^{2}+{\cal O}(\phi^{3})\,. (1.3)

The 𝒪(ϕ0){\cal O}(\phi^{0}) term in (1.3) is a negative cosmological constant, setting the radius of the asymptotically AdS5 bulk geometry to LL; the mass term,

L2m2Δ(Δ4),L^{2}m^{2}\equiv\Delta(\Delta-4)\,, (1.4)

represents the standard encoding of the dimension of the order parameter 𝒪Δ{\cal O}_{\Delta} [13, 14]. In what follows we set L=1L=1. It is vital that in real holographic models (contrary to the phenomenological toys), the full scalar potential 𝒫[ϕ]{\cal P}[\phi] is nonlinear. Unfortunately, holography is not understood at the level where given a boundary CFT, with a spectrum of gauge invariant operators, we can engineer/compute the scalar potential. In specific holographic examples, like the 𝒩=2{\cal N}=2^{*} correspondence [15, 16, 17], the cascading gauge theory duality [18, 19] or the Maldacena-Nunez model [20], the scalar potential is computed from the realization of the duality correspondence in type IIB supergravity. The construction of the holographic conformal order proposed in [5] relies on (and is applicable to) models where the leading nonlinear correction is unbounded from below444This feature is ubiquitous in top-down holographic examples. along certain directions on the scalar manifold. To be specific, we assume that

𝒫[ϕ]=12+Δ(Δ4)2ϕ2s22ϕn+𝒪(ϕn+1),{\cal P}[\phi]=-{12}+\frac{\Delta(\Delta-4)}{2}\phi^{2}-\frac{s^{2}}{2}\ \phi^{n}+{\cal O}(\phi^{n+1})\,, (1.5)

for some constant ss and an integer n3n\geq 3. Note that 𝒪(ϕn){\cal O}(\phi^{n})-truncated scalar potential is unbounded as ϕ+\phi\to+\infty. Of course, the full scalar potential 𝒫[ϕ]{\cal P}[\phi] might be bounded, but this is not important for the perturbative construction of the thermal conformal order in the phenomenologically deformed model. The phenomenologically deformed model is defined as a holographic correspondence where the top-down scalar potential 𝒫[ϕ]{\cal P}[\phi] is deformed as

𝒫[ϕ]𝒫b[ϕ]12+Δ(Δ4)2ϕ2+b(𝒫[ϕ]+12Δ(Δ4)2ϕ2),{\cal P}[\phi]\ \longrightarrow\ {\cal P}^{b}[\phi]\equiv-{12}+\frac{\Delta(\Delta-4)}{2}\phi^{2}+b\biggl{(}{\cal P}[\phi]+12-\frac{\Delta(\Delta-4)}{2}\phi^{2}\biggr{)}\,, (1.6)

with the constant deformation parameter bb being positive. The claim of [5], see also appendix A, is that the thermal conformal order always exists in the limit b+b\to+\infty, when the thermal ordered phase is holographically realized as AdS-Schwarzschild black brane, with a perturbatively small “scalar hair”

ϕ(1b)1n2TΔ𝒪Δ=γ(1b)1n2.\phi\ \propto\ \left(\frac{1}{b}\right)^{\frac{1}{n-2}}\qquad\Longleftrightarrow\qquad T^{-\Delta}\langle{\cal O}_{\Delta}\rangle=\gamma\propto\left(\frac{1}{b}\right)^{\frac{1}{n-2}}\,. (1.7)

The existence of the thermal conformal order in real holography then boils to the question whether this perturbative constructions survives as bb decreases from ++\infty to 11, since

limb1+𝒫b[ϕ]=𝒫[ϕ].\lim_{b\to 1_{+}}{\cal P}^{b}[\phi]\ =\ {\cal P}[\phi]\,. (1.8)

In this paper we continue the quest for constructing the thermal conformal order in String Theory holography. Our focus is on top-down holographic dualities between regular and fractional D3-branes on a conifold, the simplest non-compact Calabi-Yau threefold [21], and 𝒩=1{\cal N}=1 supersymmetric gauge theories - the Klebanov-Witten (KW) [22] and the Klebanov-Strassler (KS) [18] models. There are two reasons for this choice:

  • Analysis of the phenomenological model [2] revealed for the first time the exotic thermal phases in a holographic system, which are associated with the spontaneous breaking of a discrete symmetry, and persist to arbitrary high temperature. As TT\to\infty, the fact that the model of [2] was non-conformal becomes irrelevant since for any fixed mass scale m/T0m/T\to 0. In this way one can obtain a holographic thermal conformal order, as emphasized in [3]. Holographic duality between type IIB supergravity on warped deformed conifold with fluxes and the KS cascading gauge theory also reveals the exotic thermal phase [23] — the deconfined phase with spontaneously broken U(1)RU(1)_{R} chiral symmetry, that exists only above certain critical temperature TχSBT_{\chi\rm{SB}}. Like the model in [2], the cascading gauge theory is non-conformal, and has a strong coupling scale Λ\Lambda. It is natural to explore this exotic phase in the limit Λ/T0\Lambda/T\to 0, and potentially obtain a top-down holographic model of the thermal conformal order. We discuss this in section 2.

  • The limit of Λ/T0\Lambda/T\to 0 in the above example effectively removes the fractional D3-branes from the holographic model. On the gravity side one ends with type IIB supergravity on warped deformed conifold with the self-dual five-form flux. The corresponding boundary gauge theory is 𝒩=1{\cal N}=1 superconformal KW model. The gravitational bulk scalars encode the gauge invariant operators (potential conformal order parameters). As we review in section 3, the resulting holographic model (when the conifold is not deformed) is a universal example of AdS/5{}_{5}/CFT4 duality on warped and squashed Sasaki-Einstein manifolds [24]. Thus, analysis of the conformal order on the conifold will cover all such cases, see sections 3.1 and 3.2. In the case of the deformed conifold, we will have potentially the first example of the holographic conformal order with spontaneously broken continuous symmetry, see section 3.3.

We summarized our results in section 4.

2 Exotic phases of the cascading gauge theory at high temperature

Thermodynamics of the Klebanov-Strassler cascading gauge theory [18] has by now a long history [25, 26, 27, 28, 29, 30, 31, 23, 32]. We refer the reader to a recent comprehensive review [32], and focus here on the results only.

Cascading gauge theory is 𝒩=1{\cal N}=1 supersymmetric SU(N+M)×SU(N)SU(N+M)\times SU(N) gauge theory with pairs of chiral multiplets AkA_{k} and BB_{\ell}, k,=1,2,k,\ell=1,2, in the bifundamental (N+M,N¯)(N+M,\overline{N}) and (N+M¯,N)(\overline{N+M},{N}) representations, and the superpotential

𝒲ϵijϵkTr(AiBkAjB).{\cal W}\ \propto\ \epsilon^{ij}\epsilon^{k\ell}\ {\rm Tr}\left(A_{i}B_{k}A_{j}B_{\ell}\right)\,. (2.1)

The theory is not conformal, and the gauge couplings g1g_{1} and g2g_{2}, of the gauge group factors SU(N+M)SU(N+M) and SU(N)SU(N) correspondingly, run with the renormalization group scale μ^\hat{\mu},

ddln(μ^/Λ)8π2g12=3(N+M)2N(1γ),ddln(μ^/Λ)8π2g22=3N2(N+M)(1γ),\begin{split}&\frac{d}{d\ln(\hat{\mu}/\Lambda)}\ \frac{8\pi^{2}}{g_{1}^{2}}=3(N+M)-2N(1-\gamma)\,,\\ &\frac{d}{d\ln(\hat{\mu}/\Lambda)}\ \frac{8\pi^{2}}{g_{2}^{2}}=3N-2(N+M)(1-\gamma)\,,\end{split} (2.2)

where γ\gamma is the anomalous dimension of operators TrAiBj\mathop{\rm Tr}A_{i}B_{j} and Λ\Lambda is the strong coupling scale of the cascading gauge theory. To leading order in M/NM/N, γ=12\gamma=-\frac{1}{2} [22, 18], so that

8π2g128π2g22=6Mlnμ^Λ×(1+𝒪(M/N)),\frac{8\pi^{2}}{g_{1}^{2}}-\frac{8\pi^{2}}{g_{2}^{2}}=6M\ln\frac{\hat{\mu}}{\Lambda}\ \times\ \biggl{(}1+{\cal O}({M}/{N})\biggr{)}\,, (2.3)

while the sum of the gauge couplings is constant along the RG flow

8π2g12+8π2g22=const.\frac{8\pi^{2}}{g_{1}^{2}}+\frac{8\pi^{2}}{g_{2}^{2}}={\rm const}\,. (2.4)
\psfrag{t}[cc][1][0]{${T}/{\Lambda}$}\psfrag{f}[bb][1][0]{$\hat{{\cal F}}$}\psfrag{k}[tt][1][0]{$\ln\frac{{\cal K}_{\chi\rm{SB}}}{{\cal K}}$}\includegraphics[width=216.81pt]{kskts.eps}
\psfrag{t}[cc][1][0]{${T}/{\Lambda}$}\psfrag{f}[bb][1][0]{$\hat{{\cal F}}$}\psfrag{k}[tt][1][0]{$\ln\frac{{\cal K}_{\chi\rm{SB}}}{{\cal K}}$}\includegraphics[width=216.81pt]{Krkskts.eps}
Figure 1: The left panel: the reduced free energy density ^\hat{\cal F} (2.7) of the deconfined chirally symmetric phase (the black curve), and the deconfined (ordered) phase with spontaneously broken chiral symmetry (the magenta curve) of the cascading gauge theory plasma. The vertical dashed red lines indicate TχSBT_{\chi\rm{SB}} (2.8). The ordered phase is exotic: it extends for T>TχSBT>T_{\chi\rm{SB}}. The right panel: the Kretschmann scalar 𝒦{\cal K} for the corresponding phases computed for the dual black branes at the horizon. We use 𝒦χSB=𝒦(T=TχSB){\cal K}_{\chi\rm{SB}}={\cal K}(T=T_{\chi\rm{SB}}).
\psfrag{t}[cc][1][0]{${T}/{\Lambda}$}\psfrag{f}[bb][1][0]{$\hat{{\cal F}}$}\psfrag{k}[tt][1][0]{$\ln\frac{{\cal K}_{\chi\rm{SB}}}{{\cal K}}$}\includegraphics[width=216.81pt]{ktukts.eps}
\psfrag{t}[cc][1][0]{${T}/{\Lambda}$}\psfrag{f}[bb][1][0]{$\hat{{\cal F}}$}\psfrag{k}[tt][1][0]{$\ln\frac{{\cal K}_{\chi\rm{SB}}}{{\cal K}}$}\includegraphics[width=216.81pt]{Krktukts.eps}
Figure 2: The left panel: the reduced free energy density ^\hat{\cal F} (2.7) of the deconfined chirally symmetric phase with the positive specific heat (the black curves), and the deconfined chirally symmetric phase with the negative specific heat (the brown curve) of the cascading gauge theory plasma. The vertical dashed brown lines indicate TuT_{u} (2.11). The negative specific heat phase is exotic: it extends for T>TuT>T_{u} with ever increasing thermal expectation values of certain gauge invariant operators. The right panel: the Kretschmann scalar 𝒦{\cal K} for the corresponding phases computed for the dual black branes at the horizon. As in fig. 1, we normalize the Kretschmann scalar to its value at TχSBT_{\chi\rm{SB}}.

The thermal phase diagram of the theory is rich:
   At large temperatures, TΛT\gg\Lambda, its thermal equation of state is that of a conformal theory with the effective temperature-dependent central charge [25], e.g, the free energy density {\cal F} takes the form

ceff(T)T4,ceffM4ln2TΛ.{\cal F}\ \propto\ -c_{eff}(T)\ T^{4}\,,\qquad c_{eff}\propto M^{4}\ln^{2}\frac{T}{\Lambda}\,. (2.5)

   At

T=Tc=0.614(1)Λ,T=T_{c}=0.614(1)\Lambda\,, (2.6)

the theory undergoes the first-order confinement/deconfinement phase transition [29]; precisely at the transition point the free energy density vanishes,

^26π435M4Λ4|T=Tc=0,\hat{{\cal F}}\equiv\frac{2^{6}\pi^{4}}{3^{5}M^{4}}\ \frac{{\cal F}}{\Lambda^{4}}\bigg{|}_{T=T_{c}}=0\,, (2.7)

where the first equality introduces dimensionless quantity ^\hat{\cal F} we use to describe cascading gauge theory thermodynamics in the canonical ensemble.
   The next critical temperature is [31]

TχSB=0.541(9)Λ,T_{\chi\rm{SB}}=0.541(9)\Lambda\,, (2.8)

represented by the vertical dashed red lines in figs. 1 and 2. For T>TχSBT>T_{{\chi\rm{SB}}} the deconfined phase, represented by the solid black curves in figs. 1 and 2, is perturbatively stable to U(1)R2U(1)_{R}\to{\mathbb{Z}}_{2} chiral symmetry breaking fluctuations. The deconfined states of the cascading gauge theory plasma unstable to chiral symmetry breaking fluctuations are represented by the dashed solid black curves in figs. 1 and 2.
   Similar to the model of [2], the deconfined phase with spontaneously broken chiral symmetry, the solid magenta curves in fig. 1, is exotic: it exists for T>TχSBT>T_{\chi\rm{SB}} and is realized holographically as the Klebanov-Strassler black brane [23]. If we could extend this phase for TΛ\frac{T}{\Lambda}\to\infty, by analogy with [3], we would have realized the conformal order on the conifold. Alas, our numerics allowed the construction of this exotic phase in the range

TΛ=TχSBΛ×{13.9466}.\frac{T}{\Lambda}=\frac{T_{\chi\rm{SB}}}{\Lambda}\ \times\biggl{\{}1\cdots 3.9466\biggr{\}}\,. (2.9)

There is a practical and a conceptual reason for this:

  • from the practical perspective, certain normalizable components of the scalar fields near the boundary become too large for a reliable numerics555In the notations of [32], e.g., it is the coefficient fc,8,0f_{c,8,0} — see (A.55) there.;

  • the conceptual reason causing the above growth is the following: as the temperature increases, the curvature of the dual Klebanov-Strassler black brane evaluated at the horizon grows, and the construction becomes unreliable in the supergravity approximation. Specifically, in the right panel of fig. 1 we present the value of the Kretschmann scalar of the KS black brane horizon, for the exotic phase (the magenta curve) and the deconfined chirally symmetric phase (the black curve), relative to the value of the Kretschmann scalar 𝒦χSB{\cal K}_{{\chi\rm{SB}}} evaluated at T=TχSBT=T_{\chi\rm{SB}}. Over the range (2.9), the Kretschmann scalar corresponding to the exotic phase changes as

    𝒦=𝒦χSB×{18890}.{\cal K}={\cal K}_{{\chi\rm{SB}}}\times\biggl{\{}1\cdots 8890\biggr{\}}\,. (2.10)

    Thus, we conclude that there can not be a reliable conformal order on the warped deformed conifold with fluxes, arising from the TT\to\infty limit of the Klebanov-Strassler black branes — such black branes are outside the validity of the supergravity approximation.

  • Lastly, there is a terminal temperature [30]

    Tu=0.537(3)Λ,T_{u}=0.537(3)\Lambda\,, (2.11)

    represented by the vertical dashed brown lines in fig. 2, for the deconfined chirally symmetric states of the cascading gauge theory plasma — these states exist only for TTuT\geq T_{u}. In the vicinity of TuT_{u}, there are two branches of states: the stable branch with respect to the energy density fluctuations (the black curves in fig. 2), and the unstable branch (the brown curves in fig. 2). On the former branch the specific heat is positive, while it is negative on the latter branch, explaining the (in)stability to the energy density fluctuations [33]. The unstable branch of the deconfined chirally symmetric states is exotic: it extends for T>TuT>T_{u} with the ever increasing values of the thermal expectation values of 𝒪Δ={4,6,8}{\cal O}_{\Delta=\{4,6,8\}} , however, much like in the case of the exotic phase with the chiral symmetry breaking, we fail to extend it as TT\to\infty. We constructed the latter states for fairly narrow temperature range:

    TΛ=TuΛ×{11.0181}.\frac{T}{\Lambda}=\frac{T_{u}}{\Lambda}\times\biggl{\{}1\cdots 1.0181\biggr{\}}\,. (2.12)

    As the temperature of the exotic deconfined chirally symmetric phase increases, the curvature of the corresponding dual Klebanov-Tseytlin black brane evaluated at the horizon rapidly grows, see the right panel of fig. 2. Over the range (2.12), the Kretschmann scalar corresponding to this exotic phase (the brown curve) changes as

    𝒦=𝒦χSB×{1.66973.1}.{\cal K}={\cal K}_{{\chi\rm{SB}}}\times\biggl{\{}1.669\cdots 73.1\biggr{\}}\,. (2.13)

    On the contrary, the deconfined chirally symmetric phase with the positive specific heat (the black curves in fig. 2) can be extended as TΛ\frac{T}{\Lambda}\to\infty — however it is not exotic, as the thermal expectation values of 𝒪Δ={4,6,8}{\cal O}_{\Delta=\{4,6,8\}} operators vanish [29] in this limit, and we end up with the log-dressed conformal equation of state (2.5). Thus, we conclude that there can not be a reliable conformal order on the warped squashed conifold with fluxes, arising from the TT\to\infty limit of the Klebanov-Tseytlin black branes with the negative specific heat — such black branes are outside the validity of the supergravity approximation.

3 Effective theories and the conformal order on the conifold

Consistent truncation in the SU(2)×SU(2)×2SU(2)\times SU(2)\times{\mathbb{Z}}_{2} invariant sector of type IIB supergravity on warped deformed conifold with fluxes to a five dimensional manifold 5{\cal M}_{5} was derived in [31]666See [32] for a recent comprehensive review.:

S5[g^μν,Ωi=13,Φ,hi=13,{P,Ω0}]=10816πG55vol^5Ω1Ω22Ω32×{R1012(^Φ)212eΦ((h1h3)22Ω12Ω22Ω32+1Ω34(^h1)2+1Ω24(^h3)2)12eΦ(2Ω22Ω32(^h2)2+1Ω12Ω24(h2P9)2+1Ω12Ω34h22)12Ω12Ω24Ω34(4Ω0+h2(h3h1)+19Ph1)2}.\begin{split}S_{5}\biggl{[}\hat{g}_{\mu\nu},&\Omega_{i=1\cdots 3},\Phi,h_{i=1\cdots 3}\,,\,\{P,\Omega_{0}\}\biggr{]}=\frac{108}{16\pi G_{5}}\int_{{\cal M}_{5}}\hat{{\rm vol}}_{{\cal M}_{5}}\ \Omega_{1}\Omega_{2}^{2}\Omega_{3}^{2}\\ &\times\biggl{\{}R_{10}-\frac{1}{2}\left(\hat{\nabla}\Phi\right)^{2}-\frac{1}{2}e^{-\Phi}\left(\frac{(h_{1}-h_{3})^{2}}{2\Omega_{1}^{2}\Omega_{2}^{2}\Omega_{3}^{2}}+\frac{1}{\Omega_{3}^{4}}\left(\hat{\nabla}h_{1}\right)^{2}+\frac{1}{\Omega_{2}^{4}}\left(\hat{\nabla}h_{3}\right)^{2}\right)\\ &-\frac{1}{2}e^{\Phi}\left(\frac{2}{\Omega_{2}^{2}\Omega_{3}^{2}}\left(\hat{\nabla}h_{2}\right)^{2}+\frac{1}{\Omega_{1}^{2}\Omega_{2}^{4}}\left(h_{2}-\frac{P}{9}\right)^{2}+\frac{1}{\Omega_{1}^{2}\Omega_{3}^{4}}h_{2}^{2}\right)\\ &-\frac{1}{2\Omega_{1}^{2}\Omega_{2}^{4}\Omega_{3}^{4}}\left(4\Omega_{0}+h_{2}\left(h_{3}-h_{1}\right)+\frac{1}{9}Ph_{1}\right)^{2}\biggr{\}}.\end{split} (3.1)

It is a functional of the a five-dimensional metric g^μν\hat{g}_{\mu\nu} on 5{\cal M}_{5},

ds52=g^μν(y)dyμdyν,ds_{5}^{2}=\hat{g}_{\mu\nu}(y)dy^{\mu}dy^{\nu}\,, (3.2)

scalars Ωi=13\Omega_{i=1\cdots 3} describing the warping and the deformation of the conifold, a dilaton Φ\Phi, scalars hi=13h_{i=1\cdots 3} parameterizing the 3-form fluxes, a constant parameter Ω0\Omega_{0} (necessary to define the self-dual 5-form flux), and a topological parameter PP,

2P9αM,\frac{2P}{9\alpha^{\prime}}\equiv M\ \in\ {\mathbb{Z}}\,, (3.3)

related to the number of fractional D3 branes on the conifold. Finally, R10R_{10} is the Ricci scalar of the ten-dimensional type IIB metric, obtained from uplifting (3.2),

R10=R^5+(12Ω12+2Ω22+2Ω32Ω224Ω12Ω32Ω324Ω12Ω22Ω12Ω22Ω32)2^ln(Ω1Ω22Ω32){(^lnΩ1)2+2(^lnΩ2)2+2(^lnΩ3)2+(^ln(Ω1Ω22Ω32))2},\begin{split}R_{10}=\hat{R}_{5}&+\left(\frac{1}{2\Omega_{1}^{2}}+\frac{2}{\Omega_{2}^{2}}+\frac{2}{\Omega_{3}^{2}}-\frac{\Omega_{2}^{2}}{4\Omega_{1}^{2}\Omega_{3}^{2}}-\frac{\Omega_{3}^{2}}{4\Omega_{1}^{2}\Omega_{2}^{2}}-\frac{\Omega_{1}^{2}}{\Omega_{2}^{2}\Omega_{3}^{2}}\right)-2\hat{\Box}\ln\left(\Omega_{1}\Omega_{2}^{2}\Omega_{3}^{2}\right)\\ &-\biggl{\{}\left(\hat{\nabla}\ln\Omega_{1}\right)^{2}+2\left(\hat{\nabla}\ln\Omega_{2}\right)^{2}+2\left(\hat{\nabla}\ln\Omega_{3}\right)^{2}+\left(\hat{\nabla}\ln\left(\Omega_{1}\Omega_{2}^{2}\Omega_{3}^{2}\right)\right)^{2}\biggr{\}}\,,\end{split} (3.4)

and R^5\hat{R}_{5} is the five-dimensional Ricci scalar of the metric (3.2).

We find it convenient to rewrite the action (3.1) in five-dimensional Einstein frame. The latter is achieved with the following rescaling

g^μνΩ2gμν,Ω3108Ω1Ω22Ω32,\hat{g}_{\mu\nu}\ \to\ \Omega^{2}g_{\mu\nu}\,,\qquad\Omega^{-3}\equiv 108\ \Omega_{1}\Omega_{2}^{2}\Omega_{3}^{2}\,, (3.5)

leading to

108g^Ω1Ω22Ω32R^5=g(R8lnΩ12(lnΩ)2).\begin{split}108\ \sqrt{-\hat{g}}\ \Omega_{1}\Omega_{2}^{2}\Omega_{3}^{2}\ \hat{R}_{5}\ =\ \sqrt{-g}\biggl{(}\ R-8\Box\ln\Omega-12\left(\nabla\ln\Omega\right)^{2}\ \biggr{)}\,.\end{split} (3.6)

Further introducing

Ω1=13ef4w,Ω2=16ef+w+λ,Ω3=16ef+wλ,\Omega_{1}=\frac{1}{3}\ e^{f-4w}\,,\qquad\Omega_{2}=\frac{1}{\sqrt{6}}\ e^{f+w+\lambda}\,,\qquad\Omega_{3}=\frac{1}{\sqrt{6}}\ e^{f+w-\lambda}\,, (3.7)

the five-dimensional effective action becomes

S5=116πG55vol5{R403(f)220(w)24(λ)212(Φ)218e4f4wΦ[e4λ(h1)2+e4λ(h3)2]36e4f4w+Φ(h2)2𝒫flux𝒫scalar},\begin{split}&S_{5}=\frac{1}{16\pi G_{5}}\int_{{\cal M}_{5}}{\rm vol}_{{\cal M}_{5}}\ \biggl{\{}R-\frac{40}{3}\left(\nabla f\right)^{2}-20\left(\nabla w\right)^{2}-4\left(\nabla\lambda\right)^{2}-\frac{1}{2}\left(\nabla\Phi\right)^{2}\\ &-18e^{-4f-4w-\Phi}\biggl{[}e^{4\lambda}\left(\nabla h_{1}\right)^{2}+e^{-4\lambda}\left(\nabla h_{3}\right)^{2}\biggr{]}-36e^{-4f-4w+\Phi}\left(\nabla h_{2}\right)^{2}-{\cal P}_{flux}-{\cal P}_{scalar}\biggr{\}}\,,\end{split} (3.8)

where

𝒫flux=81e283f+4wΦ(h1h3)2+162e283f+4w+Φ[e4λ(h219P)2+e4λh22]+72e403f[h1(P9h2)+9h2h3+36Ω0]2,\begin{split}{\cal P}_{flux}=&81e^{-\frac{28}{3}f+4w-\Phi}(h_{1}-h_{3})^{2}+162e^{-\frac{28}{3}f+4w+\Phi}\biggl{[}e^{-4\lambda}\left(h_{2}-\frac{1}{9}P\right)^{2}+e^{4\lambda}h_{2}^{2}\biggr{]}\\ &+72e^{-\frac{40}{3}f}\biggl{[}h_{1}(P-9h_{2})+9h_{2}h_{3}+36\Omega_{0}\biggr{]}^{2}\,,\end{split} (3.9)
𝒫scalar=4e163f12w24e163f2wcosh(2λ)92e163f+8w(1cosh(4λ)).{\cal P}_{scalar}=4e^{-\frac{16}{3}f-12w}-24e^{-\frac{16}{3}f-2w}\cosh(2\lambda)-\frac{9}{2}e^{-\frac{16}{3}f+8w}\biggl{(}1-\cosh(4\lambda)\biggr{)}\,. (3.10)

Consistent truncation of (3.8), i.e.,

λ=0,h1=h3=1P(K1236Ω0),h2=P18,\lambda=0\,,\qquad h_{1}=h_{3}=\frac{1}{P}\left(\frac{K}{12}-36\Omega_{0}\right)\,,\qquad h_{2}=\frac{P}{18}\,, (3.11)

produces effective action of SU(2)×SU(2)×U(1)SU(2)\times SU(2)\times U(1) invariant sector of type IIB supergravity on warped squashed conifold with fluxes derived in [34].

The boundary holographic dual represented by (3.8) is the 𝒩=1{\cal N}=1 supersymmetric SU(N+M)×SU(N)SU(N+M)\times SU(N) cascading gauge theory, often referred to as a Klebanov-Strassler gauge theory [18]. This theory is not conformal, and has a strong coupling scale Λ\Lambda,

Λ2=2P2gseK0P2gs,\Lambda^{2}=\frac{\sqrt{2}}{P^{2}g_{s}}\ e^{-\frac{K_{0}}{P^{2}g_{s}}}\,, (3.12)

where gsg_{s} is the asymptotic value of the string coupling constant777It can always be fixed to gs=1g_{s}=1., and K0K_{0} is a parameter that can always be adjusted (using the scaling symmetries of the holographic radial coordinate) to a fixed positive value. The precise definition of K0K_{0} can be found in [35, 32]. We are interested here in the conformal holographic model, thus we must take the limit Λ0\Lambda\to 0, which is equivalent to sending P0P\to 0. Finding the conformal order in the model with 7 scalars is a daunting task — so888While it is conceivable that other general constructions of the conformal order might exist in the model, finding them without any guidances is a lost cause. we will follow the framework developed in [5], and reviewed in appendix A.
   As a first step, we truncate the general potential 𝒫flux+𝒫scalar{\cal P}_{flux}+{\cal P}_{scalar} in (3.8) to scalars that enter nonlinearly, and have “unbounded” directions on the field space manifold, at least close to the origin. Notice that hih_{i}, encoding the 3-form fluxes on the conifold, enter quadratically and generate always nonnegative contribution to the scalar potential, i.e., 𝒫flux0{\cal P}_{flux}\geq 0. Thus, in addition to P=0P=0, we set hi0h_{i}\equiv 0. Then,

𝒫flux=72(36Ω0)2e403f=8e403f,{\cal P}_{flux}=72(36\Omega_{0})^{2}\ e^{-\frac{40}{3}f}=8e^{-\frac{40}{3}f}\,, (3.13)

where in the second equality we set a constant parameter Ω0=1108\Omega_{0}=\frac{1}{108} to ensure that the radius of the asymptotically AdS5 geometry is L=1L=1.
   In the absence of 3-form fluxes the dilaton becomes an exact modulus, and we set

Φ0.\Phi\equiv 0\,. (3.14)

   We now arrive at the effective action we use to analyze conformal order on the warped deformed conifold:

S5=cKW2π25vol5{R403(f)220(w)24(λ)2𝒫[f,w,λ]},\begin{split}&S_{5}=\frac{c_{KW}}{2\pi^{2}}\int_{{\cal M}_{5}}{\rm vol}_{{\cal M}_{5}}\ \biggl{\{}R-\frac{40}{3}\left(\nabla f\right)^{2}-20\left(\nabla w\right)^{2}-4\left(\nabla\lambda\right)^{2}-{\cal P}[f,w,\lambda]\biggr{\}}\,,\end{split} (3.15)

where

cKW=27N264,c_{KW}=\frac{27N^{2}}{64}\,, (3.16)

is the central charge of 𝒩=1{\cal N}=1 superconformal SU(N)×SU(N)SU(N)\times SU(N) Klebanov-Witten model [22], and

𝒫=4e163f12w24e163f2wcosh(2λ)92e163f+8w(1cosh(4λ))+8e403f.{\cal P}=4e^{-\frac{16}{3}f-12w}-24e^{-\frac{16}{3}f-2w}\cosh(2\lambda)-\frac{9}{2}e^{-\frac{16}{3}f+8w}\biggl{(}1-\cosh(4\lambda)\biggr{)}+8e^{-\frac{40}{3}f}\,. (3.17)

From (3.15) we obtain the following equations of motion

0=f380𝒫f,0=\Box f-\frac{3}{80}\ \frac{\partial{\cal P}}{\partial f}\,, (3.18)
0=w140𝒫w,0=\Box w-\frac{1}{40}\ \frac{\partial{\cal P}}{\partial w}\,, (3.19)
0=λ18𝒫λ,0=\Box\lambda-\frac{1}{8}\ \frac{\partial{\cal P}}{\partial\lambda}\,, (3.20)
Rμν=403μfνf+20μwνw+4μλνλ+13gμν𝒫.R_{\mu\nu}=\frac{40}{3}\ \partial_{\mu}f\partial_{\nu}f+20\ \partial_{\mu}w\partial_{\nu}w+4\ \partial_{\mu}\lambda\partial_{\nu}\lambda+\frac{1}{3}g_{\mu\nu}\ {\cal P}\,. (3.21)

In the rest of this section we consider three conformal models, obtained by consistent truncations of the effective action S5=S5[f,w,λ]S_{5}=S_{5}[f,w,\lambda] (3.15):

  • Model I:

    S5I[f]S5|w=λ0;S_{5}^{I}[f]\ \equiv\ S_{5}\bigg{|}_{w=\lambda\equiv 0}\,; (3.22)
  • Model II:

    S5II[f,w]S5|λ0;S_{5}^{II}[f,w]\ \equiv\ S_{5}\bigg{|}_{\lambda\equiv 0}\,; (3.23)
  • Model III:

    S5III[f,w,λ]S5.S_{5}^{III}[f,w,\lambda]\ \equiv\ S_{5}\,. (3.24)

3.1 Conformal order in Model I

Conformal Model I realizes holographic dual to the KW gauge theory with a single dimension Δ=8\Delta=8 order parameter 𝒪8{\cal O}_{8}, represented by the bulk scalar ff. Since this scalar represents the warping of the T1,1T^{1,1} base of the singular conifold [21], the model is universal in the sense that any holographic duality between a CFT4 and a type IIB supergravity on AdS×5{}_{5}\timesY5Y_{{}_{5}}, where Y5Y_{5} is a five-dimensional Sasaki-Einstein manifold, can be truncated to our Model I [24]. In the latter case, the scalar ff represents the breathing mode of Y5Y_{5}.

We follow appendix A to construct conformal order in Model I. From (3.15) we find

𝒫I[f]=𝒫[f,0,0]=12+12803f27168027f3+𝒪(f4),{\cal P}_{I}[f]={\cal P}[f,0,0]=-12+\frac{1280}{3}f^{2}-\frac{71680}{27}f^{3}+{\cal O}(f^{4})\,, (3.25)

i.e., the leading nonlinear term is n=3n=3 (see (1.5)) and the unbounded direction is along f+f\to+\infty. The bb-deformed scalar potential takes the form, compare with (1.6),

𝒫Ib[f]=12+12803f2+b(20e163f+8e403f+1212803f2).{\cal P}_{I}^{b}[f]=-12+\frac{1280}{3}f^{2}+b\biggl{(}-20e^{-\frac{16}{3}f}+8e^{-\frac{40}{3}f}+12-\frac{1280}{3}f^{2}\biggr{)}\,. (3.26)

Using the radial coordinate as in (A.10), and introducing

c2(x)=g(x)(2xx2)1/4,c_{2}(x)=\frac{g(x)}{(2x-x^{2})^{1/4}}\,, (3.27)

we obtain the following equations of motion (ddx{}^{\prime}\equiv\frac{d}{dx})

0=f′′+fx1+(12(f)2+9h2209(x22x+2)h40x(1x)(2x)+980x2(2x)2)𝒫Ibf𝒫Ib,\begin{split}&0=f^{\prime\prime}+\frac{f^{\prime}}{x-1}+\biggl{(}-\frac{1}{2}(f^{\prime})^{2}+\frac{9h^{2}}{20}-\frac{9(x^{2}-2x+2)h}{40x(1-x)(2-x)}+\frac{9}{80x^{2}(2-x)^{2}}\biggr{)}\frac{\frac{\partial{\cal P}_{I}^{b}}{\partial f}}{{\cal P}_{I}^{b}}\,,\end{split} (3.28)
0=h4h2+(3x26x+4)hx(1x)(2x)+409(f)2,\begin{split}&0=h^{\prime}-4h^{2}+\frac{(3x^{2}-6x+4)h}{x(1-x)(2-x)}+\frac{40}{9}(f^{\prime})^{2}\,,\end{split} (3.29)
0=ghg.\begin{split}&0=g^{\prime}-h\ g\,.\end{split} (3.30)

Eqs. (3.28)-(3.29) are solved subject to the boundary conditions:
   as x0+x\to 0_{+}, i.e., at the AdS5 boundary

f=f2x2+𝒪(x3),h=329f22x3+𝒪(x4),g=A(189f22x4+𝒪(x5)),f=f_{2}x^{2}+{\cal O}(x^{3})\,,\qquad h=-\frac{32}{9}f_{2}^{2}x^{3}+{\cal O}(x^{4})\,,\qquad g=A\left(1-\frac{8}{9}f_{2}^{2}x^{4}+{\cal O}(x^{5})\right)\,, (3.31)

   as y(1x)0+y\equiv(1-x)\to 0_{+}, i.e., at the horizon

f=f0h+𝒪(y2),h=h1hy+𝒪(y3),g=A(g0h+𝒪(y2)).f=f_{0}^{h}+{\cal O}(y^{2})\,,\qquad h=h_{1}^{h}y+{\cal O}(y^{3})\,,\qquad g=A\biggl{(}g_{0}^{h}+{\cal O}(y^{2})\biggr{)}\,. (3.32)

Note that in total we have four parameters {f2,f0h,h1h,g0h}\{f_{2},f_{0}^{h},h_{1}^{h},g_{0}^{h}\}, along with an arbitrary constant AA. They determine the thermal conformal order of Model I, specifically,

s=2cKWπA3(g0h)3,=sT4,T8𝒪8=f2,T2=A2(g0h)29π2(12h1h)((320(f0h)29)(b1)+3b(5e163f0h2e403f0h)),\begin{split}&s=\frac{2c_{KW}}{\pi}\ A^{3}\left(g_{0}^{h}\right)^{3}\,,\qquad{\cal F}=-\frac{sT}{4}\,,\qquad T^{-8}\langle{\cal O}_{8}\rangle=f_{2}\,,\\ &T^{2}=\frac{A^{2}\left(g_{0}^{h}\right)^{2}}{9\pi^{2}(1-2h_{1}^{h})}\biggl{(}\left(320(f^{h}_{0})^{2}-9\right)(b-1)+3b\left(5e^{-\frac{16}{3}f^{h}_{0}}-2e^{-\frac{40}{3}f^{h}_{0}}\right)\biggr{)}\,,\end{split} (3.33)

for the entropy density ss, the free energy density {\cal F}, the temperature TT, and the thermal order parameter 𝒪8{\cal O}_{8}. Thus, see (1.1),

κI=27[(320(f0h)29)(b1)+3b(5e163f0h2e403f0h)12h11]3/2,γI,8=f2,\kappa_{I}=27\left[\ \frac{\left(320(f^{h}_{0})^{2}-9\right)(b-1)+3b\left(5e^{-\frac{16}{3}f^{h}_{0}}-2e^{-\frac{40}{3}f^{h}_{0}}\right)}{1-2h_{1}^{1}}\ \right]^{-3/2}\,,\ \ \gamma_{I,8}=f_{2}\,, (3.34)

where the subscript I refers to Model I.

\psfrag{g}[cc][1][0]{$1/\gamma_{I,8}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,8}$}\includegraphics[width=216.81pt]{mod1f2l.eps}
\psfrag{g}[cc][1][0]{$1/\gamma_{I,8}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,8}$}\includegraphics[width=216.81pt]{mod1f2s.eps}
Figure 3: Thermal conformal order parameter γ8T8𝒪8\gamma_{8}\equiv T^{-8}\langle{\cal O}_{8}\rangle of the deformed Model I diverges as bbcrit,I+0>1b\to b_{crit,I}+0>1, represented by the dashed vertical red line.
\psfrag{k}[cc][1][0]{$\kappa_{I}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$\kappa_{I}$}\includegraphics[width=216.81pt]{mod1kl.eps}
\psfrag{k}[cc][1][0]{$\kappa_{I}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$\kappa_{I}$}\includegraphics[width=216.81pt]{mod1ks.eps}
Figure 4: The coefficient κI\kappa_{I}, see (1.1), of the thermal order equation of state in the deformed Model I. The vertical dashed red line denotes b=bcrit,Ib=b_{crit,I}.

In the limit b+b\to+\infty we find, see appendix A,

f=9x2(2x)228(1+(1x)3)1b+𝒪(b2),h=𝒪(b2),g=1+𝒪(b2).f=\frac{9x^{2}(2-x)^{2}}{28(1+(1-x)^{3})}\ \frac{1}{b}+{\cal O}\left(b^{-2}\right)\,,\qquad h={\cal O}(b^{-2})\,,\qquad g=1+{\cal O}(b^{-2})\,. (3.35)

Finite bb results are obtained solving (3.28)-(3.30), subject to the asymptotics (3.31) and (3.32). In fig. 3 we plot the expectation value of the order parameter, see (1.1) and (3.33),(3.34). In fig. 4 we plot the coefficient κ\kappa of the thermal order equation of state in Model I as a function of the deformation parameter bb, see (1.1) and (3.34). Notice that the order parameter diverges as bbcrit,Ib\to b_{crit,I},

bcrit,I=1.8370(4),b_{crit,I}=1.8370(4)\,, (3.36)

from above. Since bcrit,I>1b_{crit,I}>1, there is no conformal order in universal top-down holographic Model I, representing a CFT dual to type IIB supergravity with the self-dual five-form flux on a Sasaki-Einstein manifold with a breathing mode, dual to an order parameter of dimension Δ=8\Delta=8. Furthermore, since κI(b)<1\kappa_{I}(b)<1, the conformal order in deformed Model I is subdominant in canonical and microcanonical ensembles. Following [4] we expect that this conformal order is perturbatively unstable.

3.2 Conformal order in Model II

Conformal Model II realizes holographic dual to the KW gauge theory with a pair of order parameters: a dimension Δ=8\Delta=8 operator 𝒪8{\cal O}_{8} and a dimension Δ=6\Delta=6 operator 𝒪6{\cal O}_{6}, represented by the bulk scalars ff and ww correspondingly. Since these scalars represent the warping of the T1,1T^{1,1} base, along with the squashing of the U(1)U(1) fibration of the Kähler-Einstein base of T1,1T^{1,1}, this model is universal as well: any holographic duality between a CFT4 and a type IIB supergravity on AdS×5{}_{5}\timesY~5\tilde{Y}_{{}_{5}}, where Y~5\tilde{Y}_{5} is a five-dimensional squashed Sasaki-Einstein manifold, can be truncated to our Model II [24]. In the latter case, the scalars ff and ww represent the breathing and the squashing modes of Y~5\tilde{Y}_{5}.

We follow appendix A to construct conformal order in Model II. From (3.15) we find

𝒫II[f,w]=𝒫[f,w,0]=12+(12803f2+240w2)+(7168027f31280fw21120w3)+𝒪({f,w}4),\begin{split}{\cal P}_{II}[f,w]={\cal P}[f,w,0]=&-12+\left(\frac{1280}{3}f^{2}+240w^{2}\right)+\biggl{(}-\frac{71680}{27}f^{3}-1280fw^{2}\\ &-1120w^{3}\biggr{)}+{\cal O}(\{f,w\}^{4})\,,\end{split} (3.37)

i.e., the leading nonlinear terms are n=3n=3 (see (1.5)) and the unbounded direction is along {f,w}+\{f,w\}\to+\infty. The bb-deformed scalar potential takes the form,

𝒫IIb[f,w]=12+(12803f2+240w2)+b(𝒫[f,w,0]+1212803f2240w2).{\cal P}_{II}^{b}[f,w]=-12+\biggl{(}\frac{1280}{3}f^{2}+240w^{2}\biggr{)}+b\biggl{(}{\cal P}[f,w,0]+12-\frac{1280}{3}f^{2}-240w^{2}\biggr{)}\,. (3.38)
\psfrag{u}[cc][1][0]{$f_{1,4}$}\psfrag{s}[bb][1][0]{$s$}\psfrag{i}[tt][1][0]{$f_{1,0}^{h}$}\includegraphics[width=216.81pt]{mod2fu.eps}
\psfrag{u}[cc][1][0]{$f_{1,4}$}\psfrag{s}[bb][1][0]{$s$}\psfrag{i}[tt][1][0]{$f_{1,0}^{h}$}\includegraphics[width=216.81pt]{mod2fi.eps}
Figure 5: Coefficients {f1,4,f1,0h}\{f_{1,4},f_{1,0}^{h}\}, see (3.43) and (3.44), of the perturbative conformal order in b+b\to+\infty deformation of Model II, with a parameter ss characterizing the cubic coupling between the scalars ff and gg, see (3.42). The black dot indicates s=0s=0 results reported in appendix A. There are two branches of the conformal order for s<smaxs<s_{max}, represented by the vertical dashed red lines. Since smax<1s_{max}<1, there is no perturbative conformal order in Model II.
\psfrag{u}[cc][1][0]{$w_{1,3}$}\psfrag{s}[bb][1][0]{$s$}\psfrag{i}[tt][1][0]{$w_{1,0}^{h}$}\includegraphics[width=216.81pt]{mod2wu.eps}
\psfrag{u}[cc][1][0]{$w_{1,3}$}\psfrag{s}[bb][1][0]{$s$}\psfrag{i}[tt][1][0]{$w_{1,0}^{h}$}\includegraphics[width=216.81pt]{mod2wi.eps}
Figure 6: Coefficients {w1,3,w1,0h}\{w_{1,3},w_{1,0}^{h}\}, see (3.43) and (3.44), of the perturbative conformal order in b+b\to+\infty deformation of Model II, with a parameter ss characterizing the cubic coupling between the scalars ff and gg, see (3.42). The black dot indicates s=0s=0 results reported in appendix A. There are two branches of the conformal order for s<smaxs<s_{max}, represented by the vertical dashed red lines. Since smax<1s_{max}<1, there is no perturbative conformal order in Model II.

Repeating the analysis as in section 3.1, in the limit b+b\to+\infty, we find

f=1bf1+𝒪(b2),w=1bw1+𝒪(b2),h=𝒪(b2),g=1+𝒪(b2),f=\frac{1}{b}\ f_{1}+{\cal O}(b^{-2})\,,\qquad w=\frac{1}{b}\ w_{1}+{\cal O}(b^{-2})\,,\qquad h={\cal O}(b^{-2})\,,\qquad g=1+{\cal O}(b^{-2})\,, (3.39)

where

0=f1′′+f1x1+4(56f12+9sw126f1)3x2(x2)2,0=f_{1}^{\prime\prime}+\frac{f_{1}^{\prime}}{x-1}+\frac{4\ (56f_{1}^{2}+9s\ w_{1}^{2}-6f_{1})}{3\ x^{2}(x-2)^{2}}\,, (3.40)
0=w1′′+w1x1+w1(16sf1+21w13)x2(x2)2.0=w_{1}^{\prime\prime}+\frac{w_{1}^{\prime}}{x-1}+\frac{w_{1}\ (16s\ f_{1}+21w_{1}-3)}{x^{2}(x-2)^{2}}\,. (3.41)

Notice that we modified the leading nonlinear interactions in (3.37) as

(7168027f31280fw21120w3)(7168027f31280sfw21120w3),\biggl{(}-\frac{71680}{27}f^{3}-1280fw^{2}-1120w^{3}\biggr{)}\ \to\ \biggl{(}-\frac{71680}{27}f^{3}-1280s\ fw^{2}-1120w^{3}\biggr{)}\,, (3.42)

where a constant parameter ss dials the strength of the leading nonlinear coupling between f1f_{1} and w1w_{1} scalars. Ultimately, we need to set s=1s=1, however, it is convenient to start at s=0s=0, so that the scalars f1f_{1} and w1w_{1} are decoupled and we can use the results of appendix A, and then increase s1s\to 1 . Eqs. (3.40)-(3.41) are solved subject to the boundary conditions:
   as x0+x\to 0_{+}, i.e., at the AdS5 boundary

f1=f1,4x2+(f1,434sw1,32)x3+𝒪(x4),w1=w1,3x3/2+34w1,3x5/2+𝒪(x3),f_{1}=f_{1,4}\ x^{2}+\left(f_{1,4}-\frac{3}{4}s\ w_{1,3}^{2}\right)x^{3}+{\cal O}(x^{4})\,,\qquad w_{1}=w_{1,3}\ x^{3/2}+\frac{3}{4}w_{1,3}\ x^{5/2}+{\cal O}(x^{3})\,, (3.43)

   as y(1x)0+y\equiv(1-x)\to 0_{+}, i.e., at the horizon

f1=f1,0h+𝒪(y2),w1=w1,0h+𝒪(y2).f_{1}=f_{1,0}^{h}+{\cal O}(y^{2})\,,\qquad w_{1}=w_{1,0}^{h}+{\cal O}(y^{2})\,. (3.44)

In figs. 5 and 6 we present results for {f1,4,w1,3,f1,0h,w1,0h}\{f_{1,4},w_{1,3},f_{1,0}^{h},w_{1,0}^{h}\} as a function of ss. The black dots represent s=0s=0 results999With appropriate rescaling due to different normalization of the scalar kinetic terms in (3.15) and (A.1). from appendix A for order parameters of dimensions Δ=8\Delta=8 and Δ=6\Delta=6. For any 0<s<smax0<s<s_{max},

smax=0.3956(7),s_{max}=0.3956(7)\,, (3.45)

represented by the vertical red dashed lines, we find two branched of the perturbative conformal order in deformed Model II. However, since smax<1s_{max}<1, we can not construct even a perturbative conformal order in b+b\to+\infty deformation of Model II — recall that the latter requires s=1s=1. We conclude that there is no conformal order (at least within the deformation framework of [5]) in universal top-down holographic Model II, representing a CFT dual to type IIB supergravity with the self-dual five-form flux on a squashed Sasaki-Einstein manifold with a breathing and a squashing modes, dual to the order parameters of dimensions Δ=8\Delta=8 and Δ=6\Delta=6 correspondingly.

3.3 Conformal order in Model III

Conformal Model III realizes holographic dual to the KW gauge theory with a triplet of order parameters: a dimension Δ=8\Delta=8 operator 𝒪8{\cal O}_{8}, a dimension Δ=6\Delta=6 operator 𝒪6{\cal O}_{6} and a dimension Δ=3\Delta=3 operator 𝒪3{\cal O}_{3}, represented by the bulk scalars ff, ww and λ\lambda correspondingly. A noticeable different of Model III from the models discussed in sections 3.1 and 3.2 is that the former ones represent holographic duals to the boundary CFTs with an unbroken U(1)RU(1)_{R} symmetry; while 𝒪30\langle{\cal O}_{3}\rangle\neq 0 would spontaneously break this continuous symmetry as U(1)2U(1)\to{\mathbb{Z}}_{2}.

We follow appendix A to construct conformal order in Model III. From (3.15) we find

𝒫III[f,w,λ]=𝒫=12+(12803f2+240w212λ2)+(7168027f31280fw21120w3+64λ2(f+6w))+𝒪({f,w,λ}4),\begin{split}{\cal P}_{III}[f,w,\lambda]={\cal P}=&-12+\left(\frac{1280}{3}f^{2}+240w^{2}-12\lambda^{2}\right)+\biggl{(}-\frac{71680}{27}f^{3}-1280fw^{2}\\ &-1120w^{3}+64\lambda^{2}(f+6w)\biggr{)}+{\cal O}(\{f,w,\lambda\}^{4})\,,\end{split} (3.46)

i.e., the leading nonlinear terms are n=3n=3 (see (1.5)). The relevant unbounded directions of the truncated potential are less obvious to identify: as in Model II, one would expect for {f,w}\{f,w\} scalars to become large-positive to drive their effective mass (A.19) below the BF bound; while they would be expected to become large-negative for the effective mass of λ\lambda to dip below its effective BF bound near the horizon. As we shortly demonstrate, conformal order in deformed Model III exists when both scalars {f,w}\{f,w\} are negative at the horizon. The bb-deformed scalar potential takes the form,

𝒫IIIb[f,w,λ]=12+(12803f2+240w212λ2)+b(𝒫+1212803f2240w2+12λ2).{\cal P}_{III}^{b}[f,w,\lambda]=-12+\biggl{(}\frac{1280}{3}f^{2}+240w^{2}-12\lambda^{2}\biggr{)}+b\biggl{(}{\cal P}+12-\frac{1280}{3}f^{2}-240w^{2}+12\lambda^{2}\biggr{)}\,. (3.47)

Using the radial coordinate as in (A.10), and introducing

c2(x)=g(x)(2xx2)1/4,c_{2}(x)=\frac{g(x)}{(2x-x^{2})^{1/4}}\,, (3.48)

we obtain the following equations of motion (ddx{}^{\prime}\equiv\frac{d}{dx})

0=f′′+fx1+(12(f)234(w)2320(λ)2+9h2209(x22x+2)h40x(1x)(2x)+980x2(2x)2)𝒫IIIbf𝒫IIIb,\begin{split}0=&f^{\prime\prime}+\frac{f^{\prime}}{x-1}+\biggl{(}-\frac{1}{2}(f^{\prime})^{2}-\frac{3}{4}(w^{\prime})^{2}-\frac{3}{20}(\lambda^{\prime})^{2}+\frac{9h^{2}}{20}-\frac{9(x^{2}-2x+2)h}{40x(1-x)(2-x)}\\ &+\frac{9}{80x^{2}(2-x)^{2}}\biggr{)}\frac{\frac{\partial{\cal P}_{III}^{b}}{\partial f}}{{\cal P}_{III}^{b}}\,,\end{split} (3.49)
0=w′′+wx1+(13(f)212(w)2110(λ)2+3h2103(x22x+2)h20x(1x)(2x)+340x2(2x)2)𝒫IIIbw𝒫IIIb,\begin{split}0=&w^{\prime\prime}+\frac{w^{\prime}}{x-1}+\biggl{(}-\frac{1}{3}(f^{\prime})^{2}-\frac{1}{2}(w^{\prime})^{2}-\frac{1}{10}(\lambda^{\prime})^{2}+\frac{3h^{2}}{10}-\frac{3(x^{2}-2x+2)h}{20x(1-x)(2-x)}\\ &+\frac{3}{40x^{2}(2-x)^{2}}\biggr{)}\frac{\frac{\partial{\cal P}_{III}^{b}}{\partial w}}{{\cal P}_{III}^{b}}\,,\end{split} (3.50)
0=λ′′+λx1+(53(f)252(w)212(λ)2+3h223(x22x+2)h4x(1x)(2x)+38x2(2x)2)𝒫IIIbλ𝒫IIIb,\begin{split}0=&\lambda^{\prime\prime}+\frac{\lambda^{\prime}}{x-1}+\biggl{(}-\frac{5}{3}(f^{\prime})^{2}-\frac{5}{2}(w^{\prime})^{2}-\frac{1}{2}(\lambda^{\prime})^{2}+\frac{3h^{2}}{2}-\frac{3(x^{2}-2x+2)h}{4x(1-x)(2-x)}\\ &+\frac{3}{8x^{2}(2-x)^{2}}\biggr{)}\frac{\frac{\partial{\cal P}_{III}^{b}}{\partial\lambda}}{{\cal P}_{III}^{b}}\,,\end{split} (3.51)
0=h4h2+(3x26x+4)hx(1x)(2x)+409(f)2+203(w)2+43(λ)2,\begin{split}&0=h^{\prime}-4h^{2}+\frac{(3x^{2}-6x+4)h}{x(1-x)(2-x)}+\frac{40}{9}(f^{\prime})^{2}+\frac{20}{3}(w^{\prime})^{2}+\frac{4}{3}(\lambda^{\prime})^{2}\,,\end{split} (3.52)
0=ghg.\begin{split}&0=g^{\prime}-h\ g\,.\end{split} (3.53)

Eqs. (3.49)-(3.52) are solved subject to the boundary conditions:
   as x0+x\to 0_{+}, i.e., at the AdS5 boundary

f=325bλ12x3/2+f4x2+𝒪(x5/2),w=(310bλ12lnx+w3)x3/2+𝒪(x5/2lnx),λ=x1/4(λ1x1/2+38λ1x3/2+𝒪(x2lnx)),h=310λ12x1/238λ12x3/2+𝒪(x2ln2x),g=A(115λ12x3/2320λ12x5/2+𝒪(x3ln2x)),\begin{split}&f=-\frac{3}{25}b\lambda_{1}^{2}\ x^{3/2}+f_{4}\ x^{2}+{\cal O}(x^{5/2})\,,\qquad w=\left(\frac{3}{10}b\lambda_{1}^{2}\ \ln x+w_{3}\right)x^{3/2}+{\cal O}\left(x^{5/2}\ln x\right)\,,\\ &\lambda=x^{1/4}\left(\lambda_{1}\ x^{1/2}+\frac{3}{8}\lambda_{1}\ x^{3/2}+{\cal O}(x^{2}\ln x)\right)\,,h=-\frac{3}{10}\lambda_{1}^{2}\ x^{1/2}-\frac{3}{8}\lambda_{1}^{2}\ x^{3/2}+{\cal O}(x^{2}\ln^{2}x)\,,\\ &g=A\left(1-\frac{1}{5}\lambda_{1}^{2}\ x^{3/2}-\frac{3}{20}\lambda_{1}^{2}\ x^{5/2}+{\cal O}(x^{3}\ln^{2}x)\right)\,,\end{split} (3.54)

   as y(1x)0+y\equiv(1-x)\to 0_{+}, i.e., at the horizon

f=f0h+𝒪(y2),w=w0h+𝒪(y2),λ=λ0h+𝒪(y2),h=h1hy+𝒪(y3),g=A(g0h+𝒪(y2)).\begin{split}&f=f_{0}^{h}+{\cal O}(y^{2})\,,\qquad w=w_{0}^{h}+{\cal O}(y^{2})\,,\qquad\lambda=\lambda_{0}^{h}+{\cal O}(y^{2})\,,\\ &h=h_{1}^{h}y+{\cal O}(y^{3})\,,\qquad g=A\biggl{(}g_{0}^{h}+{\cal O}(y^{2})\biggr{)}\,.\end{split} (3.55)

Note that in total we have eight parameters {f4,f0h,w3,w0h,λ1,λ0h,h1h,g0h}\{f_{4},f_{0}^{h},w_{3},w_{0}^{h},\lambda_{1},\lambda_{0}^{h},h_{1}^{h},g_{0}^{h}\}, along with an arbitrary constant AA. They determine the thermal conformal order of Model III, specifically,

s=2cKWπA3(g0h)3,=sT4,T8𝒪8=f4,T6𝒪6=w3,T3𝒪3=λ1,T2=A2(g0h)2π2(12h1h)[(3209(f0h)2(λ0h)2+20(w0h)21)(b1)+(2e163f0h2w0hcosh(2λ0h)+38e163f0h+8w0h(1cosh(4λ0h))23e403f0h13e163f0h12w0h)×b],\begin{split}&s=\frac{2c_{KW}}{\pi}\ A^{3}\left(g_{0}^{h}\right)^{3}\,,\qquad{\cal F}=-\frac{sT}{4}\,,\\ &T^{-8}\langle{\cal O}_{8}\rangle=f_{4}\,,\qquad T^{-6}\langle{\cal O}_{6}\rangle=w_{3}\,,\qquad T^{-3}\langle{\cal O}_{3}\rangle=\lambda_{1}\,,\\ &T^{2}=\frac{A^{2}\left(g_{0}^{h}\right)^{2}}{\pi^{2}(1-2h^{h}_{1})}\biggl{[}\left(\frac{320}{9}\left(f_{0}^{h}\right)^{2}-\left(\lambda_{0}^{h}\right)^{2}+20\left(w_{0}^{h}\right)^{2}-1\right)(b-1)\\ &+\biggl{(}2e^{-\frac{16}{3}f_{0}^{h}-2w_{0}^{h}}\cosh(2\lambda_{0}^{h})+\frac{3}{8}e^{-\frac{16}{3}f_{0}^{h}+8w_{0}^{h}}\left(1-\cosh(4\lambda_{0}^{h})\right)-\frac{2}{3}e^{-\frac{40}{3}f_{0}^{h}}-\frac{1}{3}e^{-\frac{16}{3}f_{0}^{h}-12w_{0}^{h}}\biggr{)}\\ &\times b\biggr{]}\,,\end{split} (3.56)

for the entropy density ss, the free energy density {\cal F}, the temperature TT, and the thermal order parameters {𝒪8,𝒪6,𝒪3}\{{\cal O}_{8},{\cal O}_{6},{\cal O}_{3}\}. Thus, see (1.1),

κIII=[12h11]3/2×[(3209(f0h)2(λ0h)2+20(w0h)21)(b1)+(2e163f0h2w0hcosh(2λ0h)+38e163f0h+8w0h(1cosh(4λ0h))23e403f0h13e163f0h12w0h)×b]3/2,γIII,{8,6,3}={f4,w3,λ1},\begin{split}&\kappa_{III}=\biggl{[}1-2h_{1}^{1}\biggr{]}^{3/2}\times\biggl{[}\left(\frac{320}{9}\left(f_{0}^{h}\right)^{2}-\left(\lambda_{0}^{h}\right)^{2}+20\left(w_{0}^{h}\right)^{2}-1\right)(b-1)\\ &+\biggl{(}2e^{-\frac{16}{3}f_{0}^{h}-2w_{0}^{h}}\cosh(2\lambda_{0}^{h})+\frac{3}{8}e^{-\frac{16}{3}f_{0}^{h}+8w_{0}^{h}}\left(1-\cosh(4\lambda_{0}^{h})\right)-\frac{2}{3}e^{-\frac{40}{3}f_{0}^{h}}-\frac{1}{3}e^{-\frac{16}{3}f_{0}^{h}-12w_{0}^{h}}\biggr{)}\\ &\times b\biggr{]}^{-3/2}\,,\ \ \gamma_{III,\{8,6,3\}}=\{f_{4},w_{3},\lambda_{1}\}\,,\end{split} (3.57)

where the subscript III refers to Model III.

\psfrag{f}[cc][1][0]{$f_{0}^{h}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{w}[tt][1][0]{$w_{0}^{h}$}\includegraphics[width=216.81pt]{mod3f0h.eps}
\psfrag{f}[cc][1][0]{$f_{0}^{h}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{w}[tt][1][0]{$w_{0}^{h}$}\includegraphics[width=216.81pt]{mod3w0h.eps}
Figure 7: Conformal order in the deformed Model III arises due to negative values (3.55) of ff and ww scalars at the Schwarzschild horizon (blue curves), resulting in the effective mass of λ\lambda-scalar, see (3.60), below the effective BF bound. The green curves represent large-bb approximations to {f0h,w0h}\{f_{0}^{h},w_{0}^{h}\}, see (3.59).
\psfrag{m}[cc][1][0]{$m_{eff,\lambda}^{2}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{w}[tt][1][0]{$w_{0}^{h}$}\includegraphics[width=289.07999pt]{mod3m2.eps}
Figure 8: Effective mass of the λ\lambda-scalar evaluated at the Schwarzschild horizon (blue curve), see (3.60). The red line represents AdS5 BF bound.

In the limit b+b\to+\infty we find, see appendix A,

{f,w,λ}={f1,w1,1}1b+𝒪(b2),h=𝒪(b2),g=1+𝒪(b2),\{f,w,\lambda\}=\{f_{1},w_{1},\ell_{1}\}\ \frac{1}{b}+{\cal O}\left(b^{-2}\right)\,,\qquad h={\cal O}(b^{-2})\,,\qquad g=1+{\cal O}(b^{-2})\,, (3.58)

with

f4=0.27(3)1b+𝒪(b2),f0h=0.008(6)1b+𝒪(b2),w3=0.16(5)1b+𝒪(b2),w0h=0.094(0)1b+𝒪(b2),λ1=0.95(2)1b+𝒪(b2),λ0h=0.59(3)1b+𝒪(b2).\begin{split}&f_{4}=0.27(3)\frac{1}{b}+{\cal O}(b^{-2})\,,\ f_{0}^{h}=-0.008(6)\frac{1}{b}+{\cal O}(b^{-2})\,,\ w_{3}=0.16(5)\frac{1}{b}+{\cal O}(b^{-2})\,,\\ &w_{0}^{h}=-0.094(0)\frac{1}{b}+{\cal O}(b^{-2})\,,\ \lambda_{1}=0.95(2)\frac{1}{b}+{\cal O}(b^{-2})\,,\ \lambda_{0}^{h}=0.59(3)\frac{1}{b}+{\cal O}(b^{-2})\,.\end{split} (3.59)

Finite bb results are obtained solving (3.49)-(3.53), subject to the asymptotics (3.54) and (3.55). Notice that both f0hf_{0}^{h} and w0hw_{0}^{h} are negative as b1b\gg 1; as fig. 7 shows they continue to stay negative for finite bb, driving the effective mass of λ\lambda at the Schwarzschild horizon below the effective BF bound,

meff,λ2|horizon=Δ(Δ4)+16b(f+6w)|horizonΔ=3=3+16b(f0h+6w0h)=12.16(1)+𝒪(b1),\begin{split}m_{eff,\lambda}^{2}\bigg{|}_{horizon}=&\Delta(\Delta-4)+16b(f+6w)\bigg{|}_{horizon}^{\Delta=3}=-3+16b\left(f_{0}^{h}+6w_{0}^{h}\right)\\ =&-12.16(1)+{\cal O}(b^{-1})\,,\end{split} (3.60)

and causing the condensation of all the scalars, see fig. 8.

\psfrag{g}[cc][1][0]{$1/\gamma_{I,8}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,8}$}\includegraphics[width=216.81pt]{mod3fl.eps}
\psfrag{g}[cc][1][0]{$1/\gamma_{I,8}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,8}$}\includegraphics[width=216.81pt]{mod3fs.eps}
Figure 9: Thermal conformal order parameter γ8T8𝒪8\gamma_{8}\equiv T^{-8}\langle{\cal O}_{8}\rangle of the deformed Model III diverges as bbcrit,III+0b\to b_{crit,III}+0, represented by the dashed vertical red line.
\psfrag{g}[cc][1][0]{$1/\gamma_{I,6}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,6}$}\includegraphics[width=216.81pt]{mod3wl.eps}
\psfrag{g}[cc][1][0]{$1/\gamma_{I,6}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,6}$}\includegraphics[width=216.81pt]{mod3ws.eps}
Figure 10: Thermal conformal order parameter γ6T6𝒪6\gamma_{6}\equiv T^{-6}\langle{\cal O}_{6}\rangle of the deformed Model III diverges as bbcrit,III+0b\to b_{crit,III}+0, represented by the dashed vertical red line.
\psfrag{g}[cc][1][0]{$1/\gamma_{I,3}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,3}$}\includegraphics[width=216.81pt]{mod3ll.eps}
\psfrag{g}[cc][1][0]{$1/\gamma_{I,3}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$1/\gamma_{I,3}$}\includegraphics[width=216.81pt]{mod3ls.eps}
Figure 11: Thermal conformal order parameter γ3T3𝒪3\gamma_{3}\equiv T^{-3}\langle{\cal O}_{3}\rangle of the deformed Model III diverges as bbcrit,III+0b\to b_{crit,III}+0, represented by the dashed vertical red line.
\psfrag{k}[cc][1][0]{$\kappa_{III}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$\kappa_{III}$}\includegraphics[width=216.81pt]{mod3kl.eps}
\psfrag{k}[cc][1][0]{$\kappa_{III}$}\psfrag{b}[bb][1][0]{$b$}\psfrag{i}[tt][1][0]{$\kappa_{III}$}\includegraphics[width=216.81pt]{mod3ks.eps}
Figure 12: The coefficient κIII\kappa_{III}, see (1.1), of the thermal order equation of state in the deformed Model III. The vertical dashed red line denotes b=bcrit,IIIb=b_{crit,III}.

In figs. 9-11 we plot the expectation values of the order parameters γIII,{8,6,3}\gamma_{III,\{8,6,3\}}, see (1.1) and (3.33),(3.57). In fig. 12 we plot the coefficient κ\kappa of the thermal order equation of state in Model III as a function of the deformation parameter bb, see (1.1) and (3.57). We obtained reliable numerical results for b1.03b\gtrsim 1.03; their extrapolation to smaller values of bb predicts that all the order parameters diverge as bbcrit,IIIb\to b_{crit,III},

bcrit,III=0.99(7)|extrapolation,b_{crit,III}=0.99(7)\bigg{|}_{\rm extrapolation}\,, (3.61)

from above. We take closeness of bcrit,IIIb_{crit,III} to 11 as the strong indication that the actual value of bcrit,IIIb_{crit,III} is precisely 11. We thus conclude that while there is a conformal order in bb-deformed Model III, it disappears once this model becomes a realization of the top-down holography, much like in top-down holographic models discussed in [5, 9]. Since κIII(b)<1\kappa_{III}(b)<1, the conformal order in deformed Model III is subdominant in canonical and microcanonical ensembles. Following [4] we expect that this conformal order is perturbatively unstable.

4 Conclusion

In this paper we searched for the holographic conformal order [3] on warped deformed conifold with fluxes [32] — representing the supergravity dual to Klebanov-Strassler cascading gauge theory [18]. We focused on two exotic phases [2] of the black branes on the conifold:

  • the Klebanov-Strassler black branes [23], realizing the spontaneous breaking of the continuous RR-symmetry;

  • the branch of Klebanov-Tseytlin black branes with the negative specific heat [30].

The two classes of the black branes have an intrinsic scale, the holographic dual to the strong coupling scale Λ\Lambda of the boundary cascading gauge theory. Thus, to find the conformal order one needs to construct these black branes in the high temperature limit TΛ0\frac{T}{\Lambda}\to 0. We showed that as one increases the temperature, these exotic black branes on the conifold become ever more stringy. As a result, the conformal order can not be reached within the controllable supergravity approximation.

In the second part of the paper we changed the strategy: instead of starting with the non-conformal exotic black branes and taking the high-temperature limit, we took the conformal limit in the effective action of type IIB supergravity on warped deformed conifold with fluxes [31]. The latter limit effectively removes the fractional D3 branes from the model. The resulting five-dimensional effective action has seven bulk scalars; but only three of these scalars have the nonlinear potential. The nonlinearity of the scalar potential appears to be crucial in holographic examples of the conformal order constructed so far [2, 3, 5, 9]. Taking this fact as a hint, we consistently truncated the “conformal” effective action to these three scalars. Within this action, we further identified two additional universal consistent truncations:

  • type IIB supergravity on warped Sasaki-Einstein manifolds with the self-dual five-form flux [30, 24];

  • type IIB supergravity on warped and squashed Sasaki-Einstein manifolds with the self-dual five-form flux [30, 24].

We showed that in both truncations one can construct the conformal order, provided one deforms the bulk scalar potential obtained from the Kaluza-Klein reduction on the Sasaki-Einstein manifold. However, before the deformation parameter is removed, the conformal order disappears. The same story repeats in the full effective action on the warped deformed conifold with the three bulk scalars — the conformal order exists, as long as the scalar potential is (arbitrarily small) deformed from the one obtained from the type IIB supergravity Kaluza-Klein reduction. Along with the previous results [5, 9] we believe this sends a powerful message: there is a ’real’ holography and a ’toy’ one; conformal order is possible only in the latter. This begs a question: how much one can trust phenomenological holographic models to predict String Theory Universe phenomena?

Acknowledgments

This research is supported in part by Perimeter Institute for Theoretical Physics. Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Colleges and Universities. This work was further supported by NSERC through the Discovery Grants program.

Appendix A Perturbative holographic thermal conformal order

Consider thermal conformal order in the phenomenological model

S5=c2π25vol5{R12(ϕ)2𝒫b[ϕ]},S_{5}=\frac{c}{2\pi^{2}}\int_{{\cal M}_{5}}{\rm vol}_{{\cal M}_{5}}\biggl{\{}R-\frac{1}{2}(\nabla\phi)^{2}-{\cal P}^{b}[\phi]\biggr{\}}\,, (A.1)

with the scalar potential 𝒫b[ϕ]{\cal P}^{b}[\phi] given by (1.6), i.e.,

𝒫b=12+12Δ(Δ4)ϕ2bk=n3pkϕk,{\cal P}^{b}=-12+\frac{1}{2}\Delta(\Delta-4)\phi^{2}-b\sum_{k=n\geq 3}^{\infty}p_{k}\ \phi^{k}\,, (A.2)

where pn>0p_{n}>0, and the following next nonzero coefficient pkp_{k} is for k=n+pn+1k=n+p\geq n+1. We assume Δ2\Delta\geq 2. To this end, we assume the black brane metric ansatz

ds5=2=c12dt2+c22d𝒙2+c32dr2,ci=ci(r),ds_{5}=^{2}=-c_{1}^{2}\ dt^{2}+c_{2}^{2}\ d\bm{x}^{2}+c_{3}^{2}\ dr^{2}\,,\qquad c_{i}=c_{i}(r)\,, (A.3)

along with ϕ=ϕ(r)\phi=\phi(r), all depending on the radial coordinate r[r0,+)r\in[r_{0},+\infty). There is the smooth Schwarzschild horizon as rr0+r\to r_{0_{+}}, i.e.,

limrr0+c12=0,\lim_{r\to r_{0_{+}}}c_{1}^{2}=0\,, (A.4)

and asymptotic AdS5 geometry, i.e., as rr\to\infty,

c12r2,c2r2,c32r2,ϕ𝒪ΔrΔ,c_{1}^{2}\to r^{2}\,,\quad c_{2}\to r^{2}\,,\quad c_{3}^{2}\to r^{-2}\,,\quad\phi\to\frac{\langle{\cal O}_{\Delta}\rangle}{r^{\Delta}}\,, (A.5)

where 𝒪Δ0\langle{\cal O}_{\Delta}\rangle\neq 0 is the order parameter, see (1.1).

From the effective action (A.1) we derive the following equations of motion

0=c1′′+c1[lnc23c3]+13c1c32𝒫b,0=c_{1}^{\prime\prime}+c_{1}^{\prime}\ \left[\ln\frac{c_{2}^{3}}{c_{3}}\right]^{\prime}+\frac{1}{3}c_{1}c_{3}^{2}\ {\cal P}^{b}\,, (A.6)
0=c2′′+c2[lnc1c22c3]+13c2c32𝒫b,0=c_{2}^{\prime\prime}+c_{2}^{\prime}\ \left[\ln\frac{c_{1}c_{2}^{2}}{c_{3}}\right]^{\prime}+\frac{1}{3}c_{2}c_{3}^{2}\ {\cal P}^{b}\,, (A.7)
0=(ϕ)212[lnc2][ln(c1c2)]2c32𝒫b,0=(\phi^{\prime})^{2}-12\left[\ln c_{2}\right]^{\prime}\left[\ln(c_{1}c_{2})\right]^{\prime}-2c_{3}^{2}\ {\cal P}^{b}\,, (A.8)
0=ϕ′′+ϕ[lnc1c23c3]c32𝒫bϕ,0=\phi^{\prime\prime}+\phi^{\prime}\ \left[\ln\frac{c_{1}c_{2}^{3}}{c_{3}}\right]^{\prime}-c_{3}^{2}\ \frac{\partial{\cal P}^{b}}{\partial\phi}\,, (A.9)

where ddr{}^{\prime}\equiv\frac{d}{dr}. Equations (A.6)-(A.9) can be systematically analyzed perturbatively as b+b\to+\infty. Introducing a new radial coordinate [36],

1x=c1c2,x(0,1],1-x=\frac{c_{1}}{c_{2}}\,,\qquad x\in(0,1]\,, (A.10)

with x1x\to 1_{-} representing the regular Schwarzschild horizon (A.4), and x0+x\to 0_{+} the boundary asymptotes (A.5), we find

c1=(1x)c2,c2=πT(2xx2)1/4[1+𝒪((1b)2n2)],c32dr2=dx24(2xx2)2[1+𝒪((1b)2n2)],ϕ=(12pnb)1n2[ϕ1(x)+𝒪((1b)2n2,(1b)pn2)],\begin{split}&c_{1}=(1-x)c_{2}\,,\qquad c_{2}=\frac{\pi T}{(2x-x^{2})^{1/4}}\biggl{[}1+{\cal O}\left(\left(\frac{1}{b}\right)^{\frac{2}{n-2}}\right)\biggr{]}\,,\\ &c_{3}^{2}\ dr^{2}=\frac{dx^{2}}{4(2x-x^{2})^{2}}\biggl{[}1+{\cal O}\left(\left(\frac{1}{b}\right)^{\frac{2}{n-2}}\right)\biggr{]}\,,\\ &\phi=\left(\frac{1}{2p_{n}b}\right)^{\frac{1}{n-2}}\biggl{[}\phi_{1}(x)+{\cal O}\left(\left(\frac{1}{b}\right)^{\frac{2}{n-2}},\left(\frac{1}{b}\right)^{\frac{p}{n-2}}\right)\biggr{]}\,,\end{split} (A.11)

where TT is the Hawking temperature of the horizon, and the leading as b+b\to+\infty scalar hair ϕ1\phi_{1} satisfies

0=ϕ1′′+ϕ1x1+ϕ1(nϕ12Δ(Δ4))8x2(2x)2,0=\phi_{1}^{\prime\prime}+\frac{\phi_{1}^{\prime}}{x-1}+\frac{\phi_{1}(n\phi_{1}-2\Delta(\Delta-4))}{8x^{2}(2-x)^{2}}\,, (A.12)

subject to the boundary conditions:
   as x0+x\to 0_{+}, i.e., at the AdS5 boundary

ϕ1f1,0xΔ/4TΔ𝒪Δf1,0(1b)1n2,\phi_{1}\to f_{1,0}\ x^{\Delta/4}\qquad\Longrightarrow\qquad T^{-\Delta}\langle{\cal O}_{\Delta}\rangle\ \propto\ f_{1,0}\ \left(\frac{1}{b}\right)^{\frac{1}{n-2}}\,, (A.13)

   as x1x\to 1_{-}, i.e., at the horizon

ϕ1f1,h.\phi_{1}\to f_{1,h}\,. (A.14)

Once the perturbative (as b+b\to+\infty) solution (A.11) is constructed, the finite-bb conformal order can be analyzed numerically, incrementally decreasing the deformation parameter bb [3, 5, 9].

The reason why we expect the specific large-bb scaling of the conformal order as in (A.11) was given in [5]:

  • Recall the story of the holographic superconductor [37, 38]. A scalar field in asymptotically AdS geometry must have a mass above the (space-time dependent) Breitenlohner-Freedman (BF) bound to avoid the condensation. In the vicinity of the Schwarzschild horizon, the BF bound can be modified either by changing the effective dimensionality of the space-time (as in the extremal limit of a Reissner-Nordstrom black brane) [38], or via nonlinear scalar interactions, leading to large negative contribution of the effective mass [37]. Thus, it is possible for a scalar field to be above the BF bound close to the AdS boundary, and below the effective BF bound close to the horizon. This scenario triggers the condensation of the scalar - the black brane develops the scalar hair. Close to the transition point, the condensation can be studied in the probe approximation — neglecting the backreaction.

  • The conformal order realizes the black brane horizon scalarization mechanism of [37]. It becomes a probe approximation in the limit b+b\to+\infty. Indeed, with the scaling

    ϕ(1b)1n2,n3,\phi\ \propto\ \left(\frac{1}{b}\right)^{\frac{1}{n-2}}\,,\qquad n\geq 3\,, (A.15)

    the scalar backreaction on the AdS5-Schwarzschild geometry vanishes as in (A.11). Furthermore, an order knk\geq n monomial in the scalar potential (A.2) scales as

    bpkϕkpk(1b)2n2+knn2,b\ p_{k}\ \phi^{k}\ \propto\ p_{k}\ \left(\frac{1}{b}\right)^{\frac{2}{n-2}+\frac{k-n}{n-2}}\,, (A.16)

    while the mass term in the scalar potential (A.2), as well as the kinetic term, scale as

    (ϕ)2Δ(Δ4)ϕ2(1b)2n2.(\nabla\phi)^{2}\ \sim\ \Delta(\Delta-4)\ \phi^{2}\ \propto\ \left(\frac{1}{b}\right)^{\frac{2}{n-2}}\,. (A.17)

    Thus, any nonlinear term in the scalar potential with k>nk>n is subleading in the b+b\to+\infty limit compare to the quadratic scalar terms in the effective action (A.1), while the leading nonlinear k=nk=n term scales precisely as the former.

  • To summarize, given (A.15), the holographic model (A.1) represents the probe approximation of the scalar ϕ\phi on AdS5-Schwarzschild geometry with the effective action:

    Sscalar=5=AdS5Schwarzschildvol5{12(ϕ)212meff2ϕ2}5=AdS5Schwarzschildvol5{12(ϕ1)212meff2ϕ12},\begin{split}S_{scalar}=&\int_{{\cal M}_{5}={\rm AdS}_{5}-{\rm Schwarzschild}}{{\rm vol}}_{{\cal M}_{5}}\biggl{\{}-\frac{1}{2}(\nabla\phi)^{2}-\frac{1}{2}m_{eff}^{2}\ \phi^{2}\biggr{\}}\\ \propto&\int_{{\cal M}_{5}={\rm AdS}_{5}-{\rm Schwarzschild}}{\rm vol}_{{\cal M}_{5}}\biggl{\{}-\frac{1}{2}(\nabla\phi_{1})^{2}-\frac{1}{2}m_{eff}^{2}\ \phi_{1}^{2}\biggr{\}}\,,\end{split} (A.18)

where

meff2=Δ(Δ4)2pnbϕn2=Δ(Δ4)ϕ1n2,\begin{split}m_{eff}^{2}=\Delta(\Delta-4)-2p_{n}b\ \phi^{n-2}=\Delta(\Delta-4)-\phi_{1}^{n-2}\,,\end{split} (A.19)

and the scalar field ϕ1\phi_{1} satisfies (A.12). Notice from (A.19) that a large positive value of ϕ1\phi_{1} at the horizon, see (A.14), can make meff2m_{eff}^{2} sufficiently negative, and trigger the scalarization as in [37].

\psfrag{d}[cc][1][0]{$\Delta$}\psfrag{u}[bb][1][0]{$f_{1,0}$}\psfrag{i}[tt][1][0]{$f_{1,h}$}\includegraphics[width=216.81pt]{uvdelta.eps}
\psfrag{d}[cc][1][0]{$\Delta$}\psfrag{u}[bb][1][0]{$f_{1,0}$}\psfrag{i}[tt][1][0]{$f_{1,h}$}\includegraphics[width=216.81pt]{irdelta.eps}
Figure 13: Normalizable coefficients, (A.13) and (A.14), of the bulk scalar field dual to the conformal order parameter 𝒪Δ{\cal O}_{\Delta} in the limit b+b\to+\infty. The red squares represent the analytical results obtained for Δ=4\Delta=4 and Δ=8\Delta=8 correspondingly.

We now return to the analysis of (A.12)-(A.14). We consider101010Generalizations to other nn and Δ\Delta are straightforward. n=3n=3, as this case will be of relevance for the discussion in section 3, and consider integer values of Δ={2,10}\Delta=\{2,\cdots 10\}. Numerical results for {f1,0,f1,h}\{f_{1,0},f_{1,h}\} are presented in fig. 13; for Δ=4\Delta=4 and Δ=8\Delta=8 (A.12)-(A.14) can be solved analytically:

Δ=4:ϕ1=323x(2x){f1,0,f1,h}={643,323},Δ=8:ϕ1=64x2(2x)2(1+(1x)2)2{f1,0,f1,h}={64,64}.\begin{split}&\Delta=4:\qquad\phi_{1}=\frac{32}{3}\ x(2-x)\qquad\Longrightarrow\qquad\{f_{1,0},f_{1,h}\}=\{\frac{64}{3},\frac{32}{3}\}\,,\\ &\Delta=8:\qquad\phi_{1}=64\ \frac{x^{2}(2-x)^{2}}{(1+(1-x)^{2})^{2}}\qquad\Longrightarrow\qquad\{f_{1,0},f_{1,h}\}=\{{64},{64}\}\,.\end{split} (A.20)

References