Abstract:Diagrammatic expansions are a paradigmatic and powerful tool of quantum many-body theory. Their evaluation to high order, e.g., by the Diagrammatic Monte Carlo technique, can provide unbiased results in strongly correlated and challenging regimes. However, calculating a factorial number of terms to acceptable precision remains very costly even for state-of-the-art methods. We achieve a dramatic acceleration of evaluating Feynman's diagrammatic series by use of specialised hardware architecture within the recently introduced combinatorial summation (CoS) framework. We present how exploiting the massive parallelism and concurrency available from GPUs leads to orders of magnitude improvement in computation time even on consumer-grade hardware. This provides a platform for making probes of novel phenomena of strong correlations much more accessible.

What problem does this paper attempt to address?

### What problem does this paper attempt to solve? This paper aims to solve the problem of efficient computation of Feynman diagram expansion in quantum many - body theory. Specifically, the author attempts to significantly accelerate the computation of Feynman diagram series by taking advantage of the large - scale parallelism and concurrency of dedicated hardware architectures such as GPUs. #### Background and problem description Feynman diagram expansion is an important tool in quantum many - body theory, which can provide unbiased results, especially in strongly - correlated and challenging physical systems. However, as the order of expansion increases, the computational complexity grows factorially, and the computational cost remains very high even with the most advanced methods. #### Main contributions of the paper 1. **Acceleration with GPU**: The author shows how to achieve a significant acceleration in the computation of Feynman diagram series by applying the recently introduced combinatorial sum (CoS) framework to the GPU. They demonstrate that on consumer - level hardware, this method can improve the computation time by several orders of magnitude. 2. **Large - scale parallel computation**: By fully exploiting the large - scale parallel processing power of the GPU, the author achieves a faster computation speed than existing CPU methods. This makes it more feasible to study new physics of strongly - correlated phenomena. 3. **Algorithm optimization**: The paper details how to map the CoS algorithm onto the GPU architecture and optimize the memory access pattern through the "flattening" transformation of the graph structure, thereby further improving the computational efficiency. #### Specific problems - **Computational cost of high - order Feynman diagrams**: In traditional computational methods, when expanding high - order Feynman diagrams, the amount of computation grows factorially, resulting in extremely high computational costs. - **Numerical stability**: In high - order computations, numerical instability and loss of precision are important issues. - **Utilization of hardware resources**: How to maximize the use of the parallel processing power of the GPU to minimize the computation time and resource consumption. #### Solutions - **Combinatorial sum (CoS) algorithm**: By factorizing the graph structure, the computational complexity is reduced from O(n!) to O(n^2 * 2^n). - **GPU parallelization**: Utilize the large - scale parallel processing power of the GPU to assign the operations of each graph edge to different threads, thereby achieving efficient parallel computation. - **Graph structure optimization**: Optimize the memory access pattern and improve the computational efficiency by performing a "flattening" transformation on the graph structure. #### Results Through these improvements, the author shows that on several hardware platforms (including CPUs and different models of NVIDIA GPUs), compared to traditional CPU methods, the new method can achieve an acceleration of up to three orders of magnitude. This not only significantly improves the computational efficiency but also provides the possibility of exploring new strongly - correlated physical phenomena. ### Summary This paper solves the problem of high computational cost of high - order terms in Feynman diagram expansion computations by introducing the GPU - based combinatorial sum algorithm, achieving significant computational acceleration and providing new tools and methods for studying strongly - correlated physical systems.

Exploiting Parallelism for Fast Feynman Diagrammatics

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

Feynman Diagrams as Computational Graphs

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

The development of computational methods for Feynman diagrams

Combinatorial summation of Feynman diagrams: Equation of state of the 2D SU(N) Hubbard model

Decomposing imaginary time Feynman diagrams using separable basis functions: Anderson impurity model strong coupling expansion

Parallelizing quantum simulation with decision diagrams

Efficient reduction of Feynman integrals on supercomputers

ALGORITHMIC CALCULATION OF TWO-LOOP FEYNMAN DIAGRAMS

Calculation of Feynman Diagrams from their Small Momentum Expansion

Decomposing Imaginary-Time Feynman Diagrams Using Separable Basis Functions: Anderson Impurity Model Strong-Coupling Expansion

High-frequency asymptotics of the vertex function: diagrammatic parametrization and algorithmic implementation

A fast, dense Chebyshev solver for electronic structure on GPUs

Parallel-in-Time Integration of the Landau-Lifshitz-Gilbert Equation with the Parallel Full Approximation Scheme in Space and Time

Distributed Time, Conservative Parallel Logic Simulation on GPUs

Data-Parallel Algorithms for String Diagrams

Parallel Transport Time-Dependent Density Functional Theory Calculations with Hybrid Functional on Summit

Parallel Computation of Combinatorial Symmetries

Semi-Deterministic and Stochastic Sampling of Feynman Diagrams with 1/N$_f$ Expansions