We consider 1/2 BPS supersymmetric circular Wilson loops in four-dimensional \( \mathcal{N} \) = 2 SU(N) SYM theories with massless matter content and non-vanishing β-function. Following Pestun’s approach, we can use supersymmetric localization on the sphere \( \mathbbm{S} \) 4 to map these observables into a matrix model, provided that the one-loop determinants are consistently regularized. Employing a suitable procedure, we construct the regularized matrix model for these theories and show that, at order g4, the predictions for the 1/2 BPS Wilson loop match standard perturbative renormalization based on the direct evaluation of Feynman diagrams on \( \mathbbm{S} \) 4. Despite conformal symmetry begin broken at the quantum level, we also demonstrate that the matrix model approaches perfectly captures the expression of the renormalized observable in flat space at this perturbative order. Moreover, we revisit in detail the difference theory approach, showing that when the β-function is non-vanishing, this method does not account for evanescent terms which are made finite by the renormalization procedure and participate to the corrections at order g6.
