Highly nonperturbative nature of the Mott metal-insulator transition: Two-particle vertex divergences in the coexistence region

We thoroughly analyze the divergences of the irreducible vertex functions occurring in the charge channel of the half-filled Hubbard model in close proximity to the Mott metal-insulator transition (MIT). In particular, by systematically performing dynamical mean-field theory (DMFT) calculations on the two-particle level, we determine the location and the number of the vertex divergences across the whole coexistence region adjacent to the first-order metal-to-insulator transition. We find that the lines in the parameter space, along which the vertex divergences occur, display a qualitatively different shape in the coexisting metallic and insulating phase, which is also associated to an abrupt jump of the number of divergences across the MIT. Physically, the systematically larger number of divergences on the insulating side of the transition reflects the sudden suppression of local charge fluctuation at the MIT. Further, a systematic analysis of the results demonstrates that the number of divergence lines increases as a function of the inverse temperature beta = (kBT )-1 by approaching the Mott transition in the zero-temperature limit. This makes it possible to identify the zero-temperature MIT as an accumulation point of an infinite number of vertex divergence lines, unveiling the highly nonperturbative nature of the underlying transition.