We present fully differential next-to-leading order results for Higgs production in association with a Z boson in gluon fusion. Our two-loop virtual contributions are evaluated numerically using sector decomposition, including full top-quark mass effects, and supplemented at high pT by an analytic high-energy expansion to order $$ {m}_Z^4,{m}_H^4,{m}_t^{32} $$
m
Z
4
,
m
H
4
,
m
t
32
. Using the expanded results we also present a study of the top-quark mass scheme uncertainty at large pT.