The functional renormalization group has become a widely used tool for the analysis of the leading low-temperature correlations in weakly to moderately coupled many-fermion lattice systems. A bottleneck for quantitatively more precise results is that the treatment of the frequency dependence of the flowing interactions is numerically quite demanding. Yet the frequency dependence is needed to compute relevant self-energies and hence for controlled results on the energy scales for ordering or for the quasiparticle properties. Here we explore an approximate parametrization of the frequency dependence of the interaction vertex that is inspired by established simplifications in the theory of superconductivity and that keeps the numerical effort bounded. We demonstrate the validity of the approximation for Cooper pairing problems and apply it to the two-dimensional Hubbard model.
The functional renormalization group has become a widely used tool for the analysis of the leading low-temperature correlations in weakly to moderately coupled many-fermion lattice systems. A bottleneck for quantitatively more precise results is that the treatment of the frequency dependence of the flowing interactions is numerically quite demanding. Yet the frequency dependence is needed to compute relevant self-energies and hence for controlled results on the energy scales for ordering or for the quasiparticle properties. Here we explore an approximate parametrization of the frequency dependence of the interaction vertex that is inspired by established simplifications in the theory of superconductivity and that keeps the numerical effort bounded. We demonstrate the validity of the approximation for Cooper pairing problems and apply it to the two-dimensional Hubbard model.
The treatment of strongly interacting two-dimensional Fermi systems constitutes one of the most challenging problems in the field of condensed matter physics. Many theoretical works focused on the strongly coupled Hubbard model, since it is expected to capture the most important physics of cuprate superconductors. Due to our methodological improvements and understanding of the frequency dependence of the two-particle vertex function, the functional renormalization group combined with the dynamical mean field theory can now be used to study competing correlations in the strongly interacting regime. While limited to moderate interaction strengths, the functional renormalization group describes efficiently systems with a hierarchy of different energy scales and competing correlations. For instance, it provides definite evidence for d-wave superconductivity in the two-dimensional Hubbard model at moderate coupling. In a first project, we study the frequency dependence of the one-particle irreducible vertex function generated by the functional renormalization group flow. The frequency dependence, which becomes singular for strong interactions, appears to be important already for moderate couplings, and it cannot be represented by separate channels each depending on a single linear combination of frequencies. For strongly interacting systems, the dynamical mean field theory captures strong local correlation effects nonperturbatively. With this approximation, we study the impact of local correlations on the magnetic susceptibility. The local dynamics strongly affects the spin response function and its momentum dependence. In contrast to the widely used random-phase approximation with self-energy corrections that predicts Néel antiferromagnetic order, the local vertex corrections favor an incommensurate order similar to the ordering instability predicted by the Fermi surface geometry, as for weakly interacting systems. The dynamical mean field theory is also used as starting point for the functional renormalization group flow. We demonstrate that, due to our improvements for the parametrization of the vertex function, this approach is actually able to access the strong coupling physics. Moreover, we derive a flow scheme that conserves the local contributions and reduces the truncation error of the flow equations. In the strongly interacting regime, we capture strong d-wave pairing correlations driven by magnetic fluctuations with a mechanism similar to the one observed in the weakly interacting system.
Using a leading algorithmic implementation of the functional renormalization group (fRG) for interacting fermions on two-dimensional lattices, we provide a detailed analysis of its quantitative reliability for the Hubbard model. In particular, we show that the recently introduced multiloop extension of the fRG flow equations for the self-energy and two-particle vertex allows for a precise match with the parquet approximation also for two-dimensional lattice problems. The refinement with respect to previous fRG-based computation schemes relies on an accurate treatment of the frequency and momentum dependences of the two-particle vertex, which combines a proper inclusion of the high-frequency asymptotics with the so-called truncated unity fRG for the momentum dependence. The adoption of the latter scheme requires, as an essential step, a consistent modification of the flow equation of the self-energy. We quantitatively compare our fRG results for the self-energy and momentum-dependent susceptibilities and the corresponding solution of the parquet approximation to determinant quantum Monte Carlo data, demonstrating that the fRG is remarkably accurate up to moderate interaction strengths. The presented methodological improvements illustrate how fRG flows can be brought to a quantitative level for two-dimensional problems, providing a solid basis for the application to more general systems.
Using a leading algorithmic implementation of the functional renormalization group (fRG) for interacting fermions on two-dimensional lattices, we provide a detailed analysis of its quantitative reliability for the Hubbard model.In particular, we show that the recently introduced multiloop extension of the fRG flow equations for the self-energy and two-particle vertex allows for a precise match with the parquet approximation also for twodimensional lattice problems.The refinement with respect to previous fRG-based computation schemes relies on an accurate treatment of the frequency and momentum dependences of the two-particle vertex, which combines a proper inclusion of the high-frequency asymptotics with the so-called truncated unity fRG for the momentum dependence.The adoption of the latter scheme requires, as an essential step, a consistent modification of the flow equation of the self-energy.We quantitatively compare our fRG results for the self-energy and momentumdependent susceptibilities and the corresponding solution of the parquet approximation to determinant quantum Monte Carlo data, demonstrating that the fRG is remarkably accurate up to moderate interaction strengths.The presented methodological improvements illustrate how fRG flows can be brought to a quantitative level for two-dimensional problems, providing a solid basis for the application to more general systems.
We present a functional renormalization group (fRG) study of the two dimensional Hubbard model, performed with an algorithmic implementation which lifts some of the common approximations made in fRG calculations. In particular, in our fRG flow; (i) we take explicitly into account the momentum and the frequency dependence of the vertex functions; (ii) we include the feedback effect of the self-energy; (iii) we implement the recently introduced multiloop extension which allows us to sum up all the diagrams of the parquet approximation with their exact weight. Due to its iterative structure based on successive one-loop computations, the loop convergence of the fRG results can be obtained with an affordable numerical effort. In particular, focusing on the analysis of the physical response functions, we show that the results become independent from the chosen cutoff scheme and from the way the fRG susceptibilities are computed, i.e., either through flowing couplings to external fields, or through a “post-processing” contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRG-based computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.
We present a functional renormalization group calculation for the two-dimensional $(t,{t}^{\ensuremath{'}})$-Hubbard model at Van Hove filling. Using a channel decomposition we describe the momentum and frequency dependence of the vertex function in the normal phase. Compared to previous studies that neglect frequency dependencies we find higher pseudocritical scales and a smaller region of d-wave superconductivity. A large contribution to the effective interaction is given by a forward-scattering process with finite frequency exchange. We test different frequency parametrizations and in a second calculation include the frequency dependence of the imaginary self-energy. We also generalize the channel decomposition to frequency-dependent fermion-boson vertex functions.
We study the impact of including self-energy feedback and frequency-dependent interactions on functional renormalization group flows for the two-dimensional Hubbard model on the square lattice at weak to moderate coupling strength. Previous studies using the functional renormalization group had ignored these two ingredients to a large extent, and the question is how much the flows to strong coupling analyzed by this method depend on these approximations. Here we include the imaginary part of the self-energy on the imaginary axis and the frequency dependence of the running interactions on a frequency mesh of 10 frequencies on the Matsubara axis. We find that (i) the critical scales for the flows to strong coupling are shifted downward by a factor that is usually of order 1 but can get larger in specific parameter regions, and (ii) that the leading channel in this flow does not depend strongly on whether self-energies and frequency dependence is included or not. We also discuss the main features of the self-energies developing during the flows.