We describe a Wick ordered functional renormalization group method for interacting Fermi systems, where the complete flow from the bare action of the microscopic model to the effective low-energy action is obtained from a differential flow equation. We apply this renormalization group approach to a prototypical two-dimensional lattice electron system, the Hubbard model on a square lattice. The flow equation for the effective interactions is evaluated numerically on 1-loop level. The effective interactions diverge at a finite energy scale which is exponentially small for small bare interactions. To analyze the nature of the instabilities signalled by the diverging interactions we compute the flow of the singlet superconducting susceptibilities for various pairing symmetries and also charge and spin density susceptibilities. Depending on the choice of the model parameters (hopping amplitudes, interaction strength and band-filling) we find antiferromagnetic order or d-wave superconductivity as leading symmetry breaking instability.
Renormalization-group analysis of the two-dimensional Hubbard model
