Partial bosonization of the two-dimensional Hubbard model focuses the functional renormalization flow on channels in which interactions become strong and local order sets in. We compare the momentum structure of the four-fermion vertex, obtained on the basis of a patching approximation, to an effective bosonic description. For parameters in the antiferromagnetic phase near the onset of local antiferromagnetic order, the interaction of the electrons is indeed well described by the exchange of collective bosonic degrees of freedom. The residual four-fermion vertex after the subtraction of the bosonic-exchange contribution is small. We propose that similar partial bosonization techniques can improve the accuracy of renormalization flow studies also for the case of competing order.
Partial bosonization of the two-dimensional Hubbard model focuses the functional renormalization flow on channels in which interactions become strong and local order sets in. We compare the momentum structure of the four-fermion vertex, obtained on the basis of a patching approximation, to an effective bosonic description. For parameters in the antiferromagnetic phase near the onset of local antiferromagnetic order, the interaction of the electrons is indeed well described by the exchange of collective bosonic degrees of freedom. The residual four-fermion vertex after the subtraction of the bosonic-exchange contribution is small. We propose that similar partial bosonization techniques can improve the accuracy of renormalization flow studies also for the case of competing order.
Despite its apparent simplicity, the two-dimensional Hubbard model for locally interacting fermions on a square lattice is widely considered as a promising approach for the understanding of Cooper pair formation in the quasi two-dimensional high-Tc cuprate materials. In the present work this model is investigated by means of the functional renormalization group, based on an exact flow equation for the effective average action. In addition to the fermionic degrees of freedom of the Hubbard Hamiltonian, bosonic fields are introduced which correspond to the different possible collective orders of the system, for example magnetism and superconductivity. The interactions between bosons and fermions are determined by means of the method of “rebosonization” (or “flowing bosonization”), which can be described as a continuous, scale-dependent Hubbard-Stratonovich transformation. This method allows an efficient parameterization of the momentumdependent effective two-particle interaction between fermions (four-point vertex), and it makes it possible to follow the flow of the running couplings into the regimes exhibiting spontaneous symmetry breaking, where bosonic fluctuations determine the types of order which are present on large length scales. Numerical results for the phase diagram are presented, which include the mutual influence of different, competing types of order.
Magnetic and superconducting instabilities in the two-dimensional $t\text{\ensuremath{-}}{t}^{\ensuremath{'}}$-Hubbard model are discussed within a functional renormalization group approach. The fermionic four-point vertex is efficiently parametrized by means of partial bosonization. The exchange of composite bosons in the magnetic, charge-density, and superconducting channels accounts for the increase in the effective couplings with increasing length scale. We compute the pseudocritical temperature for the onset of local order in various channels.
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.
The authors suggest an approximative solution of the two-dimensional Hubbard model close to half-filling. It is based on partial bosonization, supplemented by an investigation of the functional renormalization group flow. The inclusion of the fermionic and bosonic fluctuations leads in lowest order to agreement with Hartree-Fock result or Schwinger-Dyson equation and cures the mean field ambiguity. We compute the temperature dependence of the antiferromagnetic order parameter and the gap below the critical temperature. We argue that the Mermin-Wagner theorem is not practically applicable for the spontaneous breaking of the continuous spin symmetry in the antiferromagnetic state of the Hubbard model. The long distance behavior close to and below the critical temperature is governed by the renormalization flow for the effective interactions of composite Goldstone bosons and deviates strongly from the Hartree-Fock result.
The functional renormalization group (fRG) provides an unbiased framework to analyze competing instabilities in two-dimensional electron systems and has been used extensively over the past decade [1]. In order to obtain an equally unbiased tool to interprete the flow, we investigate the combination of a many-patch, one-loop calculation with higher dimensional bosonization [2] of the resulting low-energy action. Subsequently a semiclassical approximation [3] can be used to describe the resulting phases. The spinless Hubbard model on a square lattice with nearest neighbor repulsion is investigated as a test case. [1] M. Salmhofer and C. Honerkamp, Prog. Theor. Phys. 105, 1 (2001). [2] A. Houghton, H.-J. Kwon, J. B. Marston, Adv.Phys. 49, 141 (2000); P. Kopietz, Bosonization of interacting fermions in arbitrary dimensions, (Springer, Berlin, 1997). [3] H.-H. Lin, L. Balents, M. P. A. Fisher, Phys. Rev. B 56, 65696593 (1997); J. O. Fjaerestad, J. B. Marston, U. Schollwoeck, Ann. Phys. (N.Y.) 321, 894 (2006).
The phases with spontaneously broken symmetries corresponding to antiferromagnetic and d-wave superconducting order in the two-dimensional $t\ensuremath{-}{t}^{'}$-Hubbard model are investigated by means of the functional renormalization group. The introduction of composite boson fields in the magnetic, charge density, and superconducting channels allows an efficient parametrization of the four-fermion vertex and the study of regimes where either the antiferromagnetic or superconducting order parameter, or both, are nonzero. We compute the phase diagram and the temperature dependence of the order parameter below the critical temperature, where antiferromagnetic and superconducting order show a tendency of mutual exclusion.