---
title: Celestial Amplitudes in Conformal Scattering
url: https://www.emergentmind.com/topics/celestial-amplitudes
type: topic
---

# Celestial Amplitudes in Conformal Scattering

Celestial amplitudes are scattering amplitudes written in a conformal basis on the celestial sphere, obtained by Mellin-transforming external on-shell energies. In four-dimensional massless kinematics, this recasts the \(S\)-matrix into correlation functions of conformal primary operators inserted at points \((z,\bar z)\) on the two-dimensional celestial sphere, with boost weights \(\Delta_i\) replacing energies and with the Lorentz group acting as the global conformal group \(SL(2,\mathbb C)\) [2108.04801]. Within this reformulation, infrared data, asymptotic symmetries, collinear singularities, conformal block expansions, massive states, and loop effects are reorganized in terms that are natural for a two-dimensional conformal description [1904.10940][2404.01920][2307.03551].

## 1. Kinematics, Mellin transform, and conformal covariance

The basic massless kinematic parametrization writes each null momentum as
\[
p_i^\mu=\omega_i\,q^\mu(z_i,\bar z_i),\qquad 
q^\mu(z,\bar z)=(1+z\bar z,\;z+\bar z,\;-i(z-\bar z),\;1-z\bar z),
\]
with \(\omega_i>0\) and \((z_i,\bar z_i)\) stereographic coordinates on the celestial sphere [2108.04801]. The proper orthochronous Lorentz group \(SO^+(1,3)\simeq SL(2,\mathbb C)/\mathbb Z_2\) acts on \((z,\bar z)\) by Möbius transformations,
\[
z\to \frac{az+b}{cz+d},\qquad \bar z\to \frac{\bar a\bar z+\bar b}{\bar c\bar z+\bar d},
\]
so the celestial sphere carries the standard global conformal action of \(SL(2,\mathbb C)\) [2108.04801].

The celestial transform Mellin-transforms each external energy,
\[
\widetilde{\mathcal A}_n(\{\Delta_i,J_i,z_i,\bar z_i\})
=
\prod_{i=1}^n\int_0^\infty d\omega_i\,\omega_i^{\Delta_i-1}
\mathcal A_n(\{\omega_i,z_i,\bar z_i,J_i\}),
\]
where \(\Delta_i\) is the boost weight and \(J_i\) is the four-dimensional helicity [2108.04801]. In the radiative sector one uses the principal series \(\Delta_i\in 1+i\mathbb R\), with analytic continuation allowed. The associated two-dimensional weights are
\[
h_i=\frac{\Delta_i+J_i}{2},\qquad \bar h_i=\frac{\Delta_i-J_i}{2}.
\]
The conformal basis therefore diagonalizes boosts rather than energies, and global Lorentz generators act as the standard \(SL(2,\mathbb C)\) generators on celestial correlators [2108.04801].

Poincaré symmetry is more subtle in this basis. Lorentz transformations are manifest, but translations act by dimension-shifting operators. In particular, momentum operators raise \(\Delta\) by one unit on the leg they act upon, so momentum conservation becomes a set of functional constraints rather than a purely algebraic \(\delta^{(4)}(\sum p_i)\) statement [1812.01080]. For low multiplicity, these constraints force singular support. In four-point massless kinematics, momentum conservation enforces reality of the cross-ratio, so celestial correlators are supported on \(z=\bar z\) together with channel-dependent ranges such as \(z>1\) in the \(12\to 34\) channel [1904.10940].

## 2. Conformal primary wavefunctions and the celestial dictionary

The external conformal basis is built from conformal primary wavefunctions \(\Phi_{\Delta,J}^s(X;z,\bar z)\) that transform simultaneously as four-dimensional spin-\(s\) fields and two-dimensional conformal primaries. Their defining covariance is
\[
\Phi^s_{\Delta,J}(\Lambda X;z',\bar z')
=
(cz+d)^{\Delta+J}(\bar c\bar z+\bar d)^{\Delta-J}
D_s(\Lambda)\,\Phi^s_{\Delta,J}(X;z,\bar z),
\]
with \(D_s(\Lambda)\) the spin-\(s\) Lorentz representation [2108.04801].

A convenient construction uses a null tetrad built from the null direction \(q^\mu\), the spacetime point \(X^\mu\), and polarization vectors \(\epsilon_\pm^\mu\). The basic massless scalar primary is
\[
\phi^\Delta=(-q\!\cdot\! X)^{-\Delta},
\]
while gauge and graviton primaries are obtained by dressing \(\phi^\Delta\) with tetrad vectors:
\[
A_{\Delta,+1}=m\,\phi^\Delta,\qquad A_{\Delta,-1}=\bar m\,\phi^\Delta,
\]
\[
h_{\Delta,+2}=m\,m\,\phi^\Delta,\qquad h_{\Delta,-2}=\bar m\,\bar m\,\phi^\Delta.
\]
Incoming and outgoing solutions are distinguished by the \(i\epsilon\) prescription \(X_\pm^\mu=X^\mu\pm i\{1,0,0,0\}\) [2108.04801].

The shadow transform is intrinsic to the celestial basis. For a primary \(\mathcal O_{\Delta,J}\), its shadow \(\mathcal O_{2-\Delta,-J}\) is obtained by an integral transform on the sphere with kernel
\[
\mathcal O_{2-\Delta,-J}(z,\bar z)
=
\frac{k_{\Delta,J}}{2\pi}
\int d^2w'\,
\frac{\mathcal O_{\Delta,J}(w',\bar w')}
{(z-w')^{2-\Delta-J}(\bar z-\bar w')^{2-\Delta+J}},
\]
with \(k_{\Delta,J}=\Delta-1+|J|\) chosen so that the transform squares to \((-1)^{2J}\) [2108.04801]. This shadow structure enters both the operator dictionary and the treatment of low-point amplitudes.

Celestial operators are defined by pairing bulk fields with conformal primary wavefunctions through an inner product,
\[
O^{s,\pm}_{\Delta,J}(z,\bar z)
\equiv
i\big(\widehat O^s(X),\Phi^s_{\Delta^*,-J}(X_\mp;z,\bar z)\big)_\Sigma,
\]
so incoming and outgoing states become operator insertions labeled by \((\Delta,J;z,\bar z)\) [2108.04801]. This realizes the basic holographic dictionary: four-dimensional external states map to two-dimensional operator insertions; helicity maps to celestial spin via \(h-\bar h=J\); and energies Mellin-transform to boost weights \(\Delta\) [2108.04801].

Massive external states require a different conformal basis. Their momenta are parametrized on the unit hyperboloid by
\[
p^\mu=\frac{m}{2y}\big(1+y^2+z\bar z,\;z+\bar z,\;i(-z+\bar z),\;1-y^2-z\bar z\big),
\]
and the corresponding conformal basis is constructed by integrating over \((y,z,\bar z)\) rather than Mellin-transforming a single energy variable [2108.04801]. Later work computed scalar three-point celestial amplitudes with two massive scalars in terms of a hypergeometric function and with three massive scalars as a triple Mellin-Barnes integral, extending the dictionary beyond the massless sector [2404.01920].

## 3. Soft theorems, asymptotic symmetries, and celestial currents

A defining structural feature of celestial amplitudes is that infrared theorems become operator relations. In Bondi coordinates \((u,r,z,\bar z)\), asymptotically flat gravity has free data \(\{m_B,N_z,C_{zz}\}\) obeying the standard Bondi constraints, and residual diffeomorphisms preserving the falloffs are supertranslations \(f(z,\bar z)\) and superrotations \(Y^z(z)\) [2108.04801]. The associated charges satisfy a diagonal Ward identity across \(\mathcal I^\pm\),
\[
\langle\text{out}\,|\,Q^+ S-SQ^-\,|\text{in}\rangle=0,
\]
with the charges splitting into soft and hard parts \(Q^\pm=Q_S^\pm+Q_H^\pm\) [2108.04801].

In momentum space, the leading soft photon and graviton theorems read
\[
A_{n+1}^{YM}(q;\{p_i\})
\to
\Big[\sum_{i=1}^n \frac{Q_i\,\epsilon\!\cdot\! p_i}{q\!\cdot\! p_i}\Big]A_n(\{p_i\}),
\]
\[
M_{n+1}^{grav}(q;\{p_i\})
\to
\Big[\sum_{i=1}^n \kappa\,\frac{\epsilon_{\mu\nu}p_i^\mu p_i^\nu}{q\!\cdot\! p_i}\Big]M_n(\{p_i\}),
\]
and, after Mellin transformation, these soft insertions become two-dimensional currents [2108.04801]. Leading soft photons generate a Kac–Moody current \(J^a(z)\), while leading soft gravitons generate a stress tensor \(T(z)\). In gravity the stress tensor can be written as
\[
T_{zz}
=
-i\,\frac{3!}{8\pi G}
\int d^2w\,\frac{1}{(z-w)^4}\int du\,u\,\partial_u C^w_{\ \bar z},
\]
and its Ward identity has the standard conformal form
\[
\langle T(z)\,\mathcal O_k(z_k,\bar z_k)\rangle
\sim
\Big[\frac{h_k}{(z-z_k)^2}+\frac{\partial_{z_k}}{z-z_k}\Big]
\langle \mathcal O_k\rangle
\]
[2108.04801].

For Yang–Mills, the conformally soft limit \(\Delta_s\to 1\) of a positive-helicity gluon produces a holomorphic current \(j^a(z)\equiv \mathcal O^a_{0+}(z,\bar z)\), and field-theoretic analysis shows that this limit is dominated by the soft-energy region \(\omega_s\to 0\) [1903.01676]. The corresponding current Ward identity is
\[
j^a(z)\,\prod_{i=1}^n \mathcal O^{b_i}_{\lambda_i J_i}(z_i,\bar z_i)
=
\sum_{i=1}^n
\frac{f^{ab_i}{}_{c}}{z-z_i}\,
\mathcal O^{b_1}_{\lambda_1J_1}\cdots
\mathcal O^{c}_{\lambda_iJ_i}\cdots
\mathcal O^{b_n}_{\lambda_nJ_n},
\]
and the current–primary OPE is
\[
J^a(z)\,\mathcal O^{\,b}_{\Delta,h}(w,\bar w)
\sim
\frac{f^{ab}{}_{c}\,\mathcal O^{\,c}_{\Delta,h}(w,\bar w)}{z-w}
\]
[1903.01676].

The same logic extends to conformally soft gravitons. In Lorentzian three-point celestial amplitudes, a soft graviton at \(\Delta=1\) gives a primary of weights \((3/2,-1/2)\), and its level-one descendant is identified with the supertranslation current \(P_z\) [2212.07025]. This sharpens the link between soft graviton theorems, BMS supertranslations, and celestial current algebra [2212.07025].

## 4. Collinear limits, OPEs, and conformal block structure

A central operational statement is that collinear limits of four-dimensional amplitudes become operator product expansions on the celestial sphere. For gravitons, taking the holomorphic collinear limit \(z_{12}\to 0\) with \(\bar z_{12}\) fixed and Mellin-transforming the two energies produces Euler Beta functions. The resulting celestial OPEs include
\[
\mathcal O_{\Delta_1,+2}(z_1,\bar z_1)\,
\mathcal O_{\Delta_2,+2}(z_2,\bar z_2)
\sim
-\frac{\kappa}{2}\,\frac{\bar z_{12}}{z_{12}}\,
B(\Delta_1-1,\Delta_2-1)\,
\mathcal O_{\Delta_1+\Delta_2,+2}(z_2,\bar z_2),
\]
and
\[
\mathcal O_{\Delta_1,+2}(z_1,\bar z_1)\,
\mathcal O_{\Delta_2,-2}(z_2,\bar z_2)
\sim
-\frac{\kappa}{2}\,\frac{\bar z_{12}}{z_{12}}\,
B(\Delta_1-1,\Delta_2+3)\,
\mathcal O_{\Delta_1+\Delta_2,-2}(z_2,\bar z_2)
\]
[2108.04801]. In Yang–Mills, analogous OPE coefficients arise from Beta functions multiplied by Mellin normalization ratios. For example, the same-helicity gluon OPE coefficient is
\[
C^{(++)+}(\lambda_1,\lambda_2)
=
i\,\frac{g(\Delta_1)g(\Delta_2)}{g(\Delta_1+\Delta_2)}\,B(i\lambda_1,i\lambda_2)
=
1+\frac{\lambda_1\lambda_2}{(1+i\lambda_1)(1+i\lambda_2)},
\]
with corresponding formulas for \((--)\) and mixed-helicity channels [1903.01676].

Low-point celestial amplitudes exhibit both the promise and the singularity structure of the formalism. Four-point massless scalar and gluon celestial amplitudes admit conformal partial wave decompositions in a basis of global \(SL(2,\mathbb C)\) conformal partial waves \(\Psi_{\Delta,J}(z,\bar z)\), and the four-point correlator can be expanded as
\[
\widetilde{\mathcal A}_4(z,\bar z)
=
\sum_{J=0}^{\infty}{}'
\int_C d\Delta\,\rho(\Delta,J)\,\Psi_{\Delta,J}(z,\bar z)
\]
[1904.10940]. In the scalar exchange example, only \(J=0\) partial waves appear, whereas the four-gluon MHV amplitude contains all integer spins \(J\ge 0\) in its partial-wave decomposition [1904.10940]. Crossing is implemented by Möbius transformations \(z\to 1-z\), \(z\to 1/z\), and \(z\to z/(z-1)\), with the appropriate conformal weight factors [1904.10940].

Three-point amplitudes are more delicate. In plane-wave Minkowski kinematics, real massless three-point amplitudes vanish for finite momenta because momentum conservation cannot be satisfied away from soft and collinear regions. Nevertheless, celestial three-point amplitudes are nonzero because the momentum-conservation delta function has support on soft and collinear configurations, and Mellin or shadow integrals detect that support [2212.07025]. In a conformal-in / shadow-out basis, massless scalar, photon, gluon, and graviton three-point celestial amplitudes take the standard two-dimensional CFT form, and the three-point gluon amplitudes realize a spin-one current with two spin-one primaries [2212.07025].

Beyond four points, conformal block expansions reveal multi-particle operator content. In the comb channel, five-point scalar celestial amplitudes exhibit two-particle operators, while six-point scalar amplitudes exhibit new three-particle operators [2404.01920]. A plausible implication, explicitly conjectured in that work, is that \(n\)-particle exchanges appear in the comb-channel conformal block expansion of \((n+3)\)-point massless celestial amplitudes [2404.01920].

## 5. Massive states, loop amplitudes, and lower-dimensional variants

Massive celestial amplitudes retain conformal covariance but replace the one-dimensional Mellin integral over energy by hyperbolic integrals over on-shell massive momenta. This enables genuinely new three-point structures and higher-point conformal block expansions. In particular, the three-point coefficient with two massive scalars contains a hypergeometric function, and the coefficient with three massive scalars admits a triple Mellin-Barnes representation [2404.01920]. Turning on the mass of one incoming particle in five-point scalar amplitudes produces a new operator series whose leading member can be interpreted as a two-particle exchange in the OPE of one massive and one massless scalar [2404.01920].

Loop amplitudes sharpen the ultraviolet and infrared sensitivity of the celestial basis because Mellin integrals probe all energies. The 2021 lecture notes emphasize that tree-level four-graviton celestial amplitudes are formally divergent on the principal series and that stringy form factors regulate the ultraviolet behavior [2108.04801]. A concrete string-theoretic analysis of genus-one open-string four-gluon scattering shows that the Mellin transform commutes with the appropriate field-theory degeneration limit in worldsheet moduli space, and that the one-loop celestial form factor carries an overall \((\alpha')^{\beta-3}\) instead of the tree-level \((\alpha')^\beta\), making \(g_{10}^2/\alpha'^3\) the natural dimensionless loop-expansion parameter in the celestial basis [2307.03551].

On the field-theory side, one-loop rational all-plus and single-minus gluon and graviton amplitudes provide explicit finite celestial correlators beyond tree level [2007.09338]. Their Mellin transforms yield closed-form four- and five-point celestial amplitudes with the expected \(SL(2,\mathbb C)\) covariance and with the familiar four-point support on real cross-ratios [2007.09338]. A different development recasts MHV celestial amplitudes in terms of Liouville-theoretic partial differential equations; the resulting \(\mathcal O(b^2)\) corrections are logarithmic for both gluons and gravitons, and in the gluon case the deformed celestial OPE is isomorphic to the one-loop correction of the celestial OPE in pure Yang–Mills theory [2409.05936].

Lower-dimensional analogues modify the transform itself. In two-dimensional massive theories, Lorentz boosts act additively on the rapidity \(\theta\), so the celestial transform becomes a Fourier transform in rapidity rather than a Mellin transform in energies,
\[
\mathcal A(\omega)=\int_{-\infty}^{\infty} d\theta\,e^{i\omega\theta}S(\theta),
\]
or, with an \(i\epsilon\) prescription, retarded and advanced transforms \(\mathcal A^\pm(\omega)\) [2209.02776]. In the \(O(N)\) nonlinear sigma model, exact two-dimensional celestial amplitudes are Meijer \(G\)-functions, and crossing becomes
\[
\mathcal A_2(w)=e^{-2\pi w}\mathcal A_2(-w),\qquad
\mathcal A_1(w)=e^{-2\pi w}\mathcal A_3(-w)
\]
[2312.02083]. These two-dimensional constructions isolate analyticity, crossing, and unitarity in a particularly transparent form [2209.02776][2312.02083].

## 6. Reformulations, conceptual tensions, and current directions

Several later developments concern the status of distributional support and the correct notion of a celestial correlator. In \(\mathbb R^{2,2}\), an ambidextrous basis built by composing the Mellin map with improved light transforms yields low-point massless celestial amplitudes that are non-distributional in the celestial coordinates and transform covariantly in all scattering channels [2212.00962]. For the four-gluon tree amplitude, this basis gives a streamlined Grassmannian representation and an alpha-space decomposition that exposes towers of exchanged operators and selection rules between symmetric and antisymmetric \(\mathbb Z_2\) sectors [2212.00962].

A different reorganization decomposes celestial amplitudes into integrals of \(AdS_3\)-Witten diagrams associated with leaves of a hyperbolic foliation of spacetime. For the Kleinian three-point MHV amplitude, each leaf subamplitude is smooth except for the expected light-cone singularities, and the full translationally invariant celestial amplitude is the residue of the pole at \(\sum_i\Delta_i=3\) [2312.07820]. This suggests a precise relation between celestial amplitudes and \(AdS_3\) contact diagrams, with translation invariance emerging only after summing over leaves [2312.07820].

The most explicit recent critique of the conventional Mellin definition is the claim that conventional massless celestial amplitudes are distributional and fail to realize the celestial OPE, most sharply in the non-MHV paradox: OPE arguments predict nonzero tree-level celestial amplitudes with helicities \(-+++\), whereas the conventional Mellin transform of the momentum-space amplitude vanishes [2512.05882]. The proposed resolution is a class of regular celestial amplitudes obtained by smearing slightly off shell and combining massive and tachyonic conformal bases. At tree level these amplitudes are non-distributional and consistent with celestial OPEs, leading to the proposed revised dictionary that CCFT correlators are the regular, not conventional, celestial amplitudes [2512.05882]. This is a substantive conceptual shift rather than a settled consensus.

A complementary geometric reformulation appears in positive-geometry constructions. For arbitrary multiplicity \(n\) of massless states in sufficiently high spacetime dimension \(D\), all Mellin integrations of tree-level \(\phi^3\) amplitudes can be carried out and the result identified with the canonical form of a celestial associahedron [2512.15852]. In that construction, distributional support on the celestial sphere is not imposed separately but arises geometrically from the positivity domain, and the same Mellin machinery extends through scalar scaffolding to gluons and gravitons [2512.15852]. This suggests a geometric interpretation of celestial amplitudes in terms of boundary kinematics, factorization facets, and scattering-equation variables [2512.15852].

Other structural extensions include a celestial double copy, in which momentum-space kinematic numerators are promoted to generalized differential operators acting on conformal primary wavefunctions and then squared, reproducing three- and four-point gravity celestial amplitudes from Yang–Mills data [2007.15027]. The 2021 lecture notes also emphasize unresolved issues: constructing complete celestial CFTs, determining central charges and Kac–Moody levels, understanding loops and IR/UV regularization, incorporating massive spectra systematically, and identifying a nonperturbative flat-space holographic dual [2108.04801]. Taken together, these developments indicate that celestial amplitudes are not merely Mellin transforms of known \(S\)-matrix elements; they are a rapidly evolving framework in which the correct basis, operator algebra, and even the definition of the correlator remain active subjects of research.

Source: https://www.emergentmind.com/topics/celestial-amplitudes