In this talk, we discuss top-quark contributions to $ZZ$ production through gluon fusion at two loops. We use syzygies and modular arithmetic to compute the reductions to master integrals.
In order to numerically evaluate the amplitude, we express it in terms of finite integrals using a new method to construct integrable linear combinations of divergent integrands.