Functional matching and renormalization group equations at two-loop order

We present a systematic method for determining the two-loop effective Lagrangian resulting from integrating out a set of heavy particles in an ultraviolet scalar theory. We prove that the matching coefficients are entirely determined from the (double-)hard region of the loop integrals and present a master formula for matching, applicable to both diagrammatic and functional approaches. We further employ functional methods to determine compact expressions for the effective Lagrangian that do not rely on any previous knowledge of its structure or symmetries. The same methods are also applicable to the computation of renormalization group equations. We demonstrate the application of the functional approach by computing the two-loop matching coefficients and renormalization group equations in a scalar toy model.