Determinant Monte Carlo for irreducible Feynman diagrams in the strongly correlated regime

We develop a numerically exact method for the summation of irreducible Feynman diagrams for fermionic self-energy in the thermodynamic limit. The technique, based on the diagrammatic determinant Monte Carlo and its recent extension to connected diagrams, allows us to reach high (similar to 10) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a nontrivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the nonperturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices.