In order to investigate the validity of the Brueckner and the simple sum of the ladder diagrams approximations to the energy of an infinite system of fermions, we have calculated as a function of density all the terms in the perturbation expansion through the fourth order in the strength of the potential. We have done this for a repulsive, square-well, two-body potential spin-\textonehalf{} fermions. We are able to construct rigorous upper and lower bounds from the coefficients for the value of the ladder approximation and find that the standard solution procedures give reasonably accurate results (within a few percent). We find that the error in the solution as obtained in practice to the equations of the Brueckner approximation is large compared to the size of its departure from the ladder approximation. We further find for low-to-moderate densities and for low-to-moderate potential strengths that the Brueckner approximation both as a sum of a certain class of diagrams and as computed in practice lies above the ladder approximation while the complete perturbation theory lies below it. This result arises from the neglect of the ring diagrams by the Brueckner approximation.
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.
Planar parity-preserving tensor networks (ppTNs) admit a fermionic formulation in which the Gaussian, efficiently contractible limit is deformed by local interactions. We investigate how many-body methods can be used to analyze such contractions in a minimal two-parameter ppTN. Its contraction is simultaneously an interacting fermionic partition function, a loop gas, a deformed toric-code norm, and a quartic Ising model, allowing the same phase diagram to be approached with complementary tools. At its center lies a `topological island', characterized beyond the Gaussian limit by a boundary-twist
Z2 indicator that reduces to Chern-number parity and admits a loop-winding interpretation under duality. Fermionic perturbation theory predicts the interacting phase boundaries, tensor-network numerics establish their critical behavior, and the loop and spin descriptions reveal a self-dual line with a
c=1 four-state-Potts multicritical fixed point.