Novel approach for computing gradients of physical observables

Simone Bacchio
2023-11-27
Abstract:We show that an infinitesimal step of gradient flow can be used for defining a novel approach for computing gradients of physical observables with respect to action parameters. Compared to the commonly used perturbative expansion, this approach does not require calculating any disconnected contribution or vacuum expectation value and can provide results up to three orders of magnitudes more precise. On the other hand, it requires a non-trivial condition to be satisfied by the flow action, the calculation of its force and its Laplacian, and the force of the observable, whose gradient needs to be measured. As a proof of concept, we measure gradients in $\beta$ of Wilson loops in a four-dimensional SU(3) Yang-Mills theory simulated on a $16^4$ lattice using the Wilson action.
High Energy Physics - Lattice
What problem does this paper attempt to address?
The paper attempts to address the challenge of calculating the gradient of observables with respect to action parameters in lattice gauge theory. Specifically, the authors propose a new method to compute the gradient of physical observables with respect to action parameters, which has the following advantages over traditional perturbative expansion methods: 1. **No need to calculate disconnected contributions or vacuum expectation values**: Traditional methods require the calculation of these terms, which are often noisy. The new method avoids these issues. 2. **Improved accuracy**: The new method can improve accuracy by 3 orders of magnitude in certain cases. The authors validated the effectiveness of this method through numerical experiments and applied it to the gradient calculation of Wilson loops in 4-dimensional SU(3) Yang-Mills theory. The experimental results show that at small β values, the error of the new method is significantly reduced, especially at β=0, where the signal-to-noise ratio is improved by more than 3 orders of magnitude. For larger β values, although the new method remains effective, its advantages are not as pronounced as at small β values. Additionally, the authors explored the convergence of the β expansion and pointed out that it converges well for β<6; however, for larger β values, more terms may be needed to achieve good convergence. This conclusion poses a challenge to flow-based methods, particularly at high β values. In summary, the paper provides a new computational method that can significantly improve calculation accuracy in certain cases and offers new insights for future applications.