Testing the operator product expansion Wilson coefficient for A2 from lattice QCD with a dynamical charm

Gluon and ghost propagators, obtained in Landau gauge from lattice simulations with two light two heavy dynamical quark flavors (N_f=2+1+1), are described here with a running formula including a four-loop perturbative expression and a nonperturbative operator product expansion correction dominated by the local operator A². The Wilson coefficients and their variation as a function of the coupling constant are extracted from the numerical data and compared with the theoretical expressions that, after being properly renormalized, are known at O(α⁴). As Λ_MS¯ is also rather well known for N_f=2+1+1, this allows for a precise consistency test of the operator product expansion approach in the joint description of different observables.