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 (~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.