We derive an expansion of the functional renormalization (fRG) equations that treats the frequency and momentum dependencies of the vertices in a systematic manner. The scheme extends the channel-decomposed fRG equations to the frequency domain and reformulates them as a series of linear integral equations in the particle-particle, particle-hole and particle-hole exchange channels. We show that the linearity of the equations offers numerous computational advantages and leads to converged, stable solutions for a variety of Hamiltonians. As the expansion is in the coupling between channels, the truncations that are necessary to making the scheme computationally viable still lead to equations that treat contributions from all channels equally. As a first benchmark we apply the two-loop fRG equations to the single impurity Anderson model. We consider the sources of error within the fRG, the computational cost associated with each, and how the choice of regulator affects the flow of the fRG. We then use the optimal truncation scheme to study the Extended Hubbard Hamiltonian in one and two dimensions. We find that in many cases of interest the fRG flow converges to a stable vertex and self-energy from which we can extract the various correlation functions and susceptibilities of interest.
We derive an expansion of the functional renormalization (fRG) equations that treats the frequency and momentum dependencies of the vertices in a systematic manner. The scheme extends the channel-decomposed fRG equations to the frequency domain and reformulates them as a series of linear integral equations in the particle-particle, particle-hole and particle-hole exchange channels. We show that the linearity of the equations offers numerous computational advantages and leads to converged, stable solutions for a variety of Hamiltonians. As the expansion is in the coupling between channels, the truncations that are necessary to making the scheme computationally viable still lead to equations that treat contributions from all channels equally. As a first benchmark we apply the two-loop fRG equations to the single impurity Anderson model. We consider the sources of error within the fRG, the computational cost associated with each, and how the choice of regulator affects the flow of the fRG. We then use the optimal truncation scheme to study the Extended Hubbard Hamiltonian in one and two dimensions. We find that in many cases of interest the fRG flow converges to a stable vertex and self-energy from which we can extract the various correlation functions and susceptibilities of interest.
The functional renormalization group (fRG) is an established tool in the treatment of correlated electron systems, notably for the description of competing instabilities. In recent years, methodological advancements led to the multiloop extension of the fRG, which systematically includes loop corrections beyond the conventional one-loop truncation and yields a quantitatively accurate description of two-dimensional lattice systems. At the same time, the single-boson exchange (SBE) decomposition of the two-particle vertex has been shown to offer both computational and interpretative advantages paving the way to more affordable approximation schemes. We here apply their combination coined as multiloop SBE fRG to the two-dimensional Hubbard model at weak coupling. After providing a detailed account of the underlying formalism in physical channels, we analyze the results for the frequency- and momentum-dependent vertex functions. We find that the SBE approximation, i.e., neglecting the flow of the multi-boson exchange contributions, accurately reproduces the parquet approximation at loop convergence. The presented algorithmic improvement opens the route for the treatment of more challenging parameter regimes and more realistic models.
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.
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.
We present a reformulation of the functional renormalization group (fRG) for many-electron systems, which relies on the recently introduced single-boson exchange (SBE) representation of the parquet equations [F. Krien, A. Valli, and M. Capone, Phys. Rev. B 100, 155149 (2019)]. The latter exploits a diagrammatic decomposition, which classifies the contributions to the full scattering amplitude in terms of their reducibility with respect to cutting one interaction line, naturally distinguishing the processes mediated by the exchange of a single boson in the different channels. We apply this idea to the fRG by splitting the one-loop fRG flow equations for the vertex function into SBE contributions and a residual four-point fermionic vertex. Similarly as in the case of parquet solvers, recasting the fRG algorithm in the SBE representation offers both computational and interpretative advantages: The SBE decomposition not only significantly reduces the numerical effort of treating the high-frequency asymptotics of the flowing vertices, but it also allows for a clear physical identification of the collective degrees of freedom at play. We illustrate the advantages of an SBE formulation of fRG-based schemes, by computing through the merger of dynamical mean-field theory and fRG the susceptibilities and the Yukawa couplings of the two-dimensional Hubbard model from weak to strong coupling, for which we also present an intuitive physical explanation of the results. The SBE formulation of the one-loop flow equations paves a promising route for future multiboson and multiloop extensions of fRG-based algorithms.
This thesis aims at developing new schemes for the treatment of correlation effects in condensed matter systems using quantum field theoretical approaches. In particular, our goal is to extend the description of correlation physics at the two-particle level. This is necessary for an unbiased treatment of condensed matter systems that exhibit electronic correlations and competing ordering tendencies. In this respect, the functional renormalization group (fRG) approaches have surely contributed substantially over the last years, as they account for all scattering channels and their mutual feedback effects in an unbiased way. In spite of its flexibility, the application of the fRG is limited by its inherent perturbative nature. To go beyond the conventional weak-coupling implementations, we discuss the general idea to extend fRG based computational schemes by using an exactly solvable interacting reference problem as starting point for the RG flow. The systematic expansion around this solution accounts for a non-perturbative inclusion of correlations at both, the one-particle (self-energy) and two-particle (vertex functions) level. The full treatment of the two-particle vertex functions, however, poses a huge limitation to the numerical performance, not only in the fRG, but in several forefront many-body algorithms. In this perspective, we provide a detailed diagrammatic analysis of the frequency and momentum structures of the vertex functions, together with their physical interpretation. This constitutes the basis for sophisticated parametrization schemes. We then explain the technical details necessary for cutting-edge numerical implementations, and further benchmark our ideas using refined implementations of both, the fRG and the parquet approximation (PA).
We derive functional renormalization group (fRG) schemes for Fermi systems which are based on the two-particle irreducible approach to the quantum many-body problem. In a first step, the cutoff is introduced in the non-interacting propagator as it is commonly done in fRG based on one-particle irreducible vertex functions. The most natural truncation of the resulting infinite hierarchy of flow equations is shown to be fully equivalent to self-consistent perturbation theory. An earlier suggested alternative truncation strategy is considered as well. In a second step, the cutoff is introduced in the two-particle interaction. Again two truncation procedures are investigated, one of which was derived before. In the latter, the mean-field solution of the many-body problem is considered as the starting point of the renormalization group flow. We compare the performance and the required numerical resources for solving the coupled flow equations for all the approximate schemes by applying them to the problem of the quantum anharmonic oscillator. In a functional integral representation, this model has a formal similarity to the quantum many-body problem. The perspectives for applying the derived two-particle irreducible fRG approaches to zero- and one-dimensional systems of correlated fermions are discussed.