Subleading soft dressings of asymptotic states in QED and perturbative quantum gravity

We construct Faddeev-Kulish states in QED and perturbative quantum gravity to subleading order in the soft momentum expansion and to first order in the coupling constant, using the charge conservation formula of asymptotic symmetries associated with the tree-level subleading soft theorems. We demonstrate that the emission and absorption of soft photons/gravitons in dressed amplitudes vanish. The fact that no additional soft radiation may be added to a dressed amplitude supports the claim that, in the dressed state formalism, the soft and hard sectors of scattering processes are correlated. We also show that the dressed virtual amplitudes are equivalent to the infrared-finite part of the traditional amplitudes constructed using Fock states. Since there is no real soft radiation in the asymptotic Hilbert space, the dressed state formalism gives the same cross sections as the Bloch-Nordsieck method.


Introduction
It is well known that in the conventional formulation of quantum field theories, S-matrix elements are plagued by infrared divergences and are therefore ill-defined. This is due to the fact that interactions in such theories cannot be turned off at large times, which invalidates the use of Fock states as asymptotic states. In 1970, Kulish and Faddeev [1] showed for QED that the correct asymptotic Hilbert space which accounts for the non-vanishing interactions is the set of "dressed" states, which are essentially coherent states of soft photons. In 2013, an analogous construction has been made for perturbative quantum gravity by [2]. Such coherent states are referred to as Faddeev-Kulish (FK) states, or dressed states. Indeed, it was shown by Chung [3] for QED and Ware et al. [2] for gravity that scattering amplitudes in the FK state basis do not exhibit infrared divergence. Throughout the paper, we will refer to such amplitudes as FK amplitudes, as opposed to the traditional amplitudes constructed using Fock states.
The FK states are closely related to asymptotic symmetries of gauge theories. In recent years, asymptotic symmetries have received a considerable amount of attention due to their subtle connection with soft theorems . The soft charges associated with the asymptotic symmetries induce a degeneracy in the vacuum state of the theory, and charge conservation leads to a transition among the vacua in a scattering process. It has been shown that the FK states implement this vacuum transition properly, while Fock states do not [34][35][36][37]; FK dressings are operators that carry a definite leading soft charge. In this context, vanishing S-matrix elements in Fock state basis due to infrared divergence can be regarded as a consequence of a superselection rule of the degenerate vacua. This relationship between asymptotic symmetries and FK dressings have been shown to be generalizable to other asymptotic boundaries of the spacetime, such as Rindler and black hole horizons [38][39][40].
Although the FK state approach has proven to be useful, there are still unanswered questions. For example, in the traditional construction of amplitudes, one deals with infrared divergences by summing over inclusive cross sections [41], leaving only the contribution from the infrared-finite part of the amplitude [42]. The FK amplitude is free of infrared divergence, but does it agree with this infrared-finite part? For another example, given a set of dressed states as asymptotic states, may one add arbitrary soft radiation to the incoming state? Also, do the soft particles carry information about the hard particles? To answer these questions, the known FK states, which only contain terms at the leading order in the soft expansion, are not enough; one needs to extend the construction to the subleading order.
Motivated by these questions, in this paper we derive the FK dressings of QED and gravity to first order in the coupling and subleading order in the soft expansion. The known leading soft FK states have been shown to be charge eigenstates of asymptotic symmetry associated with the leading soft theorem [35,37]. Therefore, one would expect that a similar construction using asymptotic symmetries associated with the subleading soft theorem will correctly lead to the subleading dressing. However, this presents several challenges. The subleading soft charges do not commute with each other as well as with the leading soft charges, so one cannot construct simultaneous eigenstates that carry a definite leading and subleading charge. But as discussed in section 2, to first order in the coupling we may finesse this problem as far as S-matrix elements are concerned, and show that the construction is self consistent. The subleading soft dressing operator so obtained agrees with the Wilson line dressing first derived by Mandelstam [43] to first order in the coupling in the exponent. Thus to this order in the coupling and in the soft expansion, there is again the equivalence between Wilson line dressings and FK dressings discussed in [40,44]. With the subleading dressings at hand, we proceed to derive two main results, which we outline in the next two paragraphs.
First, we show that an external photon/graviton insertion to an FK amplitude vanishes in the soft limit. It follows that once Fock states are dressed to yield a set of FK states, there is no extra soft radiation in the asymptotic Hilbert space. That is, any soft photons/gravitons may only appear through FK dressings. We anticipate that this result may shed some light on the discussions regarding factorization of the soft and hard sectors, and whether the soft particles carry information of the hard degrees of freedom; see [45][46][47][48] for example. We can say that the soft and hard sectors are correlated for the FK states since specifying one fixes the other. Also, since a soft photon/graviton creation operator acting on a vacuum makes it orthogonal to all states in the asymptotic Hilbert space, one may assume that it annihilates the vacuum as far as scattering processes are concerned. This removes some complications regarding Fock vacuum and soft operators, see [34] for example.
Second, we show that the virtual FK amplitude is equivalent to the infrared-finite part of the traditional amplitude in the Fock basis, up to power-law type corrections of the soft energy scale which is negligible by definition. This is rather reassuring; if this were not the case, FK amplitudes would be in conflict with experimentally measurable cross sections, which involve the infrared-finite part of traditional amplitudes. Since there is no real soft radiation in the asymptotic Hilbert space, the FK state formalism yields the same cross sections as the Bloch-Nordsieck method [41]. This serves as a sanity check for the argument that the asymptotic states are FK states instead of Fock states.
Two important comments are in order. First, it has been shown that the subleading soft theorems receive one-loop exact quantum corrections [49] that are logarithmically divergent in the soft energy. Since we construct the subleading dressing only to first order in the coupling constant, we will not consider such corrections and work with the tree-level subleading soft theorems. While we leave details of the loop-corrected dressings for future investigation, we argue that our second result, that FK amplitudes are in agreement with infrared-finite part of traditional amplitudes, should remain unaffected by loop corrections. This follows from the fact that logarithmic divergence is an integrable singularity, and consequently the additional contributions of loop corrections to FK amplitudes integrate to negligible quantities. Secondly, there is a piece of the subleading soft contributions to scattering amplitudes which does not exponentiate. This is just the contribution usually assigned to Low's theorem [50][51][52]. This complication, again, may be ignored at leading order in the coupling constant in the exponent of the dressing operator.
The paper is organized as follows. In section 2 we begin with a brief review of superrotation and the associated charge, followed by the construction of the gravitational subleading FK dressing. Section 3 presents an analogous construction for the subleading FK dressing in QED. Using the dressings, in section 4 we show that the insertion of a soft graviton makes the amplitude vanish, and in section 5 we demonstrate the equivalence between FK amplitudes and the finite part of traditional amplitudes. The dressings of gravity and QED are very similar; although sections 4 and 5 work explicitly with gravitons, we expect the conclusions to be valid for QED as well. In appendix A we review the correspondence between subleading soft theorems and the Ward identities of the associated asymptotic symmetries, both for scalars and for particles with spin.

Review of superrotation in asymptotically flat spacetime
We start by establishing our notation regarding asymptotically flat spacetimes and reviewing the materials associated with the superrotation on I ± . We follow the construction of [7] closely. For the sake of simplicity we will assume that all matter particles are massless scalars; for particles with spin see appendix A.1.2.

Metric and mode expansions
In Bondi coordinates, the metric for an asymptotically flat spacetime near the future null infinity I + reads [7] where the first line corresponds to the flat metric. Here z = e iφ tan θ 2 is the stereographic coordinate, γ zz = 2 (1+zz) 2 is the metric on the 2-sphere and m + B is the Bondi mass aspect. The uz-component of the metric has the expansion where D z is the covariant derivative on S 2 and N + z is the angular momentum aspect. In general, m + B , N + z and C zz are functions of u, z andz. Let us define the graviton field h µν through where η µν = diag(−1, 1, 1, 1). Near I + , the graviton field can be approximated by the on-shell mode expansion where ω k = k 0 = |k|, and ± µν = ± µ ± ν are the spin-2 polarization tensors. By parametrizing the graviton momentum k µ by (ω k , z,z), we may write the polarization tensors as The out-operators satisfy the standard commutation relation, Near the past null infinity I − , the asymptotically flat metric reads We have the mode expansion for the incoming graviton field where the in-operators satisfy

Superrotation charge
Superrotation near I ± are generated by the following vector fields respectively [7], (2.13) parametrized by the same vector Y z on the sphere. We will drop the constraint that Y z is a conformal Killing vector, following [20]. The conserved charges associated with superrotation have the expressions [53] where I + − (I − + ) is the past (future) boundary of the future (past) null infinity. The charges can be decomposed into soft and hard parts, where the soft charges are given by 19) with N zz = ∂ u C zz and M zz = ∂ v D zz ; we refer to [7] for details. Here x z denotes a unit 3-vector whose direction is given by (z,z), The hard charges Q ± H are defined by their actions on the Fock states, where the momentum of the massless scalar p i is written as

Distinction between in and out operators
Notice that in (2.17) and (2.19) we follow the notation of [4] and others to distinguish between the "out" operators on I + and the "in" operators on I − . They are related by the boundary condition such that the subleading soft contribution to the S-matrix element from insertions of a in † s in the incoming state and −a out −s in the outgoing state are equivalent, see [7] for a discussion. An alternative approach, which was taken in [37], is to make this relation explicit by taking in (2.19) and (2.26), Then we can remove the superscript in Q ± S and write Either approach will yield the same result; with the distinction intact, one just has to be cautious when contracting an "out" and an "in" operator. We will employ the latter convention (2.24) and use (2.26) for it removes some complications in the amplitude computation.

Subleading soft dressing
Let us consider the scattering from some incoming state |in to some outgoing state out|the states can be either dressed states or Fock states. Superrotation symmetry states that the associated charge must be conserved in a scattering process, i.e.
where S is the scattering matrix. This can be written as where we used (2.21) and (2.22). Let us choose the vector field for which (2.28) becomes [7,20] (see appendix A.1 for a derivation) where k µ ≡ (ω, ωx z ), and S (1)− g is the subleading soft factor for negative-helicity graviton, with η i = +1 for incoming particles and η i = −1 for outgoing particles. 1 Using the identity We now claim that under certain circumstances stated below, we may write for a ket vacuum |0 , Strictly speaking, the subleading soft charge does not annihilate the vacuum state (and hence the symbol ≈), but rather creates a state containing a soft graviton. In section 4 we will show that no state may scatter to such a state and vice versa in the dressed state formalism. Therefore, Q S may be taken to annihilate |0 insofar as scattering processes are concerned.
As was done in [37] for the leading soft dressings, we aim to construct the subleading soft dressed state by using superrotation charge conservation. Since Q S (Y g ) is of the form a − a † , we want to consider a coherent state of the form we deviate from the convention used in [7]. The signs ηi derive from the different momentum-space representations of the action of angular momentum on bras and kets. Alternatively, one could define J λν i as in [7], for which case the statement of angular momentum conservation becomes i Ji = 0. We will adopt (2.32), in order for angular momentum conservation to take the more natural form i∈in Ji = j∈out Jj. The two approaches are equivalent.
where N µν in is a tensor whose components are to be determined by charge conservation, φ(ω k ) is an infrared function [1,2] that has support only in a small neighborhood of ω k = 0 satisfying φ(0) = 1. Here |p 1 , · · · , p m is an m-particle Fock state, where b † (p) is the creation operator of the scalar field. However, since we are using the tree-level subleading soft theorem, we can only construct a dressing that can be trusted to order κ in the exponent. In this spirit, let us define the incoming state as the linearized version of (2.35), By direct computation, where in the last line we used the notation N in · − ≡ N µν in − µν . Similarly, we can construct a bra state, With these states, we may write where in the second equality we employ a convenient abuse of notation to write N · − both as an operator and as its action on the amplitude in the momentum-basis. In section 5 we will see that, due to the presence of φ(ω k ), adding or removing subleading dressings do not change the value of the amplitude, that is, This shows that the non-commutativity of subleading charges and the nonexistence of simultaneous eigenstates do not cause difficulties to S-matrix calculations at this order. Using (2.47), one may write (2.46) as Thus, the superrotation charge conservation (2.31) reads One can derive a similar relation associated with If we treat supertranslation (associated with simple poles) separately as in [37], we may assume that N µν in,out do not possess poles. Then it follows that where W g (p 1 , · · · , p n ) is the gravitational n-particle dressing, which to the subleading order in soft momentum expansion and leading order in κ is given by which, keeping in mind that only order κ terms can be trusted, may be conveniently written as 2 (2.56) Here we employed the notation a µν (k) = s s * µν (k)a s (k), where s spans all polarizations. This includes unphysical polarizations, since the projection to physical polarizations in (2.32) is a consequence of our choice (2.29) of Y ; superrotation charge should be conserved for a generic vector field. Unphysical polarizations are also required for canceling out infrared divergence, see [2,36] for example. The expression (2.56) expanded to first order in κ agrees with the gravitational Wilson line dressing of Mandelstam [43]. Thus to this order in the coupling, one observes the equivalence between Wilson lines and FK dressings as discussed in [40,44].
For explicit calculations, it is convenient to define the infrared function as where λ is the infrared cutoff which we take to be zero at the very end of the calculation, and Λ is a very small energy scale below which particles are considered to be soft. With this definition, the dressing (2.56) becomes (2.58) The dressing W g acting on an n-particle Fock state will be understood as an n-particle dressing with the corresponding momenta of the hard particles, unless explicitly stated otherwise. The dressed states automatically implement conservation of supertranslation charge [36] and superrotation charge, as shown above. An important point concerning the validity of (2.58) should be emphasized here. In any scattering process there are contributions from real emissions and from virtual diagrams. The applicability of the subleading soft graviton and soft photon theorems to a 2 → 2 scattering process has been studied in [54]. There it was shown that the subleading soft photon theorem correctly reproduces the scattering amplitude to subleading order both for real and virtual photon emissions. However, for the case of soft gravitons, the subleading soft graviton theorem correctly reproduces the real external emissions, but there is a violation of the theorem for virtual gravitons. Thus, for the choice of vector field in (2.29), we expect our dressing to be correct for the case of real emissions which we discuss in section 4. For scattering involving virtual gravitons we would need to generalize (2.58).

Dressing in QED
We present here the construction of subleading soft dressing in QED, which is fairly parallel to the case of gravity in section 2. In this section we will mostly follow the notation and conventions of [12].

Mode expansions and conserved charges
The photon field near I + becomes nearly free and can be approximated by where the out-operators satisfy the commutator Likewise, near I − we have the incoming photon field, with the standard commutator The QED analog of superrotation is the asymptotic symmetry on I ± associated with Low's theorem. In terms of the mode expansions (3.1) and (3.3), the corresponding conserved charges are 3 [12] where the soft parts are given by Here x z is a unit 3-vector whose direction is given by (z,z); its Cartesian components are given in (2.20). The hard parts are defined by their actions on the Fock states, where A ∈ {z,z}, Q i is the electric charge of the i-th particle, and L Y is the Lie derivative on S 2 , see appendix A.2 for details. Scattering processes conserve Q Y , which implies that for any asymptotic states |in and out|.
As was the case in gravity, the contribution to the subleading soft matrix element from a soft photon insertion a in † α in the incoming state is equivalent to the contribution from an insertion −a out −β in the outgoing state. Therefore, we can follow the procedure of section 2.2 and define (3.11) where the soft operators are related by Q + We emphasize that this procedure is done to avoid defining separate rules for contractions between operators on I + and I − ; one may obtain the same result with the distinction intact.

Subleading soft dressing
Now we construct the subleading soft Faddeev-Kulish dressing in QED to leading order in the coupling constant e as linearized coherent states that respect Q Y charge conservation (3.10). Using (3.5)-(3.9), (3.10) can be written as (3.14) Then, (3.13) takes the form [12] (see appendix A.2 for details) where k µ ≡ (ω, ωx z ), and S (1)− e is the subleading soft factor for negative-helicity photon, with η i = +1 (−1) for outgoing (incoming) particle (by introducing η i we deviate from the convention of [12], see footnote 1). An analogous expression involving + can be derived by choosing , one obtains the following expression for the soft charge, Let us begin by considering a vacuum |0 such that, As noted in section 2.3, formally the subleading soft charge does not annihilate the vacuum, but rather adds to it a soft photon. As will be shown in section 4, in scattering processes such a state completely factors out, and therefore in S-matrix computations one may act as if Q S annihilates the vacuum (hence the symbol ≈). Now we consider states which are dressed to first order in e that take the form 20) where N µ in,out are operators to be determined, and φ(ω k ) is the infrared function that restricts the momentum integrals to soft modes. Then, and Along the same line of reasoning as in (2.48), this leads to Assume that the simple poles (associated with large gauge symmetry) have been treated separately, as in [35]. Then lim ω→0 (1 + ω∂ ω )N in,out · − = N in,out · − , and the equation of charge conservation (3.15) becomes if the matrix element out|S|in is to not vanish. A natural splitting for the dressing is Combining this with the known leading soft dressing, one deduces the dressed asymptotic state for QED to subleading order in the soft expansion and to first order in e to be p 1 , · · · , p n = W e |p 1 , · · · , p n , (3.29) where the dressing W e is which, keeping in mind that only terms to first order in e may be trusted, can be conveniently put as (3.31) The term O(e 2 ) emphasizes that the subleading dressing is valid only to order e. The photon operator is defined as a µ (k) = α α * µ (k)a α (k), where α spans all polarizations including unphysical ones, since the projection to physical polarizations is due to our choice (3.14) of Y ; the charge should be conserved for a generic vector field. The unphysical polarizations are also required to cancel out Weinberg's infrared-divergent factor. With the explicit implementation (2.57) of φ(ω k ), we obtain where λ is the infrared cutoff and Λ is the separation scale below which we consider particles to be soft. Notice that the structure is very similar to the gravitational dressing (2.58).
One can obtain the QED dressing from gravity by the replacement κ 2 (p µ i s µν ) → eQ i s ν for each particle. The comments about non-commutativity of subleading and leading soft charges with all its complications discussed towards the end of section (2.3) also apply here. In the subsequent sections, we will be working with the gravitational dressing with the understanding that same results can be shown for QED with minimal modifications.

External soft gravitons and photons
With the leading soft Faddeev-Kulish states, it is known that adding an external soft graviton does not induce infrared divergence; the divergent soft factors from the dressings cancel those from external legs [34,36]. Now that we have constructed dressings to subleading order in the soft expansion, we are in a position to investigate what happens to the O(ω 0 ) subleading soft factors. Although we will be working with gravitons, the derivation for QED is very similar and the final result is also valid for external soft photons. Every step in the calculation can be changed to the corresponding expression for QED by replacing κ 2 (p µ i s µν ) with eQ i s ν for each hard particle. As in (2.57), let Λ be the soft energy scale, below which particles are considered to be soft, and let λ be the infrared cutoff which we take to be zero at the end of calculations. Let us consider a scattering amplitude from the dressed m-particle state in to the dressed n-particle state out , with a soft graviton insertion of polarization s and momentum k µ = (ω, ωx) where ω is soft, The dressed amplitude M has the small-ω expansion where each M (n) is independent of ω. The different contributions to M(k, s; {p}) are illustrated in figure 1. It is known that the first term involving the infrared-divergent amplitude M (−1) vanishes [36]. To see this, note that M (−1) receives contribution from diagrams 1(a), 1(b), 1(d) and 1(e). Using the notation η i = +1 (−1) if i is outgoing (incoming), we may write where M ≡ out|S|in = out|W † g SW g |in is the dressed amplitude without graviton insertion. In the second equation we used (2.57) to write φ(ω) = 1. The first sum in the square brackets comes from graviton emission from external legs (figures 1(a) and 1(d)); the second sum comes from graviton emission from dressings (figures 1(b) and 1(e)).
We can determine M (0) by collecting the O(ω 0 ) terms in the amplitude. To do so, let us first decompose the gravitational dressing W g into leading and subleading parts, Then, we obtain the commutators The first and second terms in the summands correspond respectively to the leading and subleading soft contributions from figures 1(b) and 1(e). The second terms contribute to M (0) , along with emissions from internal propagators ( figure 1(c)) and external legs (figures 1(a) and 1(d)). There is one subtlety here -the second terms in the summands are missing the subleading dressing factors W (1) g . However, as we will see in section 5, insertion of such factors only add O(Λ) corrections to the amplitude, which is negligible by the definition of the soft energy scale Λ. Therefore, within the amplitudes one may replace W (0) g of (4.7) and (4.8) with W g and write (4.9) since φ(ω) = 1 for ω < Λ. We remind the reader that the sign η i in the subleading soft factor comes from the different momentum-space representations of the action of J µν i on bras and kets: At this point one may remove the infrared regulator λ → 0, and conclude that the soft emission amplitude is negligible since ω is by definition less than the soft energy scale Λ, which in turn is by definition much less than any energy scales of our interest. As the emission amplitude vanishes in the soft limit, the state containing a zero-energy graviton can be treated as null as far as scattering processes are concerned: In summary, the use of leading and subleading Faddeev-Kulish states do not allow absorption and emission of on-shell soft gravitons at tree level.

Equivalence of Faddeev-Kulish amplitudes and traditional amplitudes
In this section, we show that the Faddeev-Kulish amplitude is equivalent to the infraredfinite part of traditional amplitudes constructed using Fock states, up to power-law type corrections in the soft energy scale Λ which is negligible by definition. Keeping both leading and subleading terms to first order in κ in the exponent of the dressing function W g , we explicitly show that this equivalence is up to order Λ for radiation-less amplitudes in the case of scattering of a scalar from an external potential. In reference [36] only the infraredfiniteness was shown, keeping just the leading order term in the exponent of the dressing. Again, although we will only derive the result explicitly for gravity, the derivation for QED is similar and the result also holds for QED amplitudes. For simplicity we will consider a 1 → 1 gravitational potential scattering between the dressed states p i and p f at one-loop order. Let us define the shorthand notation and similarly P µν f and Q µν f corresponding to p f . We also use the notation where λ is the infrared cutoff and Λ is the soft energy scale. We can use these and (2.54) to write the dressed states as 3) where concatenation implies contraction, for example, First consider the contribution to the matrix element due to graviton exchange between dressings. There are two self-interactions of the dressings, each coming from W g (p i ) and W † g (p f ) (figures 2(b) and 2(c)), and one cross-interaction between W g (p i ) and W † g (p f ) ( figure 2(a)). Using the Baker-Campbell-Hausdorff (BCH) formula, we may write, Figure 3: Contributions to the amplitude due to graviton exchange between a dressing and either an external leg or the internal diagram.
and similar for the incoming dressing, where we used (see [2,55] for example) and employed the notation P f IP f ≡ P µν f I µνρσ P ρσ f . Then to first order in κ, 4 The last terms on the RHS of (5.10) and (5.11) are the self-interaction contributions of the dressings to the matrix element, The cross-interaction between the two dressings (figure 2(a)) come from the contraction between the graviton operators of (5.10) and (5.11), which introduces the term Now, let us consider the remaining contributions, namely the ones due to a dressing exchanging gravitons with either an external leg or an internal propagator, see figure 3. At the leading soft (divergent) order, graviton exchange with an internal propagator does not contribute, while each exchange with an external leg (figures 3(a)-(d)) induces a soft factor ±ηP µν , where the + (−) sign corresponds to emission (absorption) of the graviton, and η = 1 (−1) if the external leg is outgoing (incoming). At the next order, graviton exchanges with the internal diagram (figures 3(e) and 3(f)) induce one subleading soft Figure 4: Graviton exchanges that do not involve dressings. The soft graviton loops exponentiate and factor out [42].
factor ηQ µν +Q µν for each external leg, whereQ µν = O(ω 0 ) is a subleading soft factor due to the graviton being off-shell; see [54] for example. Since there are two dressings and two external legs, at one-loop level this introduces eight terms, four from W † (p f ) which correspond to diagrams 3(c), 3(d) and 3(f), and four from W (p i ) corresponding to figures 3(a), 3(b) and 3(e), Expressions (5.12), (5.13), (5.14) and (5.15) comprise the full contribution of the graviton interaction to the matrix element that involves Faddeev-Kulish dressings at one-loop. The net contribution to subleading order in soft momentum is Now, the traditional amplitude p f |S|p i is infrared divergent, owing to the presence of soft virtual graviton loops. The divergence can be factored out, such that at one loop [36,42], where p f |S|p i is the infrared-finite matrix element where all virtual graviton loop momenta below the soft energy scale Λ has been removed. In figure 4, only diagrams 4(a)-(c) contribute to this term. We will not concern ourselves with subleading corrections to this factor coming from diagrams 4(a)-(e), since such corrections do not alter our conclusion.
Putting (5.16) and (5.17) together, at one-loop we observe that Recall that the integral spans only the soft sector λ < ω k < Λ. Since the integrand is of order O(ω 0 k ), after removing the infrared regulator λ the second term on the RHS of (5.18) becomes By definition, Λ is the soft scale and is therefore negligible compared to all other energy scales of significance. Therefore, That is, the Faddeev-Kulish amplitude dressed to subleading soft order is equivalent to the infrared-finite traditional matrix element. This extends the results of [2,36] where the analyses were done only at the level of infrared-divergent terms. Now, recall that our construction of dressings does not account for loop-corrections to the subleading soft theorems. As one can see from (2.31) and (3.15), such corrections will introduce ln ω terms to the integrand in the exponent of the dressings. Then, from (5.19) we observe that the corresponding corrections to dressed amplitudes will involve infrared-finite integrals of the form lim λ→0 Λ λ dω k ln ω k = O(Λ ln Λ), (5.21) which becomes arbitrarily small as we decrease Λ, and therefore is negligible compared to other energy scales of our interest. It follows that the dressed amplitudes remain unaffected by loop corrections to the subleading soft theorem. One should check that the equivalence (5.20) reproduces the inclusive cross section obtained by the Bloch-Nordsieck method [41], which sums over emissions of real gravitons with energies below the detector resolution E res . It was shown in [42] that the cross section for a transition α → β accompanied by any number of real gravitons with total energy below E res can be written as where Γ α→β (ω 1 , · · · , ω N ) denotes the cross section for emission of N gravitons with energies ω 1 , · · · , ω N . The cross sections are now the norm-squared of the dressed amplitudes. In the previous section, we have shown in (4.12) that scattering amplitudes, and therefore cross sections, for processes that emit/absorb real gravitons with energy below the soft scale Λ are negligible. This has the effect of replacing the infrared cutoff λ in (5.22) with the soft energy scale Λ, which results in [42] Γ Γ α→β is the cross section for α → β without the undetectable real gravitons. In the original construction of [42] with Fock states, one factors out the soft loop contribution from Γ α→β . In our construction with the dressed amplitude (5.20), we have no soft graviton loops (they have been canceled on account of the dressings) and thus may write where Γ 0 α→β is the cross section computed excluding virtual graviton loop momenta below the soft scale Λ. Therefore, which agrees with the inclusive cross section computed in [42] using Fock states via the Bloch-Nordsieck method.

Discussions
Using ideas similar to the one presented in [37], at leading order in the coupling we constructed the Faddeev-Kulish dressing for gravity and QED to subleading order in the soft energy expansion. We have shown that the dressed amplitudes are equivalent to the infrared-finite part of the traditional matrix elements, up to negligible power-law type corrections of the soft energy scale. We have also shown that, to first-order in the coupling constant, the FK state formalism does not allow soft radiation of photon and graviton. This supports the proposition that for the FK states soft particles carry information about the hard particles and vice versa [47,48]. An important question that we have not addressed in this paper is the characterization of a general dressed state up to and including all subleading soft corrections. Since the subleading charges do not commute with each other or with the leading charge, it is not clear how to construct dressed states beyond the order κ term in the exponent of the dressing operator. Another manifestation of this in perturbation theory is that the subleading corrections do not exponentiate [54,56]. This connection between the non-Abelian nature of the asymptotic symmetry charges and the violation of factorization of soft quanta has not been noticed before. Indeed, in non-Abelian gauge theories there is a breakdown of exponentiation of soft gluon effects for certain processes even at the leading soft order. Further studies of their connection may enhance our understanding of both of these problems. See reference [57] for a recent work on the leading soft FK dressings and asymptotic symmetries of perturbative QCD.
We end with some comments about the loop corrections. The subleading theorem is known to receive ln ω corrections at one-loop level, see [31,49,[58][59][60] for example. Since we have not considered such corrections in our construction, this is another reason why the subleading soft dressings we derived may be trusted only to first order in the coupling constants. We have argued that the agreement of FK amplitudes with traditional amplitudes is not affected by the correction. However, the non-existence of soft radiation presented in section 4 is a tree-level statement. It would be interesting to study how loop corrections affect this result.
For the vector field Y z , let us choose We now show that the factor in the parentheses on the RHS is the subleading soft factor, Since all hard particles are assumed to be scalars, the angular momentum consists of only the orbital part, where η i = +1 (−1) if i-th particle is outgoing (incoming). Our definition of J µν is different from that of [7] by a sign, see footnote 1. Now let us parametrize and −µν = −µ −ν . The quantities (E, z,z) are related to p µ by Let us write (A.10) out as (A.14) One can show using (A.12) that Now, noting that Using these with (A.15) and (A.16), we may write (A.14) as This is exactly the expression appearing on the RHS of (A.9), which was to be shown.

A.1.2 Spin correction
Now suppose the hard particles have non-zero spin. The angular momentum J µν appearing in the subleading soft factor now contains the spin contribution S µν , Let us define helicity in terms of the Pauli-Lubansky pseudovector, where µνρσ is the Levi-Civita tensor with 0123 = 1, and in the second equation the orbital part drops out due to antisymmetry in p. In this basis, the spin angular momentum has components [61] where in the second equation we used the parametrization (A.12). For particles with spin, the action of the hard charge is [7] out|[Q H , S]|in = i where L Y is the Lie derivative on the 2-sphere with respect to Y , where S AB is the pullback of (A.25) to the 2-sphere. By coordinate transformation from x µ to (z,z), one finds that (A. 28) and similarly S zz = ihγ zz , S zz = Szz = 0. Thus, and (A.26) can be written as In the previous section, we saw that with the choice Y = (z − w) 2 (z −w) −1 ∂ z , the first two terms in the square brackets correspond to the subleading soft factor coming from the orbital angular momentum L µν . The third term then must correspond to the factor coming from spin angular momentum S µν , which is Combining the results, we obtain But from (A.8) we observe that this can be written as which is exactly the third term in the square brackets of (A. 30), showing that the formalism extends to particles with spin.

A.2 QED
The results of appendix A.1.2 is directly relevant to QED since charged particles have spin. The construction is very similar -we will start with charged scalars and move on to spin corrections.

A.2.1 Massless scalars
From the action of hard and soft charges for charged scalars, we have [12] out|[Q H , S]|in = i (1 + ω∂ ω ) out|a − (ωx z )S|in We now show that the factor on the RHS is the subleading soft factor for scalars, Comparing this to gravity (A.10), we observe that this is just S which exactly agrees with the RHS of (A.39).

A.2.2 Spin correction
In the presence of spin, the subleading factor (A.41) gains an additional spin contribution, which can be obtained using the results of appendix A.1.2. We shall show that this contribution is exactly the spin corrected term in the action of hard charge obtained by replacing Y z ∂ z with L Y in (A.35).
Just as the orbital angular momentum piece of S (1)− e , multiplying eQ i (p i · − ) −1 to (A.33) should give us the spin angular momentum part, i.e.
Thus the full subleading soft factor becomes The spin-corrected action of hard charge is obtained by replacing Y z ∂ z with L Y in (A.35), Accordingly, the charge conservation (A.39) becomes The new term is exactly the spin contribution (A.42) to the subleading soft factor S (1)− e .