Thanks for the good question. To my knowledge, some specialized techniques (other than the standard modern multi-loop techniques) have been developed for the state-of-art calculations of g-2. Nevertheless, for lower-loop calculations, AmpRed remains a highly suitable tool. It allows us to decouple the loop integrals from the spinor chain from the very beginning, and it provides tools to automatically convert Dirac spinors and gamma matrices to Weyl spinors and Pauli matrices. Furthermore, the decoupled loop integrals can be evaluated with a single command. To illustrate this in detail, I’ll prepare a short notebook in about one month.

Hi Wen! Wonderful work!

I've been looking into various HEP Mathematica packages (like FeynArts, FeynCalc, FeynHelpers, LoopTools, Package-X etc.) recently and got frustrated with how hard it is to get some results without digging very hard into theory and specific methods. I wonder how easier it is to compute a g-factor using AmpRed, at least for a single loop like it is done here for example? We need to calculate the form factor somehow by turning them into divergent single-loop integrals, simplifying amplitudes, and substituting all the constants. All of which require some steep learning curve actually to execute in practice.

I would really like to see these tools to have a more human design for particle physics in the future, similar to how it is easier to compute with molecules, bio-sequences, astronomical objects, and other physical entities now with the Wolfram Language.