Skip to content

Liouville Theory and the Virasoro Bootstrap

Liouville theory is the canonical exactly solvable CFT with a continuous Virasoro spectrum. Its bootstrap data combine a reflection identification, delta-normalized states, special-function three-point constants, and Virasoro blocks integrated along a prescribed contour. The formulas below describe spacelike Liouville theory with real b>0b>0 and generic external momenta; analytic continuation outside that domain requires an explicit pole prescription.

Evidence cutoff. This account reflects the sources consulted through 2026-08-09. Statements about analytic continuation and crossing are conditional on the spectrum, meromorphy, contour, and normalization hypotheses stated below; they are not claims about every theory called “Liouville” or every continuation to c1c\leq1.

Required background. Noncompact CFTs and continuous spectra provide direct-integral normalization and contour language. Highest-weight modules, null states, and the Kac determinant provide the degenerate Virasoro representations used to obtain shift equations.

Helpful background. Partial waves and the shadow formalism give a higher-dimensional comparison for continuum decompositions and normalization choices.

On a Euclidean surface with reference metric gg, take

SL[ϕ]=14πd2xg[gabaϕbϕ+QRgϕ+4πμe2bϕ],S_L[\phi] =\frac{1}{4\pi}\int d^2x\sqrt g\, \left[g^{ab}\partial_a\phi\partial_b\phi +Q R_g\phi+4\pi\mu e^{2b\phi}\right],

with

Q=b+b1,cL=1+6Q2,μ>0.Q=b+b^{-1}, \qquad c_L=1+6Q^2, \qquad \mu>0.

The exponential primary

Vα(z,zˉ)=e2αϕ(z,zˉ)V_\alpha(z,\bar z)=e^{2\alpha\phi(z,\bar z)}

has weights

hα=hˉα=α(Qα).h_\alpha=\bar h_\alpha=\alpha(Q-\alpha).

Normalizable scattering states lie on

α=Q2+iP,PR,hP=Q24+P2.\alpha=\frac Q2+iP, \qquad P\in\mathbb R, \qquad h_P=\frac{Q^2}{4}+P^2.

The reflection relation

Vα=R(α)VQα,R(α)R(Qα)=1,V_\alpha=R(\alpha)V_{Q-\alpha}, \qquad R(\alpha)R(Q-\alpha)=1,

identifies PP and P-P. We therefore choose a reflected basis V^P\widehat V_P with P0P\geq0 and declare

PP=2πδ(PP),1=0dP2πPP.\langle P|P'\rangle=2\pi\delta(P-P'), \qquad \mathbf1=\int_0^\infty\frac{dP}{2\pi}|P\rangle\langle P|.

The hat absorbs the reflection amplitude and a normalization factor. Raw VαV_\alpha correlators instead contain both δ(PP)\delta(P-P') and the reflected δ(P+P)\delta(P+P'); mixing these two conventions double-counts the continuum. The spectrum, two-point normalization, and degenerate analytic continuations are developed in Ribault 2018, §§3.1.1–3.1.2, pp. 69–73.

On the sphere, write ϕ=ϕ0+ϕ~\phi=\phi_0+\widetilde\phi and let

A[ϕ~]=d2xge2bϕ~.\mathcal A[\widetilde\phi] =\int d^2x\sqrt g\,e^{2b\widetilde\phi}.

The constant-mode part of an nn-point path integral is proportional to

dϕ0exp ⁣[2(iαiQ)ϕ0μAe2bϕ0].\int_{-\infty}^{\infty}d\phi_0\, \exp\!\left[ 2\left(\sum_i\alpha_i-Q\right)\phi_0 -\mu\mathcal A e^{2b\phi_0} \right].

For Re(iαiQ)>0\operatorname{Re}(\sum_i\alpha_i-Q)>0, the substitution y=μAe2bϕ0y=\mu\mathcal A e^{2b\phi_0} gives

12b(μA)sΓ(s),s=iαiQb.\frac{1}{2b} (\mu\mathcal A)^{-s}\Gamma(s), \qquad s=\frac{\sum_i\alpha_i-Q}{b}.

Outside this convergence half-plane, the expression is defined by meromorphic continuation rather than by the original absolutely convergent zero-mode integral. Its gamma-function poles foreshadow the discrete residues that appear when an OPE contour is continued. This is an interacting analogue of the free-boson zero-mode issue, but it does not produce an ordinary momentum-conservation delta function because the exponential potential breaks translations of ϕ\phi.

Let

γ(x)=Γ(x)Γ(1x),A=πμγ(b2)b22b2.\gamma(x)=\frac{\Gamma(x)}{\Gamma(1-x)}, \qquad A=\pi\mu\gamma(b^2)b^{2-2b^2}.

In the standard exponential-field normalization, the Dorn–Otto–Zamolodchikov–Zamolodchikov three-point constant is

C(α1,α2,α3)=A(Qα1α2α3)/b×Υb(0)i=13Υb(2αi)Υb(α1+α2+α3Q)Υb(α1+α2α3)×1Υb(α1+α3α2)Υb(α2+α3α1).\begin{aligned} C(\alpha_1,\alpha_2,\alpha_3) ={}&A^{(Q-\alpha_1-\alpha_2-\alpha_3)/b}\\ &\times \frac{\Upsilon_b'(0) \prod_{i=1}^3\Upsilon_b(2\alpha_i)} {\Upsilon_b(\alpha_1+\alpha_2+\alpha_3-Q) \Upsilon_b(\alpha_1+\alpha_2-\alpha_3)}\\ &\times \frac{1} {\Upsilon_b(\alpha_1+\alpha_3-\alpha_2) \Upsilon_b(\alpha_2+\alpha_3-\alpha_1)}. \end{aligned}

The entire function Υb\Upsilon_b is fixed, up to normalization, by its zero set, the duality bb1b\leftrightarrow b^{-1}, and the shift equations

Υb(x+b)=γ(bx)b12bxΥb(x),\Upsilon_b(x+b) =\gamma(bx)b^{1-2bx}\Upsilon_b(x), Υb(x+b1)=γ(x/b)b2x/b1Υb(x).\Upsilon_b(x+b^{-1}) =\gamma(x/b)b^{2x/b-1}\Upsilon_b(x).

The original derivation, reflection amplitude, and four-point construction appear in Zamolodchikov and Zamolodchikov 1996, §§2–5, pp. 579–592. The formula is meromorphic in the αi\alpha_i; its zeros and poles are part of the contour data, not removable numerical nuisances.

The degenerate field Vb/2V_{-b/2} is not a normalizable state on the principal-series contour; it is introduced as a meromorphic continuation of correlation functions. It has a level-two null descendant,

(L12+b2L2)Vb/2=0,\left(L_{-1}^2+b^2L_{-2}\right)V_{-b/2}=0,

and its OPE with a generic exponential contains only two Virasoro families:

Vb/2VαC+(α)Vαb/2+C(α)Vα+b/2.V_{-b/2}V_\alpha \sim C_+(\alpha)V_{\alpha-b/2} +C_-(\alpha)V_{\alpha+b/2}.

The null relation turns a four-point Ward identity into a second-order BPZ equation. Comparing its two finite channel decompositions relates the coefficients with intermediate momenta shifted by ±b/2\pm b/2. After eliminating the degenerate OPE normalization, one obtains a bb-shift equation for the generic three-point constant; using the dual field V1/(2b)V_{-1/(2b)} gives the b1b^{-1} shift. This derivation and the normalization-independent form of the shift equations are given in Ribault 2018, §3.1.2, pp. 71–75.

The DOZZ formula passes this test directly. Define

x0=α1+α2+α3Q,x1=α1+α2α3,x2=α1+α3α2,x3=α2+α3α1.\begin{gathered} x_0=\alpha_1+\alpha_2+\alpha_3-Q, \qquad x_1=\alpha_1+\alpha_2-\alpha_3,\\ x_2=\alpha_1+\alpha_3-\alpha_2, \qquad x_3=\alpha_2+\alpha_3-\alpha_1. \end{gathered}

Taking the ratio at α1+b\alpha_1+b and α1\alpha_1 gives

C(α1+b,α2,α3)C(α1,α2,α3)=A1Υb(2α1+2b)Υb(2α1)×Υb(x0)Υb(x1)Υb(x0+b)Υb(x1+b)×Υb(x2)Υb(x3)Υb(x2+b)Υb(x3b).\begin{aligned} \frac{C(\alpha_1+b,\alpha_2,\alpha_3)} {C(\alpha_1,\alpha_2,\alpha_3)} ={}&A^{-1} \frac{\Upsilon_b(2\alpha_1+2b)} {\Upsilon_b(2\alpha_1)}\\ &\times \frac{\Upsilon_b(x_0)\Upsilon_b(x_1)} {\Upsilon_b(x_0+b)\Upsilon_b(x_1+b)}\\ &\times \frac{\Upsilon_b(x_2)\Upsilon_b(x_3)} {\Upsilon_b(x_2+b)\Upsilon_b(x_3-b)}. \end{aligned}

Applying the Υb\Upsilon_b shift identity converts this expression into the gamma-function BPZ shift relation. The dual identity supplies the independent b1b^{-1} shift. The two shifts, reflection symmetry, meromorphy, and growth assumptions remove periodic ambiguities for generic irrational b2b^2; omitting those assumptions overstates what the finite BPZ equation alone proves.

Place four operators at 0,z,1,0,z,1,\infty. In the reflected P0P\geq0 normalization declared above, an ss-channel decomposition has the form

G(z,zˉ)=0dP2πC12(P)CP34FP(s)(z)2.\mathcal G(z,\bar z) =\int_0^\infty\frac{dP}{2\pi}\, \mathcal C_{12}(P)\mathcal C_{P34} \left|\mathcal F_P^{(s)}(z)\right|^2.

The coefficients C\mathcal C are the DOZZ constants multiplied by the conversion factors from raw VαV_\alpha fields to V^P\widehat V_P states. In the original full-line convention, the equivalent formula is one half of an integral over PRP\in\mathbb R. The choice must be made once and used in the two-point function, completeness relation, and four-point measure.

For real external αi\alpha_i in the strip where no DOZZ pole crosses α=Q/2+iR\alpha=Q/2+i\mathbb R, the contour is literal. When the αi\alpha_i are analytically continued, moving poles can cross it. The continued correlator is then

G=Ccontinueddμ(P)ρ(P)FP2+2πiPResP=P[dμ(P)ρ(P)FP2],\mathcal G =\int_{\mathcal C_{\mathrm{continued}}}d\mu(P)\, \rho(P)|\mathcal F_P|^2 +2\pi i\sum_{P_*} \operatorname*{Res}_{P=P_*} \bigl[d\mu(P)\rho(P)|\mathcal F_P|^2\bigr],

with signs fixed by the contour orientation. Zamolodchikov and Zamolodchikov 1996, §2, pp. 581–583 state the initial domain and the need for additional discrete terms after pole crossings. Crossing means that the resulting analytically continued function agrees with the corresponding tt-channel construction, including the same residues—not that the unmodified real-PP integral works for every external momentum.

  • Reflection: verify that integrating over P0P\geq0 gives the same result as the half-weighted full-line integral. A factor-of-two discrepancy signals double counting.
  • Zero mode: reproduce the power of μ\mu and the gamma-function poles from the constant-mode integral.
  • Degenerate limit: let one external momentum approach b/2-b/2; the continuum contour must pinch or reduce to the two allowed BPZ channels with the correct residues.
  • Duality: repeat the shift test with b1b^{-1}. A proposed three-point constant that satisfies only one shift equation can retain an undetected periodic factor.
  • Pole crossing: vary an external αi\alpha_i continuously and record each pole crossing before comparing channels.
  • Block normalization: check the leading OPE behavior FP(s)(z)zhPh1h2(1+O(z))\mathcal F_P^{(s)}(z)\sim z^{h_P-h_1-h_2}(1+O(z)) in the same convention used by the structure constants.
  • Numerics: if the qq-series for a Virasoro block is truncated, report the radial domain, truncation order, precision, and tail estimate. Numerical agreement supports a specified integral; it does not remove analytic-continuation hypotheses.

Spacelike b>0b>0 gives cL25c_L\geq25. Timelike or c1c\leq1 continuations involve different contours, spectra, and sometimes different structure constants; the formulas on this page cannot be transferred to them merely by substituting an imaginary bb.

The channel-changing transform of the Virasoro blocks themselves is the subject of Irrational CFT, fusion kernels, and crossing.

Show that hαh_\alpha is real and bounded below on α=Q/2+iP\alpha=Q/2+iP with real PP.

Solution hα=(Q2+iP)(Q2iP)=Q24+P2.h_\alpha =\left(\frac Q2+iP\right) \left(\frac Q2-iP\right) =\frac{Q^2}{4}+P^2.

It is real and has minimum Q2/4Q^2/4 at P=0P=0. The states remain delta-normalized because PP is continuous.

For Res>0\operatorname{Re}s>0, evaluate

I=dϕ0e2bsϕ0μAe2bϕ0.I=\int_{-\infty}^{\infty}d\phi_0\, e^{2bs\phi_0-\mu\mathcal A e^{2b\phi_0}}.
Solution

Set y=μAe2bϕ0y=\mu\mathcal A e^{2b\phi_0}, so dϕ0=dy/(2by)d\phi_0=dy/(2by) and e2bsϕ0=(y/μA)se^{2bs\phi_0}=(y/\mu\mathcal A)^s. Then

I=(μA)s2b0dyys1ey=(μA)s2bΓ(s).I=\frac{(\mu\mathcal A)^{-s}}{2b} \int_0^\infty dy\,y^{s-1}e^{-y} =\frac{(\mu\mathcal A)^{-s}}{2b}\Gamma(s).

The poles after meromorphic continuation occur at s=0,1,2,s=0,-1,-2,\ldots.

  • Ribault, Sylvain. “Conformal Field Theory on the Plane.” SciPost Physics Lecture Notes 1 (2018). DOI.
  • Zamolodchikov, Alexander B., and Al. B. Zamolodchikov. “Structure Constants and Conformal Bootstrap in Liouville Field Theory.” Nuclear Physics B 477 (1996): 577–605. DOI.