Renormalization automated by Hopf algebra
  It was recently shown that the renormalization of quantum field theory is organized by the Hopf algebra of decorated rooted trees, whose coproduct identifies the divergences requiring subtraction and whose antipode achieves this. We automate this process in a few lines of recursive symbolic code, which deliver a finite renormalized expression for any Feynman diagram. We thus verify a representation of the operator product expansion, which generalizes Chen's lemma for iterated integrals. The subset of diagrams whose forest structure entails a unique primitive subdivergence provides a representation of the Hopf algebra ${\cal H}_R$ of undecorated rooted trees. Our undecorated Hopf algebra program is designed to process the 24,213,878 BPHZ contributions to the renormalization of 7,813 diagrams, with up to 12 loops. We consider 10 models, each in 9 renormalization schemes. The two simplest models reveal a notable feature of the subalgebra of Connes and Moscovici, corresponding to the commutative part of the Hopf algebra ${\cal H}_T$ of the diffeomorphism group: it assigns to Feynman diagrams those weights which remove zeta values from the counterterms of the minimal subtraction scheme. We devise a fast algorithm for these weights, whose squares are summed with a permutation factor, to give rational counterterms. 
D.J.Broadhurst and D.Kreimer (Erwin Schr\"odinger Institute, Vienna)