Pseudogap opening in the two-dimensional Hubbard model: A functional renormalization group analysis

Using the recently introduced multiloop extension of the functional renormalization group, we compute the frequency- and momentum-dependent self-energy of the two-dimensional Hubbard model at half filling and weak coupling. We show that, in the truncated-unity approach for the vertex, it is essential to adopt the Schwinger-Dyson form of the self-energy flow equation in order to capture the pseudogap opening. We provide an analytic understanding of the key role played by the flow scheme in correctly accounting for the impact of the antiferromagnetic fluctuations. For the resulting pseudogap, we present a detailed numerical analysis of its evolution with temperature, interaction strength, and loop order.