Scalar-graviton amplitudes and celestial holography

We compute scattering amplitudes involving one massive scalar and two, three, or four gravitons. We show that when the conformal dimension of the massive scalar is set to zero, the resulting celestial correlators depend only on the coordinates of the gravitons. Such correlators of gravitons are well-defined and do not suffer from divergences associated with the Mellin transform of usual graviton amplitudes. Moreover, they are non-distributional and take the form of standard CFT correlators. We show that they are consistent with the usual OPEs but the statement of the soft theorem is modified.