Top-Down Holographic $G$-Structure Glueball Spectroscopy at (N)LO in $N$ and Finite Gauge Coupling

TL;DR

This paper computes glueball spectra in a top-down holographic QCD model at finite gauge coupling, including NLO corrections, filling a gap in the literature and validating semi-universality of corrections.

Contribution

It provides the first detailed calculation of glueball masses at finite gauge coupling with NLO corrections in a top-down holographic QCD setup.

Findings

Glueball masses computed at finite gauge coupling.

NLO corrections show a specific suppression pattern.

Validation of semi-universality of NLO corrections.

Abstract

The top-down type IIB holographic dual of large-$N$ thermal QCD as constructed in arXiv:0902.1540 involving a fluxed resolved warped deformed conifold, its delocalized type IIA SYZ mirror as well as its M-theory uplift constructed in arXiv:1306.4339 - both in the finite gauge coupling ($g_s\stackrel{<}{\sim}1$)/`MQGP' limit of arXiv:1306.4339 - were shown explicitly to possess a local $SU(3)/G_2$-structure in arXiv:1507.02692 Glueball spectroscopy at finite gauge coupling has thus far been missing in the literature. In this paper, we fill this gap by calculating the masses of the $0^{++}, 0^{-+},0^{--}, 1^{++}, 2^{++}$ (`glueball') states (which correspond to fluctuations in the dilaton or complexified two-forms or appropriate metric components) in the aforementioned backgrounds of $G$-structure in the `MQGP' limit of arXiv:1306.4339, using WKB quantization conditions on one hand and imposing Neumann/Dirichlet boundary conditions at an IR cut-off/horizon radius $r_h$ on the solutions to the equations of motion on the other. We also discuss $r_h=0$-limits of all calculations; in this context we also calculate the $0^{++}, 0^{--},1^{++}, 2^{++}$ glueball masses up to NLO in $N$ and find a $\frac{g_sM^2}{N}(g_sN_f)$-suppression similar to and further validating a similar semi-universality of NLO corrections to transport coefficients, observed in arXiv:1606.04949.

Tables

Table 1. Table 1: Summary of Glueball Spectra: m = m ~ ​ r h L 2 𝑚 ~ 𝑚 subscript 𝑟 ℎ superscript 𝐿 2 m=\tilde{m}\frac{r_{h}}{L^{2}} from Type IIB, IIA and M Theory using WKB quantization condition for r h ≠ 0 subscript 𝑟 ℎ 0 r_{h}\neq 0 , and m = m ~ ​ r 0 L 2 𝑚 ~ 𝑚 subscript 𝑟 0 superscript 𝐿 2 m=\tilde{m}\frac{r_{0}}{L^{2}} for r h = 0 subscript 𝑟 ℎ 0 r_{h}=0 (equalities in the r h = 0 subscript 𝑟 ℎ 0 r_{h}=0 column, are valid up to NLO in N 𝑁 N ); the colored triangles/square in the third column correspond to the colored triangles/square that appear in Fig. 2 in the combined plot of r h ≠ 0 subscript 𝑟 ℎ 0 r_{h}\neq 0 supergravity calculations of glueballs

S. No.	Glueball	$\tilde{m}$ using WKB $r_h\ne 0$	$\tilde{m}$ using WKB $r_h=0$
		(units of $\pi T$, up to LO in $N$)	(units of $\frac{r_0}{L^2}$, up to NLO in $N$)
		(large-$\tilde{m}$ limit)	(large-$\tilde{m}$ limit)
1	$0^{++}$	(M theory)	(M theory)
	(Fluctuations: $h_{00,r{x}_1,rr}$	$\frac{\sqrt{35+70n}}{\sqrt{𝒫ℒ\left(13.15+26.31n\right)}}$	No turning points
	in M-theory metric)
	$0^{++}$	(Type IIB)	(Type IIB)
	(Dilaton Fluctuations)	$9.18+8.08n$	$\frac{\left(4.64+6.28n\right)}{\left({\delta }^2-1\right)}\left[1-0.01\frac{g_s{M}^2}{N}\left(g_s{N}_f\right)\log N\log r_0\right]$
2	$0^{-+}$	(Type IIA)	(Type IIA)
	(1-form fluctuation $a_{{\theta }_2}$)	$11.12\left(n+\frac{1}{2}\right),n=0$	$\frac{3.72+4.36n}{\left({\delta }^2-1\right)},n=0$
		$\left(6.22+4.80n\right),n\in {\mathbb{Z}}^{+}$	$4.8\left(n+\frac{1}{2}\right),n\in {\mathbb{Z}}^{+}$
3	$0^{--}$	(Type IIB)	(Type IIB)
	2-form fluctuation $A_{23}$	$=m_n^{0^{++}}(\mathrm{dilaton},T)$	$\frac{6.28n+4.71}{\left({\delta }^2-1\right)}\left(1+\frac{0.01g_s^{}{}_{}{}^2\log N{M}^2{N}_f\log \left(r_0\right)}{N}\right),n=0$
			$\left(7.87+6.93n\right),n\in {\mathbb{Z}}^{+}$
4	$1^{++}$	(M theory)	(M theory)
	(Fluctuations: $h_{it}=h_{ti},i=x_{2,3}$	$8.08\left(n+\frac{1}{2}\right)$	$m_n^{1^{++}}\left(n=0,r_h=0\right)=m_n^{0^{--}}\left(n=0,r_h=0\right)$
	in M-theory metric)		$\left(3.46+6.93n\right),n\in {\mathbb{Z}}^{+}$
5	$2^{++}$	(M theory)	(M theory)
	(Fluctuations: $h_{x_2{x}_3}=h_{x_3{x}_2},$	$8.08\left(n+\frac{1}{2}\right)=m_n^{1^{++}}\left(T\right)$	$=m_n^{1^{++}}\left(r_h=0\right)$
	$h_{x_2{x}_2}=h_{x_3{x}_3}$ in M-theory metric)
	$2^{++}$	(Type IIB)	(Type IIB)
	(Fluctuation $h_{x_2{x}_3}=h_{x_3{x}_2}$	$9.18+8.08n=m_n^{0^{++}}(IIB,T)$	$=m_n^{1^{++}}\left(r_h=0\right)$
	in type IIB metric)

Table 2. Table 2: Comparison of [ 57 ] ’s N → ∞ → 𝑁 N\rightarrow\infty lattice results for 0 + + superscript 0 absent 0^{++} glueball with our supergravity results obtained using WKB quantization condition and redefinitions of [ 34 ] for M theory scalar metric fluctuations

State	$N\to \infty$ Entry in Table 34 of [57]	M-theory scalar metric perturbations	Type IIB Dilaton fluctuations of [61]
	in units of square root of	(6.1.2 - in units of	in units of reciprocal of
	string tension	$\frac{r_h}{L^2}$)	temporal circle’s diameter
$0^{++}$	$4.065\pm 0.055$	4.267	4.07 (normalized to match lattice)
$0^{++\ast }$	$6.18\pm 0.13$	6.251	7.02
$0^{++\ast \ast }$	$7.99\pm 0.22$	7.555	9.92
$0^{++\ast \ast }$	-	8.588	12.80
$0^{++\ast \ast \ast }$	-	9.464	15.67

Table 3. Table 3: Summary of Glueball Spectra from Type IIB, IIA and M Theory for r h ≠ 0 / r h = 0 subscript 𝑟 ℎ 0 subscript 𝑟 ℎ 0 r_{h}\neq 0/r_{h}=0 using Neumann/Dirichlet boundary conditions at the horizon r h subscript 𝑟 ℎ r_{h} /IR cut-off r 0 subscript 𝑟 0 r_{0}

S. No.	Glueball	Spectrum Using	Spectrum Using
		N(eumann)/D(irichlet)	N(eumann)/D(irichlet)
		b.c., $r=r_h$(units of $\pi T$)	b.c., $r=r_0$(units of $\frac{r_0}{L^2}$)
1	$0^{++}$	(M theory)	(M theory)
		(N) $12.25\sqrt{2+n}$	(N) 4.1
		(D) $12.25\sqrt{1+n}$
2	$0^{-+}$	(Type IIA)	(Type IIA)
		(N/D) $\frac{3.1}{\pi }\sqrt{n}$	(N) $m_{n=0}^{0^{-+}}=0,m_{n=1}^{0^{-+}}\approx 3.4,m_{n=2}^{0^{-+}}\approx 4.35$
			(D) $m_{n=0}^{0^{-+}}=0,m_{n=1}^{0^{-+}}\approx 3.06,m_{n=2}^{0^{-+}}\approx 4.81$
3	$0^{--}$	(Type IIB)	(Type IIB)
		(N/D) $m_{n=0}^{0^{--}}\left(T\right)=0,m_{n=1}^{0^{--}}\left(T\right)=\frac{32.46}{\pi },$	(large $n$)
		$m_{n=2}^{0^{--}}\left(T\right)=\frac{32.88}{\pi }$	(N/D)
			$\frac{1}{2}{5}^{3/4}\sqrt[4]{\frac{2\left(\sqrt{6}\sqrt{{\pi }^2\left(16n^2+22n+7\right)+6}+6\right)+3{\pi }^2\left(2n+1\right)}{32-3{\pi }^2}}$
4	$1^{++}$	(M theory)	(M theory)
		(N/D) $m_{n=0}^{1^{++}}\left(T\right)=2.6956,m_{n=1}^{1^{++}}\left(T\right)=2.8995$	(N) $m_{n=0}^{1^{++}}\left(r_h=0\right)\approx 1.137$
		$m_{n=2}^{1^{++}}\left(T\right)=2.9346$	(D) $m_{n=0}^{1^{++}}\left(r_h=0\right)\approx 0.665$
5	$2^{++}$	(M theory)	(M theory)
		(N) $m_{n=0}^{2^{++}}\left(T\right)=\frac{5.086}{\pi },m_{n=1}^{2^{++}}\left(T\right)=\frac{5.269}{\pi }$	$=m_n^{1^{++}}\left(r_h=0\right)$
		$m_{n=2}^{2^{++}}\left(T\right)=\frac{5.318}{\pi }$
		$m_{n=0}^{2^{++}}(D,T)=0,m_{n+1}^{2^{++}}(D,T)=m_n^{2^{++}}(N,T)$

Table 4. Table 4: Comparison of ratios of 0 − − superscript 0 absent 0^{--} glueball masses obtained from Neumann/Dirichlet boundary conditions at the horizon, with [ 61 ]

Ratio	Our Results	[61]’s Results
$\frac{m_{0^{--}}^{\ast }}{m_{0^{--}}}$	1.0129	1.5311
$\frac{m_{0^{--}}^{\ast \ast }}{m_{0^{--}}^{\ast }}$	1.0033	1.3244
$\frac{m_{0^{--}}^{\ast \ast \ast }}{m_{0^{--}}^{\ast \ast }}$	1.0013	1.2393
$\frac{m_{0^{--}}^{\ast \ast \ast \ast }}{m_{0^{--}}^{\ast \ast \ast }}$	1.0007	1.1588

Equations

$$4{\pi }^2\left(\frac{1}{g_{SU(N+M)}^2}+\frac{1}{g_{SU(N)}^2}\right)e^{\varphi }\sim \pi ;4{\pi }^2\left(\frac{1}{g_{SU(N+M)}^2}-\frac{1}{g_{SU(N)}^2}\right)e^{\varphi }\sim \frac{1}{2\pi {\alpha }^{\prime }}{\int }_{S^2}{B}_2.$$

$$N_{\mathrm{eff}}(r)=N\left[1+\frac{3g_s{M}_{\mathrm{eff}}^2}{2\pi N}\left(\log r+\frac{3g_s{N}_f^{\mathrm{eff}}}{2\pi }{\left(\log r\right)}^2\right)\right],$$

$$M_{\mathrm{eff}}(r)=M+\frac{3g_s{N}_{f}M}{2\pi }\log r+\sum _{m\ge 1}\sum _{n\ge 1}{N}_f^m{M}^n{f}_{mn}(r),$$

$$N_f^{\mathrm{eff}}(r)=N_f+\sum _{m\ge 1}\sum _{n\ge 0}{N}_f^m{M}^n{g}_{mn}(r).$$

$$L_{\mathrm{OKS}-\mathrm{BH}}\sim \sqrt{M}{N}_f^{\frac{3}{4}}\sqrt{\left(\sum _{m\ge 0}\sum _{n\ge 0}{N}_f^m{M}^n{f}_{mn}(\Lambda )\right)}{\left(\sum _{l\ge 0}\sum _{p\ge 0}{N}_f^l{M}^p{g}_{lp}(\Lambda )\right)}^{\frac{1}{4}}{g}_s^{\frac{1}{4}}\sqrt{{\alpha }^{\prime }}$$

$$\equiv N_f^{\frac{3}{4}}{\sqrt{\left(\sum _{m\ge 0}\sum _{n\ge 0}{N}_f^m{M}^n{f}_{mn}(\Lambda )\right)}{\left(\sum _{l\ge 0}\sum _{p\ge 0}{N}_f^l{M}^p{g}_{lp}(\Lambda )\right)}^{\frac{1}{4}}{L}_{\mathrm{KS}}|}_{\Lambda :\log \Lambda <\frac{2\pi }{3g_s{N}_f}},$$

$$ds^{2}=\frac{1}{\sqrt{h}}\left(-g_{1}dt^{2}+dx_{1}^{2}+dx_{2}^{2}+dx_{3}^{2}\right)+\sqrt{h}\biggl{[}g_{2}^{-1}dr^{2}+r^{2}d{\cal M}_{5}^{2}\biggr{]}.$$

$$d{\mathcal{M}}_5^2=h_1(d\psi +\cos {\theta }_{1}d{\varphi }_1+\cos {\theta }_{2}d{\varphi }_2{)}^2+h_2(d{\theta }_1^2+{\sin }^2{\theta }_{1}d{\varphi }_1^2)+$$

$$+h_4(h_{3}d{\theta }_2^2+{\sin }^2{\theta }_{2}d{\varphi }_2^2)+h_5\cos \psi \left(d{\theta }_{1}d{\theta }_2-\sin {\theta }_1\sin {\theta }_{2}d{\varphi }_{1}d{\varphi }_2\right)+$$

$$+h_5\sin \psi \left(\sin {\theta }_{1}d{\theta }_{2}d{\varphi }_1+\sin {\theta }_{2}d{\theta }_{1}d{\varphi }_2\right),$$

$$h_1=\frac{1}{9}+\mathcal{O}\left(\frac{g_s{M}^2}{N}\right),h_2=\frac{1}{6}+\mathcal{O}\left(\frac{g_s{M}^2}{N}\right),h_4=h_2+\frac{a^2}{r^2},$$

$$h_3=1+\mathcal{O}\left(\frac{g_s{M}^2}{N}\right),h_5\ne 0,L={\left(4\pi g_{s}N\right)}^{\frac{1}{4}}.$$

$$\displaystyle\hskip-32.52127pth=\frac{L^{4}}{r^{4}}\Bigg{[}1+\frac{3g_{s}M_{\rm eff}^{2}}{2\pi N}{\rm log}r\left\{1+\frac{3g_{s}N^{\rm eff}_{f}}{2\pi}\left({\rm log}r+\frac{1}{2}\right)+\frac{g_{s}N^{\rm eff}_{f}}{4\pi}{\rm log}\left({\rm sin}\frac{\theta_{1}}{2}{\rm sin}\frac{\theta_{2}}{2}\right)\right\}\Biggr{]},$$

$$h=\frac{L^4}{r^4}\left[1+\sum _{i=1}\frac{{\mathcal{H}}_i\left({\varphi }_{1,2},{\theta }_{1,2},\psi \right)}{r^i}\right],\mathrm{large}r;$$

$$h=\frac{L^4}{r^4}\left[1+\sum _{i,j;(i,j)\ne (0,0)}\frac{h_{ij}\left({\varphi }_{1,2},{\theta }_{1,2},\psi \right){\log }^{i}r}{r^j}\right],\mathrm{small}r.$$

$$(a){\tilde{F}}_3=2M{A}_1\left(1+\frac{3g_s{N}_f}{2\pi }\log r\right)e_{\psi }\wedge \frac{1}{2}\left(\sin {\theta }_{1}d{\theta }_1\wedge d{\varphi }_1-B_1\sin {\theta }_{2}d{\theta }_2\wedge d{\varphi }_2\right)$$

$$-\frac{3g_{s}M{N}_f}{4\pi }{A}_2\frac{dr}{r}\wedge e_{\psi }\wedge \left(\cot \frac{{\theta }_2}{2}\sin {\theta }_{2}d{\varphi }_2-B_2\cot \frac{{\theta }_1}{2}\sin {\theta }_{1}d{\varphi }_1\right)$$

$$-\frac{3g_{s}M{N}_f}{8\pi }{A}_3\sin {\theta }_1\sin {\theta }_2\left(\cot \frac{{\theta }_2}{2}d{\theta }_1+B_3\cot \frac{{\theta }_1}{2}d{\theta }_2\right)\wedge d{\varphi }_1\wedge d{\varphi }_2,$$

$$\displaystyle\hskip-28.90755pt(b)H_{3}={6g_{s}{A_{4}}M}\Biggl{(}1+\frac{9g_{s}N_{f}}{4\pi}~{}{\rm log}~{}r+\frac{g_{s}N_{f}}{2\pi}~{}{\rm log}~{}{\rm sin}\frac{\theta_{1}}{2}~{}{\rm sin}\frac{\theta_{2}}{2}\Biggr{)}\frac{dr}{r}$$

$$\displaystyle\hskip-21.68121pt\wedge\frac{1}{2}\Biggl{(}{\rm sin}~{}\theta_{1}~{}d\theta_{1}\wedge d\phi_{1}-{B_{4}}~{}{\rm sin}~{}\theta_{2}~{}d\theta_{2}\wedge d\phi_{2}\Biggr{)}+\frac{3g^{2}_{s}MN_{f}}{8\pi}{A_{5}}\Biggl{(}\frac{dr}{r}\wedge e_{\psi}-\frac{1}{2}de_{\psi}\Biggr{)}$$

$$\displaystyle\wedge\Biggl{(}{\rm cot}~{}\frac{\theta_{2}}{2}~{}d\theta_{2}-{B_{5}}~{}{\rm cot}~{}\frac{\theta_{1}}{2}~{}d\theta_{1}\Biggr{)}.$$

$$(i)\mathrm{weak}(g_s)\mathrm{coupling}-\mathrm{large}{\mathrm{t}}^{\prime }\mathrm{Hooft}\mathrm{coupling}\mathrm{limit}:$$

$$g_s\ll 1,g_s{N}_f\ll 1,\frac{g_s{M}^2}{N}\ll 1,g_{s}M\gg 1,g_{s}N\gg 1$$

$$\mathrm{effected}\mathrm{by}:g_s\sim {ϵ}^d,M\sim {\left(\mathcal{O}(1)ϵ\right)}^{-\frac{3d}{2}},N\sim {\left(\mathcal{O}(1)ϵ\right)}^{-19d},ϵ\ll 1,d>0$$

$$\displaystyle\hskip-12.28577pt(ii){\rm MQGP\ limit}:\frac{g_{s}M^{2}}{N}\ll 1,g_{s}N\gg 1,{\rm finite}\ g_{s},M\$$

$$\mathrm{effected}\mathrm{by}:g_s\sim {ϵ}^d,M\sim {\left(\mathcal{O}(1)ϵ\right)}^{-\frac{3d}{2}},N\sim {\left(\mathcal{O}(1)ϵ\right)}^{-39d},ϵ\lesssim 1,d>0.$$

$$x=\sqrt{h_2}{h}^{\frac{1}{4}}sin⟨{\theta }_1⟩⟨r⟩{\varphi }_1,y=\sqrt{h_4}{h}^{\frac{1}{4}}sin⟨{\theta }_2⟩⟨r⟩{\varphi }_2,z=\sqrt{h_1}⟨r⟩h^{\frac{1}{4}}\psi ,$$

$${i^{\ast }J|}_{RC/DC}\approx 0,$$

$${\Im m\left(i^{\ast }\Omega \right)|}_{RC/DC}\approx 0,$$

$${\Re e\left(i^{\ast }\Omega \right)|}_{RC/DC}\sim \mathrm{volume}\mathrm{form}\left(T^3(x,y,z)\right),$$

$$\left(\begin{array}[]{c}sin\theta_{2}d\phi_{2}\\ d\theta_{2}\end{array}\right)\rightarrow\left(\begin{array}[]{cc}cos\langle\psi\rangle&sin\langle\psi\rangle\\ -sin\langle\psi\rangle&cos\langle\psi\rangle\end{array}\right)\left(\begin{array}[]{c}sin\theta_{2}d\phi_{2}\\ d\theta_{2}\end{array}\right),$$

$$d\psi \to d\psi +f_1({\theta }_1)\cos {\theta }_{1}d{\theta }_1+f_2({\theta }_2)\cos {\theta }_{2}d{\theta }_2,$$

$$d{\varphi }_{1,2}\to d{\varphi }_{1,2}-f_{1,2}({\theta }_{1,2})d{\theta }_{1,2},$$

$${\Lambda }^1\otimes su(3{)}^{\perp }=(3\oplus \bar{3})\otimes (1\oplus 3\oplus \overline{3)}$$

$$=(1\oplus 1)\oplus (8\oplus 8)\oplus (6\oplus \bar{6})\oplus (3\oplus \bar{3})\oplus (3\oplus \bar{3}{)}^{\prime }\equiv W_1\oplus W_2\oplus W_3\oplus W_4\oplus W_5.$$

$$J=e^1\wedge e^2+e^3\wedge e^4+e^5\wedge e^6,$$

$$\Omega =(e^1+i{e}^2)\wedge (e^3+i{e}^4)\wedge (e^5+i{e}^6).$$

$$T_{SU(3)}^{\mathrm{IIB}}\in W_1\oplus W_2\oplus W_3\oplus W_4\oplus W_5\sim \frac{e^{-3\tau }}{\sqrt{g_{s}N}}\oplus {\left(g_{s}N\right)}^{\frac{1}{4}}{e}^{-3\tau }\oplus \sqrt{g_{s}N}{e}^{-3\tau }\oplus -\frac{2}{3}\oplus -\frac{1}{2}$$

$$\frac{2}{3}{W}_5^{\bar{3}}=W_4^{\bar{3}}$$

$$T_{SU(3)}^{\mathrm{IIA}}$$

Full text

**Top-Down Holographic $G$-Structure Glueball Spectroscopy at (N)LO in $N$ and Finite Coupling

**Karunava Sil111e-mail: krusldph@iitr.ac.in, Vikas Yadav222email: viitr.dph2015@iitr.ac.in and Aalok Misra333e-mail: aalokfph@iitr.ac.in

Department of Physics, Indian Institute of Technology, Roorkee - 247 667, Uttaranchal, India

The top-down type IIB holographic dual of large-$N$ thermal QCD as constructed in [1] involving a fluxed resolved warped deformed conifold, its delocalized type IIA SYZ mirror as well as its M-theory uplift constructed in [2] - both in the finite coupling ($g_{s}\stackrel{{\scriptstyle<}}{{\sim}}1$)/‘MQGP’ limit of [2] - were shown explicitly to possess a local $SU(3)/G_{2}$-structure in [3]. Glueballs spectra in the finite-gauge-coupling limit (and not just large-t’Hooft coupling limit) - a limit expected to be directly relevant to strongly coupled systems at finite temperature such as QGP [4] - has thus far been missing in the literature. In this paper, we fill this gap by calculating the masses of the $0^{++},0^{-+},0^{--},1^{++},2^{++}$ (‘glueball’) states (which correspond to fluctuations in the dilaton or complexified two-forms or appropriate metric components) in the aforementioned backgrounds of $G$-structure in the ‘MQGP’ limit of [2]. We use WKB quantization conditions on one hand and impose Neumann/Dirichlet boundary conditions at an IR cut-off (‘$r_{0}$’)/horizon radius (‘$r_{h}$’) on the solutions to the equations of motion on the other. We find that the former technique produces results closer to the lattice results. We also discuss $r_{h}=0$-limits of all calculations. In this context we also calculate the $0^{++},0^{--},1^{++},2^{++}$ glueball masses up to NLO in $N$ and find a $\frac{g_{s}M^{2}}{N}(g_{s}N_{f})$-suppression similar to and further validating a similar semi-universality of NLO corrections to transport coefficients, observed in [5].

1 Introduction

The AdS/CFT correspondence [6] remarkably establishes an equivalence between the partition functions of a five dimensional gravitational theory (bulk theory) and a four dimensional supersymmetric and scale invariant gauge theory (boundary theory). A generalization of the AdS/CFT correspondence is necessary to explore more realistic gauge theories (less supersymmetric and non-conformal) such as QCD with $SU(3)$ gauge group. The top-down model that we have considered in this work is motivated to capture QCD-like gauge theories from a suitable gravitational background. QCD is a strongly coupled theory at low energies. The low energy dynamics of QCD involves the color neutral bound states of gluons, known as glueballs. Hence, the non-perturbative aspects of QCD can be largely understood from the glueball sector of the theory. Moreover, the plasma phase of QCD (QGP) occurs at high temperatures $T>T_{c}$. In QGP medium the quarks and the gluons stay in a deconfined state due to Debye screening. However, the recent RHIC experiments indicate strongly that non-perturbative effects of QCD are present in the plasma phase. In fact the lattice results of [7] conclude that QGP must be non-perturbative in the temperature regime $T_{c}\leq T\leq 5T_{c}$. This is precisely the reason why we concentrate on the glueball spectra in the finite-gauge-coupling limit (and not just large-t’Hooft coupling limit) - a limit expected to be directly relevant to strongly coupled systems at finite temperature such as QGP [4].

QCD is a non-abelian gauge theory, in which gauge fields play the role of dynamical degrees of freedom. Non-abelian nature of QCD allows the gauge bosons to form color-neutral bound states of gluons known as glueballs (gg, ggg, etc.). Therefore, the study of glueballs and their spectra enables us to gain a better understanding of the non-perturbative regime of QCD. Glueball state is represented by quantum numbers $\text{J}^{\text{PC}}$, where $J,P$ and $C$ correspond to total angular momentum, parity and charge conjugation respectively

Different generalized versions of the AdS/CFT correspondence has thus far been proposed to study non-supersymmetric field theories with a running gauge coupling constant. The original proposal was given by Witten to obtain a gravity dual for non-supersymmetric field theories. As per Witten’s formalism, non-supersymmetric Yang-Mills theory can be obtained by compactifying one of the spatial direction on a circle and imposing antiperiodic boundary conditions on the fermions around this circle. This makes the fermions and scalars massive and they get decoupled leaving only gauge fields as degrees of freedom. The gravity dual of this compactified theory was asymptotically AdS. In a particular case of ${\cal N}=4$ SU(N) super-Yang-Mills theory dual to type IIB string theory on $AdS_{5}\times S^{5}$, the procedure described above gives an effective model of three dimensional Yang-Mills theory, i.e., $\text{QCD}_{3}$.

The gravity dual of non-supersymmetric theories in the low energy limit is typically given by supergravity backgrounds involving $\text{AdS}_{p}\times\text{M}_{q}$ where AdS is the anti de Sitter space with dimension $p$ and $\text{M}_{q}$ is the internal manifold with dimension $q$. In the supergravity theory the Klauza-Klein modes on $\text{M}_{q}$ can be classified according to the spherical harmonics of the $\text{M}_{q}$, which forms representations of the isometry group of $\text{M}_{q}$. The states carrying the non-trivial isometry group quantum numbers are heavier and do not couple to the pure gluonic operators on the boundary. Thus the glueballs are identified with singlet states of the isometry group.

In the past decade, glueballs have been studied extensively to gain new insight into the non-perturbative regime of QCD. Various holographic setups such as soft-wall model, hard wall model, modified soft wall model, etc. have been used to obtain the glueball’s spectra. In [8, 9] a soft wall holographic model was used to study the glueball spectrum. In [8] glueballs and scalar mesons were studied at finite temperature. It was found that the masses of the hadronic states decreases and the widths become broader as T increases. But for temperature range of the order of 40-60 MeV, states disappear from the glueball and meson spectral function. Both hard wall and soft wall holographic models were considered [10, 11] to obtain the glueball correlation functions to study the dynamics of QCD. Decay rates were obtained for glueballs in both models. Dynamical content of the correlators was investigated[11] by obtaining their spectral density and relating it with various other quantities to obtain the estimates for three lowest dimensional gluon condensate. In [12] a two-flavor quenched dynamical holographic QCD model was considered with two different forms for the dilaton field given as: $\Phi=\mu_{G}^{2}z^{2}$ and $\Phi=\mu_{G}^{2}z^{2}\tanh(\mu_{G}^{4}z^{2}/\mu_{G}^{2}),$ ($z$ being a radial coordinate).In [13] an AdS5 mass renormalization was implemented in a modified holographic softwall model to obtain the spectrum of scalar and higher even glueball spin states with $P=C=+1$. In [14, 15] a bottom-up approach was used to obtain the mass spectra of the scalar and vector glueballs. In this case, the Vector glueball masses were found to be heavier than thatof the scalar glueballs while higher values for both were reportedin other approaches. In [16] holographic description was used for supersymmetric and non-supersymmetric, non-commutative dipole gauge theory in 4D. WKB approximation was used to obtain the mass by solving the dilaton and antisymmetric tensor field equations in the bulk. For the supersymmetric theory, dipole length plays the role of an intrinsic scale while for non-supersymmetric theory the same role is played by the temperature Two different phases for baryons were found, a big baryon dual to the static string and a small baryon dual to a moving string. In [17] spectrum for scalar, vector and tensor two-gluon and trigluon glueballs were obtained in 5-D holographic QCD model with a metric structure deformed by the dilaton field. The spectrum was compared with the results obtained from both soft-wall and hard-wall holograpic QCD models. Here, the spectra of the two-gluon glueball was found to be in agreement with the lattice data. For trigluon glueballs, the masses for $1^{\pm-}$,$2^{--}$ were matched while masses for $0^{--},0^{+-}$ and $2^{+-}$ were lighter than lattice data which indicates that the latter glueballs are dominated by three-gluon condensate. In [18] holographic glueball spectrum was obtained in the singlet sector of ${\cal N}$=1 supersymmetric Klebanov-Strassler model. States containing the bifundamental $A_{i}$ and $B_{i}$ fields were not considered. Comparison with the lattice data showed a nice agreement for $1^{+-}$ and $1^{--}$ states while $0^{+-}$ results were different because of its fermionic component.

Glueballs appear in the meson spectra of QCD and difficulty in their identification in the meson spectra is largely due to lack of information about their coupling with mesons in strongly coupled QCD. Lattice QCD gives an estimate for the masses but it does not give any information about the glueball couplings and their decay widths both of which are required for identification of glueballs. Holographic approach gives a better understanding of glueball decay rates than lattice QCD. Various holographic models such as Witten-Sakai-Sugimoto model, Soft wall model and supersymmetric Klebanov-Strassler model, etc. have been used to obtain the coupling between mesons and glueballs to obtain expressions for the glueball decay widths.

In [19, 20, 21, 22, 23, 24] top-down Witten-Sakai-Sugimoto model was used as a holographic setup for low energy QCD to obtain the coupling of scalar glueballs to mesons and subsequently obtain their decay widths. Results obtained were compared with the experimental data available for lattice counterparts f0(1500) and f0(1710) of scalar glueballs. In [19] results obtained for decay widths and branching ratios for scalar glueball decays were found to be consistent with experimental data for f0(1500) state in [19] while in [21] results favored f0(1710) as scalar glueball candidate instead of f0(1500). Decay patterns were obtained for scalar glueball candidate f0(1710) in top-down holographic Witten-Sakai-Sugimoto model for low energy QCD in [20]. It was shown that there exists a narrow pseudoscalar glueball heavier than the scalar glueball whose decay pattern involves $\eta$ and $\eta\prime$ mesons. In [22, 23, 24] Witten-Sakai-Sugimoto model was used to study the phenomenology of scalar glueball states. A dilaton and an exotic mode were obtained as two sets of scalar glueball states in [22, 23]. Calculation of mass spectra showed that out of two modes, dilaton mass is quite close to both f0(1710) and f0(1500) scalar glueball candidates while calculation of decay width showed that f0(1710) is the favored glueball candidate corresponding to dilaton mode. In [25] the holographic top-down Witten-Sakai-Sugimoto model was used to study the tensor $2^{++}$ glueball mass spectrum and decay width . Decay width was found to be above 1 GeV for glueball mass MT=2400 MeV while for MT=2000 MeV it was reduced to 640 MeV. In [26] modified holographic soft-wall model was used to calculate the mass spectrum and Regge trajectories of lightest scalar glueball and higher spin glueball states. Results were obtained for both even and odd spins glueball states.

In this paper, we use a large-$N$ top-down holographic dual of QCD to obtain the spin $2^{++},1^{++},$ $0^{++},0^{--},0^{-+}$ glueball spectrum explicitly for QCD3 from type IIB, type IIA and M theory perspectives. Now for the computation of the glueball masses, we need to introduce a scale in our theory. In other words, the conformal invariance has to be broken. This can be done in two different ways. The first approach, after Witten [27], corresponds to the compactification of the time direction on a circle of finite radius, forming a black hole in the background. In this case the masses are determined in units of the horizon radius $r_{h}$ of the black hole. The other approach is to consider a cut-off at $r=r_{0}$ in the gravitational background ($r$ being the non-compact radial direction)[28]. This forbids the arbitrary low energy excitations of the boundary field theory and hence breaks the conformal invariance. So, in this case the required scale to address strong interaction is introduced by the IR cut-off $r_{0}$. From a top-down perspective this IR cut-off will in fact be proportional to two-third power of the Ouyang embedding parameter obtained from the minimum radial distance (corresponding to the lightest quarks) requiring one to be at the South Poles in the $\theta_{1,2}$ coordinates, in the holomorphic Ouyang embedding of flavor $D7$-branes. In the spirit of [27], the time direction for both cases will be compact with fermions obeying anti-periodic boundary conditions along this compact direction, and hence we will be evaluating three-dimensional glueball masses.

Glueball masses can be obtained by evaluating the correlation functions of gauge invariant local operator. The first step to obtain the glueball spectrum in QCD3 is to identify the operators in the gauge theory that have quantum number corresponding to the glueballs of interest. According to the gauge/gravity duality each supergravity mode corresponds to a gauge theory operator. This operator couples to the supergravity mode at the boundary of the AdS space, for example, the lowest dimension operator with quantum numbers $J^{PC}=0^{++}$ is Tr$F^{2}=$Tr$F_{\mu\nu}F^{\mu\nu}$ and this operator couples to the dilaton mode on the boundary. To calculate $0^{++}$ glueball mass we need to evaluate the correlator $\left<TrF^{2}(x)TrF^{2}(y)\right>$ =$\Sigma_{i}c_{i}e^{-m{i}|x-y|}$, where $m_{i}$ give the value for glueball mass. However the masses can also be obtained by solving the wave equations for supergravity modes which couples to the gauge theory operators on the boundary. The latter approach is used in this paper.

The 11D metric obtained as the uplift of the delocalized SYZ type IIA metric, up to LO in $N$, can be interpreted as a black $M5$-brane wrapping a two-cycle, i.e. a black $M3$-brane [29, 3]. Taking this as the starting point, compactifying again along the M-theory circle, we land up at the type IIA metric and then compactifying again along the periodic temporal circle (with the radius given by the reciprocal of the temperature), one obtains ${\rm QCD}_{3}$ corresponding to the three non-compact directions of the black $M3$-brane world volume. The Type IIB background of [1], in principle, involves $M_{4}\times$ RWDC($\equiv$ Resolved Warped Deformed Conifold); asymptotically the same becomes $AdS_{5}\times T^{1,1}$. To determine the gauge theory fields that would couple to appropriate supergravity fields a la gauge-gravity duality, ideally one should work the same out for the $M_{4}\times$ RWDC background (which would also involve solving for the Laplace equation for the internal RWDC). We do not attempt to do the same here. Motivated however by, e.g.,

(a) asymptotically the type IIB background of [1] and its delocalized type IIA mirror of [2] consist of $AdS_{5}$

and

(b) terms of the type $Tr(F^{2}(AB)^{k}),(F^{4}(AB)^{k})$ where $F^{2}=F_{\mu\nu}F^{\mu\nu},F^{4}=F_{\mu_{1}}^{\ \mu_{2}}F_{\mu_{2}}^{\ \mu_{3}}F_{\mu_{3}}^{\ \mu_{4}}F_{\mu_{4}}^{\ \mu_{1}}-\frac{1}{4}\left(F_{\mu_{1}}^{\ \mu_{2}}F_{\mu_{2}}^{\ \mu_{1}}\right)^{2}$, $A,B$ being the bifundamental fields that appear in the gauge theory superpotential corresponding to $AdS_{5}\times T^{1,1}$ in [30], form part of the gauge theory operators corresponding to the solution to the Laplace equation on $T^{1,1}$ [31] (the operator $TrF^{2}$ which shares the quantum numbers of the $0^{++}$ glueball couples to the dilaton and $TrF^{4}$ which also shares the quantum numbers of the $0^{++}$ glueball couples to trace of metric fluctuations and the four-form potential, both in the internal angular directions),

we calculate in this paper:

$\bullet$ type IIB dilaton fluctuations, which we refer to as $0^{++}$ glueball

$\bullet$ type IIB complexified two-form fluctuations that couple to $d^{abc}Tr(F_{\mu\rho}^{a}F^{b\ \rho}_{\lambda}F^{c\ \lambda}_{\ \ \ \ \ \nu})$, which we refer to as $0^{--}$ glueball

$\bullet$ type IIA one-form fluctuations that couple to $Tr(F\wedge F)$, which we refer to as $0^{-+}$ glueball

$\bullet$ M-theory metric’s scalar fluctuations which we refer to as another (lighter) $0^{++}$ glueball

$\bullet$ M-theory metric’s vector fluctuations which we refer to as $1^{++}$ glueball,

and

$\bullet$ M-theory metric’s tensor fluctuations which we refer to as $2^{++}$ glueball.

All holographic glueball spectra calculations done thus far, have only considered a large t’Hooft coupling limit: $g_{\rm YM}^{2}N\gg 1,\ N\gg 1$. However, holographic duals of thermal QCD laboratories like sQGP also require a finite gauge coupling [4]. This was addressed as part of the ‘MQGP limit’ in [2]. It is in this regard that results of this paper - which discusses supergravity glueball spectra at finite string coupling - are particularly significant. Also, the recent observation - see, e.g., [32] - that the non-perturbative properties of quark gluon plasma can be related to the change of properties of scalar and pseudoscalar glueballs, makes the study of glueballs quite important.

The rest of the paper is organized as follows. In Sec. 2, via five sub-sections, we summarize the top-down type IIB holographic dual of large-$N$ thermal QCD of [1], its delocalized SYZ type IIA mirror and its M theory uplift of [2, 3]. In Sec. 3, we discuss a supergravity calculation of the spectrum of $0^{++}$ glueball at finite horizon radius $r_{h}$ (3.1) and setting $r_{h}=0$ (3.2). The $r_{h}\neq 0$ computations are given in 3.1, corresponding to use of WKB quantization conditions using coordinate/field redefinitions of [34]. The $r_{h}=0$ calculations are subdivided into sub-section 3.2.1 corresponding to solving the $0^{++}$ equation of motion up to LO in $N$ and imposing Neumann/Dirichlet boundary condition at the horizon, and sub-section 3.2.2 corresponding to WKB quantization conditions inclusive of non-conformal/NLO-in-$N$ corrections using the redefinitions of [34]. Sec. 4 has to do with the $0^{-+}$ glueball spectrum. Further therein, 4.1.1 and 4.1.2 respectively are on obtaining the $r_{h}\neq 0$ spectrum and its $r_{h}=0$ limit using Neumann/Dirichlet boundary conditions on the solutions up to LO in $N$ respectively at the horizon and the IR cutoff. Then subsections 4.2.1 and 4.2.2 respectively are on WKB quantization at finite and zero $r_{h}$ up to LO in $N$, using redefinitions of [34]. Sec. 5 is on $0^{--}$ glueball spectrum. Therein, subsections 5.1 and 5.2 are on getting the spectrum by imposing Neumann/Dirichlet boundary condition on the solutions up to LO in $N$ to the EOM respectively at the horizon and the IR cut-off. Subsection 5.3 has to do with obtaining the spectrum using WKB quantization up to LO in $N$ at $r_{h}\neq 0$ using the redefinitions of [34]; 5.4 has to do with a similar calculation in the $r_{h}=0$ limit at LO in $N$ in 5.4.1 and up to NLO in $N$ in 5̱.4.2. Section 6 has to do with M-theory calculations of $0^{++},1^{++},2^{++}$ glueballs arising from appropriate metric fluctuations. Subsection 6.1 is on such a $0^{++}$ glueball spectrum (imposing Neumann/Dirichlet boundary conditions at the horizon/IR cut-off in 6.1.1, and using WKB quantization conditions and the redefinitions of [34] at finite/zero horizon radius in 6.1.2). Subsection 6.2 is on such a $2^{++}$ glueball spectrum (via imposing Neumann/Dirichlet boundary conditions at the horizon/IR cut-off in 6.2.1/6.2.2, and via WKB quantization conditions using redefinitions of [34] at LO in $N$ at $r_{h}\neq 0$ in 6.2.3 and zero $r_{h}$ in 6.2.4 as well as up to NLO in $N$ in the $r_{h}=0$-limit, in 6.2.5). Subsection 6.3 is on such a $1^{++}$ glueball spectrum (via WKB quantization conditions using the redefinitions of [34] up to LO in $N$ in the finite/zero horizon radius limit in 6.3.1/6.3.2, and up to NLO in $N$ in the zero-horizon limit in 6.3.3). From the point of view of comparing string theory and M theory glueball spectrum calculations, we obtain the $2^{++}$ glueball spectrum arising from tensor mode of metric fluctuations in the type IIB background of [1] in Sec. 7. Subsection 7.1 has to do with a supergravity calculation (via Neumann/Dirichlet boundary conditions at the horizon in 7.1.1 and WKB quantization condition using redefinitions of [34] in 7.1.2), and 7.2 has to do with zero-horizon radius limit calculation (WKB quantization condition using redefinitions of [34] at LO in $N$ in 7.2.1 and up to NLO in $N$ in 7.2.2). Section 8 contains a summary of and discussion on the results obtained in this paper. There a appendix A on the square of different fluxes that appear in EOM relevant to spin-two perturbations of the type IIB metric.

2 Background: A Top-Down Type IIB Holographic Large-$N$ Thermal QCD and its M-Theory Uplift in the ‘MQGP’ Limit

In this section, via five sub-sections we will:

• •

provide a short review of the type IIB background of [1] which is supposed to provide a UV complete holographic dual of large-$N$ thermal QCD, as well as their precursors in subsection 2.1,

• •

discuss the ’MQGP’ limit of [2] and the motivation for considering the same in subsection 2.2,

• •

briefly review issues as discussed in [2] pertaining to construction of delocalized S(trominger) Y(au) Z(aslow) mirror and approximate supersymmetry, in subsection 2.3,

• •

briefly review the new results of [3] and [5] pertaining to construction of explicit $SU(3)$ and $G_{2}$ structures respectively of type IIB/IIA and M-theory uplift,

• •

briefly discuss the new Physics-related results of [3] and [5], in subsection 2.4

2.1 Type IIB Dual of Large-$N$ Thermal QCD

In this subsection, we will discuss a UV complete holographic dual of large-$N$ thermal QCD as given in Dasgupta-Mia et al [1]. As partly mentioned in Sec. 1, this was inspired by the zero-temperature Klebanov-Witten model [30], the non-conformal Klebanov-Tseytlin model [35], its IR completion as given in the Klebanov-Strassler model [36] and Ouyang’s inclusion [37] of flavor in the same 444See [70] for earlier attempts at studying back-reacted $D3/D7$ geometry at zero temperature; we thank L. Zayas for bringing [70, 40] to our attention., as well as the non-zero temperature/non-extremal version of [38] (the solution however was not regular as the non-extremality/black hole function and the ten-dimensional warp factor vanished simultaneously at the horizon radius), [39] (valid only at large temperatures) of the Klebanov-Tseytlin model and [40] (addressing the IR), in the absence of flavors.

(a) Brane construction

In order to include fundamental quarks at non-zero temperature in the context of type IIB string theory, to the best of our knowledge, the following model proposed in [1] is the closest to a UV complete holographic dual of large-$N$ thermal QCD. The KS model (after a duality cascade) and QCD have similar IR behavior: $SU(M)$ gauge group and IR confinement. However, they differ drastically in the UV as the former yields a logarithmically divergent gauge coupling (in the UV) - Landau pole. This necessitates modification of the UV sector of the KS model apart from inclusion of non-extremality factors. With this in mind and building up on all of the above, the type IIB holographic dual of [1] was constructed. The setup of [1] is summarized below.

• •

From a gauge-theory perspective, the authors of [1] considered $N$ black $D3$-branes placed at the tip of six-dimensional conifold, $M\ D5$-branes wrapping the vanishing two-cycle and $M\ \overline{D5}$-branes distributed along the resolved two-cycle and placed at the outer boundary of the IR-UV interpolating region/inner boundary of the UV region.

• •

More specifically, the $M\ \overline{D5}$ are distributed around the antipodal point relative to the location of $M\ D5$ branes on the blown-up $S^{2}$. If the $D5/\overline{D5}$ separation is given by ${\cal R}_{D5/\overline{D5}}$, then this provides the boundary common to the outer UV-IR interpolating region and the inner UV region. The region $r>{\cal R}_{D5/\overline{D5}}$ is the UV. In other words, the radial space, in [1] is divided into the IR, the IR-UV interpolating region and the UV. To summarize the above:

• –

$r_{0}/r_{h}<r<|\mu_{\rm Ouyang}|^{\frac{2}{3}}(r_{h}=0)/{\cal R}_{D5/\overline{D5}}(r_{h}\neq 0)$: the IR/IR-UV interpolating regions with $r\sim\Lambda$: deep IR where the $SU(M)$ gauge theory confines

• –

$r>|\mu_{\rm Ouyang}|^{\frac{2}{3}}(r_{h}=0)/{\cal R}_{D5/\overline{D5}}(r_{h}\neq 0)$: the UV region.

• •

$N_{f}\ D7$-branes, via Ouyang embedding, are holomorphically embedded in the UV (asymptotically $AdS_{5}\times T^{1,1}$), the IR-UV interpolating region and dipping into the (confining) IR (up to a certain minimum value of $r$ corresponding to the lightest quark) and $N_{f}\ \overline{D7}$-branes present in the UV and the UV-IR interpolating (not the confining IR). This is to ensure turning off of three-form fluxes, constancy of the axion-dilaton modulus and hence conformality and absence of Landau poles in the UV.

• •

The resultant ten-dimensional geometry hence involves a resolved warped deformed conifold. Back-reactions are included, e.g., in the ten-dimensional warp factor. Of course, the gravity dual, as in the Klebanov-Strassler construct, at the end of the Seiberg-duality cascade will have no $D3$-branes and the $D5$-branes are smeared/dissolved over the blown-up $S^{3}$ and thus replaced by fluxes in the IR.

The delocalized S(trominger) Y(au) Z(aslow) type IIA mirror of the aforementioned type IIB background of [1] and its M-theory uplift had been obtained in [2, 3, 5].

(b) Seiberg duality cascade, IR confining $SU(M)$ gauge theory at finite temperature and

$N_{c}=N_{\rm eff}(r)+M_{\rm eff}(r)$

IR Confinement after Seiberg Duality Cascade: Footnote numbered 3 shows that one effectively adds on to the number of $D3$-branes in the UV and hence, one has $SU(N+M)\times SU(N+M)$ color gauge group (implying an asymptotic $AdS_{5}$) and $SU(N_{f})\times SU(N_{f})$ flavor gauge group, in the UV: $r\geq{\cal R}_{D5/\overline{D5}}$. It is expected that there will be a partial Higgsing of $SU(N+M)\times SU(N+M)$ to $SU(N+M)\times SU(N)$ at $r={\cal R}_{D5/\overline{D5}}$ [41]. The two gauge couplings, $g_{SU(N+M)}$ and $g_{SU(N)}$ flow logarithmically and oppositely in the IR:

[TABLE]

Had it not been for $\int_{S^{2}}B_{2}$, in the UV, one could have set $g_{SU(M+N)}^{2}=g_{SU(N)}^{2}=g_{YM}^{2}\sim g_{s}\equiv$ constant (implying conformality) which is the reason for inclusion of $M$ $\overline{D5}$-branes at the common boundary of the UV-IR interpolating and the UV regions, to annul this contribution. In fact, the running also receives a contribution from the $N_{f}$ flavor $D7$-branes which needs to be annulled via $N_{f}\ \overline{D7}$-branes. The gauge coupling $g_{SU(N+M)}$ flows towards strong coupling and the $SU(N)$ gauge coupling flows towards weak coupling. Upon application of Seiberg duality, $SU(N+M)_{\rm strong}\stackrel{{\scriptstyle\rm Seiberg\ Dual}}{{\longrightarrow}}SU(N-(M-N_{f}))_{\rm weak}$ in the IR; assuming after repeated Seiberg dualities or duality cascade, $N$ decreases to 0 and there is a finite $M$, one will be left with $SU(M)$ gauge theory with $N_{f}$ flavors that confines in the IR - the finite temperature version of the same is what was looked at by [1]. 2. 2.

Obtaining $N_{c}=3$, and Color-Flavor Enhancement of Length Scale in the IR: So, in the IR, at the end of the duality cascade, what gets identified with the number of colors $N_{c}$ is $M$, which in the ‘MQGP limit’ to be discussed below, can be tuned to equal 3. One can identify $N_{c}$ with $N_{\rm eff}(r)+M_{\rm eff}(r)$, where $N_{\rm eff}(r)=\int_{\rm Base\ of\ Resolved\ Warped\ Deformed\ Conifold}F_{5}$ and $M_{\rm eff}=\int_{S^{3}}\tilde{F}_{3}$ (the $S^{3}$ being dual to $\ e_{\psi}\wedge\left(\sin\theta_{1}d\theta_{1}\wedge d\phi_{1}-B_{1}\sin\theta_{2}\wedge d\phi_{2}\right)$, wherein $B_{1}$ is an asymmetry factor defined in [1], and $e_{\psi}\equiv d\psi+{\rm cos}~{}\theta_{1}~{}d\phi_{1}+{\rm cos}~{}\theta_{2}~{}d\phi_{2}$) where $\tilde{F}_{3}(\equiv F_{3}-\tau H_{3})\propto M(r)\equiv 1-\frac{e^{\alpha(r-{\cal R}_{D5/\overline{D5}})}}{1+e^{\alpha(r-{\cal R}_{D5/\overline{D5}})}},\alpha\gg 1$ [42]. The effective number $N_{\rm eff}$ of $D3$-branes varies between $N\gg 1$ in the UV and 0 in the deep IR, and the effective number $M_{\rm eff}$ of $D5$-branes varies between 0 in the UV and $M$ in the deep IR (i.e., at the end of the duality cacade in the IR). Hence, the number of colors $N_{c}$ varies between $M$ in the deep IR and a large value [even in the MQGP limit of (2.2) (for a large value of $N$)] in the UV. Hence, at very low energies, the number of colors $N_{c}$ can be approximated by $M$, which in the MQGP limit is taken to be finite and can hence be taken to be equal to three. However, in this discussion, the low energy or the IR is relative to the string scale. But these energies which are much less than the string scale, can still be much larger than $T_{c}$. Therefore, for all practical purposes, as regard the energy scales relevant to QCD, the number of colors can be tuned to three.

In the IR in the MQGP limit, with the inclusion of terms higher order in $g_{s}N_{f}$ in the RR and NS-NS three-form fluxes and the NLO terms in the angular part of the metric, there occurs an IR color-flavor enhancement of the length scale as compared to a Planckian length scale in KS for ${\cal O}(1)$ $M$, thereby showing that quantum corrections will be suppressed. Using [1]:

[TABLE]

it was argued in [3] that the length scale of the OKS-BH metric in the IR will be given by:

[TABLE]

which implies that in the IR, relative to KS, there is a color-flavor enhancement of the length scale in the OKS-BH metric. Hence, in the IR, even for $N_{c}^{\rm IR}=M=3$ and $N_{f}=2$ (light flavors) upon inclusion of of $n,m>1$ terms in $M_{\rm eff}$ and $N_{f}^{\rm eff}$ in (2), $L_{\rm OKS-BH}\gg L_{\rm KS}(\sim L_{\rm Planck})$ in the MQGP limit involving $g_{s}\stackrel{{\scriptstyle\sim}}{{<}}1$, implying that the stringy corrections are suppressed and one can trust supergravity calculations. As a reminder one will generate higher powers of $M$ and $N_{f}$ in the double summation in $M_{\rm eff}$ in (2), e.g., from the terms higher order in $g_{s}N_{f}$ in the RR and NS-NS three-form fluxes that become relevant for the aforementioned values of $g_{s},N_{f}$. 3. 3.

Further, the global flavor group in the UV-IR interpolating and UV regions, due to presence of $N_{f}$ $D7$ and $N_{f}\ \overline{D7}$-branes, is $SU(N_{f})\times SU(N_{f})$, which is broken in the IR to $SU(N_{f})$ as the IR has only $N_{f}$ $D7$-branes.

Hence, the following features of the type IIB model of [1] make it an ideal holographic dual of thermal QCD:

• •

the theory having quarks transforming in the fundamental representation, is UV conformal and IR confining with the required chiral symmetry breaking in the IR and restoration at high temperatures

• •

the theory is UV complete with the gauge coupling remaining finite in the UV (absence of Landau poles)

• •

the theory is not just defined for large temperatures but for low and high temperatures

• •

(as will become evident in Sec. 3) with the inclusion of a finite baryon chemical potential, the theory provides a lattice-compatible QCD confinement-deconfinement temperature $T_{c}$ for the right number of light quark flavors and masses, and is also thermodynamically stable; given the IR proximity of the value of the lattice-compatible $T_{c}$, after the end of the Seiberg duality cascade, the number of quark flavors approximately equals $M$ which in the ‘MQGP’ limit of (2.2) can be tuned to equal 3

• •

in the MQGP limit (2.2) which requires considering a finite gauge coupling and hence string coupling, the theory was shown in [2] to be holographically renormalizable from an M-theory perspective with the M-theory uplift also being thermodynamically stable.

(d) Supergravity solution on resolved warped deformed conifold

The working metric is given by :

[TABLE]

$g_{i}$’s are black hole functions in modified OKS(Ouyang-Klebanov-Strassler)-BH (Black Hole) background and are assumed to be: $g_{1,2}(r,\theta_{1},\theta_{2})=1-\frac{r_{h}^{4}}{r^{4}}+{\cal O}\left(\frac{g_{s}M^{2}}{N}\right)$ where $r_{h}$ is the horizon, and the ($\theta_{1},\theta_{2}$) dependence come from the ${\cal O}\left(\frac{g_{s}M^{2}}{N}\right)$ corrections. The $h_{i}$’s are expected to receive corrections of ${\cal O}\left(\frac{g_{s}M^{2}}{N}\right)$ [41]. We assume the same to also be true of the ‘black hole functions’ $g_{1,2}$. The compact five dimensional metric in (4), is given as:

[TABLE]

$r\gg a,h_{5}\sim\frac{({\rm deformation\ parameter})^{2}}{r^{3}}\ll 1$ for $r\gg({\rm deformation\ parameter})^{\frac{2}{3}}$, i.e. in the UV/IR-UV interpolating region. The $h_{i}$’s appearing in internal metric as well as $M,N_{f}$ are not constant and up to linear order depend on $g_{s},M,N_{f}$ are given as below:

[TABLE]

One sees from (2.1) and (2.1) that one has a non-extremal resolved warped deformed conifold involving an $S^{2}$-blowup (as $h_{4}-h_{2}=\frac{a^{2}}{r^{2}}$), an $S^{3}$-blowup (as $h_{5}\neq 0$) and squashing of an $S^{2}$ (as $h_{3}$ is not strictly unity). The horizon (being at a finite $r=r_{h}$) is warped squashed $S^{2}\times S^{3}$. In the deep IR, in principle one ends up with a warped squashed $S^{2}(a)\times S^{3}(\epsilon),\ \epsilon$ being the deformation parameter. Assuming $\epsilon^{\frac{2}{3}}>a$ and given that $a={\cal O}\left(\frac{g_{s}M^{2}}{N}\right)r_{h}$ [41], in the IR and in the MQGP limit, $N_{\rm eff}(r\in{\rm IR})=\int_{{\rm warped\ squashed}\ S^{2}(a)\times S^{3}(\epsilon)}F_{5}(r\in{\rm IR})\ll M=\int_{S^{3}(\epsilon)}F_{3}(r\in{\rm IR})$; we have a confining $SU(M)$ gauge theory in the IR.

The warp factor that includes the back-reaction, in the IR is given as:

[TABLE]

where, in principle, $M_{\rm eff}/N_{f}^{\rm eff}$ are not necessarily the same as $M/N_{f}$; we however will assume that up to ${\cal O}\left(\frac{g_{s}M^{2}}{N}\right)$, they are. Proper UV behavior requires [41]:

[TABLE]

In the IR, up to ${\cal O}(g_{s}N_{f})$ and setting $h_{5}=0$, the three-forms are as given in [1]:

[TABLE]

The asymmetry factors in (2.1) are given by: $A_{i}=1+{\cal O}\left(\frac{a^{2}}{r^{2}}\ {\rm or}\ \frac{a^{2}\log r}{r}\ {\rm or}\ \frac{a^{2}\log r}{r^{2}}\right)+{\cal O}\left(\frac{{\rm deformation\ parameter}^{2}}{r^{3}}\right),$ $B_{i}=1+{\cal O}\left(\frac{a^{2}\log r}{r}\ {\rm or}\ \frac{a^{2}\log r}{r^{2}}\ {\rm or}\ \frac{a^{2}\log r}{r^{3}}\right)+{\cal O}\left(\frac{({\rm deformation\ parameter})^{2}}{r^{3}}\right)$. As in the UV, $\frac{({\rm deformation\ parameter})^{2}}{r^{3}}\ll\frac{({\rm resolution\ parameter})^{2}}{r^{2}}$, we will assume the same three-form fluxes for $h_{5}\neq 0$.

Further, to ensure UV conformality, it is important to ensure that the axion-dilaton modulus approaches a constant implying a vanishing beta function in the UV. This was discussed in detail in appendix B of [3], wherein in particular, assuming an F-theory uplift involving, locally, an elliptically fibered $K3$, it was shown that UV conformality and the Ouyang embedding are mutually consistent.

2.2 The ‘MQGP Limit’

In [2], we had considered the following two limits:

[TABLE]

(the limit in the first line though not its realization in the second line, considered in [1]);

[TABLE]

Let us enumerate the motivation for considering the MQGP limit which was discussed in detail in [3]. There are principally two.

Unlike the AdS/CFT limit wherein $g_{\rm YM}\rightarrow 0,N\rightarrow\infty$ such that $g_{\rm YM}^{2}N$ is large, for strongly coupled thermal systems like sQGP, what is relevant is $g_{\rm YM}\sim{\cal O}(1)$ and $N_{c}=3$. From the discussion in the previous paragraphs specially the one in point (c) of sub-section 2.1, one sees that in the IR after the Seiberg duality cascade, effectively $N_{c}=M$ which in the MQGP limit of (2.2) can be tuned to 3. Further, in the same limit, the string coupling $g_{s}\stackrel{{\scriptstyle<}}{{\sim}}1$. The finiteness of the string coupling necessitates addressing the same from an M theory perspective. This is the reason for coining the name: ‘MQGP limit’. In fact this is the reason why one is required to first construct a type IIA mirror, which was done in [2] a la delocalized Strominger-Yau-Zaslow mirror symmetry, and then take its M-theory uplift. 2. 2.

From the perspective of calculational simplification in supergravity, the following are examples of the same and constitute therefore the second set of reasons for looking at the MQGP limit of (2.2):

• •

In the UV-IR interpolating region and the UV, $(M_{\rm eff},N_{\rm eff},N_{f}^{\rm eff})\stackrel{{\scriptstyle\rm MQGP}}{{\approx}}(M,N,N_{f})$

• •

Asymmetry Factors $A_{i},B_{j}$(in three-form fluxes)$\stackrel{{\scriptstyle MQGP}}{{\rightarrow}}1$ in the UV-IR interpolating region and the UV.

• •

Simplification of ten-dimensional warp factor and non-extremality function in MQGP limit

With ${\cal R}_{D5/\overline{D5}}$ denoting the boundary common to the UV-IR interpolating region and the UV region, $\tilde{F}_{lmn},H_{lmn}=0$ for $r\geq{\cal R}_{D5/\overline{D5}}$ is required to ensure conformality in the UV. Near the $\theta_{1}=\theta_{2}=0$-branch, assuming: $\theta_{1,2}\rightarrow 0$ as $\epsilon^{\gamma_{\theta}>0}$ and $r\rightarrow{\cal R}_{\rm UV}\rightarrow\infty$ as $\epsilon^{-\gamma_{r}<0},\lim_{r\rightarrow\infty}\tilde{F}_{lmn}=0$ and $\lim_{r\rightarrow\infty}H_{lmn}=0$ for all components except $H_{\theta_{1}\theta_{2}\phi_{1,2}}$; in the MQGP limit and near $\theta_{1,2}=\pi/0$-branch, $H_{\theta_{1}\theta_{2}\phi_{1,2}}=0/\left.\frac{3g_{s}^{2}MN_{f}}{8\pi}\right|_{N_{f}=2,g_{s}=0.6,M=\left({\cal O}(1)g_{s}\right)^{-\frac{3}{2}}}\ll 1.$ So, the UV nature too is captured near $\theta_{1,2}=0$-branch in the MQGP limit. This mimics addition of $\overline{D5}$-branes in [1] to ensure cancellation of $\tilde{F}_{3}$.

2.3 Approximate Supersymmetry, Construction of the Delocalized SYZ IIA Mirror and Its M-Theory Uplift in the MQGP Limit

A central issue to [2, 29] has been implementation of delocalized mirror symmetry via the Strominger Yau Zaslow prescription according to which the mirror of a Calabi-Yau can be constructed via three T dualities along a special Lagrangian $T^{3}$ fibered over a large base in the Calabi-Yau. This sub-section is a quick review of precisely this.

To implement the quantum mirror symmetry a la S(trominger)Y(au)Z(aslow) [43], one needs a special Lagrangian (sLag) $T^{3}$ fibered over a large base (to nullify contributions from open-string disc instantons with boundaries as non-contractible one-cycles in the sLag). Defining delocalized T-duality coordinates, $(\phi_{1},\phi_{2},\psi)\rightarrow(x,y,z)$ valued in $T^{3}(x,y,z)$ [2]:

[TABLE]

using the results of [44] it was shown in [29, 5] that the following conditions are satisfied:

[TABLE]

for the $T^{2}$-invariant sLag of [44] for a deformed conifold. Hence, if the resolved warped deformed conifold is predominantly either resolved or deformed, the local $T^{3}$ of (12) is the required sLag to effect SYZ mirror construction.

Interestingly, in the ‘delocalized limit’ [45] $\psi=\langle\psi\rangle$, under the coordinate transformation:

[TABLE]

and $\psi\rightarrow\psi-\cos\langle{\bar{\theta}}_{2}\rangle\phi_{2}+\cos\langle\theta_{2}\rangle\phi_{2}-\tan\langle\psi\rangle ln\sin{\bar{\theta}}_{2}$, the $h_{5}$ term becomes $h_{5}\left[d\theta_{1}d\theta_{2}-sin\theta_{1}sin\theta_{2}d\phi_{1}d\phi_{2}\right]$, $e_{\psi}\rightarrow e_{\psi}$, i.e., one introduces an isometry along $\psi$ in addition to the isometries along $\phi_{1,2}$. This clearly is not valid globally - the deformed conifold does not possess a third global isometry.

To enable use of SYZ-mirror duality via three T dualities, one also needs to ensure a large base (implying large complex structures of the aforementioned two two-tori) of the $T^{3}(x,y,z)$ fibration. This is effected via [46]:

[TABLE]

for appropriately chosen large values of $f_{1,2}(\theta_{1,2})$. The three-form fluxes remain invariant. The fact that one can choose such large values of $f_{1,2}(\theta_{1,2})$, was justified in [2]. The guiding principle is that one requires the metric obtained after SYZ-mirror transformation applied to the non-Kähler resolved warped deformed conifold is like a non-Kähler warped resolved conifold at least locally. Then $G^{IIA}_{\theta_{1}\theta_{2}}$ needs to vanish [2]. This was explicitly shown in [3].

As in Klebanov-Strassler construction, a single T-duality along a direction orthogonal to the $D3$-brane world volume, e.g., $z$ of (12), yields $D4$ branes straddling a pair of $NS5$-branes consisting of world-volume coordinates $(\theta_{1},x)$ and $(\theta_{2},y)$. Further, T-dualizing along $x$ and then $y$ would yield a Taub-NUT space from each of the two $NS5$-branes [47]. The $D7$-branes yield $D6$-branes which get uplifted to Kaluza-Klein monopoles in M-theory [48] which too involve Taub-NUT spaces. Globally, probably the eleven-dimensional uplift would involve a seven-fold of $G_{2}$-structure, analogous to the uplift of $D5$-branes wrapping a two-cycle in a resolved warped conifold [49].

2.4 $G$-Structures

The mirror type IIA metric after performing three T-dualities, first along $x$, then along $y$ and finally along $z$, utilizing the results of [45] was worked out in [2]. The type IIA metric components were worked out in [2].

Now, any metric-compatible connection can be written in terms of the Levi-Civita connection and the contorsion tensor $\kappa$. Metric compatibility requires $\kappa\in\Lambda^{1}\otimes\Lambda^{2}$, $\Lambda^{n}$ being the space of $n$-forms. Alternatively, in $d$ complex dimensions, since $\Lambda^{2}\cong so(d)$, $\kappa$ also be thought of as $\Lambda^{1}\otimes so(d)$. Given the existence of a $G$-structure, we can decompose $so(d)$ into a part in the Lie algebra $g$ of $G\subset SO(d)$ and its orthogonal complement $g^{\perp}=so(d)/g$. The contorsion $\kappa$ splits accordingly into $\kappa=\kappa^{0}+\kappa^{g}$, where $\kappa^{0}$ - the intrinsinc torsion - is the part in $\Lambda^{1}\otimes g^{\perp}$. One can decompose $\kappa^{0}$ into irreducible $G$ representations providing a classification of $G$-structures in terms of which representations appear in the decomposition. Let us consider the decomposition of $T^{0}$ in the case of $SU(3)$-structure. The relevant representations are $\Lambda^{1}\sim 3\oplus\bar{3},g\sim 8,g^{\perp}\sim 1\oplus 3\oplus\bar{3}.$ Thus the intrinsic torsion, an element of $\Lambda^{1}\oplus su(3)^{\perp}$, can be decomposed into the following $SU(3)$ representations:

[TABLE]

The $SU(3)$ structure torsion classes [50] can be defined in terms of J, $\Omega$, dJ, $d{\Omega}$ and the contraction operator $\lrcorner:{\Lambda}^{k}T^{\star}\otimes{\Lambda}^{n}T^{\star}\rightarrow{\Lambda}^{n-k}T^{\star}$, $J$ being given by:

[TABLE]

and the (3,0)-form $\Omega$ being given by

[TABLE]

The torsion classes are defined in the following way:

• •

$W_{1}\leftrightarrow[dJ]^{(3,0)}$, given by real numbers $W_{1}=W_{1}^{+}+W_{1}^{-}$ with $d{\Omega}_{+}\wedge J={\Omega}_{+}\wedge dJ=W_{1}^{+}J\wedge J\wedge J$ and $d{\Omega}_{-}\wedge J={\Omega}_{-}\wedge dJ=W_{1}^{-}J\wedge J\wedge J$;

• •

$W_{2}\leftrightarrow[d\Omega]_{0}^{(2,2)}$ : $(d{\Omega}_{+})^{(2,2)}=W_{1}^{+}J\wedge J+W_{2}^{+}\wedge J$ and $(d{\Omega}_{-})^{(2,2)}=W_{1}^{-}J\wedge J+W_{2}^{-}\wedge J$;

• •

$W_{3}\leftrightarrow[dJ]_{0}^{(2,1)}$ is defined as $W_{3}=dJ^{(2,1)}-[J\wedge W_{4}]^{(2,1)}$;

• •

$W_{4}\leftrightarrow J\wedge dJ$ : $W_{4}=\frac{1}{2}J\lrcorner dJ$;

• •

$W_{5}\leftrightarrow[d\Omega]_{0}^{(3,1)}$: $W_{5}=\frac{1}{2}{\Omega}_{+}\lrcorner d{\Omega}_{+}$ (the subscript 0 indicative of the primitivity of the respective forms).

In [29], we saw that the five $SU(3)$ structure torsion classes, in the MQGP limit, satisfied (schematically):

[TABLE]

$(r\sim e^{\frac{\tau}{3}})$, such that

[TABLE]

in the UV-IR interpolating region/UV, implying a Klebanov-Strassler-like supersymmetry [51]. Locally around $\theta_{1}\sim\frac{1}{N^{\frac{1}{5}}},\theta_{2}\sim\frac{1}{N^{\frac{3}{10}}}$, the type IIA torsion classes of the delocalized SYZ type IIA mirror metric (LABEL:Type_IIA_metric), were shown in [3] to be:

[TABLE]

Further,

[TABLE]

indicative of supersymmetry after constructing the delocalized SYZ mirror.

Apart from quantifying the departure from $SU(3)$ holonomy due to intrinsic contorsion supplied by the NS-NS three-form $H$, via the evaluation of the $SU(3)$ structure torsion classes, to our knowledge for the first time in the context of holographic thermal QCD at finite gauge coupling in [3]:

(i) the existence of approximate supersymmetry of the type IIB holographic dual of [1] in the MQGP limit near the coordinate branch $\theta_{1}=\theta_{2}=0$ was demonstrated, which apart from the existence of a special Lagrangian three-cycle (as shown in [29, 3]) is essential for construction of the local SYZ type IIA mirror;

(ii) it was demonstrated that the large-$N$ suppression of the deviation of the type IIB resolved warped deformed conifold from being a complex manifold, is lost on being duality-chased to type IIA - it was also shown that one further fine tuning $\gamma_{2}=0$ in $W_{2}^{\rm IIA}$ can ensure that the local type IIA mirror is complex;

(iii) for the local type IIA $SU(3)$ mirror, the possibility of surviving approximate supersymmetry was demonstrated which is essential from the point of view of the end result of application of the SYZ mirror prescription.

We can get a one-form type IIA potential from the triple T-dual (along $x,y,z$) of the type IIB $F_{1,3,5}$ in [2] and using which the following $D=11$ metric was obtained in [2] ($u\equiv\frac{r_{h}}{r}$):

[TABLE]

If $V$ is a seven-dimensional real vector space, then a three-form $\varphi$ is said to be positive if it lies in the $GL\left(7,\mathbb{R}\right)$ orbit of $\varphi_{0}$, where $\varphi_{0}$ is a three-form on $\mathbb{R}^{7}$ which is preserved by $G_{2}$-subgroup of $GL(7,\mathbb{R})$. The pair $\left(\varphi,g\right)$ for a positive $3$-form $\varphi$ and corresponding metric $g$ constitute a $G_{2}$-structure. The space of $p$-forms decompose as following irreps of $G_{2}$ [52]:

[TABLE]

The subscripts denote the dimension of representation and components of same representation/dimensionality, are isomorphic to each other. Let $M$ be a $7$-manifold with a $G_{2}$-structure $\left(\varphi,g\right)$. Then the components of spaces of $2$-, $3$-, $4$-, and $5$-forms are given in [52, 53]. The metric $g$ defines a reduction of the frame bundle F to a principal $SO\left(7\right)$-sub-bundle $Q$, that is, a sub-bundle of oriented orthonormal frames. Now, $g$ also defines a Levi-Civita connection $\nabla$ on the tangent bundle $TM$, and hence on $F$. However, the $G_{2}$-invariant $3$-form $\varphi$ reduces the orthonormal bundle further to a principal $G_{2}$-subbundle $Q$. The Levi-Civita connection can be pulled back to $Q$. On $Q$, $\nabla$ can be uniquely decomposed as

[TABLE]

where $\bar{\nabla}$ is a $G_{2}$-compatible canonical connection on $P$, taking values in the sub-algebra $\mathfrak{g}_{2}\subset\mathfrak{so}\left(7\right)$, while $\mathcal{T}$ is a $1$-form taking values in $\mathfrak{g}_{2}^{\perp}\subset\mathfrak{so}\left(7\right)$; $\mathcal{T}$ is known as the intrinsic torsion of the $G_{2}$-structure - the obstruction to the Levi-Civita connection being $G_{2}$-compatible. Now $\mathfrak{so}\left(7\right)$ splits under $G_{2}$ as

[TABLE]

But $\Lambda_{14}^{2}\cong\mathfrak{g}_{2}$, so the orthogonal complement $\mathfrak{g}_{2}^{\perp}\cong\Lambda_{7}^{2}\cong V$. Hence $\mathcal{T}$ can be represented by a tensor $T_{ab}$ which lies in $W\cong V\otimes V$. Now, since $\varphi$ is $G_{2}$-invariant, it is $\bar{\nabla}$-parallel. So, the torsion is determined by $\nabla\varphi$. from Lemma 2.24 of [54]:

[TABLE]

Due to the isomorphism between the $\Lambda^{a=1,...,5}_{7}$s, $\nabla\varphi$ lies in the same space as $T_{AB}$ and thus completely determines it. Equation (25) is equivalent to:

[TABLE]

where $T_{AB}$ is the full torsion tensor. Equation (26) can be inverted to yield:

[TABLE]

The tensor $T_{A}^{\ M}$, like the space W, possesses 49 components and hence fully defines $\nabla\varphi$. In general $T_{AB}$ cab be split into torsion components as

[TABLE]

where $T_{1}$ is a function and gives the $\mathbf{1}$ component of $T$ . We also have $T_{7}$, which is a $1$-form and hence gives the $\mathbf{7}$ component, and, $T_{14}\in\Lambda_{14}^{2}$ gives the $\mathbf{14}$ component. Further, $T_{27}$ is traceless symmetric, and gives the $\mathbf{27}$ component. Writing $T_{i}$ as $W_{i}$, we can split $W$ as

[TABLE]

From [55], we see that a $G_{2}$ structure can be defined as:

[TABLE]

where $A,B,C=1,...,6,10;a,b,c,=1,...,6$ and ${f}_{ABC}$ are the structure constants of the imaginary octonions. Using the same, the $G_{2}$-structure torsion classes were worked out around $\theta_{1}\sim\frac{1}{N^{\frac{1}{5}}},\theta_{2}\sim\frac{1}{N^{\frac{3}{10}}}$ in [3] to:

[TABLE]

Hence, the approach of the seven-fold, locally, to having a $G_{2}$ holonomy ($W_{1}^{G_{2}}=W_{2}^{G_{2}}=W_{3}^{G_{2}}=W_{4}^{G_{2}}=0$) is accelerated in the MQGP limit.

As stated earlier, the global uplift to M-theory of the type $IIB$ background of [1] is expected to involve a seven-fold of $G_{2}$ structure (not $G_{2}$-holonomy due to non-zero $G_{4}$). It is hence extremely important to be able to see this, at least locally. It is in this sense that the results of [2] are of great significance as one explicitly sees, for the first time, in the context of holographic thermal QCD at finite gauge coupling, though locally, the aforementioned $G_{2}$ structure having worked out the non-trivial $G_{2}$-structure torsion classes.

3 $0^{++(***...}$ Glueball spectrum from type $IIB$ supergravity background

In this section we discuss the $0^{++}$ glueball spectrum by solving the dilaton wave equation in the type $IIB$ background discussed in section $2$. The type $IIB$ metric as given in equation 4 with the warp factor $h$ given in 7 can be simplified by working around a particular value of $\theta_{1}$ and $\theta_{2}$: $\{\theta_{1}=\frac{1}{N^{1/5}},\theta_{2}=\frac{1}{N^{3/10}}\}$, then keeping terms upto (N)ext to (L)eading (O)rder in $N$ in the large $N$ limit. The simplified type $IIB$ metric is given as:

[TABLE]

with the components $g_{tt},g_{x_{1}x_{1}}(=g_{x_{2}x_{2}}=g_{x_{3}x_{3}}),g_{rr}$ as given below:

[TABLE]

where $B(r)=\frac{3M^{2}g_{s}\log(r)\left(12N_{f}g_{s}\log(r)+6N_{f}g_{s}+N_{f}g_{s}(-\log(N))-2\log(4)N_{f}g_{s}+8\pi\right)}{32\pi^{2}N}$, and $a$ being the resolution parameter is proportional to the horizon radius $r_{h}$. Hence while computing the spectrum with a cut-off in the radial direction and no horizon, we must put both $r_{h}$ and $a$ to zero in the above equation.

Moreover the dilaton profiles with/without the black-hole are given below as:

[TABLE]

Again working around the particular choices of $\theta_{1}$ and $\theta_{2}$ the above profile can be simplified upto NLO in $N$. The dilaton equation that has to be solved is given as:

[TABLE]

To solve the above dilaton equation we assume $\phi$ in (35) to be of the form $\phi=e^{ik.x}\tilde{\phi}(r)$. Now with this, we adopt the WKB method to get to the final result. The first step towards the WKB method is to convert the glueball equation of motion into a Schrödinger-like equation. Then the WKB quantization condition can be applied on the potential term obtained from the Schrödinger-like equation. For $(0^{++})$ glueball spectrum with no horizon $(r_{h}=0)$, one of the solution was obtained by imposing Neumann boundary condition at the cut-off.

3.1 $r_{h}\neq 0$ using WKB quantization method

In the following we discuss the spectrum of $0^{++}$ glueball in the type $IIB$ background with a black hole, implying a horizon of radius $r_{h}$ in the geometry. The results corresponding to the coordinate and field redefinitions of [34], are discussed below.

Using the redefinitions of [34] with $r=\sqrt{y},r_{h}=\sqrt{y_{h}}$ and finally $y=y_{h}\left(1+e^{z}\right)$, the $0^{++}$ EOM (35), with $k^{2}=-m^{2}$ can be written as:

[TABLE]

where at leading order in $N$, $E_{z}$ and $F_{z}$ are given with $L=(4\pi g_{s}N)^{1/4}$ as:

[TABLE]

[TABLE]

Now, redefining the wave function $\tilde{\phi}$ as $\psi(z)=\sqrt{E_{z}}\tilde{\phi}(z)$ equation (36) reduces to a Schrödinger-like equation

[TABLE]

where the potential $V(z)$ is a rather cumbersome expression which we will not explicitly write out. The WKB quantization condition becomes: $\int_{z_{1}}^{z_{2}}\sqrt{V(z)}=\left(n+\frac{1}{2}\right)\pi$ where $z_{1,2}$ are the turning points of $V(z)$. We will work below with a dimensionless glueball mass $\tilde{m}$ assumed to be large and defined via: $m=\tilde{m}\frac{r_{h}}{L^{2}}$. To determine the turning points of the potential $V(z)$, we consider two limits of the same - $r\in[r_{h},\sqrt{3}a(r_{h})\approx\sqrt{3}br_{h}]\cup[\sqrt{3}br_{h},\infty)\equiv(IR,IR/UV\ {\rm interpolating}\ +UV)$. In the IR, we have to take the limit $z\rightarrow-\infty$. Now in the large $\tilde{m}$ and large $\log{N}$ limit this potential at small $z$ can be shown to be given as:

[TABLE]

Hence, there are no turning points in the IR.

Now, in the UV, apart from taking the large $z$ limit we also have to take $N_{f}=M=0$, to get:

[TABLE]

The turning points of (41) are:

$z_{1}=\frac{1}{8}\left(3b^{2}-\sqrt{9b^{4}-6b^{2}\left(7{y_{h}}\tilde{m}^{2}+18\right)+{y_{h}}\tilde{m}^{4}-36{y_{h}}\tilde{m}^{2}-108}+{y_{h}}\tilde{m}^{2}+6\right)$,

$z_{2}=\frac{1}{8}\left(3b^{2}+\sqrt{9b^{4}-6b^{2}\left(7{y_{h}}\tilde{m}^{2}+18\right)+{y_{h}}\tilde{m}^{4}-36{y_{h}}\tilde{m}^{2}-108}+{y_{h}}\tilde{m}^{2}+6\right)$ which in the large $\tilde{m}$ limit is given as: \Biggl{\{}z_{1}=(3+3b^{2})+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right),z_{2}=\frac{\tilde{m}^{2}}{4}-\frac{3(2+3b^{2})}{4}+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\Biggr{\}}.

To obtain a real spectrum, one first notes:

[TABLE]

and

[TABLE]

Therefore

[TABLE]

yielding:

[TABLE]

3.2 Glueball Mass for $r_{h}=0$ and an IR cut-off $r_{0}$

In this background the type $IIB$ metric and the dilaton profile has to be modified by the limit $r_{h}\rightarrow 0$ and hence with $a\rightarrow 0$. This time also we have provided two solutions to the dilaton equation. The first solutions was obtained by first following the redefinition of variables in [33] and then imposing a neumann boundary condition at the radial cut-off. For the other solution we again consider the WKB method after a redefinition of variables as given in [34].

3.2.1 Neumann boundary condition at $r_{0}$

Following [33], we redefine the radial coordinate as $z=\frac{1}{r}$. With this change of variable, the radial cut-off now maps to $z=z_{0}$, with $z_{0}=\frac{1}{r_{0}}$. The dilaton equation using the metric and the dilaton background in the limit $(r_{h},a)\rightarrow 0$ is given as:

[TABLE]

where upto NLO in $N$ we have,

[TABLE]

Now to convert the above equation in a one-dimensional schrodinger like form we introduce a new field variable $\psi(z)$ as: $\psi(z)=e^{-U}\tilde{\phi}(z)$.

With this one can write the equation in the following schrodinger like form,

[TABLE]

The potential $V(z)$, in the large-$N$ large-$\log{N}$ limit is given as:

[TABLE]

Hence, the Schrödinger equation becomes:

[TABLE]

whose solution is given as under:

[TABLE]

Requiring finiteness of $\psi(z)$ at $z=0$ requires setting $c_{2}=0$. Then imposing Neumann boundary condition on $\tilde{\phi}(z)$ at $z=z_{0}$ implies:

[TABLE]

implying in the large-$N$ large-$z$ (as the Neumann boundary condition will be implemented in the IR) limit:

[TABLE]

where $x_{0}\equiv 2\sqrt{g_{s}N\pi}mz_{0}$. The graphical solution points out that the ground state has a zero mass and the lightest (first excited state) glueball mass is approximately given by $3.71\frac{r_{0}}{L^{2}}$.

3.2.2 WKB method: Including the Non-Conformal/NLO (in $N$) Corrections

Again following the redefinition of variables as given in [34]: $r=\sqrt{y},y=y_{0}\left(1+e^{z}\right)$, where $r_{0}=\sqrt{y_{0}}$ is the radial cut-off, and using the type $IIB$ metric as well as the dilaton profile in the limit $(rh,a)\rightarrow 0$, the dilaton equation (35) can be written as:

[TABLE]

where $C_{z}$ and $D_{z}$ are given up to NLO in $N$ by,

[TABLE]

[TABLE]

Defining a new variable $\psi(z)$ such that: $\psi(z)=\sqrt{C_{z}}\tilde{\phi}(z)$, the above equation can be converted into a Schrödinger-like equation,

[TABLE]

where the potential is given up to NLO in $N$ as:

[TABLE]

The domain of integration over which $V(z)>0$ can be shown to be:

$\{\log\left(\frac{3{g_{s}}^{2}M^{2}{N_{f}}\log(N)\log({y_{0}})}{64\pi^{2}\tilde{m}N}+\frac{1}{\tilde{m}}\right),\log\left(\delta^{2}-1\right)\}$, where $\mu^{\frac{2}{3}}=\delta\sqrt{y_{0}}$. Note that the IR cut-off $r_{0}$ or $y_{0}$ is not put in by hand but is proportional to the Ouyang embedding parameter (corresponding to the embedding of the flavor $D7$-branes ) raised to two-third power. The proportionality constant $\delta$ could be determined, as discussed in point number 5 in Section 8, by matching with lattice calculations and turns out to be ${\cal O}(1)$.

Expanding $\sqrt{V}$ first in $N$ and then in $\tilde{m}$ and then integrating over the above mentioned domain, one gets the following quantization condition,

[TABLE]

yielding:

[TABLE]

4 Scalar Glueball ($0^{-+(***...)}$) Masses

As $Tr(F\tilde{F})$ has $P=-,C=+$ and it couples to $A_{1}$ in the Wess-Zumino term for the type IIA $D4-$brane: $\int_{\Sigma_{4,1}}A\wedge F\wedge F$, one considers $A_{1}$’s EOM:

[TABLE]

where $\mu,\nu,...=a(\equiv 0,1,2,3),r,\alpha(\equiv 5,...,9)$. Like [56], assume $A_{\mu}=\delta^{\theta_{2}}_{\mu}a_{\theta_{2}}(r)e^{ik\cdot x},k^{2}=-m^{2}$ as the fluctuation about the type IIA $A_{1}$ that was worked out in [2]. The $0^{-+}$ EOM then reduces to:

[TABLE]

4.1 Neumann/Dirichlet Boundary Conditions

4.1.1 $r_{h}\neq 0$

Then taking the large-N limit followed by a small-$\theta_{1}$ small-$\theta_{2}$ limit one can show that the equation of motion (62) yields::

[TABLE]

Working near $r=r_{h}$, approximating (64) by:

[TABLE]

whose solution is given as under:

[TABLE]

Imposing Neumann boundary condition: $a_{\theta_{2}}^{\prime}(r=r_{h})=0$, utilizing:

[TABLE]

$r_{h}=T\sqrt{4\pi g_{s}N}$ - up to LO in $N$ - requires $c_{2}=0$ and

[TABLE]

implying:

[TABLE]

One can show that imposing Dirichlet boundary condition $a_{\theta_{2}}(r=r_{h})=0$, yields the same spectrum as (68). If the temperature $T$ gets identified with $a$ of [56], then the ground state, unlike [56], is massless; the excited states for lower $n$’s are closer to $a=0$ and the higher excited states are closer to $a\rightarrow\infty$ in [56].

4.1.2 $r_{h}=0$ Limit of (63)

The $r_{h}=0$ limit of (63) gives:

[TABLE]

which near $r=r_{0}$ yields:

[TABLE]

Imposing Neumann boundary condition on (70) yields:

[TABLE]

One can similarly show that imposing Dirichlet boundary condition on (70) for $c_{2}=0$ yields:

[TABLE]

4.1.3 WKB Quantization for $r_{h}\neq 0$

The potential corresponding to the Schrödinger-like equation a la [34], substituting $m=\tilde{m}\frac{\sqrt{y_{h}}}{L^{2}}$, is given by:

[TABLE]

Therefore, in the IR:

[TABLE]

the turning points being given by $-\infty$ and $\log\left(\frac{2\ 3^{2/3}\tilde{m}^{2}+24}{3\ 3^{2/3}\tilde{m}^{2}+88}\right)\approx-0.405$. But only $z\in(-\infty,-2.526]$ corresponds to the IR in our calculations. So,

[TABLE]

which obtains:

[TABLE]

Similarly, in the UV:

[TABLE]

whose turning points are: $\log\left(\frac{1}{25}\left(2\ 3^{2/3}\tilde{m}^{2}\pm\sqrt{6}\sqrt{2\sqrt[3]{3}\tilde{m}^{4}-36\ 3^{2/3}\tilde{m}^{2}+111}+21\right)\right)\\ =\log(3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)),\log\left(0.33\tilde{m}^{2}-1.32+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right)$. Now: $\sqrt{V(UV,T)}=\frac{1}{2}\sqrt[3]{3}\tilde{m}e^{-z}\sqrt{e^{z}-3}+{\cal O}\left(\frac{1}{\tilde{m}}\right)$. Therefore,

[TABLE]

which obtains:

[TABLE]

4.1.4 WKB Quantization at $r_{h}=0$

The ‘potential’ is given by:

[TABLE]

Therefore, in the IR:

[TABLE]

and in the IR, the domain of integration becomes: $[\log\left(\frac{\sqrt{2}}{3^{\frac{1}{3}}\tilde{m}}\right),\log(\delta^{2}-1)]$: $\int_{\log\left(\frac{\sqrt{2}}{3^{\frac{1}{3}}\tilde{m}}\right)}^{\log(\delta^{2}-1)}\sqrt{V(IR,r_{h}=0)}=\frac{3^{\frac{1}{3}}(\delta^{2}-1)}{2}\tilde{m}-1.1126$, yielding:

[TABLE]

Also, in the UV:

[TABLE]

whose turning points are: $\frac{1}{25}\left(23^{2/3}\tilde{m}^{2}\pm\sqrt{12\sqrt[3]{3}\tilde{m}^{4}-2163^{2/3}\tilde{m}^{2}-2134}+21\right)=\left(3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right),\frac{4}{25}3^{\frac{2}{3}}\tilde{m}^{2}-\frac{33}{25}\right)$, yielding: $\int_{3}^{\frac{4}{25}3^{2/3}\tilde{m}^{2}-\frac{33}{25}}\sqrt{V(UV,r_{h}=0)}=\frac{\pi}{43^{\frac{1}{6}}}=\left(n+\frac{1}{2}\right)\pi$ which obtains:

[TABLE]

5 Glueball ($0^{--(***...)}$) Masses

5.1 $r_{h}\neq 0$ and Neumann/Dirichlet Boundary Conditions at $r=r_{h}$

Given the Weiss-Zumino term $A^{\mu\nu}d^{abc}{\rm Tr}\left(F_{\mu\rho}^{a}F^{b\ \rho}_{\lambda}F^{c\ \lambda}_{\nu}\right)$ and the two-form potential $A_{\mu\nu}$ is dual to a pseudo-scalar, for $r_{h}\neq 0$, corresponding to ${\rm QCD}_{3}$, one writes down the EOM for the fluctuation $\delta A^{23}$. The $B_{MN},C_{MN}$ EOMs are:

[TABLE]

or defining $A_{MN}=B_{MN}+iC_{MN}$, (5.1) can be rewritten as:

[TABLE]

Now, $A_{MN}\rightarrow A^{(0)}_{MN}+\delta A_{MN}$ with $\delta A^{MN}=\delta^{M}_{2}\delta^{N}_{3}\delta A_{23}$, the EOM satisfied by $\delta A_{23}(x^{0,1,2,3},r)=\int\frac{d^{4}k}{\left(2\pi\right)^{4}}e^{ik\cdot x}g_{22}G(r)$ reduces to:

[TABLE]

Assuming $a=\left(\alpha+\beta\frac{g_{s}M^{2}}{N}+\gamma\frac{g_{s}M^{2}}{N}\log r_{h}\right)r_{h}$ [41, 5] for $(\alpha,\beta,\gamma)=(0.6,4,4)$ [5] and $k^{\mu}=(\omega,k_{1},0,0):k^{2}=-m^{2}$, and defining ${\cal G}(r)\equiv g_{22}G(r)$ the EOM for $G(r)$ is:

[TABLE]

The EOM (5.1), near $r=r_{h}$ can be approximated as:

[TABLE]

where

[TABLE]

The solution to (89) is given by:

[TABLE]

implying:

[TABLE]

Assuming $c_{2}=0$, the Neumann boundary condition at $r=r_{h}$ can be satisfied by setting the argument of the gamma function to a negative integer $n$. It runs out setting $\frac{-2{a_{2}}+{b_{1}}+\sqrt{{b_{1}}^{2}-4{b_{2}}}}{2\sqrt{{b_{1}}^{2}-4{b_{2}}}}=-n\in\mathbb{Z}^{-}\cup\left\{0\right\}$ produces a negligible ground state $0^{--}$ mass. Hence, we consider $\frac{-2{a_{2}}+{b_{1}}+3\sqrt{{b_{1}}^{2}-4{b_{2}}}}{2\sqrt{{b_{1}}^{2}-4{b_{2}}}}=-n\in\mathbb{Z}^{-}\cup\left\{0\right\}$, which gives a finite ground state mass. This condition up to LO in $N$ yields:

[TABLE]

the solution to which are given below:

[TABLE]

One can show that one obtains the same spectrum as in (5.1) after imposing Dirichlet boundary condition $G(r=r_{h})=0$.

5.2 $r_{h}=0$ limit of (5.1)

[TABLE]

The EOM (5.2) near $r=r_{0}$ - IR cut-off at $r_{h}=0$, reduces to:

[TABLE]

where

[TABLE]

The solution to (96) is given by:

[TABLE]

One can then work out $G^{\prime}(r=r_{0})=\left.\left(\frac{{\cal G}(r)}{g_{22}}\right)^{\prime}\right|_{r=r_{0}}$. Now, setting $c_{2}=0$, defining $\tilde{m}$ via: $m=\tilde{m}\frac{r_{0}}{L^{2}}$, and using the large $\tilde{m}$-limit of Hermite functions:

[TABLE]

and

[TABLE]

one can show that the Neumann(/Dirichlet: $G(r=r_{0})=0$) boundary condition at $r=r_{0}$ is equivalent to the condition:

[TABLE]

yielding:

[TABLE]

5.3 $0^{--}$ Glueball Spectrum from WKB Method for $r_{h}\neq 0$

Using the variables of [34], the potential in the IR is given as:

[TABLE]

where in the ‘large’ $\tilde{m}$-limit, $V(IR,T)>0$ for (i)$z:e^{z}>\frac{2.31799\times 10^{31}-3.86332\times 10^{28}\tilde{m}^{2}}{6.25374\times 10^{31}-5.79497\times 10^{29}\tilde{m}^{2}}=0.067+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)$ if $\tilde{m}>24.495$ and (ii) $0<e^{z}<\frac{2.31799\times 10^{31}-3.86332\times 10^{28}\tilde{m}^{2}}{6.25374\times 10^{31}-5.79497\times 10^{29}\tilde{m}^{2}}=0.067+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)$ if $10.388<\tilde{m}<24.495$. One can show that: $\int_{\log(0.067)}^{-2.526}\sqrt{V(IR,T)}=\int_{-\infty}^{\log(0.067)}\sqrt{V(IR,T)}\approx 0$, implying there is no contribution to the WKB quantization condition in the IR.

Now, consider:

[TABLE]

For $\tilde{m}>4.29$ the turning points of $V(UV,T)$ are $0.125\tilde{m}^{2}-0.025\sqrt{25.\tilde{m}^{4}+18.\tilde{m}^{2}-8775.}+4.125=4.08+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)<e^{z}<0.125\tilde{m}^{2}+0.025\sqrt{25.\tilde{m}^{4}+18.\tilde{m}^{2}-8775.}+4.125=0.25\tilde{m}^{2}+4.17+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)$. Hence,

[TABLE]

Hence one obtains isospectrality with $0^{++}$; for large $n$, there is also isospectrality with $2^{++}$.

5.4 WKB Method at $r_{h}=0$

In this section we will discuss obtaining the spectrum at $r_{h}=0$ using WKB quantization condition at LO in $N$ in 5.3.1 and up to NLO in $N$ in 5.3.2.

5.4.1 LO in $N$

In the IR, the WKB ‘potential’ can be shown to be given by:

[TABLE]

with turning points: $(\log\left(\frac{1}{m_{0}}+{\cal O}\left(\frac{1}{m_{0}^{3}}\right)\right)\approx-\log(m_{0}),\log\left(\delta^{2}-1\right))$. Further dropping ${\cal O}\left(\frac{1}{m_{0}^{3}}\right)$ terms, $\int_{-\log m_{0}}^{-2.526}\sqrt{V(IR,r_{h}=0)}=\frac{\left(\delta^{2}-1\right)}{2}\tilde{m}-0.785=\left(n+\frac{1}{2}\right)\pi$ yielding:

[TABLE]

In the UV,

[TABLE]

with turning points: $\left(\log\left(3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right),\log\left(\frac{\tilde{m}^{2}}{4}-\frac{3}{2}+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right)\right)$, and $\sqrt{V(UV,r_{h}=0)}=\frac{e^{-z}}{2}\sqrt{e^{z}-3}\tilde{m}+{\cal O}\left(\frac{1}{\tilde{m}}\right)$:

[TABLE]

implying:

[TABLE]

5.4.2 NLO (in $N$)/Non-Conformal Corrections

Up to NLO in $N$, in the IR, the potential ‘$V(IR,r_{h}=0)$’ is given by:

[TABLE]

The turning points of (5.4.2) up NLO in $N$ are given by:

\Biggl{[}\log\left(\frac{1}{\tilde{m}}\left[1-\frac{{g_{s}}M^{2}\log({y_{0}})(-{g_{s}}{N_{f}}(6{\log N}-36+\log(16777216))+36{g_{s}}{N_{f}}\log({y_{0}})+48\pi)}{128\pi^{2}N}\right]\right),\log(\delta^{2}-1)\Biggr{]}. After evaluation of the integral of $\sqrt{V(IR,r_{h}=0)}$ between the aforementioned turning points, in the large-$\tilde{m}$-limit, one obtains the following quantization condition:

[TABLE]

which yields:

[TABLE]

6 Glueball Masses from M theory

The glueball spectrum for spin $0^{++},1^{++}$ and $2^{++}$ is calculated in this section from the M-theory perspective both by considering $r_{h}\neq 0$ and an IR radial cut-off (for $r_{h}=0$) in the background. The $11$ dimensional M-theory action is given as:

[TABLE]

where $G_{4}=dC_{3}+A_{1}\wedge dB_{2}+dx_{10}\wedge dB_{2}$, and $C_{\mu\nu 10}^{M}=B_{\mu\nu}^{IIA},C_{\mu\nu\rho}^{M}=C_{\mu\nu\rho}^{IIA}$. Now, as shown in [2], no $F_{4}^{IIA}$ (to be obtained via a triple T-dual of type IIB $F_{1,3,5}$ where $F_{1}\sim F_{x/y/z},F_{3}\sim F_{xyr/\theta_{1}/\theta_{2}},F_{xzr/\theta_{1}/\theta_{2}},F_{yzr/\theta_{1}/\theta_{2}}$ and $F_{5}\sim F_{xyz\beta_{1}\beta_{2}}$ where $\beta_{i}=r/\theta_{i}$) can be generated. 555Consider $T_{x}$ followed by $T_{y}$ followed by $T_{z}$ where $T_{i}$ means T-dualizing along i-th direction. As an example, $T_{x}F_{x}^{IIB}\rightarrow{\rm non-dynamical\ 0-form\ field\ strength}^{IIA}\cite[cite]{[\@@bibref{}{kiritsis-book}{}{}]},T_{y}T_{x}F_{x}^{IIB}\rightarrow F_{y}^{IIB}$,

$T_{z}F_{y}^{IIB}\rightarrow F_{yz}^{IIA}$ implying one can never generate $F_{4}^{IIA}$ from $F_{1}^{IIB}$. As also an example consider $T_{x}F_{xy\beta_{i}}^{IIB}\rightarrow F_{y\beta_{i}}^{IIA},T_{y}F_{y\beta_{i}}^{IIA}\rightarrow F_{\beta_{i}}^{IIB},T_{z}F_{\beta_{i}}^{IIB}\rightarrow F_{\beta_{i}z}^{IIA}$ again not generating $F_{4}^{IIA}$; $T_{x}F_{xyz\beta_{1}\beta_{2}}^{IIB}\rightarrow F_{yz\beta_{1}\beta_{2}}^{IIA}$, $T_{y}F_{yz\beta_{1}\beta_{2}}^{IIA}\rightarrow F_{z\beta_{1}\beta_{2}}^{IIB},T_{z}F_{z\beta_{1}\beta_{2}}^{IIB}\rightarrow F_{\beta_{1}\beta_{2}}^{IIB}$; thus one can not generate $F_{4}^{IIA}$. Thus, the four-form flux $G_{4}=d\left(C_{\mu\nu 10}dx^{\mu}\wedge dx^{\nu}\wedge dx_{10}\right)+\left(A^{F_{1}}_{1}+A^{F_{3}}_{1}+A^{F_{5}}_{1}\right)\wedge H_{3}=H_{3}\wedge dx_{10}+A\wedge H_{3}$, where $C_{\mu\nu 10}\equiv B_{\mu\nu}$.

[TABLE]

Now, $H_{3}\wedge dx_{10}\wedge*_{11}\left(H_{3}\wedge A\right)=0$ as neither $H_{3}$ nor $A$ has support along $x_{10}$. Hence,

[TABLE]

where $H_{\mu\nu\rho 10}=H_{\mu\nu\rho}$, and

[TABLE]

with $H_{[\mu_{1}\mu_{2}\mu_{3}}A_{\mu_{4}]}\equiv H_{\mu_{1}\mu_{2}\mu_{3}}A_{\mu_{4}}-\left(H_{\mu_{2}\mu_{3}\mu_{4}}A_{\mu_{1}}-H_{\mu_{3}\mu_{4}\mu_{1}}A_{\mu_{2}}+H_{\mu_{4}\mu_{1}\mu_{2}}A_{\mu_{3}}\right)$.

It was shown in [2] that in the MQGP limit (2.2), in equation (LABEL:flux_action_D=11-iii), contribution of $H_{3}\wedge dx^{11}\wedge*_{11}(H_{3}\wedge dx^{11})$ is always dominated by $\sqrt{G}H_{\mu\nu\omega}G^{\mu\sigma}G^{\nu\psi}G^{\omega\alpha}G^{10\ 10}H_{\sigma\psi\alpha}$ term and in equation (117), contribution of $(H\wedge A)\wedge*_{11}(H\wedge A)$ is dominated by $\sqrt{G}H_{\theta_{1}\theta_{2}y}A_{\psi}G^{\theta_{1}\alpha}G^{\theta_{2}\rho}G^{y\sigma}G^{\psi\beta}H_{\alpha\rho\beta}$ term. Therefore, for simplicity in calculations, we assume that leading contribution in equations (LABEL:flux_action_D=11-iii) and (117) are governed by aforementioned terms.

The relevant inverse components of the 11-dimensional metric, in the MQGP limit (2.2), as worked out in [2] using which the most dominant contribution near $\theta_{1,2}=0$ in the equations (LABEL:flux_action_D=11-iii) and (117), as shown in [2], are given by the following analytical expressions:

[TABLE]

and

[TABLE]

where ${\cal F}(\alpha_{\theta_{1}},\alpha_{\theta_{2}};a;g_{s},M,N_{f})$ is a well-defined function of the parameters indicated. One hence notes that $\lim_{r_{\Lambda}\rightarrow\infty}\int_{0}^{\frac{\sqrt{4\pi g_{s}N}}{r_{h}}}\int_{r_{h}}^{r_{\Lambda}}\left(H\wedge A\right)\wedge*_{11}\left(H\wedge A\right)\sim\frac{r_{\Lambda}}{r_{h}}N^{\frac{27}{20}}$, and is UV-divergent. Also, this yields a large cosmological constant in the IR because : $\frac{G_{4}\wedge*G_{4}}{\sqrt{G}}\ni\frac{|H\wedge A|^{2}}{\sqrt{G}}\sim\frac{N^{\frac{17}{20}}}{r^{3}N^{\frac{17}{20}}}=\frac{1}{r^{3}}$. To take care of both these issues, from the discussion on holographic renormalizability of the $D=11$ supergravity action in [2], one sees that this term can be cancelled by a boundary counter term: $\int_{r=r_{\Lambda}}\sqrt{G}|G_{4}|^{2}$.

Now, using:

[TABLE]

and (118), one sees that one obtains a large-$N$ suppressed cosmological constant from the second term in $\frac{|H\wedge dx^{10}|^{2}}{\sqrt{G}}$ that remains small $\forall r>r_{h}$. To ensure one does not generate an $N$-enhanced cosmological constant from the first term in (118), one imposes the condition: $-4\alpha_{\theta_{2}}^{2}+27\alpha_{\theta_{1}}^{6}=0$, i.e., $\alpha_{\theta_{2}}=\frac{3^{\frac{3}{4}}\alpha_{\theta_{1}}^{\frac{3}{2}}}{\sqrt{2}}$.

One hence obtains the following flux-generated cosmological constant (with a slight abuse of notation):

[TABLE]

Metric Fluctuations:

The background metric $g^{(0)}_{\mu\nu}$ is linearly perturbed as $g_{\mu\nu}=g^{(0)}_{\mu\nu}+h_{\mu\nu}$. With this perturbation the equation of motion follows from the action (114) as:

[TABLE]

Now we assume the perturbation to have the following form: $h_{\mu\nu}=\epsilon_{\mu\nu}(r)e^{ikx_{1}}$. Clearly there is a $SO(2)$ rotational symmetry in the $x_{2}-x_{3}$ plane which allow us to classify different perturbations into three categories, namely tensor, vector and scalar type of metric perturbations.

The mass spectrum was obtained by (i) solving equation (122) and applying Neumann/Diriclet boundary condition near $r_{h}/r_{0}$, (ii) following the redefinition of variables in [34] and then considering the WKB quantization condition.

6.1 $0^{++}$ Glueball spectrum

The $0^{++}$ glueball in M-theory corresponds to scalar metric perturbations [73]:

[TABLE]

where $g_{tt}$, $g_{x_{i}x_{i}}$ and $g_{rr}$ are the metric components of the M-theory and is given in equation (2.4).

6.1.1 M-theory background with $r_{h}\neq 0$

Considering these components at leading order in $N$, and taking into account the above perturbation, we get the following differential equation for $q_{3}(r)$ with $q^{2}=-m^{2}$ from (122):

[TABLE]

(a) Spectrum from Neumann/Dirichlet Boundary Condition:

Equation (6.1.1), for $a=0.6r_{h}$ near $r=r_{h}$ (writing $m=\tilde{m}\frac{r_{h}}{L^{2}}$) simplifies to:

[TABLE]

Lets write the above equation of the following form,

[TABLE]

where we have, $p=\left(2-0.00666667\tilde{m}^{2}\right)$,  $s=\frac{\left(0.76-0.103333\tilde{m}^{2}\right)}{{r_{h}}}$.

The solution to (126) is given by:

[TABLE]

Setting $c_{2}=0$, one can verify that one satisfy the Neumann boundary condition: $h^{\prime}(r=r_{h})=0$ provided:

[TABLE]

implying:

[TABLE]

One can similarly show that by imposing Dirichlet boundary condition: $h(r=r_{h})=0$:

[TABLE]

(b) Spectrum using WKB method:

Following the redefinition of ([34]), equation (6.1.1) can be rewritten as a Schrodinger like form, where for $a=0.6r_{h}$ and setting $g_{s}=0.9,N\sim(g_{s})^{-39}\sim 100$ - in the MQGP limit of [2] - the ‘potential’ in the IR can be shown to be given by:

[TABLE]

The potential (131) is positive for $z\in(-\infty,\log(0.349)]$, but to remain within the IR, one truncates this domain to $z\in(-\infty,-2.526]$ and the same yields:

[TABLE]

or

[TABLE]

In the UV, one can show that:

[TABLE]

implying no turning points in the UV.

6.1.2 M-theory background with IR Cut-Off $r_{0}$

In this case, equation (6.1.1) is modified only by the limit $(r_{h},a)\rightarrow 0$. The equation in this limit is given as,

[TABLE]

(a) Spectrum from Neumann/Dirichlet Boundary Condition:

Near the cut-off at $r=r_{0}$, equation (135) is given by:

[TABLE]

whose solution is given by:

[TABLE]

Thus:

[TABLE]

Numerically/graphically we see that for $c_{1}=-0.509c_{2}$, one gets $q_{3}(r=r_{0},\tilde{m}\approx 4.1)=0$. We hence estimate the ground state of $0^{++}$ from metric fluctuations in M-theory to be $4.1\frac{r_{0}}{L^{2}}$.

(b) Spectrum using WKB method:

Defining $m=\tilde{m}\frac{\sqrt{y_{0}}}{L^{2}}$ and following [34], the potential term in the Schrödinger-like equation from (135) is given as,

[TABLE]

We hence see that in the large-$N$ limit, $V(z)<0$ and hence has no turning points. The WKB method a la [34] does not work in this case.

6.2 $2^{++}$ Glueball spectrum

To study the spectrum of spin $2^{++}$ glueball, we consider the tensor type of metric perturbations where the non-zero perturbations are given as:

[TABLE]

where $g_{x_{1}x_{1}}$ is as given in (2.4).

6.2.1 M-theory background with $r_{h}\neq 0$

Considering the tensor modes of metric perturbations and the M-theory metric components corrected upto NLO in $N$, as given in (2.4), we obtain a second order differential equation in $H(r)$ from (122),

[TABLE]

where we assume $k^{2}=-m^{2}$ with $m$ being the mass of the corresponding glueball.

(a) Spectrum from Neumann/Dirichlet Boundary Condition:

Near $r=r_{h}$, the solution to the above equation will be given on the same lines as 5.1 for $0^{--}$ glueballs, and the analog of (93) is:

[TABLE]

the solutions to which are given as:

[TABLE]

One can impose Dirichlet boundary condition: $H(r=r_{h})=0$, and show that,

$m_{n}^{2^{++}\ {(\rm Neumann)}}=m_{n+1}^{2^{++}\ {(\rm Dirichlet)}}$, for $n=0,1,2,..$.

(b) Spectrum from WKB method:

Using the variables of [34], The potential term in the schrodinger like equation for $2^{++}$ glueball can be obtained from (6.2.1) and in the IR region written in terms of $m=\tilde{m}\frac{r_{h}}{L^{2}}=\tilde{m}\frac{\sqrt{y_{h}}}{L^{2}}$, can be shown to given by:

[TABLE]

Now, negative $z$ implies being closer to the IR and the boundary would correspond to $r\leq\sqrt{3}a=0.6\sqrt{3}r_{h}$ corresponding to $z=-2.53$. It can be shown that for $\tilde{m}>7.211$ $V(z\in[-2.71,-2.53])>0$. Now $\int_{z=-2.71}^{z=-2.53}\sqrt{V_{IR}(z)}\approx 0$. Hence, the IR does not contribute to the $2^{++}$ glueball spectrum.

In the UV we must consider the limit $(z\rightarrow\infty)$. Moreover, in the UV $N_{f}=M=0$.

[TABLE]

whose turning points are: $\{\log\left(4.08+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right),\infty\}$, giving the WKB quantization as:

[TABLE]

implying:

[TABLE]

6.2.2 M-theory background with an IR Cut-Off $r_{0}$

Considering the limit $(r_{h},a\rightarrow 0)$ equation (6.2.1) is given by,

[TABLE]

(a) Neumann/Dirichlet Boundary Condition at $r=r_{0}$:

Up to LO in $N$ near $r=r_{0}$, the above equation is given by:

[TABLE]

The solution of (149) is given by:

[TABLE]

The Neumann boundary condition $H^{\prime}(r=r_{0})=0$, numerically yields that for $c_{1}=-0.509c_{2}$, the lightest $2^{++}$ glueball has a mass $1.137\frac{r_{0}}{L^{2}}$. Similarly, by imposing Dirichlet boundary condition: $H(r=r_{h})=0$, for $c_{1}=-0.509c_{2}$, the lightest $2^{++}$ glueball has a mass $0.665\frac{r_{0}}{L^{2}}$.

(b) Spectrum using WKB method:

Following [34], the ‘potential’ term, in the IR region, up to leading order in $N$ with $m=\tilde{m}\frac{r_{0}}{L^{2}}=\tilde{m}\frac{\sqrt{y_{0}}}{L^{2}}$ is given as:

[TABLE]

and $V_{IR}(z)>0$ for $z\in[\log\left(\frac{1}{\tilde{m}}+\frac{1}{\tilde{m}^{3}}\right),\log(\delta^{2}-1)]\approx[-\log\tilde{m},\log(\delta^{2}-1)]$. Hence WKB quantization condition gives,

[TABLE]

implying:

[TABLE]

In the UV, we have:

[TABLE]

and $V_{UV}(z)>0$ for $z>\log\left(\frac{1}{8}\left(-{\tilde{m}}^{2}+\sqrt{{\tilde{m}}^{4}+28{m_{0}}^{2}-140}+10\right)\right)=\log\left(3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right)$. Further, $\int_{\log 3}^{\infty}\sqrt{V_{UV}(z)}=\frac{\pi\tilde{m}}{8\sqrt{3}}$, implying:

[TABLE]

(c) NLO-in-$N$/Non-Conformal Corrections using WKB method:

The ‘potential’ inclusive of NLO-in-$N$ terms, in the IR region in the $r_{h}=0$ limit, is given by:

[TABLE]

whose turning points are given by:

$\left[\log\left\{\frac{1-\frac{{g_{s}}M^{2}\log({y_{0}})({g_{s}}{N_{f}}(-12{\log N}+72+\log(4096)-15\log(16))+72{g_{s}}{N_{f}}\log({y_{0}})+96\pi)}{256\pi^{2}N}}{{m_{0}}}\right\},\log\left(\delta^{2}-1\right)\right]$. The integral of

$\sqrt{V(IR,r_{h}=0)}$ between these turning points, in the large-$\tilde{m}$ limit, yields the same spectrum as $0^{--}$ up to NLO in $N$.

6.3 Spin-$1^{++}$ Glueball spectrum

Here we need to consider the vector type of metric perturbation with the non-zero components given as: $h_{ti}=h_{it}=g_{x_{1}x_{1}}G(r)e^{ikx_{1}}$, $i={x_{2},x_{3}}$.

6.3.1 M theory background with $r_{h}\neq 0$

Substituting the above ansatz for the perturbation in (122), the differential equation in $G(r)$ is given with $k^{2}=-m^{2}$ as,

[TABLE]

(a) Neumann/Diriclet boundary condition at $r=r_{h}$:

Near $r=r_{h}$, equation (6.3.1) up to LO in $N$, is given by:

[TABLE]

whose solution is given by:

[TABLE]

Imposing Neumann boundary condition at $r=r_{h}$ yields:

[TABLE]

Considering $p=\frac{0.01\tilde{m}^{2}}{\sqrt{5.0016-0.57\tilde{m}^{2}}}-\frac{1}{\sqrt{5.0016-0.57\tilde{m}^{2}}}+1$ and setting $c_{2}=0$ in (6.3.1), and then using $\lim_{z\rightarrow 0}U(p,2,z\sim 0)\sim\frac{z^{-1}\ _{1}F_{1}(p-1;0;z)}{\Gamma(p)}$, one notes that one can satisfy the Neumann boundary condition at $r=r_{h}$ provided $\lim_{z\rightarrow 0}\ {}_{1}F_{1}(p-1;0;z)=\lim_{b\rightarrow 0}\lim_{z\rightarrow 0}\ {}_{1}F_{1}(p-1;b;z)$ (i.e. first set $z$ to 0 and then $b$), $p=-n\in\mathbb{Z}^{-}$. Hence:

[TABLE]

One can show that one obtains the same spectrum as (6.3.1) even upon imposing Dirichlet boundary condition: $G(r=r_{h})=0$.

(b) Spectrum using WKB method:

Following [34], the ‘potential’ $V$ in the Schrödinger-like equation working with the dimensionless mass variable $\tilde{m}$ defined via: $m=\tilde{m}\frac{r_{h}}{L^{2}}=\tilde{m}\frac{\sqrt{y_{h}}}{L^{2}}$, in the IR, can be shown to be given by:

[TABLE]

The zeros of the potential as function of $e^{z}$, in (162) are given by:

$\frac{-1.8765\times 10^{14}\tilde{m}^{2}\pm 6.022232598554301\times 10^{-7}\sqrt{9.70917\times 10^{40}\tilde{m}^{4}+1.26219\times 10^{44}\tilde{m}^{2}-5.04877\times 10^{44}}+1.8765\times 10^{16}}{5.70456\times 10^{16}-5.6295\times 10^{15}\tilde{m}^{2}}=-\frac{25}{\tilde{m}^{2}}+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right),$ $0.07+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)$; the first not being permissible. Now, in the IR $r\in[r_{h},\sqrt{3}a\approx 0.6\sqrt{3}r_{h}]$ or in terms of $z:(-\infty,-2.526]$. Therefore the allowed domain of integration is: $[\log(0.07),-2.526]$. Thus, in the IR:

[TABLE]

implying a null contribution to the WKB quantization condition in the IR.

In the UV, the potential is given by:

[TABLE]

which for $b=0.6$ obtains:

[TABLE]

The zeros of the potential in the UV, as a function of $e^{z}$,in (165) are given by:

$-0.125\tilde{m}^{2}\pm 0.005\sqrt{625.\tilde{m}^{4}+26950.\tilde{m}^{2}-185671.}+1.385$ $=(-0.25\tilde{m}^{2}-1.31,4.08+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right))$; the former not being permissible. Hence, the allowed domain of integration over which the potential is positive, is: $([\log(4.08),\infty)$. Performing a large-$\tilde{m}$-expansion, one obtains:

[TABLE]

yielding:

[TABLE]

6.3.2 M theory background with an IR Cut-Off $r_{0}$

(a) Neumann/Diriclet boundary condition at $r=r_{0}$:

Considering the limit of $(r_{h},a)\rightarrow 0$ in equation (6.3.1) up to LO in $N$ and imposing Neumann boundary condition at the IR cut-off $r=r_{0}$, yields isospectrality with $2^{++}$ glueball spectrum at $r_{h}=0$.

(b) Spectrum using WKB method:

Using the redefinition of [34], the ‘potential’ up to leading order in $N$ is given by:

[TABLE]

In the IR region we get the potential as:

[TABLE]

giving the turning points as, $z\in[\log\left(\frac{1}{\tilde{m}}+{\cal O}\left(\frac{1}{\tilde{m}^{3}}\right)\right)\approx-\log\tilde{m},\log(\delta^{2}-1)]$ and the WKB quantization condition becomes:

[TABLE]

Therefore:

[TABLE]

Further, in the UV:

[TABLE]

implying that $V_{UV}(z)>0$ for $z\in[\log\left(\frac{1}{8}\left(-\tilde{m}^{2}+\sqrt{\tilde{m}^{4}+28\tilde{m}^{2}-140}+10\right)=\log\left(3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right),\infty\right)$. This yields the following WKB quantization condition:

[TABLE]

which obtains:

[TABLE]

(c) NLO in $N$/Non-Conformal Corrections using WKB method:

In the IR region, the ‘potential’ including NLO-in-$N$ corrections in the $r_{h}=0$ limit, is given by:

[TABLE]

whose turning points are:

$\left[\log\left(\frac{1}{\tilde{m}}\left[1-\frac{{g_{s}}M^{2}\log({y_{0}})({g_{s}}{N_{f}}(-12{\log N}+72+\log(4096)-15\log(16))+72{g_{s}}{N_{f}}\log({y_{0}})+96\pi)}{256\pi^{2}N}\right]\right),\log(\delta^{2}-1)\right]$, and the integral of $\sqrt{V(IR,r_{h}=0)}$ between these turning points yields an isospectrality with the $0^{--}$ and $2^{++}$ NLO-in-$N$ spectrum.

7 $2^{++}$ Glueball Masses from Type IIB

7.1 $r_{h}\neq 0$

The $10$-dimensional type IIB supergravity action in the low energy limit is given by,

[TABLE]

where $\phi$ is the dilaton, $G_{MN}$ is the $10$-d metric and $F_{1}$, $H_{3}$, $\widetilde{F_{3}}$, $\widetilde{F_{5}}$ are different fluxes.

The five form flux $\widetilde{F_{5}}$ and the three form flux $\widetilde{F_{3}}$ are defined as,

[TABLE]

where $F_{5}$ and $F_{3}$ are sourced by the $D_{3}$ and $D_{5}$ branes respectively. $B_{2}$ is the NS-NS two form and $C_{0}$ is the axion. The three form fluxes $\widetilde{F_{3}}$, $H_{3}$, the two form $B_{2}$ and the axion $C_{0}$ are given as [1] - see (2.1). Now varying the action in (176) with respect to the metric $g_{\mu\nu}$ one get the following equation of motion,

[TABLE]

we consider the following linear perturbation of the metric,

[TABLE]

where as before $\mu,\nu=\{t,x_{1},x_{2},x_{3},r,\theta_{1},\theta_{2},\phi_{1},\phi_{2},\psi\}$. Here the only non zero component according to the tensor mode of metric fluctuation is $h_{x_{2}x_{3}}$. Since the non zero components of $F_{1}$, $\widetilde{F_{3}}$ and $H_{3}$ has no indices as $\{x_{2},x_{3}\}$, the final equation of motion gets simplified and can be shown to be given as:

[TABLE]

Working at a particular value of $\theta_{1}$ and $\theta_{2}$ given as: $\theta_{1}=N^{-1/5}$ and $\theta_{2}=N^{-3/10}$, the square of different fluxes figuring in (180) at the lowest-order in $N$, are given in (A). Writing the perturbation $h_{x_{2}x_{3}}$ as $h_{x_{2}x_{3}}=\frac{r^{2}}{2(g_{s}\pi N)^{1/2}}H(r)e^{ikx_{1}}$, (180) reduces to the following second order differential equation in $H(r)$:

[TABLE]

7.1.1 Mass Spectrum from Neumann Boundary Condition at $r=r_{h}$

To get a sensible answer, one has to perform a small-$T$ expansion when one rewrites and solves (7.1) around $r=r_{h}$. This time around the analog of (93) becomes:

[TABLE]

which for $g_{s}=0.8,N=g_{s}^{-39}\sim 6000,M=3$ yields:

[TABLE]

So, one obtains an approximate match between the ground state $2^{++}$ mass from M theory and type IIB string theory.

7.1.2 WKB Quantization Method

Using the variables of [34], the ‘potential’, defining $m=\tilde{m}\frac{\sqrt{y_{h}}}{L^{2}}$, yields:

[TABLE]

In the IR, the potential is given by:

[TABLE]

and in the IR, for $\tilde{m}>2.986$, $V(IR<T)>0$ for $z\in[\log\left(0.067+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right)\right)\approx-2.708,-2.526]$ and: $\int_{2.708}^{-2.526}\sqrt{V(IR,T)}\approx 0$ - hence the IR provides no contribution to the WKB quantization.

In the UV:

[TABLE]

the turning points for $\tilde{m}>7.141$ are $\left(\log(4.08+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right),0.25\tilde{m}^{2}-1.77)\right)$. Hence,

[TABLE]

which obtains:

[TABLE]

Hence, the string theory $2^{++}$ glueball is isospectral with $0^{++}$; in the large $n$-limit of the spectrum, the M-theory and type IIB spectra coincide.

7.2 $r_{h}=0$

7.2.1 WKB Method

Once again defining $r=\sqrt{y},r_{h}=\sqrt{y_{0}},y=y\left(1+e^{z}\right)$ analogous to their $r_{h}\neq 0$ redefinitions of [34], the ‘potential’ in the IR, is given by:

[TABLE]

The domain in the IR over which $V(r_{h}=0)>$ is: $[-\frac{1}{2}\log(5+\tilde{m}^{2}),\log(\delta^{2}-1)]$ and:

[TABLE]

yielding:

[TABLE]

In the UV

[TABLE]

the zeros of which, as functions of $e^{z}$, are at $(\frac{1}{8}\left(-\tilde{m}^{2}\pm\sqrt{\tilde{m}^{4}+28\tilde{m}^{2}-140}+10\right)$ $=(-\frac{\tilde{m}^{2}}{4}-\frac{1}{2},3+{\cal O}\left(\frac{1}{\tilde{m}^{2}}\right))$. Hence the domain of integration over which $V(UV,r_{h}=0)>0$ is: $[\log 3,\infty)$. Therefore:

[TABLE]

yielding:

[TABLE]

7.2.2 NLO-in-$N$/Non-Conformal Corrections in the IR in the $r_{h}=0$ Limit

The ‘potential’ inclusive of the NLO-in-$N$ corrections in the IR in the $r_{h}=0$ limit, reads as:

[TABLE]

whose turning points, in the large-$\tilde{m}$ limit are: $[\log\left\{\frac{1}{\tilde{m}}+{\cal O}\left(\frac{1}{\tilde{m}^{3}}\right)\right\},\log(\delta^{2}-1)]\approx[-\log\tilde{m},\log\left(\delta^{2}-1\right)]$. Now,

[TABLE]

which yields the same LO spectrum as $0^{--},1^{++}$ and the $2^{++}$ spectrum obtained from M theory. Hence, the type IIB at $r_{h}=0$ is unable to capture the non-conforamal NLO-in-$N$ corrections in the $2^{++}$ corrections, as the same either precisely cancel out or are $\frac{1}{\tilde{m}}$-suppressed in the large-$\tilde{m}$ limit, in a IIB computation.

8 Summary and Discussion

Supergravity calculations of glueball spectra in top-down holographic duals of large-$N$ thermal QCD at finite string/gauge coupling and not just $g_{s}N\gg 1$, have thus far, been missing in the literature. Such a limit is particularly relevant to sQGP [4]. This work fills in this gap by working out the spectra of $0^{++},0^{-+},0^{--},1^{++},2^{++}$ glueballs in a type IIB/delocalized SYZ IIA mirror/(its) M-theory (uplift) model corresponding to the top-down holographic dual of [1] in the MQGP limit introduced in [2]. As discussed in 2.1, towards the end of a Seiberg duality cascade in the IR, despite setting $M(=N_{c}$ in the IR) to three in the MQGP limit, due a flavor-color enhancement of the length scale as compared to the KS model, one can trust supergravity calculations without worrying about stringy corrections. Further, in the MQGP limit, all physical quantities as seen in all calculations in [5], receive non-conformal corrections that appear at the NLO in $N$ and display a universal $\frac{g_{s}M^{2}(g_{s}N_{f})}{N}$-suppression.

It should be noted that a numerical computation like the ‘shooting computation’ used in a lot of holographic glueball spectrum computations will not be feasible to use for the following reason. In the ‘shooting method’, like [61], one can first solve the EOMs in the UV using the infinite series/Frobenius method and then numerically (via Euler’s method, etc.) obtain the solution at the horizon where one imposes a Neumann boundary condition. By matching the value obtained by numerically ‘shooting’ from the UV to the horizon in the IR and matching the radial derivative of the solution so obtained to zero, one can obtain quantized values of the glueball masses. The caveat is that one should have at hand the exact radial profile of the effective number of fractional $D3$-branes ($D5$-branes wrapping the small two-cycle) and the number of flavor branes which would correctly interpolate between $(M,N_{f})=(0,0)$ in the UV and ($M=3,N_{f}=2$) in the IR. But, we do not have this information - we know the values in the IR and the UV but not for the interpolating region. Hence, numerical methods such as the ‘shooting method’ could at best be used, to obtain only the LO-in-$N$ results, not the NLO-in-$N$ results which is one of the main objectives of our computations.

The summary of all calculations is given in tables 1 (and Fig. 2) and 3 - the former table/graph having to do with a WKB quantization calculation using the coordinate/field redefinitions of [34] and the latter table having to do with obtaining the mass spectrum by imposing Neumann/Dirichlet boundary condition at $r_{h}$/IR cut-off $r_{0}$. Some of the salient features of the results are given as separate bullets.

It should be noted that the last two columns in Tables 1 and 3 have been prepared in the same spirit as the last columns in Table 2 of [75].

The $r_{h}\neq 0$ glueball spectra is plotted in Figure 2.

Some of the salient features of Table 1 and Figure 2 are presented below:

Interestingly, via a WKB quantization condition using coordinate/field redefinitions of [34], the lightest $0^{++}$ glueball spectrum for $r_{h}\neq 0$ coming from scalar metric fluctuations in M theory compares rather well with the $N\rightarrow\infty$ lattice results of [57] - refer to Table 2. Also, similar to [58], the $0^{++}$ coming from the scalar fluctuations of the M theory metric is lighter than the $0^{++}$ coming from type IIB dilaton fluctuations. Further, interestingly, one can show that by using the coordinate and field redefinitions of [33] when applied to the EOM for dilaton fluctuation to yield a WKB quantization condition, for $a=0.6r_{h}$ - as in [5] - one obtains a match with the UV limit of the $0^{++}$ glueball spectrum as obtained in [34]. For our purpose, the method based on coordinate/field redefinitions of [33], is no good for obtaining the $0^{++}$ glueball ground state and was not used for any other glueball later on in subsequent calculations in this paper.

Also, from Table 1/Figure 2, $m_{n>0}^{2^{++}}>m_{n>0}^{0^{++}}$(scalar metric perturbations), similar to [58]. 3. 3.

The higher excited states of the type IIA $0^{-+}$ glueball, for both $r_{h}\neq 0$ and $r_{h}=0$, are isospectral. This is desirable because large-$n$ corresponds to the UV and that takes one away from the BH geometry, i.e., towards $r_{h}=0$. 4. 4.

The non-conformal corrections up to NLO in $N$, have a semi-universal behavior of $\frac{(g_{s}M^{2})(g_{s}N_{f})\log r_{0}}{N}$ and turn out to be multiplied by a numerical pre-factor of ${\cal O}(10^{-2})$; we could disregard the same in the MQGP limit. 5. 5.

As per a more recent lattice calculation [59]666We thank P.Majumdar for bringing this reference to our attention., the $0^{++}$-glueball has a mass $4.16\pm 0.11\pm 0.04$ (in units of the reciprocal of the ‘hadronic scale parameter’ of [60]), which compares rather well with $m_{n=0}^{0^{++}}=4.267$ (in units of $\frac{r_{h}}{L^{2}}$) of Table 2 coming from scalar fluctuations of the M theory metric. Similarly, the $0^{-+}$-glueball in [59] has a mass $6.25\pm 0.06\pm 0.06$ and from Table 1, which matches rather nicely with $m_{n=0}^{0^{-+}}(\delta=1.26)=6.25$ (in units of $\frac{r_{0}}{L^{2}}$) of Table 1 coming from type IIA one-form fluctuation. 6. 6.

The ground state and the $n\gg 1$ excited states of $1^{++}$ and $0^{--}$ glueballs are isospectral. 7. 7.

The higher excited $r_{h}\neq 0$ $2^{++}$ glueball states corresponding to metric fluctuations of the M-theory metric and the ones corresponding to fluctuations of the type IIB metric, are isospectral. The $r_{h}=0$ $2^{++}$ glueball states corresponding to metric fluctuations of the M-theory/type IIB string theory, are isospectral. Further, it turns out that due to internal cancellation of terms and $\frac{1}{\tilde{m}}$-suppression, a type IIB $r_{h}=0$ $2^{++}$ glueball spectrum, unlike an M-theoretic computation, is unable to capture the NLO-in-$N$ corrections to the LO-in-$N$ type IIB $2^{++}$ glueball spectrum. 8. 8.

$m^{2^{++}}_{n}({\rm NLO},r_{h}=0)=m_{n}^{1^{++}}({\rm NLO},r_{h}=0)\stackrel{{\scriptstyle n\gg 1}}{{\longrightarrow}}m_{n}^{0^{--}}({\rm NLO},r_{h}=0)$, where the ‘NLO’ implies equality with the inclusion of NLO-in-$N$ corrections.

Some salient features of Table 3 are presented below:

• •

The following is the comparison of ratios of $0^{--}$ glueball masses obtained in this work from Neumann/Dirichlet boundary conditions at the horizon, with [61]:

Hence, for higher excited states, the ratio of masses of successive excited states approaches unity faster as per our results as compared to [61].

• •

From a comparison of results in Tables 1/2 or Figures 2 with $N\rightarrow\infty$ lattice results, it appears that WKB quantization-based spectra are closer to $N\rightarrow\infty$ lattice results than the computations involving imposing Neumann/Dirichlet boundary conditions at the horizon/IR cut-off. In particular, it is pleasantly surprising that the WKB quantization method applied to the $0^{++},0^{-+}$ glueball spectra, is able to provide a good agreement (in fact for the lightest $0^{++}$ glueball spectrum, better than the classic computations of [61]) with lattice results even for the ground and the lower excited states.

Acknowledgements

KS is supported by is supported by a senior research fellowship (SRF) from the Ministry of Human Re-source and Development (MHRD), Govt. of India. VY is supported by a junior research fellowship (JRF) from the University Grants Commission, Govt. of India. We would like to thank P. Majumdar and S. Gupta for useful communications, H. B. Filho for a useful clarification, and R. Das for participating in the material in 3.2.1 as part of his Masters’ project work. One of us (AM) would like to dedicate this work to the memory of his father, the late R. K. Misra who was and (whose memory) still remains a true driving force and inspiration for him.

Appendix A $\widetilde{F^{2}_{5}},\widetilde{F^{2}_{3}},H_{3}^{2}$

The expressions of squares of various fluxes that figure in the EOM (180) are given below for ready reference:

[TABLE]