“…Recently, a new use of the gradient flow method [1,2,3,4] was proposed to calculate correctly renormalized observables from lattice simulations [5,6]: Making use of the finiteness of flowed operators, non-perturbative estimates of observables are extracted by taking a vanishing flow-time extrapolation. The new method was first applied to the energy-momentum tensor for which the explicit violation of the Poincaré invariance on the lattice has been a hard obstacle in obtaining a non-perturbative estimate [5,7,8]. However, because the method is quite general, we can apply it also to obtain correctly renormalized chiral observables from lattice simulations [6].…”