Mixed electroweak-QCD corrections to e^+e^-to HZ at Higgs factory

The prospective Higgs factories, exemplified by ILC, FCC-ee and CEPC, plan to conduct the precision Higgs measurements at the $e^+e^-$ center-of-mass energy around 250 GeV. The cross sections for the dominant Higgs production channel, the Higgsstrahlung process, can be measured to a (sub-) percent accuracy. Merely incorporating the well-known next-to-leading order (NLO) electroweak corrections appears far from sufficient to match the unprecedented experimental precision. In this work, we make an important advancement toward this direction by investigating the mixed electroweak-QCD corrections to $e^+e^-\to HZ$ at next-to-next-to-leading order (NNLO) for both unpolarized and polarized $Z$ boson. The corrections turn out to reach one percent level of the Born-order results, thereby must be incorporated in the future confrontation with the data.