Controlling quark mass determinations non-perturbatively in three-flavour QCD

The determination of quark masses from lattice QCD simulations requires a non-perturbative renormalization procedure and subsequent scale evolution to high energies, where a conversion to the commonly used MS-bar scheme can be safely established. We present our results for the non-perturbative running of renormalized quark masses in Nf=3 QCD between the electroweak and a hadronic energy scale, where lattice simulations are at our disposal. Recent theoretical advances in combination with well-established techniques allows to follow the scale evolution to very high statistical accuracy, and full control of systematic effects.