HU-EP-14/36

Factorized soft graviton theorems

at loop level

Johannes Broedel, Marius de Leeuw,

[0.3cm] Jan Plefka and Matteo Rosso

Institut für Theoretische Physik,

Eidgenössische Technische Hochschule Zürich,

Wolfgang-Pauli-Strasse 27, 8093 Zürich, Switzerland

[0.1cm] {jbroedel,mrosso}@itp.phys.ethz.ch

Niels Bohr Institute,

Copenhagen University,

Blegdamsvej 17, 2100 Copenhagen Ø, Denmark

[0.1cm] deleeuwm@nbi.ku.dk

Institut für Physik und IRIS Adlershof,

Humboldt-Universität zu Berlin,

Zum Großen Windkanal 6, 12489 Berlin, Germany

[0.1cm] plefka@physik.hu-berlin.de

Abstract

We analyze the low-energy behavior of scattering amplitudes involving gravitons at loop level in four dimensions. The single-graviton soft limit is controlled by soft operators which have been argued to separate into a factorized piece and a non-factorizing infrared divergent contribution. In this note we show that the soft operators responsible for the factorized contributions are strongly constrained by gauge and Poincaré invariance under the assumption of a local structure. We show that the leading and subleading orders in the soft-momentum expansion can not receive radiative corrections. The first radiative correction occurs for the sub-subleading soft graviton operator and is one-loop exact. It depends on only two undetermined coefficients which should reflect the field content of the theory under consideration.## 1 Introduction

Recently the low energy behavior of graviton scattering amplitudes with a single soft graviton momentum has been related [Strominger:2013jfa, He:2014laa, Kapec:2014opa] to Ward identities of the extended Bondi, van der Burg, Metzner and Sachs (BMS) symmetry [Bondi:1962px, Sachs:1962wk, Barnich:2010eb]. The role of BMS symmetry as a potential hidden symmetry of the quantum gravity S-matrix in asymptotically flat four-dimensional space-time has triggered considerable interest.

In the soft limit the -point scattering amplitude is dominated by a soft pole as was shown by Weinberg about fifty years ago [Weinberg:1964ew]. The soft amplitude factorizes on the pole into a universal soft function and the remaining hard graviton -point amplitude. This property holds in any spacetime dimension.

Revived by the work of Cachazo and Strominger [Cachazo:2014fwa] the universal factorization has been shown to extend to sub- and sub-subleading order in the soft-momentum expansion [Schwab:2014xua, Afkhami-Jeddi:2014fia, Kalousios2014, Zlotnikov2014]. In distinction to the leading Weinberg pole, which is a function of the soft and hard momenta and the soft graviton polarization, the sub- and sub-subleading soft behavior of the -point graviton amplitude can be expressed in terms of differential operators in the hard momenta and hard graviton polarizations acting on the hard -point amplitude.

Gluon amplitudes exhibit a similar universal leading and subleading soft behavior at tree level. Known as Low’s theorem [Low:1958sn, Burnett:1967km], the form of the soft operators was recently recast into the language of modern on-shell techniques [Casali:2014xpa]. The soft limit of open-string tree amplitudes was studied in refs. [Schwab:2014fia, Bianchi:2014gla], where further corrections to the field-theory soft theorems have been excluded. Interestingly, the subleading soft theorems for gravitons and gluons may also be derived from the soft limit of vertex operators in a dual ambitwistor string theory [Geyer:2014lca].

In a series of recent papers [Bern:2014vva, White:2014qia] the existence of subleading soft-graviton and soft-gluon theorems was shown to be a consequence of on-shell gauge invariance and Feynman diagrammatic reasoning. Simultaneously, an alternative point of view was put forward by the present authors in ref. [Broedel:2014fsa]: The tree-level soft-graviton and soft-gluon expressions are composed from a highly restricted class of operators compatible with on-shell gauge invariance, factorization and Poincaré symmetry. This result could be established without referring to underlying Feynman diagrammatics at all. For the graviton case, this reasoning supplements the established sub- and sub-subleading soft operators with one further candidate each.

Given these results it is natural to ask whether the soft-graviton theorems receive corrections at loop level [Bern:2014oka, He:2014bga]. In four-dimensional gravity the dimension of the coupling constant leads to strong constraints on the possible loop corrections. While the leading soft function is known to be free of radiative corrections [Weinberg:1965nx, Bern:1998sv], the subleading and sub-subleading operators do not receive corrections beyond one and two loops respectively, as argued in [Bern:2014oka, He:2014bga]. Therefore – upon controlling the soft limit of leg with momentum by the parameter – the -loop graviton amplitude should behave as

(1) | ||||

(2) |

where denotes the graviton amplitude and the leading
(), subleading () or sub-subleading () soft graviton operators
at various loop levels^{1}^{1}1In the one-loop case the last term of course
drops out, similarly at tree-level the last two terms are absent..

In ref. [Cachazo:2014dia] it was proposed that loop corrections to soft theorems could be generally suppressed by taking the soft limit prior to removing the dimensional regulator. Effectively this amounts to a study of the soft limit at the level of the loop integrand. However, such an order of limits appears physically unjustified to us. Moreover, as argued in refs. [Bern:2014oka, Bern:2014vva] and extensively shown in ref. [Bonocore:2014wua], this order of limits would invalidate the cancellation of leading IR divergences when applied to QCD.

A central point in the analysis of loop corrections to soft operators concerns the expected separation into two distinct pieces [Bern:2014vva]: The “factorizing” and the “non-factorizing” contributions. The factorizing pieces correspond to the first and the last graph in figure 1: while the leading soft factor acting on the one-loop amplitude is depicted in the first graph, the last graph shows loop-corrected soft operators acting on the hard amplitude. The remaining middle graph is associated to the non-factorizing contributions.

The separation into factorizing and non-factorizing contributions was initially laid out for collinear limits in [Bern:1995ix]. Working in dimensional regularization the poles in are to be subtracted from the factorized diagrams and attributed to the non-factorized contributions to the amplitude. Related to infrared divergences, the latter are under good control. This structure was established in gauge theory in [Bern:1995ix, Kosower:1999rx, Kosower:2003cz]. To our understanding a proof of an analogous separation into factorizing and non-factorizing pieces does not exist for gravitational amplitudes yet. Given the results for gauge theory, however, it appears highly plausible to exist at least at the one-loop level. An analysis of the graviton soft operators responsible for the leading divergent contributions to the operators was performed up to the two-loop level in [Bern:2014oka].

Two different concepts are important for the description of soft operators in the current article: “locality” and “universality”. A soft operator is said to be local, if the soft particle interacts with each hard particle separately. Universality of a soft operator implies that the operator is of the same form for any amplitude, independent of the helicity configuration in question. In the current article, this universality is implemented by expressing all soft operators in terms of differential operators in momenta and polarization vectors. Accordingly, spinor-helicity expressions appearing in a calculation of the soft behavior might differ for different helicity configurations despite originating from a single universal operator expressed in polarization and momentum vectors.

Related to the discussion of locality and universality above, an important
subtlety was discussed in ref. [Bern:2014vva]^{2}^{2}2We would like to
thank Zvi Bern, Paolo di Vecchia, Josh Nohle and Scott Davies for pointing out
and explaining this subtlety.: Contrary to the naive expectation, the
factorizing contributions in gauge theory and gravity do contain
non-local parts. While in the
gauge-theory scenario this means that the soft particle interacts with
non-neighboring legs, it implies interaction with more than one hard leg at a
time for gravitational theories. While those contributions can be determined
for every particular configuration, a general formula reproducing those parts
is not known. Accordingly, the analysis in the current article will be limited
to the universal, local and factorizing contributions.

In the present note we extend our tree-level analysis [Broedel:2014fsa] of the operators appearing in soft-graviton theorems to loop level. In order to do so, we will assume that the separation into a factorizing and a non-factorizing part extends to higher loops for gravitons in four dimensions. Starting again from on-shell gauge invariance, factorization and Poincaré symmetry, we adapt the formalism to the appropriate mass dimensions at different loop levels and identify the restricted class of operators which can appear at leading, subleading and sub-subleading order. A crucial additional assumption is again the “local” form of the soft operators explained above. As our method allows to constrain operators composed of hard polarizations and momenta only, it is blind to possible non-universal contributions to the factorized part.

Before moving on to gravity, let us comment on the situation in four-dimensional gauge theory. Here there are no novel loop-level operators arising from our analysis because of the dimensionless coupling constant. This implies in particular that the new subleading operator identified in [Broedel:2014fsa] captures the form of all radiative corrections. In fact, the form reported in eq. (34) of [Broedel:2014fsa] coincides with the one-loop contribution quoted in [Bern:2014vva]. The undetermined coefficient depends on the matter content of the gauge theory in question, whereas the tensorial form of the subleading soft-gluon operator is universal. However, it is precisely at this order where one first encounters the non-local parts mentioned above: As was analyzed carefully in section 3.2 of [He:2014bga] the single-minus one-loop -gluon amplitude with a soft leg of positive helicity develops a factorized non-local subleading soft pole, whose corresponding operator depends on next-to-nearest neighboring hard legs of the soft leg. Unfortunately our formalism is unable to capture this behavior at present.

## 2 Method and previous results

For completeness we repeat our method and the resulting soft-graviton operators at tree-level [Broedel:2014fsa]. We denote four-dimensional graviton amplitudes by where is the total momentum. The momentum of leg is expressed as , which can be made soft by sending . The remaining hard momenta are labeled with . Moreover, we write the polarization tensor of the soft leg as and similarly for the hard legs . This splitting is an identity in four dimensions. In this language, the soft-graviton theorem [Cachazo:2014fwa] reads

(3) |

where we have suppressed the polarization dependence of the amplitudes . The soft operators are in general differential operators in the hard momenta and polarizations and also depend on the soft data and .

### 2.1 Constraints on factorized soft operators

We recall here the assumptions on the form of the soft operators and the constraints they must satisfy, both of which were analyzed in detail in ref. [Broedel:2014fsa].

##### Locality.

The only assumption on the form of the soft terms is what we termed locality. We assume that the soft operators can be expressed as a sum of local operators, each of which acts on the soft leg and one single hard leg:

(4) |

Clearly this assumption is justified at tree level. At loop level, however, there are contributions which are either non-local or non-universal or even both, as shown in explicit computations of the IR divergences in refs. [Bern:2014vva, He:2014bga].

##### Distributional constraint.

The self-consistency of eq. (3) gives rise to the distributional constraint [Broedel:2014fsa]

(5) |

The expansion of in leads to relations between the operators . In particular, the terms in coupling to derivatives become constrained by lower-order soft operators .

##### Gauge invariance.

Naturally, each soft operator must be gauge invariant. For the soft leg this implies that needs to be invariant under the shift for . This can be achieved by demanding

(6) |

where the symbol indicates vanishing modulo Poincaré generators

(7) |

Notice that eq. (6) is a necessary but not sufficient condition for gauge invariance, because it corresponds to the choice . Similarly, a gauge transformation on the soft leg may be represented by the generator

(8) |

As the gauge invariance of
implies the vanishing of the commutator^{3}^{3}3Strictly speaking, this is only a sufficient condition for
to be gauge invariant. Nevertheless, we will see that
the results in the following sections vanish upon the substitution with .

(9) |

which further constrains the form of the soft operators. Here, the zero in (9) is modulo Poincaré and gauge transformations which annihilate . In ref. [Broedel:2014fsa] this constraint was satisfied by considering a particular combination of differential operators to appear in the soft theorems:

(10) |

In the current article, however, we will stay more general and refrain from using this assumption.

##### Mass dimensions and loop counting.

Finally, we need to consider the correct mass dimensions for the soft operators. In four-dimensional gravity ought to have vanishing mass dimension. This is very important to keep in mind when considering loop corrections in four dimensions, as the coupling constant is dimensionful. Accordingly, for every loop order an extra factor of appears.

### 2.2 Ansätze

Let us start by noting Weinberg’s leading soft function

(11) |

We then write down the most general ansätze for the sub- and sub-subleading soft operators compatible with the above constraints. The distributional constraint eq. (5) and the form of eq. (11) require to be a differential operator of first order in whereas should be of second order in . Hence, the building blocks for the ansätze are single- and double-derivative operators of the schematic form

(12) |

where the tensor can take the values

(13) |

Moreover, the vector takes one of the three values

(14) |

where is a “dummy” index vector waiting to be contracted with itself as in

(15) |

Obviously, only even powers of are kept in the ansätze of eqs. (2.2). Multiple occurrences of are allowed and encode all possible contractions weighted with independent coefficients. In this manner we generate all allowed terms.

The number in eq. (2.2) counts the effective power of in the square brackets. Note that for and for . The allowed terms emerging from eq. (2.2) are severely reduced by the on-shell conditions . Thus armed, it is straightforward to identify the potential single- and double-derivative contributions to the soft operators , , at tree and loop level, see table 1.

Tree | 1-loop | 2-loop | 3-loop | |
---|---|---|---|---|

, | – | – | ||

, | , | – | ||

, | , | , |

Finally let us comment on the possibility of a local soft operator without any derivative terms. It needs to have the form (cf. eq. (2.2))

(16) |

because the only admissible value for is . Imposing gauge invariance on the soft leg, i.e. acting with , immediately restricts the invariant operators to , which is the Weinberg soft factor eq. (11). Hence, we conclude that there is no other possible no-derivative structure for .

### 2.3 Tree level results

After applying the constraints, the tree level results read

(17) | ||||

(18) | ||||

(19) | ||||

with the Lorentz generator density

(20) |

and four undetermined coefficients . In ref. [Broedel:2014fsa] the constants , were claimed to a priori be different for each hard leg. However, since gravity amplitudes are invariant under any permutation of the external legs, the above soft factors are compatible with this additional symmetry only if and agree for all hard legs. Furthermore, the terms proportional to and were not displayed in [Broedel:2014fsa], because double-derivative terms for the subleading operator were not considered there. Finally, notice that this result is derived without any reference to the precise form of the amplitudes and .

The only features used are gauge and Lorentz invariance together with total momentum conservation (ensured by an overall momentum-conserving delta function). The above results apply in particular to hard amplitudes involving scalars or photons along with gravitons. While for a hard scalar the single derivative vanishes, for a hard photon only the double derivative does.

In ref. [Cachazo:2014fwa], it was shown that for tree-level pure-graviton amplitudes in four dimensions.

## 3 Soft-graviton operators at loop level

Let us now study the admissible soft operators at loop level. Considering loops requires taking the mass dimension of the four-dimensional gravitational coupling into account. Following the derivation of [Broedel:2014fsa] outlined above, we find several additional permissible terms in the ansätze. A naïve way to produce consistent higher loop corrections satisfying gauge invariance and obeying the distributional constraints consists of promoting to by multiplication with . However, it is quickly checked that the terms proportional to in eqns. (2.3) as well as the promotion of violate gauge invariance.

### 3.1 Vanishing of the highest-loop contributions

As stated above, there are no loop corrections to [Weinberg:1965nx, Bern:1998sv]. Employing our method, this can be easily rederived for the local and universal factorizing contributions by looking at table 1: The potential tree-level derivative operators and vanish by explicit computation. The double-derivative candidate for a one-loop contribution to takes the form

(21) |

Upon inserting either form of in eq. (13) we see that vanishes by virtue of the on-shell relations . Hence there are no factorizing loop corrections to the leading soft-graviton behavior. This simple argument directly applies to the potential two-loop contribution to in as well as to the potential three-loop contribution to in . Since they differ by one and two factors of from eq. (21) respectively, they vanish by the same reasoning.

### 3.2 Non-vanishing loop contributions

The absence of loop corrections to the leading soft operator has an immediate consequence for the possible form of the subleading soft operators at higher loop orders. The distributional constraint eq. (5) applied to the higher-loop contributions entails that the subleading soft operators and may not contain momentum derivatives terms but only derivatives with respect to the hard polarizations . This further restricts our ansätze and leads to the following results: For the subleading soft operator all possible loop contributions vanish because of incompatibility with the constraints.

Turning to the sub-subleading soft operator we find that there are non-vanishing local and universal factorizing contributions at one-loop order exclusively. We find four contributions to , all of which originate from the tree-level contributions to in eq. (2.3) by multiplying with as explained above. Collecting everything, we find

(22) |

As a matter of fact, the term proportional to is effectively equal to
the term proportional to because and the
remaining acts – together with the further
differential operator – as a gauge transformation. Hence we may set
^{4}^{4}4Notice that the term proportional to already originates at tree-level,
but can be reabsorbed in other terms in a way similar to what explained
before. .

Beyond these three terms all local and universal factorizing loop corrections to the soft operators vanish. In particular there cannot be any factorizing loop corrections to or . Moreover, the factorizing loop corrections to are one-loop exact as shown above.

## 4 Rewriting polarization derivatives

All loop corrections appear in the form of polarization derivatives. Using the completeness relation [Dixon:1996wi]

(23) |

these can be simplified and expressed in a more convenient form. Here is the polarization of opposite helicity compared to
and is the reference vector needed to define the
polarizations^{5}^{5}5Notice that .. Using this
language, the -graviton amplitude can be rewritten as

(24) |

Accordingly, we define the conjugate-helicity amplitude as well as the effective scalar-graviton amplitude as

(25) |

The last relation follows from inserting the completeness relation eq. (23) by virtue of . Using the above notation, we find after dropping the arguments of the amplitudes

(26) | ||||

(27) |

The particular combinations appearing in the soft operators can now be expressed in the following way

(28) | ||||

(29) |

where and are the gauge invariant quantities

(30) |

An important relation is

(31) |

which is most easily seen by choosing as the reference vector for
and . Then .
The later expression always vanishes: depending on the soft and hard
helicities either the first or the second factor is zero^{6}^{6}6
In spinor helicity language one finds
(32)
.
Using these relations one may now simplify the soft-operators
at tree and loop level obtained in the previous sections.

##### Tree level.

For the novel tree-level terms in our soft operators we find at subleading order

(33) |

The contributions proportional to are of the same form as the loop-level contribution to which will be discussed below. The remaining undetermined terms of the sub-subleading soft operator parametrized by and read

(34) | ||||

(35) |

where in the last step we used the identities and as well as .

##### Loop level.

The loop corrections in eq. (22) can be rewritten as

(36) |

which can be further simplified by choosing the reference momentum to equal . This is possible because is gauge invariant by itself. After all, we arrive at the final result

(37) |

##### Interpretation.

The loop correction couples to and exclusively. Note that the latter correction is precisely the structure that was uncovered in an explicit computation for a higher-derivative gravity theory where is a single-scalar -graviton amplitude [Bianchi:2014gla]. To our knowledge there is no general theorem stating that the polarization-stripped and traced amplitude vanishes in pure gravity at loop level. While it does so for four gravitons at tree level, which can be tested using the explicit results in ref. [Sannan:1986tz], calculations are involved already at the one-loop level. At least the trace of the four-point integrand does not vanish.

## 5 Discussion

In this paper we have determined all admissible soft factors for local and universal factorizing loop corrections in four-dimensional gravity.

We found that only the sub-subleading soft operator receives loop corrections. It turns out that these corrections are one-loop exact and have two contributions proportional to and . This is consistent with the results in ref. [Bern:2014vva], where loop corrections arising from virtual scalars in graviton amplitudes were shown to not contribute to subleading soft factors.

Hence the local and universal factorizing contributions to the soft graviton operators do not receive corrections at first and second perturbative order. This implies that the soft-graviton Ward identities of extended BMS symmetry do not suffer from anomalies caused by local universal operators in the factorized sector studied here. They will be affected, however, by the non-factorized contributions entangled with the infrared divergences. It would be interesting to understand this sector in detail as well.

While the tensorial form of the loop contributions to the sub-subleading soft graviton theorems is fixed by our analysis, the undetermined scalar prefactors will reflect the field content of the theory. For example, our coupling to exactly matches the structure found in [Bianchi:2014gla] for gravity.

Finally, let us once more comment on gauge theory in four dimensions where the coupling constant is dimensionless. Therefore, the loop corrections simply correspond to the undetermined part of the subleading soft operator of section 3 in ref. [Broedel:2014fsa], which in turn agrees with the one-loop correction to quoted in [Bern:2014vva]. The undetermined coefficient again depends on the matter content of the gauge theory in question. The tensorial form of the soft operators is, however, universal.

#### Acknowledgments

We thank Z. Bern, S. Davies, B. Schwab, J. Nohle, P. di Vecchia and C. Vergu for discussions. JP thanks the Pauli Center for Theoretical Studies Zürich and the Institute for Theoretical Physics at the ETH Zürich for hospitality and support in the framework of a visiting professorship. The work of MdL and MR is partially supported by grant no. 200021-137616 from the Swiss National Science Foundation. MdL was also supported in part by FNU through grant number DFF–1323–00082.

## References

\nc@nextpunct\@bibitemStrominger:2013jfaA. Strominger, “On BMS Invariance of Gravitational Scattering”, arxiv:1312.2229.

He:2014laaT. He, V. Lysov, P. Mitra and A. Strominger, “BMS supertranslations and Weinberg’s soft graviton theorem”, arxiv:1401.7026.

Kapec:2014opaD. Kapec, V. Lysov, S. Pasterski and A. Strominger, “Semiclassical Virasoro Symmetry of the Quantum Gravity S-Matrix”, arxiv:1406.3312.

Bondi:1962pxH. Bondi, M. van der Burg and A. Metzner, “Gravitational waves in general relativity. 7. Waves from axisymmetric isolated systems”, Proc.Roy.Soc.Lond. A269, 21 (1962).

Sachs:1962wkR. Sachs, “Gravitational waves in general relativity. 8. Waves in asymptotically flat space-times”, Proc.Roy.Soc.Lond. A270, 103 (1962).

Barnich:2010ebG. Barnich and C. Troessaert, “Aspects of the BMS/CFT correspondence”, JHEP 1005, 062 (2010), arxiv:1001.1541.

Weinberg:1964ewS. Weinberg, “Photons and Gravitons in s Matrix Theory: Derivation of Charge Conservation and Equality of Gravitational and Inertial Mass”, Phys.Rev. 135, B1049 (1964).

Cachazo:2014fwaF. Cachazo and A. Strominger, “Evidence for a New Soft Graviton Theorem”, arxiv:1404.4091.

Schwab:2014xuaB. U. W. Schwab and A. Volovich, “Subleading soft theorem in arbitrary dimension from scattering equations”, arxiv:1404.7749.

Afkhami-Jeddi:2014fiaN. Afkhami-Jeddi, “Soft Graviton Theorem in Arbitrary Dimensions”, arxiv:1405.3533.

Kalousios2014C. Kalousios and F. Rojas, “Next to subleading soft-graviton theorem in arbitrary dimensions”, arxiv:1407.5982.

Zlotnikov2014M. Zlotnikov, “Sub-sub-leading soft-graviton theorem in arbitrary dimension”, JHEP 1410, 148 (2014), arxiv:1407.5936.

Low:1958snF. Low, “Bremsstrahlung of very low-energy quanta in elementary particle collisions”, Phys.Rev. 110, 974 (1958).

Burnett:1967kmT. Burnett and N. M. Kroll, “Extension of the low soft photon theorem”, Phys.Rev.Lett. 20, 86 (1968).

Casali:2014xpaE. Casali, “Soft sub-leading divergences in Yang-Mills amplitudes”, arxiv:1404.5551.

Schwab:2014fiaB. U. Schwab, “Subleading Soft Factor for String Disk Amplitudes”, JHEP 1408, 062 (2014), arxiv:1406.4172.

Bianchi:2014glaM. Bianchi, S. He, Y.-t. Huang and C. Wen, “More on Soft Theorems: Trees, Loops and Strings”, arxiv:1406.5155.

Geyer:2014lcaY. Geyer, A. E. Lipstein and L. Mason, “Ambitwistor strings at null infinity and subleading soft limits”, arxiv:1406.1462.

Bern:2014vvaZ. Bern, S. Davies, P. Di Vecchia and J. Nohle, “Low-Energy Behavior of Gluons and Gravitons from Gauge Invariance”, Phys.Rev. D90, 084035 (2014), arxiv:1406.6987.

White:2014qiaC. White, “Diagrammatic insights into next-to-soft corrections”, arxiv:1406.7184.

Broedel:2014fsaJ. Broedel, M. de Leeuw, J. Plefka and M. Rosso, “Constraining subleading soft gluon and graviton theorems”, Phys.Rev. D90, 065024 (2014), arxiv:1406.6574.

Bern:2014okaZ. Bern, S. Davies and J. Nohle, “On Loop Corrections to Subleading Soft Behavior of Gluons and Gravitons”, Phys.Rev. D90, 085015 (2014), arxiv:1405.1015.

He:2014bgaS. He, Y.-t. Huang and C. Wen, “Loop Corrections to Soft Theorems in Gauge Theories and Gravity”, arxiv:1405.1410.

Weinberg:1965nxS. Weinberg, “Infrared photons and gravitons”, Phys.Rev. 140, B516 (1965).

Bern:1998svZ. Bern, L. J. Dixon, M. Perelstein and J. Rozowsky, “Multileg one loop gravity amplitudes from gauge theory”, Nucl.Phys. B546, 423 (1999), hep-th/9811140.

Cachazo:2014diaF. Cachazo and E. Y. Yuan, “Are Soft Theorems Renormalized?”, arxiv:1405.3413.

Bonocore:2014wuaD. Bonocore, E. Laenen, L. Magnea, L. Vernazza and C. D. White, “The method of regions and next-to-soft corrections in Drell-Yan production”, arxiv:1410.6406.

Bern:1995ixZ. Bern and G. Chalmers, “Factorization in one loop gauge theory”, Nucl.Phys. B447, 465 (1995), hep-ph/9503236.

Kosower:1999rxD. A. Kosower and P. Uwer, “One loop splitting amplitudes in gauge theory”, Nucl.Phys. B563, 477 (1999), hep-ph/9903515.

Kosower:2003czD. A. Kosower, “All orders singular emission in gauge theories”, Phys.Rev.Lett. 91, 061602 (2003), hep-ph/0301069.

Dixon:1996wiL. J. Dixon, “Calculating scattering amplitudes efficiently”, hep-ph/9601359.

Sannan:1986tzS. Sannan, “Gravity as the Limit of the Type II Superstring Theory”, Phys.Rev. D34, 1749 (1986).