Asymptotic charges in massless QED revisited: A view from Spatial Infinity

Hamada and Shiu have recently shown that tree level amplitudes in QED satisfy an infinite hierarchy of soft photon theorems, the first two of which are Weinberg and Low's theorems respectively. In this paper we propose that in tree level massless QED, this entire hierarchy is equivalent to a hierarchy of (asymptotic) conservation laws. We prove the equivalence explicitly for the case of sub-subleading soft photon theorem and give substantial evidence that the equivalence continues to hold for the entire hierarchy. Our work also brings out the (complimentary) relationship between the asymptotic charges associated to soft theorems and the well known Newman-Penrose charges.

1 Introduction and outline The relationship between soft theorems in gauge theories and gravity and asymptotic symmetries in these theories is an active area of investigation. The essential idea is simple and can be understood without referring to specific theories. Given a large gauge transformation 1 which is parametrized by a gauge parameter ǫ on a cross section S 2 of null infinity I, the soft theorems are statements regarding conservation laws: where ǫ ± parametrizes large gauge transformations at I ± and are identified via anti-podal matching conditions. Q ± [ǫ ± ] are charges which are evaluated on the celestial sphere S 2 at u = −∞ (I + − ) and v = +∞ (I − + ) respectively. Thus, a rather natural question to ask is, if there indeed are such infinity of conservation laws in theories like Electrodynamics, should we not be able to derive them in classical theory? First step in this direction was taken in [1] where it was shown that the infinity of charges associated to large U (1) gauge transformations were indeed conserved in classical theory (from previous works [2,3] these conservation laws in quantum theory were known to be equivalent to soft photon theorem). This conservation law was derived by analyzing the equations of the theory at spatial infinity. By considering a compactification scheme where spatial infinity is a hyperboloid H (a three dimensional Lorentzian de Sitter space in fact) as opposed to a point, one could relate fields at u = −∞ and v = +∞ by using the field equations.
We revisit this idea below and show that it leads to (classical) conservation law for charges associated to leading as well as subleading soft photon theorem. 2 However this derivation leads to an interesting consequence. That, if we restrict ourselves to a suitable subset in the set of all radiative data, there is in fact an entire tower of conservation laws that can be shown to be valid. A priori, it is not clear what the Ward identities corresponding to these hierarchy of conservation laws imply in quantum theory.
We thus have the following question to answer: (A) If such infinite tower of charges are conserved, we would expect (along the lines of leading and sub-leading soft photon theorems) an infinite hierarchy of soft theorems. Do such theorems exist? 1 In this paper we will use the expressions "asymptotic symmetry" and "large gauge transformation" interchangeably. 2 Subleading soft photon theorem is non-universal [4], however this most general form of subleading theorem can be understood in terms of asymptotic symmetries [5].
As it turns out, from the side of soft theorems, a complimentary puzzle already existed: (B) In [6,7], it was shown that in fact for tree level scattering amplitudes in QED, there do exist an infinity of soft theorems. From the perspective of these soft theorems, a natural question to ask would be, just as leading and sub-leading soft theorems are Ward identities for certain asymptotic symmetries. Is the same true for the higher order theorems?
In this paper we present substantial evidence that questions (A) and (B) mutually answer each other. That is, the conservation of asymptotic charges (beyond the ones associated to sub-leading soft theorem) are equivalent to sub-n (n > 1) soft theorems.
The outline of the paper is as follows. In section 2, we revisit the asymptotic analysis of Maxwell's equations at null infinity. For simplicity, we consider charged massless scalar fields coupled fo U (1) gauge fields but our analysis can be generalized to the situation when charged fields are massive and are scalars or Fermions. We then revisit the derivation of asymptotic charges associated to leading and sub-leading soft photon theorems. Our derivation of these charges is along the lines of [8,9] in that we first obtain them as integrals over S 2 at u = −∞, and then show that they can be written as fluxes over I. The reason for revisiting these charges is that our derivation is amenable to a direct generalization to an infinite hierarchy of further conservation laws. In section 3 we review the infinity of soft theorems derived in [6,7] for tree level scattering amplitudes and show how these theorems can be written in terms of Ward identities. In section 4 we argue that these Ward identities precisely correspond to the hierarchy of conservation laws proposed above. Finally in section 5, we show how this infinite hierarchy of conservation laws are indeed true in classical theory, provided one restricts attention to certain subset of radiative data that is compatible with tree-level scattering. We end with certain speculations and future directions.

Maxwell equations at I + and asymptotic charges
In this section, we review the asymptotic expansion of Maxwell fields at null infinity and study Maxwell's equations in an 1 r expansion. We will work in terms of self-dual fields since this simplifies the field equations (bringing them into a form equivalent to the Newman-Penrose formulation [10]) and because the charges associated to (positive) negative helicity soft photon theorem can be written in terms of (anti) self-dual fields. Our first aim will be to understand how the self-dual field can be determined order by order in 1 r in terms of free data at I + . At each order there will appear new 'integration constants' that will be interpreted as asymptotic charges under the assumption of certain strong |u| → ∞ falloffs. We will discuss in detail the first two set of asymptotic charges and show how they correspond to the charges associated to the leading and subleading (negative helicity) soft photon theorems. We will finally discuss the relation with Newman-Penrose charges.
We consider U (1) gauge field minimally coupled to massless scalar field. In order to analyze the behavior of the radiation fields at null infinity, we rewrite the equations of motion in terms of retarded coordinates (u, r, z,z) as, where sphere indices are raised with the unit-sphere metric γ AB . These should be supplemented with Bianchi identities ∂ [a F bc] = 0. An alternative description can be given in terms of the self-dual field strength as follows. Define, dual, self-dual and anti-self-dual fields by: Then the self-dual field satisfies the equations From the perspective of soft theorems, the self-dual and anti-self dual equations are better suited than the ordinary Maxwell's equations in terms of real fields since the positive and negative helicity soft photon theorems are related to charges constructed from the F ∓ ru fields. For definitiveness we work with self-dual field F + ab whose quantization gives single particle states associated to negative helicity photons.
The standard fall-off conditions (in 1 r ) which accommodates radiative data are 3 As shown in appendix A, Maxwell equations for self-dual fields gives rise to the following recursive equation for n F ru : (2.10)

u → −∞ fall-offs and candidates for asymptotic charges
As shown in [3], in the case of massless QED, the asymptotic charge associated to large gauge transformations is given by 4 Now, the standard u → −∞ fall-offs for the radiative data associated to generic solutions of Maxwell's equations is [14] Thus the non-vanishing and finite limit of . Whence the charge density which gives rise to the leading (negative-helicity) soft theorem is given by Note that (2.9) together with (2.10) implies n F ru = O(u n )+O(u n−ǫ ). However consider a subset of radiative fields with the following fall-off in u: where the remainder R(u) falls off faster than |u| −n ∀ n.
Then using Eq. (2.9), the u → −∞ behavior of n F ru can be specified to be One may think that this subspace is too restrictive to be physically interesting. However, as we show in appendix B, these fall-offs are equivalent to the assumption that radiative data for gauge field A A (ω,x) has a Laurent expansion in ω which begins at 1 ω . This precisely corresponds to the case of tree level scattering amplitudes where absence of infrared loop effects imply that Laurent expansion remains valid [13].
For this space of radiative data, there is a natural prescription to extract out a finite and non-trivial "moment" of n F ru as u → −∞. From Eq. (2.14) we see that n,n F ru is the u-independent term of n F ru as u → −∞. This term may be thought of as an integration constant that arises when solving Eq. (2.9). In terms of the lower order coefficients n,k F ru (x), k < n this quantity can be expressed as: As we show below, each of these n,n F ru (x) generates an infinity of asymptotic charges which are in fact conserved in classical theory so far as we restrict ourselves to the subspace of radiative data defined above. We illustrate this idea with an example of charge whose Ward identities lead to subleading soft photon theorem for tree level scattering amplitudes. That is, we will see how this charge (or equivalently fluxes which are integrals over I + ) which was defined in [15,16] is nothing but

The subleading charge
We would like to compute Now using Eq. (2.9) we have, In the u → −∞ limit the current term vanishes. Evaluating (2.18) with the expansion (2.14) one finds To evaluate this "charge density" in terms of the free data, we write (2.20) as an integral over u (here we assume there is no contribution from u = +∞, which is consistent with the absence of massive charges), where we used Eq. (2.18). The ∂ u 0 j r , being a total u-derivative of a current term will not contribute and thus we have been able to express the charge-density as a flux (per unit solid angle) at I + as We immediately see that this charge matches with the charges obtained in [16] which was shown to be associated to sub-leading soft photon (of negative helicity) insertions. In fact as shown in [16], Q 1 [ǫ + ] equals the charge obtained in [15] if we define An exactly analogous procedure yields the corresponding charge Q − 1 [ǫ − ] at I − . Along with anti-podal matching conditions on ǫ ± , the conservation of charge would follow if we could show that (2.26) We emphasize that in contrast to ideas and viewpoint employed in [16], here we never had to use covariant phase space techniques neither "divergent" gauge transformations. Our analysis only refers to radiative data and structures available at null (and as we see below, spatial) infinity.

Towards a tower of asymptotic charges
Naturally we are now tempted to consider higher order "charge densities" 5 n,n F ru (x) defined at I + − . As a first example of this kind, in this section we focus on 2,2 F ru (x) and write the corresponding charge Q + 2 [ǫ + ] as an integral over I + as a functional of the free data. We will refer to this charge as sub-subleading charge as the corresponding Ward identities will turn out to be equivalent to the ("projected" [7]) sub-subleading soft photon theorem for tree-level QED. 6 That is we would like to evaluate, Using, Eq. (2.9) we can immediately see that Evaluating the equation in the u → −∞ and using (2.14) as before one finds Of course in order to justify calling them charge densities, we need to show that the corresponding charges are indeed conserved. This will be shown in section 5. 6 The soft expansion at sub-subleading order does not factorize, and hence we have to project out the unfactorized contribution to get a sub-subleading soft photon theorem. We refer to these statements as projected soft theorems.
From the first equation we find how to cancel the u 2 divergence of 2 F ru : As before, we now write (2.34) as an integral over u, assuming no contribution arises from u = +∞. Using (2.28) and (2.18), one again finds cancellation of terms, resulting in: (2.35) The expression can be further simplified by integrating by parts the term proportional to u∂ u u j r and dropping the total derivative term ∂ u 1 j r .
We note that unlike the leading and sub-leading charge densities (i.e. is a sum of two kinds of terms. First set of terms involve the integrand which is a local function of the free data, φ(u,x), A A (u,x) and in the second case, the integrand is a non-local function of the free data involving u φ(u ′ ,x) du ′ .
Using the above charge densities, we can define the corresponding charges, which are parametrized by functions on the sphere ǫ(x) as, . In fact, one can define a tower of asymptotic charges labelled by n ≥ 0 and ǫ n ∈ C ∞ (S 2 ): . We postpone to section 4 a detailed description of these charges for arbitrary n. For now we would like to comment on a particular feature of these charges that we will discover below: The spherical harmonic decomposition of n,n F ru (x) starts at l = n. In other words, if we take ǫ n = Y l,m , the charge is non-zero only for l ≥ n. It turns out that the l = n − 1 case corresponds to the so called Newman-Penrose charges [10].

Relationship with Newman Penrose charges
We thus see we have a "doubly infinite" family of charges parametrized by n, ǫ n ∈ C ∞ (S 2 ) as in Eq. (2.39). As shown above and later in section 4, these charges (which are localised at u = −∞ or v = +∞) can be written as fluxes integrated over entire null infinity, if we assume that the radiative fields are trivial at v = −∞ and u = +∞.
In [10], Newman and Penrose showed that in a free theory with massless fields (that is pure Maxwell theory in our case), one could construct infinitely many "charges" defined by taking a celestial 2-sphere located at u = constant (or analogously, v = constant) and integrating certain densities which were constructed out of radiative data and spherical harmonics Y l,m . We note here that in the absence of sources, these charges are precisely an (infinite) subset of the asymptotic charges constructed above. 7 That is if we consider the subset of charges parametrized by (n + 1, Y n,m ), then in vacuum these are the Newman Penrose (NP) charges.
NP charges are defined as follows (our choice of normalization is for later convenience) Although the charge density is evaluated at a fixed u, by using vacuum equations of motion we can easily see that the charge is independent of u. Substituting Eq. (2.41) in Eq.
(2.40) and integrating over the sphere, we see that due to defining equation of the spherical harmonics we have We can now see how these charges are related to the asymptotic charges we are considering in this paper. As our charges are in terms of F ru instead of F rz we need to find a relation between the two. This follows from the equation of motion involving the self-dual field, Asymptotic expansion at future null infinity yields, 7 See [9,11,12] for earlier explorations between soft theorems charges and NP charges.
Using Eq.(2.44) in Eq.(2.40) we see that NP charge is given by As these charges are finite, they can be evaluated at any u and in particular at u = −∞ they can be written as which are same as the asymptotic charges We conclude by noting that our fall-off conditions imply the vanishing of the NP charges. This corresponds to the observation in [10] that the charges (at future null infinity) vanish for 'outgoing solutions' (see beginning of p. 184 in [10]). The vanishing of NP charges can also be seen as a consequence of the relation with the soft theorem charges: In section 4 we will see that the asymptotic charges can be written as where ρ z...z depends on the free data at null infinity. Since (D z ) n+1 Y n,m = 0 this implies the vanishing of the NP charges.

From Ward identities to sub-n soft theorems
We would now like to show that if we assume Q + 2 [ǫ + ] = Q − 2 [ǫ − ] then the corresponding Ward identity is equivalent to the sub-subleading soft photon theorem in tree level scattering amplitudes [6,7].
However before establishing this equivalence we take a small detour and review the hierarchy of subleading soft theorems in tree level QED. At the 1st and 2nd order of hierarchy, these are nothing but sub-leading and sub-subleading soft photon theorems. As we will see, the higher order soft theorems are far less constraining then the previous ones, but are present nonetheless. We will refer to this entire hierarchy as sub-n (n ≥ 1) soft theorems

sub-n soft theorems
Consider an un-stripped (that is, including the momentum conserving delta function) treelevel amplitude in QED consisting of N -charged particles and a photon of energy ω. 8 Regarded as a function of ω, the amplitude has an expansion of the form As is well known, the first term is given by Weinberg's soft theorem 9 where ǫ µ is the photon polarization and q = (1,q) the photon 4-momentum direction. The next term in the soft expansion also factorizes and is given by Low's subleading soft theorem: where J µν = p µ ∂ ν − p ν ∂ µ is the angular momentum of the particle with momentum p. 10 For tree level scattering amplitudes and if we restrict to minimal coupling, in fact more is true. In [7] it was shown that gauge invariance implies the n-th term in (3.1) for n ≥ 1 is given by where A µν 1 ...ν n−1 is antisymmetric under the exchange of µ and a ν index but its dependence on the hard momenta is undetermined by the requirement of gauge invariance. Thus at the onset it appears that there is no factorization theorem at sub-subleading order and beyond in tree-level QED. However as was shown in [6], the second term in Eq. (3.6) can be projected out and one obtains, what we call sub-n (n ≥ 1) soft theorem. This can be understood as follows. The n-th term in (3.1) can be extracted as where R (n) denotes the reminder term in (3.6) For definitiveness in the following we restrict attention to the case where the soft photon has negative helicity. Parametrizing the soft momentum direction q µ in terms of 9 To simplify expressions we display signs as if all particles are outgoing and have positive charge. For negative outgoing or positive incoming charges the factor comes with opposite sign. 10 The sub-leading soft photon theorem is not universal even for tree-level scattering amplitudes. In fact it was shown in [4] that there is a class of higher derivative terms which can modify the sub-leading factor. Thus the most general sub-leading factor comprises of the universal term S (0) given above and an additive non-universal term. In this paper we restrict ourselves to universal sub-leading factor.

Ward identities from soft theorems
In fact, from Eq. (3.11), we can write down the Ward identities which are equivalent to these soft theorems. Here we follow the same strategy that is used in deriving Ward identities from soft theorems in the leading and sub-leading case. We first note that the LHS of (3.11) corresponds to insertion of an operator where a − (ωq) is the Fock operator of a negative helicity photon. In analogy with other soft theorems we can construct a smeared version of (3.12) such that it takes a simple form at null infinity. Recall the 'free-data' of the Maxwell field at null infinity is given by the angular components of the vector potential, A A (u,x). These are related to the Fock operators by Aw(ω,q) = 1 4πi where Aw(ω,q) is the time-Fourier transform of Aw(u,q). We can hence see that by smearing both sides of Eq.(3.11) by d 2 w T w...w we obtain a formal 11 Ward identity with a 'soft' charge given by where is the holomorphic component of a rank n (symmetric, trace-free) sphere tensor that parametrizes the charge. In going from the first to second line we used that √ γγ ww = 1.
From (3.13) and (3.14), the soft charge can be written as the operation Formal in the sense that we need to show that this identity arises from conservation laws.
acting on (3.12). In other words, the smearing (3.16) acting on the LHS of the sub-n soft theorem (3.7) can be interpreted as arising from the insertion of a soft charge (3.14). In order to look for the existence of a hard charge, we need to perform the same smearing on the RHS of (3.11).
Up to an i/e factor, such smearing defines the following differential operator on the p variable, where Kw := 1 2πe (the minus subscript denotes the negative helicity of the soft photon under consideration). We will see this operator satisfies two key properties: The first property will be a consequence of the identity [15]: where (E, z,z) parametrize the momentum p of the hard particle. The second property will ensure the smearing (3.16) on the RHS of (3.7) can be understood as arising from a hard charge: 12 Classically, Eq. (3.19) is the condition that the infinitesimal transformation δb = Tb, δb * = (−1) n Tb * is symplectic. The charge reproducing (3.21) and (3.22) will then be given by: where (b ↔ c) is the contribution from the antiparticles.
To compare with the usual definition of hard charges which are integrals over I, we will finally need to write (3.23) in terms of the free data of the scalar field at I, φ(u,x). Equivalently we can write the hard charge as integral over The (−1) n sign arises because, in the factorization formula, an incoming 4-momentum is expressed as outgoing by reversing its sign.
From these expressions, and using the fact that under E → −E, T → (−1) n+1 T, the charge (3.23) can be written as We now calculate T and associated charges Q n [T ] = Q soft n [T ] + Q hard n [T ] for n = 1, 2. We expect these charges to be related to the subleading and sub-subleading charges Q n [ǫ], n = 1, 2 which were defined in the previous section.

n = 1 case: subleading charge from soft theorem
When n = 1 one finds It is easy to verify which corresponds to property (3.19). To compute the charge (say at I + ) we use the identities One can then verify (3.25) becomes On the other hand, the soft charge (3.14) can be written as Finally, taking T z = −2D z ǫ with ǫ(x) a function on the sphere one can verify the total (hard plus soft) charge takes the form with the charge density 3.4 n = 2 case: sub-subleading charge from soft theorem For n = 2 one finds that, One can check it verifies condition (3.19). To compute the charge we use the identities One then finds (3.25) to be given by (3.37) Thus hard charge is a sum of two terms. One term is an integral over local functionals of the radiative data and the other term involves non-local fields which are then integrated. In the above equations, these are the terms which involve uφ (or u φ).
On the other hand, the soft charge (3.14) can be written as To compare with 2,2 F ru (x) obtained in the previous section, we take As shown in subsection 3.4.1, there is a similar matching between the 'non-local' terms, so that the total charges coincide: Whence the question we would like to ask is if the charges Q n [T ] defined from the sub-n soft theorems are also the same as Q n [ǫ] defined at I when n > 2. We turn to this question in section 4.

Matching the non-local terms
We want to show the matching of the 'non-local' terms in Eq. (3.41). This corresponds to the equality We start with the integrand of the LHS of (3.42) using the expressions for the current in terms of the scalar field given in Eq. (C.8). Discarding total spherederivative terms we have: (we omit a "−(φ ↔φ)" at the end of each equation) where we used ∆f = 2D z D z f = 2D z D z f for a scalar function f . The last two terms in (3.46) can be brought to a different form by integrating by parts in u. Up to total u and sphere derivatives one can show the identity for any u-independent sphere function f (x). Using this identity for f = ∆ǫ, as well as the sphere-derivative relations one can show:

Higher order charges and a conjectured equivalence
In this section we extend our previous discussion to the higher order soft theorems. That is, we would like to argue that sub-n soft theorems for n > 2 are equivalent to Ward identities associated to Q n [ǫ]. Although we do not provide a complete proof, we give several hints in this direction and conjecture the equivalence between sub-n soft theorems and Ward identities of higher-n charges. Before embarking on a tedious analysis, let us summarize what we are able to prove below. We show that if one defines a relationship between T and ǫ as, For the hard charges, we are able to bring Q n [ǫ] hard and Q n [T ] hard into a form that allows for a direct comparison. Establishing their equality however becomes a technical problem we have not been able to solve for arbitrary n. The first step in understanding the relationship between Q n [T ] and Q n [ǫ] is to write Q n [ǫ] in terms of free data at I ± . As always we focus on analyzing the charge at I + .
In order to simplify expressions, let us introduce the following notation: Using the field equations (4.6) it is straightforward to show that Properties (4.8) and (4.9) can be used to express n,n F ru (x) in terms of data at null infinity by the same procedure as in the n ≤ 2 charges studied before. We start with Eq. (4.9) and an integration by parts in u: (4.10) In the absence of massive fields the last term vanishes (this term should otherwise give the contribution to the charge from the massive particles) and one is left with the integral over u. From Eq. (4.8) and Eq. (2.10), (4.14) We finally need to express the current coefficients in terms of the scalar field radiative data φ. For the tree-level charges of interest here the scalar field can be treated as free when computing the current. For the first term of (4.14) we have 0 j u = ieφ∂ uφ + c.c. For the second term, one can show the combination of current coefficients can be written as where k ϕ(u,x) is the 1/r k+1 coefficient in the 1/r expansion of the scalar field ϕ. These coefficient can be expressed in terms of the free data by recursively solving the free field equation where ∂ −k u denotes the k-th u-primitive. Substituting (4.17) in (4.16) and then in (4.14) gives an expression of n ρ hard (u,x) in terms of the scalar field free data φ(u,x).

Higher-n charges from soft theorems
From the discussion of section 3.2, the (smeared) sub-n soft theorem can be written as [T ] given in Eqs. (3.14) and (3.25) respectively. We now discuss how to cast these quantities in terms of free data at null infinity. For the soft charge, we simply Fourier transform (3.14) to obtain For the hard charge we need to bring the differential operator T defined in Eq. (3.17) into a simpler form. As in the n = 1, 2 cases, the integral localizes in the hard momentum direction. The general form, derived in appendix D.2, is found to be Substituting this in the expression for the hard charge (3.25) and using the Fourier transform identities

Conjectured equivalence
We wish to generalize the equivalence found in section 3 between n,n F ru and the soft theorem charges Q n [T ]. The idea is to parametrize the tensor T in terms of a sphere function ǫ(x) as T ∼ (D z )ǫ, and then identify ǫ(x) with a smearing function for n,n F ru (x). That is, we wish to find T ǫ such that (4.23) Using the identity (D z ) n (D z ) n ǫ = (−1) n n!∆(n, 1)ǫ (4.24) one can show the soft part of (4.23) is satisfied provided one sets with T ǫ given by (4.25). A strategy to compare both sides of (4.26) is to bring the RHS into a form where ǫ appears with no derivatives. Substituting (4.25) in (4.22) and using the identity 13 (D z ) n−k (D z ) n ǫ = (−1) n n! (−1) k k! (D z ) k ∆(n, k + 1)ǫ, k = 0, . . . , n (4.27) one finds, after integration by parts on the sphere, (4.28) Expression (4.28) defines a "charge density" that should be compared with n ρ hard (u,x) as given in Eqs. (4.14) to (4.17). One can readily check that the k = 0 term in (4.28) agrees with the first term in (4.14). The comparison of the remaining k ≥ 1 terms in (4.28) with the second term in (4.14) is non-trivial due to the possibility of integration by parts in u and the interchanges ϕ ↔φ. This was already seen in the n = 1, 2 cases studied earlier, where the equality of charges required the use of several non-trivial identities. Unfortunately, we have not been able to find a general form of such identities that would allow us to establish the equality for arbitrary n.

Proof of conservation laws in classical theory
In the previous sections we showed that the sub-n soft theorems can be written as Ward identities of a tower of asymptotic charges and conjectured that these charges were the charges generated from n,n F ru (x) at I + and n,n F rv (x) at I − . Our conjecture was motivated by showing this equivalence in the case of leading, sub-leading and sub-subleading soft photon theorems.
However even assuming the conjecture to be valid, a natural question arises. Why do we expect this infinite hierarchy of infinite dimensional (asymptotic) charges to be conserved in classical theory? After all, the theory is not expected to be integrable. In this section we show that classically these charges are indeed conserved. That is, by analyzing the theory at spatial infinity, we show that 14 Our proof is an extension of the analysis done in [1] for the n = 0 case. We will summarize the key idea first and then provide a detailed analysis of these conservation laws.

Key ideas
n,n F ru (x) and n,n F rv (x) are charge densities localised on S 2 at u = −∞ and v = +∞ respectively. These two celestial spheres can also be understood as (future and past) boundaries of a three dimensional de Sitter space representing space-like infinity of Minkowski (or more generally asymptotically flat) spacetime in a particular compactification [17]. More in detail, when working with hyperbolic coordinates (τ, ρ, z,z), if we take ρ → ∞ while keeping τ, z,z fixed, we reach Lorentzian three dimensional de Sitter. We will refer to this (blow up of) spatial infinity as H. The advantage of working with this definition of spatial infinity is that it dovetails nicely with I ± . Thus the boundary spheres at τ = ±∞ are mapped onto S 2 at u = −∞, v = ∞ respectively.
In order to prove Eq. (5.1), we need to relate the fields at H with fields on the respective boundaries of I ± . By analyzing the equations of motion at spatial infinity and assuming appropriate fall-offs of the fields as u → −∞ and v → ∞ we relate fields at H with n,n F ru (x) and n,n F rv (x).
14 In this section F ab will stand for the real field strength rather than its self-dual part. The Ward identity associated to the negative (positive) sub-n soft theorem corresponds to the (anti) self-dual part of Eq. (5.1). The conservation statements can be understood as Eq. (5.1) for the real field strength, together with its Hodge dual version.
We will thus first analyze the equations of motion at spatial infinity in section 5.2 and in section 5.3, show how to relate n,n F ru (x) and n,n F rv (x) with certain data at H.

Maxwell equations at spatial infinity
We introduce hyperbolic coordinates in the region r > |t|, where ρ := √ x µ x µ = √ r 2 − t 2 and y = y α are coordinates on the unit hyperboloid H. The Minkowski line element takes the form with h αβ the unit hyperboloid metric. For concreteness we take y α = (τ, (5.5) Following [18], we assume a 1/ρ expansion of the Maxwell field near spatial infinity: For the scalar field, we assume that their fall-off behavior at spatial infinity is faster than O(ρ −m ) ∀ m. This assumption is motivated by the fact that we do not consider soft limit of the charged particles. Under this assumption the field equations reduce to source-free Maxwell equations at spatial infinity. The reason we can consistently do this is due to the fact that there is a clear separation between the source and electromagnetic radiation.
In the case of non-abelian gauge theories or perturbative gravity at second order there is no such separation and hence our analysis has to be generalized. Note that for massive charged particles the assumption can be proved since in a flat background massive fields decay exponentially fast at spatial infinity. Under this expansion, the source-free Maxwell equations take the form 15 where D α is the covariant derivative on H. For n = 0 the equations and their solutions are related to the leading soft photon theorem charges [1]. The idea is that the subleading 15 We remind the reader that in this section we are working with real (rather than self-dual) field strengths.
charges will be related to the remaining values of n. For n > 1 one can decouple the equations by combining the second equations of (5.8), (5.9) to eliminate n F αβ . Using the identity [D β , D α ]V β = 2V α one arrives at [20]: Once a solution to n F ρα is found, one can obtain n F αβ from the second equation in (5.9). Our next task is to specify τ → ∞ fall-offs. Due to the divergence-free condition, the α = A and α = τ fall-offs are not independent: If we assume When such asymptotics are inserted in the wave equation (5.10) one finds h = ±n. In the next subsection we show that our prescribed fall-offs at null infinity imply h = −n. Thus, at spatial infinity the problem to solve is (5.10) with the τ → ∞ asymptotic condition

Relating field expansions at null and spatial infinity
As shown in appendix B in the limit u → −∞, the falloff for the 1/r k+2 coefficient of F ru is given by, where the dots denote terms that fall off faster than any power of 1/u. We now use this expansion to obtain the τ → ∞ fall-offs at spatial infinity. Combining (2.6) with (5.14) we obtain the double sum expansion 16 We regard this as an expansion in the two small parameters: |u/r| ≪ 1, and |1/u| ≪ 1, (5.16) 16 In this section k,l F ru denotes the coefficients of the real field strength Fru rather than its self-dual part F + ru . The two have the same fall-off behavior.
For later comparison with the spatial infinity expansion, it will be convenient to rewrite (5.15) in a way that the l-sum appears first: We now express the two small parameters (5.17) in terms of the (ρ, τ ) coordinates. From the change of coordinates, one finds: Substituting these in (5.17) yields an expansion in the small parameters τ /ρ ≪ 1, and 1/τ 2 ≪ 1, (5.21) that can be related to the spatial infinity expansion. Note that the ρ dependence appears only in the l-sum. From the relation one finds the ρ dependence is the same as in the expansion (5.6): Furthermore, for a given 1/ρ power, the dominant τ → ∞ term is determined by the k = l summand. Specifically, one finds: as anticipated in Eq. (5.12). We thus conclude that where we included a + upperscript to indicate this is a τ → +∞ coefficient 17 . The analogous analysis at τ → −∞ yields 17 Not to be confused with the self-dual field. We apologize for the overlap of notation.
A similar analysis for the other components of the field strength (in the limit τ → ±∞) revels that Eqs. (5.25), (5.27) (and the analogues for past infinity) are the key equations that will allow us to relate the future and past charge densities. Notice that relation (5.13) becomes, upon using (5.25), (5.27), equivalent to relation (2.44) discussed in section 2.4. 18

Establishing the conservation
From the perspective of spatial infinity it is natural to work with n,0 F ρA (≡ n−1,n F rA ) as the charge density, since it represents 'free data' for the differential equation (5.10). The idea is to start with data at the future asymptotic boundary of H and evolve it backwards to the asymptotic past of H. One may object that this conservation is rather trivial: One is simply casting the free field equations in hyperbolic coordinates, but for free fields the conservation should be trivial! The non-triviality arises from the strong fall-off conditions imposed at future and past null infinities. As shown in the previous subsection, these select O(|τ | −n ) fall-offs on H which need not be satisfied for generic solutions to Eq. (5.10). It turns out that this specific large τ decay is related to another consequence of the strong fall-offs at null infinity: that the (vector) spherical harmonic decomposition of ( n−1,n F rA ) n,n F ru (x) starts at l = n. As shown below, the O(|τ | −n ) fall-off together with the l ≥ n spherical harmonic property imply that solutions to Eq. (5.10) satisfy specific parity conditions on H under inversion Y µ → −Y µ . This parity property then automatically implies the conservation law.
Our strategy is to find an explicit 'boundary to bulk' Green's function that solves (5.10) and satisfies all the aforementioned properties. A first candidate can be found by generalizing the Green's functions used to extend superrotations vector fields to timeinfinity [19]. This leads to the following family of Green's functions: where J µν α and L B µν are the angular momentum vector fields on H and S 2 respectively, With the help of identities given in appendix E.1 one can readily verify that (5.28) satisfies Eq. (5.10) (with respect to the y variable). It is also not difficult to verify that solutions constructed from (5.28) have the 'wrong' O(τ n ) fall-offs. 19 18 Recall in our conventions F ab = 1 2 (F + ab + F − ab ) and F − rz = 0 = F + rz . Adding (2.44) to its complex conjugated then leads to Eq. (5.13). 19 We expect the Green's functions (5.28) can be used to define smearing fields which would allow to extend the charges Qn[ǫ] to spatial infinity, as in the n = 0 case [1]. We leave for future work the study of smeared charges at spatial infinity. See also [20] for a set of closely related charges.

Conclusions and outlook
In this paper we have shown that for tree level scattering processes, there exists an infinite hierarchy of conservation laws such that at each level n of the hierarchy there is an infinite dimensional family of conserved charges Q n [ǫ n ] labeled by functions on the sphere ǫ n (x). 22 20 The ansatz (5.31) was also inspired by the integral representation of solutions to free Maxwell equations used by Herdegen, see e.g. [21]. 21 The operator ∆(n, 1) acting on covectors is given as in the definition for scalars, Eqs. For n = 0, 1 these charges are the well known asymptotic charges associated to leading and sub-leading soft photon theorems. Our analysis divorced the derivation of asymptotic charges from the paradigm of large gauge transformations. In fact, as argued in [16], we do not expect the charges corresponding to levels n ≥ 2 to be associated to any large gauge transformations. These are intrinsically boundary charges and are purely a consequence of the fall-off behavior of the radiative data and equations of motion of the theory at infinity. A natural question is, what kind of factorization theorems arise by considering Ward identities associated to charges with n ≥ 2.
We showed that the Ward identities associated to n = 2 charges are equivalent to the sub-subleading soft photon theorem and conjectured that this equivalence holds ∀ n. That is, Ward identities associated to Q n [ǫ n ] are equivalent to sub-n soft photon theorems. Note that these sub-n soft theorems only capture the factorized piece of the ω n−1 coefficient in the ω → 0 expansion of the amplitude. For n ≥ 2 there are non-factorized terms that are not seen by these charges.
Beyond the obvious open issue of proving the aforementioned conjecture, there are several open questions which emerge out of this work and which have not been addressed in this paper. We expect that our analysis can easily be generalized to tree-level amplitudes in perturbative gravity. In fact, an interesting puzzle in this direction concerns the softexactness of tree level MHV amplitudes in pure gravity [22]. That is, in the case of Gravity MHV amplitudes, the sub-n soft limit exactly factorizes and there is no remainder R n . Whence we expect the asymptotic charges to constrain the MHV amplitudes completely as opposed to only in the infrared sector.
A related question concerns the algebra of charges Q n [ǫ n ]. If it does form an algebra then the quantization of this infinite dimensional algebra could shed interesting light in the structure of gauge theories.
Once loop corrections are taken into account, the story changes completely beyond the leading order in soft limit. As shown in [13], the sub-leading soft photon theorems receives corrections at one loop and are associated to terms proportional to ln ω as opposed to ω 0 . Such a soft expansion clearly implies that the large u fall-off behavior that we have assumed for tree level scattering data breaks down and one needs to consider a slower fall-off where the radiative field contains 1 u terms as |u| → ∞ [23]. Whether such soft theorems with logarithmic corrections can be derived from asymptotic charges remains to be seen.

Acknowledgements
We are grateful to Ashoke Sen for raising the issue of proving classical conservation law associated to subleading soft photon theorem which led us to this inquiry as well as for discussions on infinity of soft theorems. We are also grateful to Amitabh Virmani for urging us to look at the relationship between asymptotic charges and Newman-Penrose charges. AL would like to thank IISER, Pune and ICTS, Bangalore for their hospitality during course of this project. MC would like to thank the Center for Theoretical Physics at Columbia University for hospitality during the final stages of this project.
A Maxwell equations for self-dual fields at I + In retarded coordinates, the components of the 4-volume form take the form with ǫ AB the unit-sphere area form. The duality operation for the various components reads:F The idea is to use Eqs. (2.5) to arrive at an equation for F + ru in terms of the source. We then consider Eqs. (2.1), (2.2), (2.3) with F ab replaced by F + ab . From (2.3) and the self-duality condition one arrives at: Next we combine the sphere divergences of Eqs. (2.3) and (A.4) so as to keep only the j z part of the current. Since in holomorphic coordinates the area form is given by the appropriate combination is: where we used that D A D B F + AB = 0. Finally, using Eq. (2.1) one can write F + rA appearing in (A.6) in terms of j r and F + ru . The resulting equation only involves F + ru and the current components j z and j r . Such equation can then be expanded in 1/r to yield a recursive relation for the F + ru components. Assuming the standard 1/r expansion, (to simplify notation we omit the + superscript in the coefficients for F + ru ) one finds, 2(n + 1)∂ u n+1 F ru + (∆ + n(n + 1)) n F ru = −2D z n j z + 2∂ u n j r + (n + 1) where A A (u,x) is the leading term of the sphere components of the vector potential and we assumed A u = O(r −1 ) so that F uA (r, u,x) = ∂ u A A (u,x) + O(r −1 ). In fact by using the leading equations in (2.2) and (A.2) it can be easily shown that, where A A (u,x) is the leading term of the sphere components of the vector potential and we have assumed that The ω → 0 expansion of the tree-level amplitude (3.1) implies the radiative free data for the photon field has a similar expansion in frequency space: In real space this translates into an expansion of the type with θ the step function. Thus, the u → ±∞ falloffs for the null infinity free data are of the type where the dots are terms that go to zero faster than any power of 1/u. For the matter field φ(u,x) we assume a similar expansion but with vanishing u → ±∞ limit (corresponding to the fact that there are no 'soft' charge particles). Equations (2.9) and (2.10) imply then the following u → −∞ falloff for the 1/r k+2 coefficient of F ru : where the dots denote terms that fall-off faster than any power of 1/u.

C Asymptotics of charged field and current at null infinity
For the massless scalar field we assume the standard 1/r expansion at null infinity: This implies a 1/r expansion for the current j a = ieϕ∂ aφ + c.c., (n − 2k + 1) The actual current sourcing the field strength has an extra term −2e 2 A µ ϕφ. However such term is not seen in the tree-level formula we are interested here. The scalar field coefficients in the expansion (C.1). Should be determined by the field equations. For the tree-level expression of interest here, it suffices to consider the free-field equation ϕ = 0, which takes the form The lowest components of the current coefficients are where φ ≡ 0 ϕ is the scalar field 'free data' and We start with the expression for the differential operator where T ≡ T n+1 w . . . w . The conventions in this appendix are slightly different than in the rest of the paper. The differential operator in the main text, Eq. (3.17), is obtained from (D.7) by doing the replacement n → n − 1 and including a multiplicative factor (−i) n+1 2 . We will work here with the form given in (D.7) and make the above replacements at the end of the calculation.
n F ρA can be treated along similar lines. Here a difficulty stems from the need to work with tensorial indices. A way to deal with this problem is to decompose the vector field m V A into grad plus curl pieces, A ∂ B g and find scalar expressions involving f and g. This allows one to use the formula (E.23) (valid for scalars) to arrive at the desired asymptotic form as given in Eq. (5.32).