We describe the first version (v1.0.0) of the code ADG that automatically (1) generates all valid Bogoliubov many-body perturbation theory (BMBPT) diagrams and (2) evaluates their algebraic expression to be implemented for numerical applications . This is achieved at any perturbative order p for a Hamiltonian containing both two-body (four-legs) and three-body (six-legs) interactions (vertices). The automated generation of BMBPT diagrams of order p relies on elements of graph theory, i.e., it is achieved by producing all oriented adjacency matrices of size ( p + 1 ) × ( p + 1 ) satisfying topological Feynman’s rules. The automated evaluation of BMBPT diagrams of order p relies both on the application of algebraic Feynman’s rules and on the identification of a powerful diagrammatic rule providing the result of the remaining p -tuple time integral. The diagrammatic rule in question constitutes a novel finding allowing for the straight summation of large classes of time-ordered diagrams at play in the time-independent formulation of BMBPT. Correspondingly, the traditional resolvent rule employed to compute time-ordered diagrams happens to be a particular case of the general rule presently identified. The code ADG is written in Python2.7 and uses the graph manipulation package NetworkX . The code is also able to generate and evaluate Hartree–Fock-MBPT (HF-MBPT) diagrams and is made flexible enough to be expanded throughout the years to tackle the diagrammatics used in various many-body formalisms that already exist or are yet to be formulated. Program summary Program Title: ADG Program Files doi: http://dx.doi.org/10.17632/6h4xrydwfb.1 Licensing provisions: GPLv3 Programming language: Python2.7 Nature of problem: As formal and numerical developments in many-body-perturbation-theory-based ab initio methods make higher orders reachable, producing and evaluating all the diagrams become rapidly undoable on a handmade basis as both their number and complexity grow quickly, making it prone to mistakes and oversights. Solution method : BMBPT diagrams are encoded as square matrices known as oriented adjacency matrices in graph theory, and then turned into graph objects using the NetworkX package. Checks on the diagrams and evaluation of their time-integrated expression are then done on a purely diagrammatic basis. HF-MBPT diagrams are produced and evaluated as well using the same principle.We describe the first version (v1.0.0) of the code ADG that automatically (1) generates all valid Bogoliubov many-body perturbation theory (BMBPT) diagrams and (2) evaluates their algebraic expression to be implemented for numerical applications . This is achieved at any perturbative order p for a Hamiltonian containing both two-body (four-legs) and three-body (six-legs) interactions (vertices). The automated generation of BMBPT diagrams of order p relies on elements of graph theory, i.e., it is achieved by producing all oriented adjacency matrices of size ( p + 1 ) × ( p + 1 ) satisfying topological Feynman’s rules. The automated evaluation of BMBPT diagrams of order p relies both on the application of algebraic Feynman’s rules and on the identification of a powerful diagrammatic rule providing the result of the remaining p -tuple time integral. The diagrammatic rule in question constitutes a novel finding allowing for the straight summation of large classes of time-ordered diagrams at play in the time-independent formulation of BMBPT. Correspondingly, the traditional resolvent rule employed to compute time-ordered diagrams happens to be a particular case of the general rule presently identified. The code ADG is written in Python2.7 and uses the graph manipulation package NetworkX . The code is also able to generate and evaluate Hartree–Fock-MBPT (HF-MBPT) diagrams and is made flexible enough to be expanded throughout the years to tackle the diagrammatics used in various many-body formalisms that already exist or are yet to be formulated. Program summary Program Title: ADG Program Files doi: http://dx.doi.org/10.17632/6h4xrydwfb.1 Licensing provisions: GPLv3 Programming language: Python2.7 Nature of problem: As formal and numerical developments in many-body-perturbation-theory-based ab initio methods make higher orders reachable, producing and evaluating all the diagrams become rapidly undoable on a handmade basis as both their number and complexity grow quickly, making it prone to mistakes and oversights. Solution method : BMBPT diagrams are encoded as square matrices known as oriented adjacency matrices in graph theory, and then turned into graph objects using the NetworkX package. Checks on the diagrams and evaluation of their time-integrated expression are then done on a purely diagrammatic basis. HF-MBPT diagrams are produced and evaluated as well using the same principle.

arXiv:1809.01187Computer Physics CommunicationsDec 10, 201853 citations

Many-Body Perturbation Theory

Yoshio Kuramoto

Lecture notes in physicsJan 1, 20202 citations

High-order perturbative calculations of nuclear ground states: Automated evaluation of many-body diagrams

Zhen Li, A. Tichai, A. Schwenk +1 more

We advance the many-body perturbation theory (MBPT) calculations of the ground-state energy and radius of closed-shell nuclei beyond third order. Using automated diagram generation and evaluation up to fifth order, we present ground-state properties of selected closed-shell nuclei up to $^{78}$ Ni with two- and three-nucleon interactions derived from chiral effective field theory. A clear convergence trend is observed for the ground-state energy enabling calculations at improved accuracy. We further investigate in detail the decomposition of the fourth-order contributions. For the ground-state energy, the magnitude of the fourth-order contribution is typically less than half of the third order, and a typical cancellation among different classes of diagrams is observed. Finally, we perform a comprehensive comparison between MBPT and non-perturbative in-medium similarity renormalization group (IMSRG) calculations, with the goal to provide insight into many-body uncertainties associated with the IMSRG(2) truncation.We advance the many-body perturbation theory (MBPT) calculations of the ground-state energy and radius of closed-shell nuclei beyond third order. Using automated diagram generation and evaluation up to fifth order, we present ground-state properties of selected closed-shell nuclei up to $^{78}$Ni with two- and three-nucleon interactions derived from chiral effective field theory. A clear convergence trend is observed for the ground-state energy enabling calculations at improved accuracy. We further investigate in detail the decomposition of the fourth-order contributions. For the ground-state energy, the magnitude of the fourth-order contribution is typically less than half of the third order, and a typical cancellation among different classes of diagrams is observed. Finally, we perform a comprehensive comparison between MBPT and non-perturbative in-medium similarity renormalization group (IMSRG) calculations, with the goal to provide insight into many-body uncertainties associated with the IMSRG(2) truncation.

arXiv:2607.00119Jun 30, 20260 citations

^{78}

Ni with two- and three-nucleon interactions derived from chiral effective field theory. A clear convergence trend is observed for the ground-state energy enabling calculations at improved accuracy. We further investigate in detail the decomposition of the fourth-order contributions. For the ground-state energy, the magnitude of the fourth-order contribution is typically less than half of the third order, and a typical cancellation among different classes of diagrams is observed. Finally, we perform a comprehensive comparison between MBPT and non-perturbative in-medium similarity renormalization group (IMSRG) calculations, with the goal to provide insight into many-body uncertainties associated with the IMSRG(2) truncation.