Functional renormalization group beyond the perturbative regime

Nils Wentzell

Published Jul 11, 2017• 0 citations

Abstract

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).

DOI: 10.15496/publikation-18318

Related papers

Renormalization group approaches to strongly correlated electron systems

Fabian B. Kugler

Strongly correlated electron systems host a plethora of fascinating physical phenomena and pose formidable challenges in their theoretical analysis. The challenges originate from the inherent complexity of the quantum many-body problem - no classical computer will ever be able to fully simulate these systems - and the lack of an effective single-particle picture, as the strong mutual interactions make it impossible to regard the electrons as independent from each other. As a consequence, most systems of correlated electrons can only be tackled approximately and numerically. In this thesis, we develop a set of numerical methods for strongly correlated electrons, which are inspired by the renormalization group (RG) idea of including degrees of freedom successively from high to low energies. This enables an efficient organization of the diverse fluctuations and is key for an accurate treatment of interacting quantum systems, where collective behavior and composite objects emerge at energy scales far below those of the microscopic constituents. In a first part, we consider the functional renormalization group (fRG), a versatile framework to study the flow of correlation functions upon modulating the underlying action. Though widely used, it has often acted more as a qualitative rather than quantitative method, due a nontransparent approximation induced by truncating the hierarchy of flow equations. We develop an iterative multiloop fRG (mfRG) scheme, which ameliorates this approximation and eliminates many of the drawbacks of fRG experienced hitherto. In particular, it restores the independence of results on the choice of RG regulator and establishes a rigorous relation to the parquet formalism. Furthermore, we show how to derive the flow equations directly from self-consistent many-body relations. This establishes a form of diagrammatic resummations at the two-particle level which circumvents ill-behaved two-particle-irreducible vertices. An application to the prototypical two-dimensional Hubbard model illustrates how our multiloop scheme elevates the fRG approach to correlated electron systems to a quantitative level. Secondly, we employ the numerical renormalization group (NRG), based on the iterative diagonalization of impurity Hamiltonians, in conjunction with the dynamical mean-field theory (DMFT) to describe local correlations in multiorbital systems. Having access to arbitrarily low temperatures and energies, NRG is a unique, real-frequency impurity solver for DMFT. It has been pivotal to the understanding of Hund metals, where strong correlations arise from Hund's rules even at moderate Coulomb repulsion. Building on recent methodological advances, we extend the range of application of DMFT+NRG from orbital-degenerate models to more realistic setups: We first study orbital differentiation in a three-orbital Hund-metal model and unravel key effects of the orbital-selective Mott transition. In a real-materials setting, we then incorporate the bandstructure from density functional theory (DFT) and analyze the archetypal Hund-metal material Sr2RuO4. We particularly follow its RG flow to the Fermi-liquid regime at previously inaccessible low temperatures and generally present DFT+DMFT+NRG as a new computational paradigm for strongly correlated materials. As a side project of our fRG work, we develop an algorithm to count Feynman diagrams from closed many-body relations, which reveals the surprising outcome that totally irreducible contributions are responsible for the factorial growth in the number of diagrams. Additionally, we use NRG to study transport through three-level quantum dots and provide benchmark data for other RG methods, which aim at further describing these systems in nonequilibrium.

Jan 1, 20190 citations

Multiloop functional renormalization group from single bosons

Fraboulet, Kilian, Al-Eryani, Aiman, Heinzelmann, Sarah +2 more

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.

arXiv:2512.11190Dec 12, 20250 citations

Quantitative functional renormalization group description of the two-dimensional Hubbard model

Cornelia Hille, Fabian B. Kugler, Christian Eckhardt +5 more

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.

Jan 1, 202055 citations

Applications of the functional renormalization group to quantum liquids

Sascha Ledowski

The topic of this thesis is the functional renormalization group. We discuss some approximations schemes. Thereafter we apply these approximations to study different fields of condensed matter physics. Generally we have to evaluate an infinite set of vertex functions describing the scattering of particles. These vertex functions get renormalized away from their bare values governed by an infinite hierarchy of flow equations. We cannot expect to actually solve these equations but have to apply a couple of approximations. The aim is to somehow separate relevant contributions from irrelevant ones. One possible scheme opens up if we rescale fields and vertices. Here "relevance" is used in a quantitative way to describe the scaling behaviour of vertices close to a fixed point of the RG. One disadvantage of describing the system in terms of infinitely many vertices is that the majority of these vertices we have to evaluate are not of interest to us. In most cases we are just looking for the self-energy or the two-particle effective interaction. However there might be contributions to the flow of these vertices that are generated by irrelevant vertices. We generally assume that we can express irrelevant vertices in terms of the relevant and marginal ones. Then in turn it should be possible to write the contributions of these irrelevant vertices to the flow of relevant and marginal ones in terms of relevant and marginal vertices as well. We show how this can be achieved by what we term the adiabatic approximation. We now consider weakly interacting bosons at the critical point of Bose-Einstein condensation. As the transition takes place at a finite temperature this temperature defines an effective ultraviolet cut-off. For the investigation of physical properties that depend on momenta smaller than this cut-off it is therefore sufficient to describe the system by a classical field theory. Our central topic here is the self-energy of the bosons and we are able to evaluate it with the full momentum dependence. For small momenta it approaches a scaling form and as the momentum is gradually increased we observe a crossover to the perturbative regime. As a test for the reliability of our expression for the selfenergy we investigate the interaction induced shift of the critical. Our results compare quite satisfactory to the best available estimates for this shift. For the anomalous dimension our approach predicts the correct order of magnitude however with a considerable error. As an improvement we include more vertices into our calculations. Here we observe that our fixed point estimates indeed approach the best known results but this convergence is quite weak. We turn toward systems of interacting fermions. The formulation of the functional renormalization group implicitly requires knowledge of the true Fermi surface of the full interacting system. In general however we can just calculate it a-posteriori from the self-energy. The requirement to flow into a fixed point can be translated into a fine-tuning of the frequency/momentum independent part r_0 of the rescaled 2-point function. We show how this bare value is related to the momentum dependent effective interaction along the complete trajectory of the RG. On the other hand r_0 expresses the difference between the bare and the true Fermi surface. Putting both equations together results into an exact selfconsistency equation for the Fermi surface. We apply our self-consistency equation above to tackle the problem of finding the true Fermi surface of interacting fermions in low dimensions. The most simple non-trivial model with an inhomogeneous Fermi surface is a system of two coupled metallic chains. The process of interband backward scattering leads to a smoothing of the Fermi surface. Of special interest is if the Fermi momenta of the two bands collapse into just one value. We propose the term confinement transition for this behaviour. We bosonize the interband backward scattering by means of a Hubbard-Stratonovich transformation and treat our system as a single channel problem. This bosonization together with the adiabatic approximation allows us to investigate the system even at strong coupling. Within a simple one-loop treatment our method predicts a confinement transition at strong coupling. However taken vertex renormalizations into account we observe that this confinement is destroyed by fluctuations beyond one-loop. Actually we observe how the confined phase can be stabilized by the inclusion of interband umklapp scattering. Thereafter we consider the physically more relevant case of a two-dimensional system of infinitely many coupled metallic chains. Here the Fermi surface consists of two disconnected weakly curved sheets. We are able to repeat the calculations we have performed for our toy model. Within a self-consistent 2-loop calculation indeed signs for a confinement transition at finite coupling strength emerge.

Sep 3, 20080 citations

Quantitative functional renormalization group description of the two-dimensional Hubbard model

Cornelia Hille, Fabian B. Kugler, Christian J. Eckhardt +5 more

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.

arXiv:2002.02733Physical Review ResearchSep 8, 202061 citations

Multiloop functional renormalization group for the two-dimensional Hubbard model: Loop convergence of the response functions

Agnese Tagliavini, Cornelia Hille, Fabian B. Kugler +3 more

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.

arXiv:1807.02697SciPost PhysicsJan 18, 201974 citations