We analyze the interplay of antiferromagnetism and pairing in the two-dimensional Hubbard model with a moderate repulsive interaction. Coupled charge, magnetic, and pairing fluctuations above the energy scale of spontaneous symmetry breaking are treated by a functional renormalization group flow, while the formation of gaps and order below that scale is treated in mean-field theory. The full frequency dependences of the interaction vertices and gap functions are taken into account. We compute the magnetic and pairing gap functions as a function of doping
p and compare with results from a static approximation. In spite of the strong frequency dependences of the effective interactions and of the pairing gap, important physical results from previous static functional renormalization group calculations are confirmed. In particular, there is a sizable doping regime with robust pairing coexisting with N\'eel or incommensurate antiferromagnetism. The critical temperature for magnetic order is interpreted as the pseudogap crossover temperature. Computing the Kosterlitz-Thouless temperature from the superfluid phase stiffness, we obtain a superconducting dome in the
(p,T) phase diagram centered around 15% hole doping.
We describe a low-temperature approach to the two-dimensional half-filled Hubbard model which allows us to study both antiferromagnetism and single-particle properties. This approach ignores amplitude fluctuations of the antiferromagnetic (AF) order parameter and is valid below a crossover temperature
TX which marks the onset of AF short-range order. Directional fluctuations (spin waves) are described by a nonlinear sigma model $(\mathrm{NL}\ensuremath{\sigma}\mathrm{M})$ that we derive from the Hubbard model. The parameters of the $\mathrm{NL}\ensuremath{\sigma}\mathrm{M}$---the spin stiffness and spin-wave velocity---are calculated as a function of the Coulomb repulsion U. The $\mathrm{NL}\ensuremath{\sigma}\mathrm{M}$ is solved by a saddle-point approximation within the
CP1 representation where the N\'eel field is parametrized by two Schwinger bosons. At zero temperature, there is always Bose condensation of the Schwinger bosons, which signals AF long-range order for any value of the Coulomb repulsion. At finite temperature, the AF long-range order is suppressed (in agreement with the Mermin-Wagner theorem), but the AF correlation length remains exponentially large. In the
CP1 representation, the fermion field is naturally expressed as the product of a Schwinger boson and a pseudofermion whose spin is quantized along the (fluctuating) N\'eel field. This allows us to write the fermion Green's function as the product (in direct space) of the Schwinger boson propagator (which is derived from the $\mathrm{NL}\ensuremath{\sigma}\mathrm{M})$ and the pseudofermion propagator. At zero temperature and weak coupling, our results are typical of a Slater antiferromagnet. The AF gap is exponentially small; there are well-defined Bogoliubov quasiparticles (QP's) (carrying most of the spectral weight) coexisting with a high-energy incoherent excitation background. As U increases, the Slater antiferromagnet progressively becomes a Mott-Heisenberg antiferromagnet. The Bogoliubov bands evolve into Mott-Hubbard bands separated by a large AF gap. A significant fraction of spectral weight is transferred from the Bogoliubov QP's to incoherent excitations. At finite temperature, there is a metal-insulator transition between a pseudogap phase at weak coupling and a Mott-Hubbard insulator at strong coupling. Finally, we point out that our results straightforwardly translate to the half-filled attractive Hubbard model, where the $\mathbf{q}=(\ensuremath{\pi},\ensuremath{\pi})$ charge and
q=0 pairing fluctuations combine to form an order parameter with SO(3) symmetry.
Almost all known high temperature superconductors are cuprates, which can be suitably modelled by the two dimensional Hubbard model. To understand the interplay of various long range properties as antiferromagnetism and superconductivity, one can calculate the phase diagram of the Hubbard model in the charge density-temperature plane. This analysis is conveniently carried out by means of exact renormalization group equations that we apply in the formalism of the effective average action. For this purpose, we derive an equivalent version of the Hubbard model that takes the form of a Yukawa theory. From this modified model long range order in various channels can be extracted by simple calculation of bosonic expectation values. We are able to reproduce the main features of the phase diagram of high temperature superconductors. Furthermore, our analysis shows how the Mermin-Wagner theorem can be reconciled with the existence of antiferromagnetic long range order at non vanishing temperature and how the inclusion of different kinds of bosonic fluctuations affect the phase diagram.
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.
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.
We present a consistent fusion of functional renormalization group and mean-field theory which explicitly introduces a bosonic field via a Hubbard-Stratonovich transformation at the critical scale, at which the order sets in. We show that a minimal truncation of the flow equations, that neglects order-parameter fluctuations, is integrable and fulfills fundamental constraints as the Goldstone theorem and the Ward identity connected with the broken global symmetry. To introduce the bosonic field, we present a technique to factorize the most singular part of the vertex, even when the full dependence on all its arguments is retained. We test our method on the two-dimensional attractive Hubbard model at half-filling and calculate the superfluid gap as well as the Yukawa couplings and residual two-fermion interactions in the ordered phase as functions of fermionic Matsubara frequencies. Furthermore, we analyze the gap and the condensate fraction for weak and moderate couplings and compare our results with previous functional renormalization group studies, and with quantum Monte Carlo data. Our formalism constitutes a convenient starting point for the inclusion of order-parameter fluctuations by keeping a full, nonsimplified, dependence on fermionic momenta and/or frequencies.