Full and unbiased solution of the Dyson-Schwinger equation in the functional integro-differential representation

Abstract

We provide a full and unbiased solution to the Dyson-Schwinger equation illustrated for φ4 theory in 2D. It is based on an exact treatment of the functional derivative ∂Γ/∂G of the four-point vertex function Γ with respect to the two-point correlation function G within the framework of the homotopy analysis method (HAM) and the Monte Carlo sampling of rooted tree diagrams. The resulting series solution in deformations can be considered as an asymptotic series around G=0 in a HAM control parameter c0G, or even a convergent one up to the phase transition point if shifts in G can be performed (such as by summing up all ladder diagrams). These considerations are equally applicable to fermionic quantum field theories and offer a fresh approach to solving functional integro-differential equations beyond any truncation scheme.