Dyson-Schwinger equations in zero dimensions and polynomial approximations

The Dyson-Schwinger (DS) equations for a quantum field theory in $D$-dimensional space-time are an infinite sequence of coupled integro-differential equations that are satisfied exactly by the Green's functions of the field theory. This sequence of equations is underdetermined because if the infinite sequence of DS equations is truncated to a finite sequence, there are always more Green's functions than equations. An approach to this problem is to close the finite system by setting the highest Green's function(s) to zero. One can examine the accuracy of this procedure in $D=0$ because in this special case the DS equations are just a sequence of coupled polynomial equations whose roots are the Green's functions. For the closed system one can calculate the roots and compare them with the exact values of the Green's functions. This procedure raises a general mathematical question: When do the roots of a sequence of polynomial approximants to a function converge to the exact roots of that function? Some roots of the polynomial approximants may (i) converge to the exact roots of the function, or (ii) approach the exact roots at first and then veer away, or (iii) converge to limiting values that are unequal to the exact roots. In this study five field-theory models in $D=0$ are examined, Hermitian $\phi^4$ and $\phi^6$ theories and non-Hermitian $i\phi^3$, $-\phi^4$, and $-i \phi^5$ theories. In all cases the sequences of roots converge to limits that differ by a few percent from the exact answers. Sophisticated asymptotic techniques are devised that increase the accuracy to one part in $10^7$. Part of this work appears in abbreviated form in Phys.~Rev.~Lett.~{\bf 130}, 101602 (2023).