In this work, we compared the performance of three MLP methods, ACE, NequIP and GAP, to study different aspects of modestly complex (≈ 50 atoms) Diels-Alder reactions. We also introduce a freely available automated fitting code mlp-train to facilitate the generation of these potentials. While all three methods provide reasonable potentials, ACE and NequIP emerged as more efficient methods for generating accurate MLPs within 1 kcal mol -1 accuracy to the ground-truth surface, in particular for highly exothermic reactions.