Functional Renormalization Group flows as diffusive Hamilton-Jacobi-type equations
Adrian Koenigstein, Martin J. Steil, Stefan Floerchinger
TL;DR
This work reframes FRG flows for two-point functions as viscous Hamilton-Jacobi-type PDEs, enabling stable, high-resolution numerical treatment via a KT-HJ finite-volume scheme. It starts from a zero-dimensional fermion-boson toy model to derive exact and numerical flow equations, then demonstrates the method on higher-dimensional truncations including a $d=3$ Z2 model and a $1+1$-dimensional Gross-Neveu–Yukawa system, highlighting how advection and diffusion terms govern field-space dynamics and how shocks, phase transitions, and sign problems manifest in the flows. The key contributions are the identification of the functional HHJ structure, the development of a hybrid conservative/nonconservative numerical scheme that handles field-dependent couplings, and the demonstration of convergence and stability across challenging regimes, with insights into convexity, Yang-Lee zeros, and stiffness. The results open avenues for more robust nonperturbative FRG calculations and suggest connections to stochastic control and viscosity solution theory, potentially guiding future numerical and analytical developments in FRG studies. The practical impact lies in enabling reliable, scalable simulations of complex FRG flows with field-dependent interactions, aiding exploration of phase structure and critical phenomena in quantum and statistical field theories.
Abstract
In order to find reliable and efficient numerical approximation schemes, we suggest to identify the Functional Renormalization Group flow equations of one-particle irreducible two-point functions as Hamilton-Jacobi(-Bellman)-type partial differential equations. Based on this reformulation and reinterpretation we adopt a numerical scheme for the solution of field-dependent flow equations as nonlinear partial differential equations. We demonstrate this novel approach by first applying it to a simple fermion-boson system in zero spacetime dimensions - which itself presents as an interesting playground for method development. Afterwards, we show, how the gained insights can be transferred to more interesting problems: One is the bosonic $\mathbb{Z}_2$-symmetric model in three Euclidean dimensions within a truncation that involves the field-dependent effective potential and field-dependent wave-function renormalization. The other example is the $(1 + 1)$-dimensional Gross-Neveu model within a truncation that involves a field-dependent potential and a field-dependent fermion mass/Yukawa coupling at nonzero temperature, chemical potential, and finite fermion number.
