Extended BMS algebra of celestial CFT

We elaborate on the proposal of flat holography in which four-dimensional physics is encoded in two-dimensional celestial conformal field theory (CCFT). The symmetry underlying CCFT is the extended BMS symmetry of (asymptotically) flat space time. We use soft and collinear theorems of Einstein-Yang-Mills theory to derive the OPEs of BMS field operators generating superrotations and supertranslations. The energy momentum tensor, given by a shadow transform of a soft graviton operator, implements superrotations in the Virasoro subalgebra of 𝔟𝔪𝔰4. Supertranslations can be obtained from a single translation generator along the light-cone direction by commuting it with the energy-momentum tensor. This operator also originates from a soft graviton and generates a flow of conformal dimensions. All supertranslations can be assembled into a single primary conformal field operator on celestial sphere.


Introduction
Bondi-van der Burg-Metzner-Sachs (BMS) group is a symmetry of asymptotically flat fourdimensional spacetimes at null infinity [1,2]. It is a semi-direct product of the SL(2, C) subgroup of global conformal transformations of the celestial sphere CS 2 at null infinity (isomorphic to Lorentz transformations) times the abelian subgroup of so-called supertranslations. In 2009, Barnich and Troessaert [3] argued that SL(2, C) should be extended to the group of all local conformal transformations (diffeomorphisms) of celestial sphere, a.k.a. superrotations. In 2013, Strominger showed that such extended BMS is also a symmetry of the S-matrix describing elementary particles and gravitons [4]. This symmetry plays central role in Strominger's proposal of a flat spacetime hologram on CS 2 [5]. 1 In JHEP03(2020)130 particular, superrotations set the stage for two-dimensional celestial conformal field theory (CCFT) on CS 2 which, according to Strominger, should encode four-dimensional physics.
In CCFT, each particle is associated to a conformal field operator. The correlators of these operators are identified with four-dimensional S-matrix elements transformed to the basis of conformal wave packets. They can be obtained by applying Mellin transformations (with respect to the energies of external particles) to traditional, momentum space amplitudes [8][9][10][11]. Within this framework, the operator insertion points z ∈ CS 2 are identified with the asymptotic directions of four-momenta while their dimensions ∆ are identified with the dimensions of wave packets. The conformal wave packets associated to stable, helicity ℓ particles have conformal weights (h,h) = ( ∆+ℓ 2 , ∆−ℓ 2 ) with Re(∆) = 1 [12]. Two of us have recently shown that, given the proper definition of CCFT energy-momentum tensor, the operators representing gauge bosons and gravitons do indeed transform under diffeomorphisms of CS 2 (superrotations) as primary conformal fields [13].
There are two special cases of conformal dimensions that lead to important insights into CCFT. For gauge bosons, it is the "conformally soft" limit of ∆ = 1 [14] in which conformal wave packets describe pure gauge fields. They correspond to asymptotically "large" gauge configurations that have observable effects similar to Goldstone modes. The operators emerging in the ∆ = 1 (h = 1,h = 0) limit of gauge boson operators are the holomorphic currents J a carrying gauge group charges. In ref. [15], we derived OPEs of J a with other operators and showed that respective Ward identities agree with soft theorems. For gravitons, both ∆ = 0 and ∆ = 1 spin 2 wave packets represent (large) diffeomorphisms. Upon a shadow transformation, the operator associated to ∆ = 0 (h = 1,h = −1) graviton (which is outside the Re(∆) = 1 stability domain) changes its dimension to ∆ = 2 (h = 2,h = 0) and becomes the CCFT energy-momentum tensor T generating superrotations [16,17]. The graviton associated to ∆ = 1 (h = 3/2,h = −1/2) yields a primary field operator with a single-derivative descendant P that generates supertranslations. Hence BMS symmetries are controlled by ∆ = 0 and ∆ = 1 soft limits. On the other hand, collinear limits are useful for studying OPEs because identical momentum directions correspond to the operator insertion points coinciding on CS 2 . We will utilize them together with the soft limits in the following discussion of the extended BMS symmetry algebra bms 4 .
While BMS is a highly nontrivial symmetry of asymptotically flat spacetimes, its implementation in CCFT is rather straightforward. Virasoro subalgebra must emerge from the standard OPE of T T products. We also know what to expect from the OPE of the supertranslation operator P with T because P is a well-defined descendant of a primary field. What is non-trivial in the context of CCFT is to show that such OPEs follow from the properties of Einstein-Yang-Mills (EYM) theory of gauge bosons coupled to gravitons. In this work, we extract bms 4 from the collinear and soft limits of EYM amplitudes.
The paper is organized as follows. In section 2, we revisit OPEs of the operators associated to gauge bosons and gravitons previously discussed in refs. [15] and [18]. We rewrite them in uniform normalization coventions. The relevant collinear limits are collected in appendix A, where they are derived by using EYM Feynman rules. In section 3, we derive the OPEs of superrotations generated by the energy-momentum tensor and of supertranslations generated by a descendant of a soft graviton operator, with the operators associated JHEP03(2020)130 to gauge bosons and gravitons. Although the energy-momentum tensor is defined by a nonlocal shadow integral, its OPEs are always localized. In section 4, we use soft theorems [15,[19][20][21][22][23] to derive the OPEs of BMS generators. 2 In section 5 we connect OPEs with bms 4 .

Preliminaries: OPEs of spin 1 and operators
The connection between light-like four-momenta p µ of massless particles and points z ∈ CS 2 relies on the following parametrization: where ω is the light-cone energy and q µ is a null vector -the direction along which the massless state propagates, determined by z. The basis of wave functions required for transforming scattering amplitudes into CCFT correlators consist of conformal wave packets characterized by z, dimension ∆ and helicity ℓ. The starting point for constructing such packets are Mellin transforms of spin 1 and 2 plane waves: where ∂ J = ∂ z for ℓ = +1, +2 and ∂ J = ∂z for ℓ = −1, −2. The conformal (quasi-primary) wave functions can be written as with the normalization constants The presence of these normalization factors makes it clear that, as mentioned in the Introduction, fields with spin 1 become pure gauge when ∆ = 1 while fields with spin 2 become pure diffeomorphisms for ∆ = 0, 1. The CCFT correlators are identified with the S-matrix elements transformed from the plane wave basis into conformal basis (2.4) by using properly normalized Mellin transformations [8][9][10][11]: where M ℓ 1 ...ℓ N are EYM Feynman's matrix elements with helicities ℓ n and c n are the normalization constants JHEP03(2020)130 see eq. (2.5). In eq. (2.6) ǫ n = +1 or −1 depending whether the particles are incoming or outgoing, respectively. We skipped the gauge group indices and the corresponding group factors that can always be written in the basis of single or multiple trace Chan-Paton factors.
OPEs can be extracted from the correlators (2.6) by considering the limits of coinciding insertion points, which on the r.h.s. correspond to the collinear limits of scattering amplitudes. Furthermore, in the cases of operators with ∆ = 1 and 0, the zeros of normalization constants c n must be canceled by "soft" poles. Indeed, finite OPE coefficients appear from such singular soft and collinear limits.
In ref. [15], we derived the OPE coefficients of the products of spin 1 operators representing gauge bosons. They follow from the well-known collinear limits of Yang-Mills amplitudes. For two gauge bosons labeled by gauge indices a and b, with identical helicities, we obtained and a similar expression with (z −w) −1 pole and the same OPE coefficient for two −1 helicities. Two gauge bosons of opposite helicities can fuse into a single operator in the same way as in the case of identical helicities. The form of the corresponding OPE terms [15] follows from the collinear limit of gauge interactions. In the case of opposite helicities, however, gravitational interactions allow the fusion of two gauge bosons into a graviton operator [18]. As explained in appendix A, EYM amplitudes involving gravitational couplings do not blow up in the collinear limit but contain another type of singularity. It is due to a phase ambiguity reflecting an azimuthal asymmetry. In appendix A, we show that when z 1 → z 2 , Here, P is the combined momentum of the collinear pair: with and

JHEP03(2020)130
In eq. (2.10), the last two terms, which originate from gravitational couplings, 3 contain the phase factors z 12 /z 12 andz 12 /z 12 depending on the direction from which z 1 approaches z 2 or equivalently, on the azimuthal angle of the plane spanned by the collinear pair about the axis of the combined momentum P . The "regular" part does not contain either z 12 nor z 12 in denominators and is well-defined, finite in the z 1 = z 2 limit. By following the same steps as in ref. [15], we obtain , (2.16) The above result agrees with refs. [15,18]. The OPEs involving gravitons can be extracted in a similar way. For two gravitons with identical helicities, we show in appendix A that As a result, and a similar expression with the conjugated phase factor and the same OPE coefficient for two −2 helicities.

JHEP03(2020)130
The collinear limit of two gravitons with opposite helicities is also derived in appendix A. It reads: The corresponding OPE is . (2.25) We close this section by listing the OPEs of gauge bosons with gravitons. They follow from the collinear limits of EYM amplitudes discussed in appendix A: .
The OPEs of operators with + and − interchanged have the phase factors (z − w)/(z −w) conjugated and the same coefficients.

OPEs of superrotations and supertranslations with spin 1 and spin 2 operators
In this section, we discuss the OPEs of superrotations generated by the energy momentum tensor T and of supertranslations generated by the operator P , with the operators associated to gauge bosons and gravitons. We know what to expect from the OPE of T with JHEP03(2020)130 primary fields. In ref. [13], we showed that the operators O a ∆,± representing gauge bosons are indeed such primary fields. To that end, we used the energy-momentum tensor defined in [16,17] as the shadow [25] of dimension ∆ = 0 graviton operator O 0,−2 : and took the collinear limit of the shadow with gauge bosons as in eq. (2.27). The limit of ∆ → 0 was taken at the end. Actually, as we will see in this section, there is a fast way of deriving and generalizing the results of ref. [13]. We will first take the ∆ → 0 limit and then take the shadow transform in order to extract the OPEs. We will proceed in the same way in the case of P : by taking the ∆ → 1 limit first.

Superrotations
We start by inserting T (z) into the correlator (2.6): where we used the definition (3.1) and which is the Mellin-transformed amplitude with a "soft" graviton. Its ∆ 0 → 0 limit has been studied in refs. [21][22][23], with the following result: where A on the r.h.s. is the Mellin transform of the amplitude without the graviton. Here, h i is the chiral weight of the ith operator and ξ is an arbitrary reference point on CS 2 reflecting the gauge choice for the graviton. When the reference point is varied, the variation of eq. (3.5) vanishes as a consequence of conformal (Lorentz) Ward identities. In our case, it is convenient to set ξ → ∞. The shadow transform can be performed by writing and integrating by parts three times. After repeatedly using the identity

7)
JHEP03(2020)130 we obtain Note that see eq. (2.5). As a result, we obtain Following the same route for the antiholomorphicT (z), one obtains The OPEs (3.10) and (3.11), which are valid for both spin 1 and spin 2 particles, prove that the respective CCFT operators are primary fields.

Supertranslations
We start by inserting P (z 0 ) into the correlator (2.6): where we used the definition (3.2). In this case, the relevant ∆ 0 → 1 soft limit is [21][22][23]: 13) In this case, varying the reference point ξ does not change the form of eq. (3.13) because its variation vanishes a consequence of Ward identities associated to momentum conservation [26]. Here again, as in eq. (3.5), we can set ξ → ∞. From eqs. (3.12) and (3.13) it follows that Note that As a result, we obtain and similar expressions forP (z), with complex conjugate poles. 4 Note the presence of (∆ − 1) factors in the above OPE coefficients. They imply that the products P (z)J a (w), P (z)J a (w), P (z)P (w) and P (z)P (w) are regular. P (z) is not a holomorphic current: it is a dimension 2 descendant of a primary field, with SL(2, C) weights (3/2,1/2). Nevertheless, the r.h.s. of eq. (3.14) implies the same Ward identity as translational symmetry of scattering amplitudes discussed in refs. [26] and [27]. In fact, P contains the P 0 + P 3 component of the momentum operator which generates translations along the light-cone directions. The remaining components can be obtained by superrotations of P . This will be discussed in section 5.

T T
The product of energy-momentum tensors involves two shadow transforms. Let us define We take the limit ∆ 0 → 0 first, as in eq. (3.5): which leads to (4.5) In eq. (4.4), the integrals over z 0 can be performed by using the conformal integrals of ref. [28], see appendix B. With the reference point ξ → ∞, the result is The same result can be obtained by a direct application of Ward identity (3.8).
It is clear from eq. (4.6) that the second term (involving the sum over i ≥ 2) is nonsingular in the limit of w → z, therefore only the first term needs to be included in the derivation of OPE: In order to simplify notation, we introduce the following variables: We anticipate a single pole term of the following form: where in the last step we used integration by parts. Indeed, we can rewrite the r.h.s. of (4.7) as

JHEP03(2020)130
In the last line, the first term can be written as (4.11) After combining eqs. (4.10) and (4.11), we obtain: (4.12) The final step is to show that the last term vanishes due to global conformal invariance.
To that end, we take the soft limit ∆ 1 → 0: After choosing the reference point σ = z, the derivative terms become where we used conformal integrals of ref. [28], see appendix B, in particular eq. (B.3). 5 The remaining terms can be written as (4. 16) In this way, the last term of eq. (4.12) becomes The integrals are divergent as it is obvious from the presence of the Γ(0) prefactor. We can regularize the divergence by taking the soft limit on the shadow transform at the end. In other words use and take λ → 0 at the end. Our results will be the same.

JHEP03(2020)130
with both terms vanishing separately as a consequence of global conformal (Lorentz) invariance [26]: The final result is the expected T (w)T (z) OPE:

T T
By repeating the steps leading to eq. (4.7), we obtain which can be simplified to After using we obtain Note that in the soft limit therefore the derivative∂ 3 z acting on the r.h.s. of eq. (4.23) gives zero modulo delta-function terms localized at z = w and z = z n . These terms do not affect the Virasoro algebra following from eq. (4.19). We conclude that up to such delta functions, T (w)T (z) = regular .

Recall (3.2) that the supertranslation current is defined as
The graviton primary field operator O 1,+2 has chiral weights ( 3 2 , − 1 2 ). We start from (4.27) In the ∆ 0 → 0 limit, cf. eq. (4.3) with ℓ 1 = +2 and the reference point ξ → ∞, with G + defined in eq. (4.21) and h 1 = 3 2 . As in previous cases, the sum over i ≥ 2 is non-singular at z = w, therefore we are left with ∂z lim Upon inserting the above expression into the shadow integral of eq. (4.27), the first two terms give rise to delta-function singularities at z = w, but vanish for z = w. In particular, where we used conformal integrals (B.1) and (B.5) from appendix B. Similarly, In this way, we obtain O(z n ,z n ) + · · · JHEP03(2020)130 After performing the shadow integrals, we obtain T (w)P (z) = 3 2(w − z) 2 P (z) + 1 w − z ∂ z P (z) + regular , (4.33) modulo delta-function singularities at z = w. Next, we consider (4.34) After taking the ∆ 0 → 0 limit, we obtain where, as in the previous case, we ignored delta-function singularities at z = w. As a result, which is the OPE expected for a level one descendant of a primary field.

Extended BMS algebra
In this section, we discuss connections between the OPEs and extended BMS algebra bms 4 ; we also make contact with ref. [29]. The OPEs of the stress energy tensor confirm the presence of Virasoro subalgebra with zero central charge. 6 Here, we concentrate on the role of supertranslations and derive the complete form of the algebra based on the OPEs. It will be convenient to work with primary fields normalized in a slightly different way. Up to this point, we considered The correlation functions of "bare" φ fields are given directly by Mellin-transformed amplitudes: In the following discussion, the primary field operators will be identified with φ h,h (z,z), although the normalization of soft graviton operators O 0,±2 and O 1,±2 , and of the related T and P will remain unchanged.

JHEP03(2020)130
The four-momentum (translation) operators were constructed in ref. [26]. Anticipating their placement in bms 4 , we define We will use P − 1 2 ,− 1 2 as a starting point for constructing all supertranslation generators. To that end, we will utilize the extended BMS algebra [29]. In addition to the Virasoro subalgebra, bms 4 consists of [L n , P kl ] = 1 2 n − l P k,n+l + n(n 2 − 1)C k,n+l , with m, n ∈ Z and i, j, k, l ∈ Z + 1 2 . Here C andC are possible central extension operators. The subset P n− 1 2 ,− 1 2 can be generated in the following way: In order to find out how these operators act on a generic primary field, we first determine the OPE of T (w), P − 1 2 ,− 1 2 and φ h,h (z,z). We use where the factor −2 is coming from the ratio of normalization factors, f (0)/f (1) = −2, and take the ∆ = 1 soft limit as in eq. (3.13), to show that After taking the shadow transform, we obtain Finally, after combining eq. (5.6) with the above OPE, we obtain which agrees with eq. (5.4) for n = 0, 1 and is in agreement with supertranslations of primary field operators proposed in ref. [29]. The remaining bms 4 generators can be constructed as By following the same steps as before, it is easy to show that Following ref. [29], we can combine all supertranslation generators into one primary field operator P(z,z) ≡ n,m∈Z P is a scalar operator with dimension 3, (h,h) = ( 3 2 , 3 2 ). The commutators (5.12) are equivalent to the OPE At this point, we can make connection to the supertranslation operator P (z) discussed in sections 3.2 and 4.3. For bare fields, the OPEs of eqs. (3.17) and (3.18) read and compare eq. (5.15) with (5.12), we find 4 P n− 1 and a similar expression forP (z). Note that both P (z) andP (z) miss the "mixed" operators P n− 1 2 ,m− 1 2 with n = 0 and m = 0 simultaneously. This is to be contrasted with the primary conformal field operator P(z,z) of eq. (5.13) which includes all BMS supertranslations.
The operator P − 1 2 ,− 1 2 = P 0 + P 3 plays a special role in our construction of BMS generators. With a primary field expanded as shifts the modes as: Four-dimensional Minkowski space can be foliated by using Euclidean AdS 3 slices [7,17]. On the other hand it is known that string theory on AdS 3 has a spectrum generated by using a spectral flow transformation on the standard representations of SL(2; R). This transformation shifts conformal dimensions and the mode expansions of primary fields. The transformation in (5.19) is akin to the spectral flow encountered in string theory on AdS 3 [30].
To summarize, we used the correlators of primary fields, obtained from tree-level, fourdimensional scattering amplitudes, to derive the extended BMS algebra (5.5). The method employed here shows that the central charges C =C = 0. We cannot exclude though that these charges are non-zero after we take into account quantum loop corrections.

JHEP03(2020)130 6 Conclusions
Extended BMS is the symmetry of asymptotically flat spacetime and consequently, of its hologram on CS 2 . According to Strominger's proposal of flat holography, four-dimensional physics is encoded in two-dimensional CCFT. In particular, Mellin-transformed fourdimensional scattering amplitudes are equivalent to the correlators of primary field operators on CS 2 . In this work, we followed this link. We investigated the structure of OPEs of the operators generating BMS transformations by using the relation of the collinear limit to the limit of coinciding points on CS 2 and the relation of the usual soft limit to the conformally soft limit of operators with integer dimensions. We established a connection between these OPEs and the extended BMS algebra bms 4 . We also showed how the OPEs of the operators representing gauge bosons and gravitons follow directly from the Feynman diagrams of Einstein-Yang-Mills theory.
BMS superrotations are implemented in a rather straightforward way as diffeomorphisms of CS 2 . They are generated by the energy-momentum tensor T (z) given by a shadow transform of the soft graviton operator with dimension 0. Although general shadow transformations are non-local, we found that for soft gravitons, the shadow integrals become localised at operator insertions. As a result, we recovered the standard OPEs which are equivalent to the Virasoro subalgebra without central extension. 7 Our analysis was based on the tree-level approximation in Einstein-Yang-Mills theory, hence it is possible that a Virasoro central charge will appear after taking into account loop corrections.
BMS supertranslations are more subtle. We investigated the OPEs of the supertranslation operator P (z), a dimension 2 descendant of the soft graviton operator with dimension 1. Indeed, the OPEs of this operator generate shifts of primary field dimensions, as expected from translations and supertranslations. In particular, it contains the momentum component P 0 + P 3 which generates a flow of chiral weights, (h,h) → (h + 1 2 ,h + 1 2 ). By comparing with BMS algebra, we found however that P (z) andP (z) contain only a subset of supertranslations. They miss not only the non-holomorphic P 0 − P 3 but also an infinite set of supertranslations. Nevertheless, we were able to generate the full set by commuting T (z) andT (z) with P 0 + P 3 which we identified with the P − 1 2 ,− 1 2 spectral flow generator of bms 4 . All supertranslations can be assembled into a single primary field operator P(z,z) of dimension 3.
In CCFT, 2D correlators of primary fields arise from 4D scattering amplitudes expressed in the basis of conformal primary wavefunctions. As shown by Pasterski and Shao [12], the continuum of conformal wavefunctions with dimensions ∆ = 1 + iR spans the complete set of integer spin solutions of wave equations, normalizable with respect to the Klein-Gordon norm. The primary fields associated to gauge bosons and gravitons have complex dimensions, therefore the theory is not unitary with respect to the usual radial quantization. It is important to keep in mind that the notion of unitarity always relies on a choice of time direction and a norm: radial quantization (where the "time" is the Euclidean radial direction) is not natural on the celestial sphere. 8 The questions what is JHEP03(2020)130 the relation between unitarity in 2D and 4D and what is the role of central charges need to be addressed in the future.
There are many directions for future investigations of CCFT. Since the present discussion relies on the properties of tree-level Einstein-Yang-Mills theory, it would be very interesting to see how four-dimensional quantum corrections are implemented in CCFT. But first and foremost, it is very important to understand the spectrum of CCFT. It is clear that supertranslations generate states very different from conformal wave packets of stable particles. This sector, together with the role of the spectral flow, are certainly worth investigating.

Acknowledgments
We are grateful to Glen Barnich

A Collinear limits in EYM theory
In this appendix, we use the notation and conventions of ref. [31]. In order to give a precise definition of the leading and subleading collinear parts, we need to specify how the collinear limit is reached from a generic kinematic configuration. Let us specify to generic light-like momenta p 1 , p 2 and introduce two light-like vectors P and r such that the momentum spinors decompose as [32]: λ 1 = λ P cos θ − ǫλ r sin θ ,λ 1 =λ P cos θ −ǫλ r sin θ , (A.1) λ 2 = λ P sin θ + ǫλ r cos θ ,λ 2 =λ P sin θ +ǫλ r cos θ , (A.2) hence p 1 = c 2 P − sc(ǫλ rλP +ǫλ Pλr ) + ǫǫs 2 r , (A.3) We also have 12 = ǫ P r , The total momentum is: issue.

JHEP03(2020)130
The collinear configuration will be reached in the limit of ǫ =ǫ = 0. With the momenta parametrized by celestial sphere coordinates as in eq. (2.1), this limit corresponds to z 1 = z 2 , i.e. p 1 + p 2 = P ≡ ω P q P , (A. 8) with Yang-Mills amplitudes are singular in the collinear limit. They contain 1/ǫ and 1/ǫ poles while subleading terms are of order 1 [32]. Here, we are interested in gravitational couplings. These contain singular ratios 12 [12] = ǫ ǫ P r Here, ǫ ± P are the polarization vectors of a massless vector boson with momentum P with the gauge reference vector r. Note that , (A.14) .
For all polarization vectors labeled by i = 1, 2, P , Other scalar products are also easy to obtain: The singularities (A.11) appear only in Feynman diagrams involving "splitting" of a massless virtual particle into the collinear pair. They originate from the propagator poles 1/(p 1 + p 2 ) 2 ∼ (ǫǫ) −1 . The vertices describing gravitational splitting "soften" these poles by terms which are quadratic in ǫ andǫ. Diagrams without such splitting are non-singular, therefore collinear singularities have a universal form that can be extracted from three-point vertices. In our discussion, we will be using the graviton polarization tensors

B Conformal integrals
In the present paper, several integrals over celestial sphere positions have been performed by using the formulas given in [28]. They have the general form

JHEP03(2020)130
The two conditions above are necessary for I n to be covariant under the group SL(2; C) of global conformal transformations and single valued. Convergence of the integral requires q i +q i < 2 for all i, but we can always extend it by using analytic continuation. When n = 3, assuming that all z i are different, the result is given by where the normalization factor If two z i 's are coincident, the integral becomes Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.