de Soto, FelicianoBlossier, B.Boucaud, P.Brinet, M.Morenas, V.Pène, O.Petrov, K.Rodríguez-Quintero, J.2014-06-302014-06-302013-04Physical Review D - Particles, Fields, Gravitation and Cosmology, 87 (7), art. no. 0740331550-7998http://hdl.handle.net/10433/999Gluon and ghost propagators, obtained in Landau gauge from lattice simulations with two light two heavy dynamical quark flavors (Nf=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 A2. 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(¿4). As ¿MS¯¯¯¯¯¯ is also rather well known for Nf=2+1+1, this allows for a precise consistency test of the operator product expansion approach in the joint description of different observables.application/pdfenAmerican Physical SocietyLatticeGluon propagatorsGhost propagatorsjournal articleopen access