Automatic reduction of four-loop bubbles

We give technical details about the computational strategy employed in a recently completed investigation of the four-loop QCD free energy. In particular, the reduction step from generic vacuum bubbles to master integrals is described from a practical viewpoint, for fully massive as well as QED-type integrals.