Skip to content

Three-Point Functions and Cubic Couplings

The previous two pages developed the two ingredients needed for the first interacting calculation in AdS/CFT:

  1. A boundary source ϕ(0)(x)\phi_{(0)}(x) determines a classical bulk solution through a bulk-to-boundary propagator.
  2. The renormalized on-shell action generates CFT correlation functions.

Two-point functions are mostly kinematics plus normalization. For a scalar primary, conformal symmetry fixes

O(x1)O(x2)=COx122Δ,xij=xixj.\langle \mathcal O(x_1)\mathcal O(x_2)\rangle = \frac{C_{\mathcal O}}{|x_{12}|^{2\Delta}}, \qquad x_{ij}=x_i-x_j.

The first genuinely dynamical CFT data appear in three-point functions. For scalar primaries,

O1(x1)O2(x2)O3(x3)=C123x12Δ1+Δ2Δ3x13Δ1+Δ3Δ2x23Δ2+Δ3Δ1\boxed{ \langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3) \rangle = \frac{C_{123}}{ |x_{12}|^{\Delta_1+\Delta_2-\Delta_3} |x_{13}|^{\Delta_1+\Delta_3-\Delta_2} |x_{23}|^{\Delta_2+\Delta_3-\Delta_1}} }

at separated points. Conformal symmetry fixes the entire position dependence. The number C123C_{123} is dynamical. Equivalently, after normalizing two-point functions, it is the OPE coefficient.

In holography, the leading large-NN contribution to C123C_{123} is computed by a cubic interaction in the bulk:

Sint=λ123AdSdd+1Xgϕ1ϕ2ϕ3.S_{\rm int} = \lambda_{123} \int_{\rm AdS}d^{d+1}X\sqrt g\,\phi_1\phi_2\phi_3.

The corresponding diagram is the simplest contact Witten diagram: three bulk-to-boundary propagators meet at one integrated bulk point.

A three-point contact Witten diagram from a cubic bulk coupling

A cubic bulk vertex computes a scalar CFT three-point function. Each boundary source Ji(xi)J_i(x_i) creates a classical bulk field through KΔi(X;xi)K_{\Delta_i}(X;x_i); the interaction point XX is integrated over AdS. The result has the unique scalar conformal structure, while the overall coefficient is determined by the renormalized cubic coupling.

This page is conceptually important because it explains how the dynamical bulk Lagrangian becomes the CFT operator algebra. Masses determine dimensions. Kinetic terms determine two-point normalizations. Cubic couplings determine three-point coefficients.

A CFT is not defined only by its list of scaling dimensions. For local scalar primaries, the basic data include

Δi,Ci,Cijk,\Delta_i, \qquad C_i, \qquad C_{ijk}, \qquad \ldots

where

Oi(x)Oj(0)=Ciδijx2Δi\langle \mathcal O_i(x)\mathcal O_j(0)\rangle = \frac{C_i\delta_{ij}}{|x|^{2\Delta_i}}

in a diagonal basis. The operator product expansion then has the schematic form

Oi(x)Oj(0)kCijkCkOk(0)xΔi+ΔjΔk+descendants.\mathcal O_i(x)\mathcal O_j(0) \sim \sum_k \frac{C_{ijk}}{C_k} \frac{\mathcal O_k(0)}{|x|^{\Delta_i+\Delta_j-\Delta_k}} + \text{descendants}.

If we define normalized operators

O^i=OiCi,\widehat{\mathcal O}_i = \frac{\mathcal O_i}{\sqrt{C_i}},

then

O^i(x)O^j(0)=δijx2Δi,\langle \widehat{\mathcal O}_i(x)\widehat{\mathcal O}_j(0)\rangle = \frac{\delta_{ij}}{|x|^{2\Delta_i}},

and the normalized OPE coefficient is

λijkCFT=CijkCiCjCk.\boxed{ \lambda_{ijk}^{\rm CFT} = \frac{C_{ijk}}{\sqrt{C_iC_jC_k}}. }

This distinction is not pedantic. Holographic computations naturally produce numbers in a chosen bulk-field normalization. To compare with bootstrap conventions, one must divide by the square roots of the two-point coefficients.

For scalar primaries, the three-point position dependence is fixed by translations, rotations, dilatations, and special conformal transformations. It is useful to define

δ12=Δ1+Δ2Δ32,δ13=Δ1+Δ3Δ22,δ23=Δ2+Δ3Δ12.\delta_{12}=\frac{\Delta_1+\Delta_2-\Delta_3}{2}, \qquad \delta_{13}=\frac{\Delta_1+\Delta_3-\Delta_2}{2}, \qquad \delta_{23}=\frac{\Delta_2+\Delta_3-\Delta_1}{2}.

Then the scalar three-point function is

O1(x1)O2(x2)O3(x3)=C123x122δ12x132δ13x232δ23.\langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3)\rangle = \frac{C_{123}}{|x_{12}|^{2\delta_{12}}|x_{13}|^{2\delta_{13}}|x_{23}|^{2\delta_{23}}}.

The goal of the bulk calculation is to compute C123C_{123}.

Work in Euclidean Poincare AdSd+1_{d+1},

ds2=L2z2(dz2+dxμdxμ),z>0.ds^2 = \frac{L^2}{z^2}(dz^2+d x^\mu d x_\mu), \qquad z>0.

Consider three scalar fields ϕi\phi_i with masses

mi2L2=Δi(Δid),i=1,2,3.m_i^2L^2 = \Delta_i(\Delta_i-d), \qquad i=1,2,3.

A convenient schematic Euclidean action is

SE=i=13Ni2dd+1Xg(gabaϕibϕi+mi2ϕi2)+λ123dd+1Xgϕ1ϕ2ϕ3+Sct.S_E = \sum_{i=1}^3 \frac{\mathcal N_i}{2} \int d^{d+1}X\sqrt g \left( g^{ab}\partial_a\phi_i\partial_b\phi_i+m_i^2\phi_i^2 \right) + \lambda_{123} \int d^{d+1}X\sqrt g\,\phi_1\phi_2\phi_3 + S_{\rm ct}.

Here Ni\mathcal N_i fixes the two-point normalization of Oi\mathcal O_i. The coupling λ123\lambda_{123} is written in the same field normalization. In a top-down compactification, Ni\mathcal N_i and λ123\lambda_{123} are not arbitrary: they come from reducing the ten- or eleven-dimensional supergravity or string action. In a bottom-up model, they are phenomenological parameters specifying the CFT data one wants to model.

For three distinct fields the cubic term has no symmetry factor. If the interaction is instead

Sint=λ3!gϕ3,S_{\rm int}=\frac{\lambda}{3!}\int\sqrt g\,\phi^3,

the three functional derivatives with respect to the same source produce a factor 3!3!, so the final answer again contains λ\lambda rather than λ/3!\lambda/3!.

At leading order in λ123\lambda_{123}, the classical solution sourced by Ji(x)J_i(x) is just the free solution:

ϕi(X)=ddxKΔi(X;x)Ji(x)+O(J2).\phi_i(X) = \int d^d x\,K_{\Delta_i}(X;x)J_i(x)+O(J^2).

The normalized Euclidean bulk-to-boundary propagator is

KΔ(z,x;y)=CΔ(zz2+xy2)Δ,CΔ=Γ(Δ)πd/2Γ(Δd/2).K_\Delta(z,x;y) = C_\Delta \left( \frac{z}{z^2+|x-y|^2} \right)^\Delta, \qquad C_\Delta = \frac{\Gamma(\Delta)}{\pi^{d/2}\Gamma(\Delta-d/2)}.

The normalization is chosen so that

limz0zΔdKΔ(z,x;y)=δ(d)(xy).\lim_{z\to0}z^{\Delta-d}K_\Delta(z,x;y) = \delta^{(d)}(x-y).

Substituting the free solutions into the cubic term gives the cubic part of the renormalized generating functional at separated points:

Sren(3)[J1,J2,J3]=λ123AdSdd+1Xgi=13[ddxiKΔi(X;xi)Ji(xi)]+local contact terms.S_{\rm ren}^{(3)}[J_1,J_2,J_3] = \lambda_{123} \int_{\rm AdS}d^{d+1}X\sqrt g \prod_{i=1}^3 \left[ \int d^d x_i\,K_{\Delta_i}(X;x_i)J_i(x_i) \right] + \text{local contact terms}.

Using the convention of the previous pages, where functional derivatives of SrenS_{\rm ren} generate correlators, the separated-point three-point function is

O1(x1)O2(x2)O3(x3)conn=λ123AdSdd+1XgKΔ1(X;x1)KΔ2(X;x2)KΔ3(X;x3)+contact terms.\boxed{ \langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3) \rangle_{\rm conn} = \lambda_{123} \int_{\rm AdS}d^{d+1}X\sqrt g\, K_{\Delta_1}(X;x_1)K_{\Delta_2}(X;x_2)K_{\Delta_3}(X;x_3) + \text{contact terms}. }

Different Euclidean sign conventions move an overall sign between the action and the generating functional. That sign is not the interesting physics here. The physical data are the separated-point coefficient, the normalization of the operators, and the tensor structure.

Set L=1L=1 for the moment and strip off the normalization constants CΔiC_{\Delta_i}. Define

DΔ1Δ2Δ3(x1,x2,x3)=0dzzd+1ddxi=13(zz2+xxi2)Δi.D_{\Delta_1\Delta_2\Delta_3}(x_1,x_2,x_3) = \int_0^\infty\frac{dz}{z^{d+1}} \int d^d x \prod_{i=1}^3 \left( \frac{z}{z^2+|x-x_i|^2} \right)^{\Delta_i}.

Let

Σ=Δ1+Δ2+Δ3.\Sigma=\Delta_1+\Delta_2+\Delta_3.

For non-extremal dimensions in the domain where the integral converges, and elsewhere by analytic continuation plus holographic renormalization, the result is

DΔ1Δ2Δ3=πd/22Γ(Σd2)Γ(δ12)Γ(δ13)Γ(δ23)Γ(Δ1)Γ(Δ2)Γ(Δ3)1x122δ12x132δ13x232δ23.\boxed{ D_{\Delta_1\Delta_2\Delta_3} = \frac{\pi^{d/2}}{2} \frac{ \Gamma\left(\frac{\Sigma-d}{2}\right) \Gamma(\delta_{12}) \Gamma(\delta_{13}) \Gamma(\delta_{23}) }{ \Gamma(\Delta_1)\Gamma(\Delta_2)\Gamma(\Delta_3) } \frac{1}{ |x_{12}|^{2\delta_{12}} |x_{13}|^{2\delta_{13}} |x_{23}|^{2\delta_{23}} }. }

Restoring the propagator normalizations and the AdS radius gives

O1(x1)O2(x2)O3(x3)=C123x122δ12x132δ13x232δ23,\langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3)\rangle = \frac{\mathfrak C_{123}}{ |x_{12}|^{2\delta_{12}} |x_{13}|^{2\delta_{13}} |x_{23}|^{2\delta_{23}} },

with

C123=λ123Ld+1i=13CΔiπd/22Γ(Σd2)Γ(δ12)Γ(δ13)Γ(δ23)Γ(Δ1)Γ(Δ2)Γ(Δ3)\boxed{ \mathfrak C_{123} = \lambda_{123}L^{d+1} \prod_{i=1}^3 C_{\Delta_i} \frac{\pi^{d/2}}{2} \frac{ \Gamma\left(\frac{\Sigma-d}{2}\right) \Gamma(\delta_{12}) \Gamma(\delta_{13}) \Gamma(\delta_{23}) }{ \Gamma(\Delta_1)\Gamma(\Delta_2)\Gamma(\Delta_3) } }

up to the field-normalization and sign conventions already stated.

This formula is one of the cleanest examples of the AdS/CFT dictionary. The bulk integral knows about the boundary conformal group because AdS isometries act as boundary conformal transformations. The bulk-to-boundary propagator transforms like a primary of dimension Δ\Delta, and the invariant AdS measure does the rest.

There is a quick conceptual derivation of the position dependence. Under an AdS isometry corresponding to a boundary conformal transformation xxx\mapsto x', the bulk point transforms as XXX\mapsto X', the measure is invariant,

dd+1Xg=dd+1Xg,d^{d+1}X\sqrt g=d^{d+1}X'\sqrt g,

and the bulk-to-boundary propagator transforms as

KΔ(X;x)=xxΔ/dKΔ(X;x).K_\Delta(X;x) = \left|\frac{\partial x'}{\partial x}\right|^{\Delta/d} K_\Delta(X';x').

Therefore the integral transforms as a product of three scalar primaries. There is no conformal cross-ratio for three points, so the answer must be the unique scalar three-point structure. The only nontrivial information left is the coefficient.

This is a good diagnostic principle:

If a tree-level scalar contact diagram in pure AdS gives a separated-point answer that is not the scalar CFT three-point structure, the calculation has violated the boundary condition, the normalization, or conformal covariance.

A worked special case: three equal dimensions

Section titled “A worked special case: three equal dimensions”

For Δ1=Δ2=Δ3=Δ\Delta_1=\Delta_2=\Delta_3=\Delta, one has

δ12=δ13=δ23=Δ2,Σ=3Δ.\delta_{12}=\delta_{13}=\delta_{23}=\frac{\Delta}{2}, \qquad \Sigma=3\Delta.

The stripped contact integral becomes

DΔΔΔ=πd/22Γ(3Δd2)Γ(Δ2)3Γ(Δ)31x12Δx13Δx23Δ.D_{\Delta\Delta\Delta} = \frac{\pi^{d/2}}{2} \frac{ \Gamma\left(\frac{3\Delta-d}{2}\right) \Gamma\left(\frac{\Delta}{2}\right)^3 }{ \Gamma(\Delta)^3 } \frac{1}{|x_{12}|^{\Delta}|x_{13}|^{\Delta}|x_{23}|^{\Delta}}.

Including the propagator normalization gives

CΔΔΔ=λLd+1[Γ(Δ)πd/2Γ(Δd/2)]3πd/22Γ(3Δd2)Γ(Δ2)3Γ(Δ)3.\mathfrak C_{\Delta\Delta\Delta} = \lambda L^{d+1} \left[ \frac{\Gamma(\Delta)}{\pi^{d/2}\Gamma(\Delta-d/2)} \right]^3 \frac{\pi^{d/2}}{2} \frac{ \Gamma\left(\frac{3\Delta-d}{2}\right) \Gamma\left(\frac{\Delta}{2}\right)^3 }{ \Gamma(\Delta)^3 }.

For example, for three marginal scalar operators in d=4d=4, so Δ=4\Delta=4, the normalized integral part is proportional to

1x124x134x234,\frac{1}{|x_{12}|^4|x_{13}|^4|x_{23}|^4},

as conformal invariance demands. The number multiplying this structure depends on the normalization of the bulk fields and on the cubic coupling.

Real supergravity actions rarely contain only ϕ1ϕ2ϕ3\phi_1\phi_2\phi_3. They contain derivative interactions, mixing terms, gauge couplings, and terms produced by dimensional reduction. A simple example is

Sint=a123dd+1Xgϕ3aϕ1aϕ2.S_{\rm int}^{\nabla} = a_{123} \int d^{d+1}X\sqrt g\, \phi_3\nabla_a\phi_1\nabla^a\phi_2.

For separated scalar three-point functions, many derivative interactions reduce to the same conformal structure with a modified coefficient. Using the free equations of motion

2ϕi=mi2ϕi,\nabla^2\phi_i=m_i^2\phi_i,

and integrating by parts, one obtains

gϕ3aϕ1aϕ2=12(m32m12m22)gϕ1ϕ2ϕ3+boundary terms.\int\sqrt g\, \phi_3\nabla_a\phi_1\nabla^a\phi_2 = \frac{1}{2} \left(m_3^2-m_1^2-m_2^2\right) \int\sqrt g\,\phi_1\phi_2\phi_3 + \text{boundary terms}.

Thus this derivative coupling shifts the effective scalar cubic coefficient by

λ123eff=λ123+a1232(m32m12m22)+,\lambda_{123}^{\rm eff} = \lambda_{123} + \frac{a_{123}}{2} \left(m_3^2-m_1^2-m_2^2\right) + \cdots,

where the dots denote other derivative structures and possible boundary contributions.

The phrase “boundary terms” should not be skipped over. At separated, non-extremal points, many such terms are contact terms or vanish after renormalization. In extremal and near-extremal correlators, however, boundary terms can carry finite physical information. This is one reason why top-down three-point computations require the properly reduced and renormalized action, not only the bulk equations of motion.

For spinning operators the situation is richer. Derivative couplings are no longer merely coefficient shifts. They produce different conformal tensor structures. A cubic vertex such as

Aaϕ1aϕ2A^a\phi_1\nabla_a\phi_2

computes a current-scalar-scalar three-point function, while graviton couplings compute stress-tensor correlators. Ward identities then relate some three-point coefficients to two-point normalizations.

Suppose the holographic two-point functions are

Oi(x)Oi(0)=Cix2Δi,\langle \mathcal O_i(x)\mathcal O_i(0)\rangle = \frac{C_i}{|x|^{2\Delta_i}},

and the cubic computation gives

O1(x1)O2(x2)O3(x3)=C123x122δ12x132δ13x232δ23.\langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3)\rangle = \frac{\mathfrak C_{123}}{|x_{12}|^{2\delta_{12}}|x_{13}|^{2\delta_{13}}|x_{23}|^{2\delta_{23}}}.

Then the normalized OPE coefficient is

λ123CFT=C123C1C2C3.\boxed{ \lambda_{123}^{\rm CFT} = \frac{\mathfrak C_{123}}{\sqrt{C_1C_2C_3}}. }

This is often the number that should be compared with conformal-bootstrap conventions.

In a top-down compactification, the cubic coupling itself often has the schematic form

λ123(d+1)116πGDYdqygYY1(y)Y2(y)Y3(y)×(derivatives, flux factors, metric factors),\lambda_{123}^{(d+1)} \sim \frac{1}{16\pi G_D} \int_Y d^q y\sqrt{g_Y}\, \mathcal Y_1(y)\mathcal Y_2(y)\mathcal Y_3(y) \times \left(\text{derivatives, flux factors, metric factors}\right),

where D=d+1+qD=d+1+q and Yi\mathcal Y_i are internal harmonics or wavefunctions on the compact space YY. This immediately explains many selection rules. For AdS5×S5\mathrm{AdS}_5\times S^5, the internal harmonics transform under SO(6)SO(6), so a three-point function is allowed only when the tensor product of the corresponding SO(6)SO(6) representations contains a singlet.

The internal-space overlap is the bulk version of a CFT statement: global symmetry representations constrain which OPE coefficients can be nonzero.

The overall size of cubic Witten diagrams is controlled by the gravitational coupling. For a classical Einstein-like bulk dual, the effective action has the schematic form

SbulkCTdd+1Xg[12(ϕ)2+12m2ϕ2+g3ϕ3+g4ϕ4+],S_{\rm bulk} \sim C_T \int d^{d+1}X\sqrt g \left[ \frac{1}{2}(\nabla\phi)^2 + \frac{1}{2}m^2\phi^2 +g_3\phi^3 +g_4\phi^4 +\cdots \right],

where

CTLd1Gd+1C_T \sim \frac{L^{d-1}}{G_{d+1}}

is proportional to the stress-tensor two-point coefficient of the CFT. In the canonical AdS5_5/CFT4_4 example, CTN2C_T\sim N^2.

With the source-normalized fields above, connected nn-point functions of unnormalized single-trace operators scale as

O1OnconnCT\langle \mathcal O_1\cdots\mathcal O_n\rangle_{\rm conn} \sim C_T

at tree level, because the whole classical action is multiplied by CTC_T. After normalizing each operator to have a unit two-point function,

O^i=OiCT,\widehat{\mathcal O}_i=\frac{\mathcal O_i}{\sqrt{C_T}},

the connected tree-level scaling becomes

O^1O^nconntreeCT1n/2.\boxed{ \langle \widehat{\mathcal O}_1\cdots\widehat{\mathcal O}_n \rangle_{\rm conn}^{\rm tree} \sim C_T^{1-n/2}. }

Thus normalized scalar three-point coefficients scale as

λ123CFTCT1/21N\lambda_{123}^{\rm CFT} \sim C_T^{-1/2} \sim \frac{1}{N}

for matrix large-NN theories with CTN2C_T\sim N^2. This is the same large-NN factorization seen from the boundary. In canonically normalized bulk fields,

ϕcanCTϕ,\phi_{\rm can} \sim \sqrt{C_T}\,\phi,

the cubic coupling is explicitly suppressed by 1/CT1/\sqrt{C_T}. Bulk perturbation theory is therefore the same expansion as large-NN factorization.

The contact integral contains gamma functions

Γ(δ12),Γ(δ13),Γ(δ23).\Gamma(\delta_{12}), \qquad \Gamma(\delta_{13}), \qquad \Gamma(\delta_{23}).

If one dimension equals the sum of the other two, for example

Δ3=Δ1+Δ2,\Delta_3=\Delta_1+\Delta_2,

then

δ12=0,\delta_{12}=0,

and the naive bulk integral diverges. This is the extremal case. Geometrically, the divergence comes from a near-boundary region of the bulk integral. In many protected supergravity examples, the corresponding cubic coupling vanishes at the same rate as the integral diverges, so the product is finite after analytic continuation. In other cases, boundary terms are essential.

The lesson is simple but important:

Extremal three-point functions cannot be computed by blindly inserting dimensions into the non-extremal contact integral.

One must use the renormalized action, including boundary terms and the correct limiting procedure. These subtleties are especially common for protected chiral-primary correlators in supersymmetric examples.

There are also logarithmic cases. When dimensions make powers in the near-boundary expansion collide, holographic renormalization produces logarithmic counterterms. These affect contact terms and anomaly-related terms. At separated points the nonlocal conformal structure remains fixed, but the precise local terms are scheme-dependent.

The formula above is a separated-point result. Local counterterms can add terms such as

δ(d)(x1x2)1x132Δ3,x1δ(d)(x1x2)δ(d)(x1x3),\delta^{(d)}(x_1-x_2)\frac{1}{|x_{13}|^{2\Delta_3}}, \qquad \Box_{x_1}\delta^{(d)}(x_1-x_2)\delta^{(d)}(x_1-x_3),

and similar descendants. These are contact terms. They matter for Ward identities, anomalies, integrated correlators, and source-dependent one-point functions. They do not change the separated-point OPE coefficient λ123CFT\lambda_{123}^{\rm CFT}.

Bulk field redefinitions can also move terms between different cubic vertices and boundary terms. For example, redefining

ϕ3ϕ3+aϕ1ϕ2\phi_3\to \phi_3+a\phi_1\phi_2

changes the apparent cubic Lagrangian. A physical separated-point CFT three-point coefficient cannot depend on such a choice. The invariant object is the properly renormalized on-shell action as a functional of the boundary sources.

To compute a scalar three-point coefficient holographically:

  1. Identify the operators and fields. Determine mi2L2=Δi(Δid)m_i^2L^2=\Delta_i(\Delta_i-d) and the correct standard or alternate quantization.
  2. Normalize the quadratic action. Compute the two-point coefficients CiC_i in the same conventions.
  3. Find the cubic action. Include derivative terms, mixing terms, boundary terms, and compactification overlap factors.
  4. Evaluate the contact integral. Reduce derivative vertices when possible, or compute the appropriate tensor integral.
  5. Renormalize. Separate the physical separated-point coefficient from contact terms and scheme-dependent local terms.
  6. Normalize the OPE coefficient. Divide by C1C2C3\sqrt{C_1C_2C_3} if comparing to unit-normalized CFT operators.
  7. Check symmetries. Global charges, spins, parity, supersymmetry, and Ward identities should match the CFT expectations.

This recipe generalizes. Four-point contact diagrams use the same logic but are no longer fixed by conformal symmetry because four points have cross-ratios. Exchange diagrams require bulk-to-bulk propagators. Loops require the quantum bulk theory. Those are the subject of the next Witten-diagram page.

MistakeWhy it is wrongSafer statement
“The cubic coupling is the OPE coefficient.”It is only true after propagator normalizations, AdS integrals, two-point normalizations, and renormalization are accounted for.Cubic couplings determine OPE coefficients.
Forgetting two-point normalizationBootstrap OPE coefficients usually use unit-normalized operators.Always compute λ123CFT=C123/C1C2C3\lambda_{123}^{\rm CFT}=\mathfrak C_{123}/\sqrt{C_1C_2C_3}.
Ignoring derivative termsSupergravity cubic actions often contain derivatives and field mixing.Reduce derivative vertices on shell only after checking boundary terms.
Treating extremal limits as ordinary integralsGamma functions can diverge and cubic couplings may vanish simultaneously.Use analytic continuation and the full renormalized action.
Confusing contact terms with separated OPE dataCounterterms change local distributions but not separated-point scalar structures.Specify whether the result is at separated points or as a distribution.
Assuming large NN means classical Einstein gravityLarge NN suppresses loops, but locality also requires a large gap to stringy or higher-spin states.Classical supergravity requires both large CTC_T and a sparse low-dimension single-trace spectrum.

Exercise 1: Fix the scalar three-point exponents

Section titled “Exercise 1: Fix the scalar three-point exponents”

Assume the scalar three-point function has the form

G(x1,x2,x3)=Cx12ax13bx23c.G(x_1,x_2,x_3) = \frac{C}{|x_{12}|^{a}|x_{13}|^{b}|x_{23}|^{c}}.

Use scaling at each insertion to show that

a=Δ1+Δ2Δ3,b=Δ1+Δ3Δ2,c=Δ2+Δ3Δ1.a=\Delta_1+\Delta_2-\Delta_3, \qquad b=\Delta_1+\Delta_3-\Delta_2, \qquad c=\Delta_2+\Delta_3-\Delta_1.
Solution

Under a global scale transformation xiΛxix_i\to \Lambda x_i, the whole correlator must scale as

G(Λx1,Λx2,Λx3)=Λ(Δ1+Δ2+Δ3)G(x1,x2,x3),G(\Lambda x_1,\Lambda x_2,\Lambda x_3) = \Lambda^{-(\Delta_1+\Delta_2+\Delta_3)}G(x_1,x_2,x_3),

so

a+b+c=Δ1+Δ2+Δ3.a+b+c=\Delta_1+\Delta_2+\Delta_3.

A sharper way is to apply inversion or, equivalently, demand the correct conformal weight at each point. The factors containing x1x_1 are x12a|x_{12}|^{-a} and x13b|x_{13}|^{-b}, so insertion 11 has weight

a+b=2Δ1.a+b=2\Delta_1.

Similarly,

a+c=2Δ2,b+c=2Δ3.a+c=2\Delta_2, \qquad b+c=2\Delta_3.

Solving these three equations gives

a=Δ1+Δ2Δ3,b=Δ1+Δ3Δ2,c=Δ2+Δ3Δ1.a=\Delta_1+\Delta_2-\Delta_3, \qquad b=\Delta_1+\Delta_3-\Delta_2, \qquad c=\Delta_2+\Delta_3-\Delta_1.

Exercise 2: Evaluate the contact integral by Schwinger parameters

Section titled “Exercise 2: Evaluate the contact integral by Schwinger parameters”

Use

1AΔ=1Γ(Δ)0dααΔ1eαA\frac{1}{A^\Delta} = \frac{1}{\Gamma(\Delta)} \int_0^\infty d\alpha\,\alpha^{\Delta-1}e^{-\alpha A}

to derive the position dependence of

DΔ1Δ2Δ3=0dzzd+1ddxi=13(zz2+xxi2)Δi.D_{\Delta_1\Delta_2\Delta_3} = \int_0^\infty\frac{dz}{z^{d+1}} \int d^d x \prod_{i=1}^3 \left( \frac{z}{z^2+|x-x_i|^2} \right)^{\Delta_i}.

You do not need to reproduce every gamma-function coefficient, but you should obtain the powers of xij|x_{ij}|.

Solution

Write each denominator with a Schwinger parameter:

(z2+xxi2)Δi=1Γ(Δi)0dαiαiΔi1eαi(z2+xxi2).(z^2+|x-x_i|^2)^{-\Delta_i} = \frac{1}{\Gamma(\Delta_i)} \int_0^\infty d\alpha_i\, \alpha_i^{\Delta_i-1} e^{-\alpha_i(z^2+|x-x_i|^2)}.

The exponent involving xx is

iαixxi2=Axxˉ2+1Ai<jαiαjxij2,\sum_i\alpha_i|x-x_i|^2 = A|x-\bar x|^2 + \frac{1}{A}\sum_{i<j}\alpha_i\alpha_j|x_{ij}|^2,

where

A=α1+α2+α3,xˉ=1Aiαixi.A=\alpha_1+\alpha_2+\alpha_3, \qquad \bar x=\frac{1}{A}\sum_i\alpha_i x_i.

The xx integral gives a Gaussian factor Ad/2A^{-d/2}. The zz integral is also elementary after combining eAz2e^{-Az^2} with the numerator zΣd1z^{\Sigma-d-1}:

0dzzΣd1eAz2=12A(Σd)/2Γ(Σd2).\int_0^\infty dz\,z^{\Sigma-d-1}e^{-Az^2} = \frac{1}{2}A^{-(\Sigma-d)/2} \Gamma\left(\frac{\Sigma-d}{2}\right).

Now change variables to αi=Aui\alpha_i=A u_i with u1+u2+u3=1u_1+u_2+u_3=1. The remaining integral over the simplex has exponential

exp[Ai<juiujxij2].\exp\left[-A\sum_{i<j}u_iu_j|x_{ij}|^2\right].

The final AA integral and simplex integral produce the unique homogeneous conformal structure. Dimensional covariance and symmetry among the three points then fix the powers to be

DΔ1Δ2Δ31x122δ12x132δ13x232δ23,D_{\Delta_1\Delta_2\Delta_3} \propto \frac{1}{|x_{12}|^{2\delta_{12}}|x_{13}|^{2\delta_{13}}|x_{23}|^{2\delta_{23}}},

with

δ12=Δ1+Δ2Δ32,δ13=Δ1+Δ3Δ22,δ23=Δ2+Δ3Δ12.\delta_{12}=\frac{\Delta_1+\Delta_2-\Delta_3}{2}, \qquad \delta_{13}=\frac{\Delta_1+\Delta_3-\Delta_2}{2}, \qquad \delta_{23}=\frac{\Delta_2+\Delta_3-\Delta_1}{2}.

Keeping track of the beta-function integrals gives the gamma-function coefficient displayed in the main text.

Show that, for free external scalar solutions satisfying 2ϕi=mi2ϕi\nabla^2\phi_i=m_i^2\phi_i,

gϕ3aϕ1aϕ2=12(m32m12m22)gϕ1ϕ2ϕ3+boundary terms.\int\sqrt g\,\phi_3\nabla_a\phi_1\nabla^a\phi_2 = \frac{1}{2}(m_3^2-m_1^2-m_2^2) \int\sqrt g\,\phi_1\phi_2\phi_3 + \text{boundary terms}.
Solution

Use the identity

aϕ1aϕ2=12[2(ϕ1ϕ2)ϕ12ϕ2ϕ22ϕ1].\nabla_a\phi_1\nabla^a\phi_2 = \frac{1}{2} \left[ \nabla^2(\phi_1\phi_2)-\phi_1\nabla^2\phi_2-\phi_2\nabla^2\phi_1 \right].

Multiplying by ϕ3\phi_3 and integrating gives

I=12gϕ32(ϕ1ϕ2)12(m12+m22)gϕ1ϕ2ϕ3.I = \frac{1}{2} \int\sqrt g\,\phi_3\nabla^2(\phi_1\phi_2) - \frac{1}{2}(m_1^2+m_2^2) \int\sqrt g\,\phi_1\phi_2\phi_3.

Integrate the first term by parts twice:

gϕ32(ϕ1ϕ2)=g(2ϕ3)ϕ1ϕ2+boundary terms.\int\sqrt g\,\phi_3\nabla^2(\phi_1\phi_2) = \int\sqrt g\,(\nabla^2\phi_3)\phi_1\phi_2 + \text{boundary terms}.

Using 2ϕ3=m32ϕ3\nabla^2\phi_3=m_3^2\phi_3 yields

I=12(m32m12m22)gϕ1ϕ2ϕ3+boundary terms.I = \frac{1}{2}(m_3^2-m_1^2-m_2^2) \int\sqrt g\,\phi_1\phi_2\phi_3 + \text{boundary terms}.

Exercise 4: Large-NN normalization of a three-point coefficient

Section titled “Exercise 4: Large-NNN normalization of a three-point coefficient”

Assume a matrix large-NN CFT has CTN2C_T\sim N^2. A single-trace operator O\mathcal O is normalized so that

O(x)O(0)N2x2Δ,\langle \mathcal O(x)\mathcal O(0)\rangle\sim \frac{N^2}{|x|^{2\Delta}},

and the connected three-point function scales as

OOOconnN2.\langle \mathcal O\mathcal O\mathcal O\rangle_{\rm conn}\sim N^2.

Define O^=O/N\widehat{\mathcal O}=\mathcal O/N. How does

O^O^O^conn\langle \widehat{\mathcal O}\widehat{\mathcal O}\widehat{\mathcal O}\rangle_{\rm conn}

scale with NN?

Solution

Each normalized operator contributes a factor 1/N1/N. Therefore

O^O^O^conn=1N3OOOconnN2N3=1N.\langle \widehat{\mathcal O}\widehat{\mathcal O}\widehat{\mathcal O}\rangle_{\rm conn} = \frac{1}{N^3} \langle \mathcal O\mathcal O\mathcal O\rangle_{\rm conn} \sim \frac{N^2}{N^3} = \frac{1}{N}.

This agrees with the bulk statement that canonically normalized cubic couplings are suppressed by 1/N1/N in theories with CTN2C_T\sim N^2.

Let Δ3=Δ1+Δ2\Delta_3=\Delta_1+\Delta_2. Which gamma function in the contact integral diverges? What is the expected position dependence of the separated three-point function after a proper limiting procedure?

Solution

If

Δ3=Δ1+Δ2,\Delta_3=\Delta_1+\Delta_2,

then

δ12=Δ1+Δ2Δ32=0.\delta_{12} = \frac{\Delta_1+\Delta_2-\Delta_3}{2} =0.

The contact integral contains Γ(δ12)=Γ(0)\Gamma(\delta_{12})=\Gamma(0), so the naive non-extremal formula diverges. This is the extremal divergence.

The separated-point conformal structure is still fixed by conformal symmetry. Substituting Δ3=Δ1+Δ2\Delta_3=\Delta_1+\Delta_2 gives

O1(x1)O2(x2)O3(x3)1x132Δ1x232Δ2,\langle \mathcal O_1(x_1)\mathcal O_2(x_2)\mathcal O_3(x_3)\rangle \propto \frac{1}{ |x_{13}|^{2\Delta_1} |x_{23}|^{2\Delta_2} },

with no power of x12|x_{12}|. The coefficient must be computed by the correct renormalized limiting procedure, including possible vanishing cubic couplings and boundary terms.