Celestial CFT, Virasoro Algebra, and BMS
1. Celestial Amplitudes via the Mellin Transform
Standard scattering amplitudes $\mathcal{M}(p_1,\ldots,p_n)$ are functions of four-momenta. In the celestial basis we trade the energy variable $\omega$ for a conformal dimension $\Delta$ via the Mellin transform. For each external particle $k$:
$$\tilde{\mathcal{O}}_{\Delta_k}(z_k,\bar{z}_k) = \int_0^\infty d\omega_k\;\omega_k^{\Delta_k - 1}\;\mathcal{O}_{\omega_k}(z_k,\bar{z}_k)$$
The celestial amplitude is the Mellin transform of the momentum-space amplitude:
$$\widetilde{\mathcal{M}}(\Delta_1,\ldots,\Delta_n;\{z_k,\bar{z}_k\}) = \prod_{k=1}^n \int_0^\infty d\omega_k\;\omega_k^{\Delta_k - 1}\;\mathcal{M}(\{p_k(\omega_k,z_k,\bar{z}_k)\})$$
The key property is that Lorentz transformations act on $(z,\bar{z})$ as$SL(2,\mathbb{C})$ Mobius transformations, so the celestial amplitude transforms as a correlator of conformal primaries on $S^2$.
2. Massless Momentum Parametrization
A null momentum $p^2 = 0$ can be written using the spinor-helicity parametrization. Introduce stereographic coordinates $(z,\bar{z})$ on the celestial sphere. A null ray with energy $\omega$ pointing in the direction labelled by $(z,\bar{z})$ has:
$$p^\mu(\omega,z,\bar{z}) = \frac{\omega}{2}\left(1 + z\bar{z},\; z + \bar{z},\; -i(z - \bar{z}),\; 1 - z\bar{z}\right)$$
Verify nullity: $p^0 = \frac{\omega}{2}(1+z\bar{z})$, $|\vec{p}|^2 = \frac{\omega^2}{4}[(z+\bar{z})^2 + (z-\bar{z})^2 + (1-z\bar{z})^2]$. Expanding:
$$|\vec{p}|^2 = \frac{\omega^2}{4}\left[4z\bar{z} + 1 - 2z\bar{z} + z^2\bar{z}^2\right] = \frac{\omega^2}{4}(1 + z\bar{z})^2 = (p^0)^2$$
Under an $SL(2,\mathbb{C})$ Lorentz transformation $z \mapsto (az+b)/(cz+d)$, the energy scales as $\omega \mapsto |cz+d|^2\,\omega$. A celestial primary of dimension$\Delta$ therefore transforms as:
$$\tilde{\mathcal{O}}_\Delta(z,\bar{z}) \mapsto |cz+d|^{2\Delta}\;\tilde{\mathcal{O}}_\Delta\!\left(\frac{az+b}{cz+d},\frac{\bar{a}\bar{z}+\bar{b}}{\bar{c}\bar{z}+\bar{d}}\right)$$
3. Conformal Primary Weights
For a massless particle of helicity $s$ (graviton: $s = \pm 2$), the conformal weights $(h,\bar{h})$ are derived from the Lorentz transformation. Under a rotation $z \mapsto e^{i\theta}z$, the operator picks up a phase $e^{-is\theta}$from helicity and $e^{-i\Delta\theta}$ from the Mellin weight. Matching:
$$h = \frac{\Delta + s}{2}, \qquad \bar{h} = \frac{\Delta - s}{2}$$
For a positive-helicity graviton ($s = +2$):$h = \frac{\Delta + 2}{2}$, $\bar{h} = \frac{\Delta - 2}{2}$. The spin is $h - \bar{h} = s = 2$ and the scaling dimension is $h + \bar{h} = \Delta$, as required. For a negative-helicity graviton ($s = -2$):$h = \frac{\Delta - 2}{2}$, $\bar{h} = \frac{\Delta + 2}{2}$.
4. Graviton OPE and the Beta Function Coefficient
The operator product expansion of two positive-helicity celestial gravitons is derived from the collinear limit of the 3-point amplitude. When $z_1 \to z_2$:
$$\mathcal{O}^+_{h_1}(z_1)\,\mathcal{O}^+_{h_2}(z_2) \sim \frac{\kappa}{2}\,\frac{B(h_1-1,\,h_2-1)}{z_{12}}\;\mathcal{O}^+_{h_1+h_2-1}(z_2) + \cdots$$
The Euler beta function arises from the Mellin convolution. Starting from the 3-graviton amplitude $\mathcal{M}_3 \sim \kappa\,\omega_1\omega_2/\omega_3$ in the collinear limit, the Mellin transform yields:
$$\int_0^\infty d\omega_1\;\omega_1^{h_1-1}\int_0^\infty d\omega_2\;\omega_2^{h_2-1}\;\frac{\omega_1\omega_2}{\omega_1+\omega_2}\;\delta(\omega_3 - \omega_1 - \omega_2)$$
Substituting $\omega_1 = x\omega_3$, $\omega_2 = (1-x)\omega_3$:
$$= \omega_3^{h_1+h_2-2}\int_0^1 dx\;x^{h_1-1}(1-x)^{h_2-1}\;x(1-x) = \omega_3^{h_1+h_2-2}\,B(h_1,h_2)$$
The $1/z_{12}$ singularity comes from the holomorphic collinear splitting function of gravity, confirming the chiral algebra structure.
5. Virasoro Generators and Central Extension
The BMS superrotation generators are promoted from global conformal transformations to local ones. The stress tensor mode expansion gives:
$$T(z) = \sum_{n \in \mathbb{Z}} L_n\,z^{-n-2}, \qquad L_n = \oint \frac{dz}{2\pi i}\;z^{n+1}\,T(z)$$
To compute $[L_m, L_n]$, use the OPE $T(z)T(w) \sim \frac{c/2}{(z-w)^4} + \frac{2T(w)}{(z-w)^2} + \frac{\partial T(w)}{z-w}$. The double contour integral gives:
$$[L_m, L_n] = \oint \frac{dw}{2\pi i}\,w^{n+1}\oint \frac{dz}{2\pi i}\,z^{m+1}\left[\frac{c/2}{(z-w)^4} + \frac{2T(w)}{(z-w)^2} + \frac{\partial T(w)}{z-w}\right]$$
Evaluating each residue: the $(z-w)^{-4}$ pole gives $\frac{c}{12}m(m^2-1)\delta_{m+n,0}$, the $(z-w)^{-2}$ pole gives $2(m+1)L_{m+n}$, and the $(z-w)^{-1}$pole gives $-(m+n+2)L_{m+n}$. Combining:
$$[L_m, L_n] = (m - n)\,L_{m+n} + \frac{c}{12}(m^3 - m)\,\delta_{m+n,0}$$
6. Ward Identity as Spin Memory
The stress tensor Ward identity for an $n$-point correlator of primaries $\mathcal{O}_k$ with weights $h_k$ at positions $z_k$ is:
$$\langle T(z)\,\mathcal{O}_1(z_1)\cdots\mathcal{O}_n(z_n)\rangle = \sum_{k=1}^n \left[\frac{h_k}{(z-z_k)^2} + \frac{\partial_{z_k}}{z-z_k}\right]\langle\mathcal{O}_1\cdots\mathcal{O}_n\rangle$$
In BMS language, $T(z)$ is the superrotation charge density. Its insertion measures the angular momentum flux. The $L_{-1}$ Ward identity specifically encodes infinitesimal translations on the celestial sphere:
$$\langle L_{-1}\,\mathcal{O}_1\cdots\mathcal{O}_n\rangle = \sum_{k=1}^n \partial_{z_k}\langle\mathcal{O}_1\cdots\mathcal{O}_n\rangle$$
The spin memory effect is the integrated version: the total superrotation charge flux between $\mathscr{I}^+_-$ and $\mathscr{I}^+_+$ equals the sum of angular momentum shifts at each insertion point. This is exactly the celestial CFT Ward identity applied to the S-matrix as a celestial correlator.
7. Perelman via Polyakov: Effective Central Charge
The worldsheet stress tensor of the nonlinear sigma model with target metric $G_{\mu\nu}$ is:
$$T_{zz} = -\frac{1}{\alpha'}\,G_{\mu\nu}(X)\,\partial_z X^\mu\,\partial_z X^\nu$$
The beta function for the target metric is $\beta^G_{\mu\nu} = \alpha' R_{\mu\nu} + O(\alpha'^2)$. Conformal invariance ($\beta = 0$) gives the Ricci-flat condition. The Zamolodchikov c-function is constructed from the two-point function $\langle T(z)T(0)\rangle = c/(2z^4)$:
$$c_{\rm eff}(t) = c_{\rm UV} - 12\pi\int_0^t \beta^i\,G_{ij}\,\beta^j\,dt'$$
Since $\beta^i G_{ij}\beta^j \geq 0$ (the Zamolodchikov metric is positive definite),$c_{\rm eff}$ is monotonically non-increasing. Now identify $\beta^G_{\mu\nu} = \alpha' R_{\mu\nu}$ with the Ricci flow$\partial_t g_{\mu\nu} = -2R_{\mu\nu}$ (setting $\alpha' = 2$). Then:
$$\frac{dc_{\rm eff}}{dt} = -12\pi\,R_{\mu\nu}\,G^{\mu\alpha}G^{\nu\beta}\,R_{\alpha\beta} = -12\pi\,|Ric|^2 \leq 0$$
This is Perelman's monotonicity in disguise: the $\mathcal{W}$-entropy is the natural c-function for the gravitational RG flow, and its monotonicity is the c-theorem.
Celestial CFT: Amplitudes, OPE, Virasoro, and c-Theorem
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server