  With a key improvement, the auxiliary mass flow method is now able to compute
Feynman integrals encountered in cutting-edge collider processes. We have
successfully applied it to compute some integrals involved in two-loop
electroweak corrections to $e^+e^-\to HZ$, two-loop QCD corrections to $3j$,
$W/Z/H+2j$, $t\bar{t}H$ and $4j$ production at hadron colliders, and three-loop
QCD corrections to $t\bar{t}$ production at hadron colliders, all of which are
crucial for precision frontier in collider physics in the following decade. Our
results are important building blocks and benchmarks for future studies of
these processes.
