Training normalizing flows with computationally intensive target probability distributions

Abstract

Machine learning techniques, in particular the so-called normalizing flows, are becoming increasingly popular in the context of Monte Carlo simulations as they can effectively approximate target probability distributions. In the case of lattice field theories (LFT) the target distribution is given by the exponential of the action. The common loss function's gradient estimator based on the "reparametrization trick" requires the calculation of the derivative of the action with respect to the fields. This can present a significant computational cost for complicated, non-local actions like e.g. fermionic action in QCD. In this contribution, we propose an estimator for normalizing flows based on the REINFORCE algorithm that avoids this issue. We apply it to two dimensional Schwinger model with Wilson fermions at criticality and show that it is up to ten times faster in terms of the wall-clock time as well as requiring up to $30\%$ less memory than the reparameterization trick estimator. It is also more numerically stable allowing for single precision calculations and the use of half-float tensor cores. We present an in-depth analysis of the origins of those improvements. We believe that these benefits will appear also outside the realm of the LFT, in each case where the target probability distribution is computationally intensive.

arXiv comment16 pages, 5 figures, 6 tables, 3 listings. Revised version as published in CPC. Added results for a other values of hoping parameter

References in corpus (11)

• Normalizing Flows: An Introduction and Review of Current Methods
• Solving Statistical Mechanics Using Variational Autoregressive Networks
• Flow-based generative models for Markov chain Monte Carlo in lattice field theory
• Equivariant flow-based sampling for lattice gauge theory
• Estimation of Thermodynamic Observables in Lattice Field Theories with Deep Generative Models
• Asymptotically unbiased estimation of physical observables with neural samplers
• Gauge-equivariant flow models for sampling in lattice field theories with pseudofermions
• Flow-based sampling in the lattice Schwinger model at criticality
• Analysis of autocorrelation times in Neural Markov Chain Monte Carlo simulations
• Hierarchical autoregressive neural networks for statistical systems
• Simulating first-order phase transition with hierarchical autoregressive networks

Cited by in corpus (1)

• NeuMC -- a package for neural sampling for lattice field theories