Eikonal methods applied to gravitational scattering amplitudes
  We apply factorization and eikonal methods from gauge theories to scattering amplitudes in gravity. We hypothesize that these amplitudes factor into an IR-divergent soft function and an IR-finite hard function, with the former given by the expectation value of a product of gravitational Wilson line operators. Using this approach, we show that the IR-divergent part of the n-graviton scattering amplitude is given by the exponential of the one-loop IR divergence, as originally discovered by Weinberg, with no additional subleading IR-divergent contributions in dimensional regularization. 
Stephen G. Naculich and Howard J. Schnitzer
hep-th gr-qc hep-ph