  We find that scattering amplitudes in massive scalar QCD can manifest the
duality between color and kinematics at loop-level. Specifically we construct
the one-loop integrands for four-point scattering between two distinct massive
scalars, and the five-point process encoding the first correction to massive
scalar scattering with gluonic radiation. We find that factorization and the
color-kinematics duality are sufficient principles to entirely bootstrap these
calculations, allowing us to construct all contributions ultimately from the
three-point tree-level amplitudes which are themselves entirely constrained by
symmetry. Double-copy construction immediately provides the associated
predictions for massive scalars scattering in the so called N=0 supergravity
theory.
