Sphaleron rate of $N_f=2+1$ QCD

We compute the sphaleron rate of $N_f=2+1$ QCD at the physical point for a range of temperatures $200$ MeV $\lesssim T \lesssim 600$ MeV. We adopt a strategy recently applied in the quenched case, based on the extraction of the rate via a modified version of the Backus-Gilbert method from finite-lattice-spacing and finite-smoothing-radius Euclidean topological charge density correlators. The physical sphaleron rate is finally computed by performing a continuum limit at fixed physical smoothing radius, followed by a zero-smoothing extrapolation.Comment: Main text: 5 pages, 4 figures. Supplementary Material: 8 pages, 29 figure