  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.
