Vertex functions are a crucial ingredient of several forefront many-body algorithms in condensed matter physics. However, the full treatment of their frequency and momentum dependence severely restricts numerical calculations. A significant advancement requires an efficient treatment of the high-frequency asymptotic behavior of the vertex functions. In this work, we first provide a detailed diagrammatic analysis of the high-frequency structures and their physical interpretation. Based on these insights, we propose a parametrization scheme, which captures the whole high-frequency domain for arbitrary values of the local Coulomb interaction and electronic density, and we discuss the details of its algorithmic implementation in many-body solvers based on parquet-equations as well as functional renormalization group schemes. Finally, we assess its validity by comparing our results for a single impurity Anderson model with exact diagonalization calculations. The proposed parametrization is pivotal for the algorithmic development of all quantum many-body methods based on vertex functions arising from local microscopic interactions, such as the diagrammatic approaches including spatial correlations beyond dynamical mean-field theory.
Vertex functions are a crucial ingredient of several forefront many-body algorithms in condensed matter physics. However, the full treatment of their frequency and momentum dependence severely restricts numerical calculations. A significant advancement requires an efficient treatment of the high-frequency asymptotic behavior of the vertex functions. In this work, we first provide a detailed diagrammatic analysis of the high-frequency structures and their physical interpretation. Based on these insights, we propose a parametrization scheme, which captures the whole high-frequency domain for arbitrary values of the Coulomb interaction and electronic density, and we discuss the details of its algorithmic implementation in many-body solvers based on parquet equations as well as functional renormalization group schemes. Finally, we assess its validity by comparing our results for a single impurity Anderson model with exact diagonalization calculations. The proposed parametrization is pivotal for the algorithmic development of all quantum many-body methods based on vertex functions arising from both local and nonlocal static microscopic interactions as well as effective dynamic interactions which uniformly approach a static value for large frequencies. In this way, our present technique can substantially improve vertex-based diagrammatic approaches including spatial correlations beyond dynamical mean-field theory.