$ZH$ production in gluon fusion: two-loop amplitudes with full top quark mass dependence

We present results for the two-loop helicity amplitudes entering the NLO QCD corrections to the production of a Higgs boson in association with a $Z$-boson in gluon fusion. The two-loop integrals, involving massive top quarks, are calculated numerically. Results for the interference of the finite part of the two-loop amplitudes with the Born amplitude are shown as a function of the two kinematic invariants on which the amplitudes depend.