Superconductivity in the attractive Hubbard model: functional renormalization group analysis
Abstract
We present a functional renormalization group analysis of superconductivity in the ground state of the attractive Hubbard model on a square lattice. Spontaneous symmetry breaking is treated in a purely fermionic setting via anomalous propagators and anomalous effective interactions. In addition to the anomalous interactions arising already in the reduced BCS model, effective interactions with three incoming legs and one outgoing leg (and vice versa) occur. We accomplish their integration into the usual diagrammatic formalism by introducing a Nambu matrix for the effective interactions. From a randomphase approximation generalized through use of this matrix we conclude that the impact of the effective interactions is limited, especially considering the effective interactions important for the determination of the order parameter. The exact hierarchy of flow equations for oneparticle irreducible vertex functions is truncated on the twoparticle level, with higherorder selfenergy corrections included in a scheme proposed by Katanin. Using a parametrization of effective interactions by patches in momentum space, the flow equations can be integrated numerically to the lowest scales without encountering divergences. Momentumshell as well as interactionflow cutoff functions are used, including a small external field or a large external field and a counterterm, respectively. Both approaches produce momentumresolved order parameter values directly from the microscopic model. The size of the superconducting gap is in reasonable agreement with expectations from other studies.
1 Introduction
Renormalization group (RG) techniques facilitate the systematic investigation of interacting manyparticle systems even in situations where perturbation theory breaks down. A central idea underlying many renormalization group methods, pioneered by Wilson [1], is to take into account the degrees of freedom successively, ordered by their kinetic energy. Formally, this is accomplished by integrating a set of differential equations called flow equations, where some energy scale (usually a cutoff) plays the role of the flow parameter. The RG flow generates a family of effective actions which interpolates smoothly between the bare action of the system and the final effective action where all degrees of freedom have been integrated out.
The first significant applications of RG methods to interacting Fermi systems in more than one dimension were undertaken by mathematicians aiming at a rigorous control of certain systems and their properties [2, 3, 4]. The usefulness of RG concepts for Fermi systems was soon recognized also by physicists [5, 6]. RG calculations for Fermi systems are often termed ”functional” as they involve the flow of functions instead of just a few running couplings. An early computational success of the functional RG (fRG) was unambiguous evidence for dwave superconductivity in the repulsive twodimensional Hubbard model [7, 8, 9].
Of particular interest in interacting Fermi systems is the spontaneous breaking of symmetries. Formally, such a symmetry breaking occurs if a product of two fermionic fields – a bilinear – which does not appear in the original Hamiltonian acquires a nonzero expectation value. This bilinear can be replaced by a bosonic field through a HubbardStratonovich transformation, which simplifies the treatment of the order parameter and its fluctuations [10]. However, these simplifications come at the price of introducing a bias towards breaking the symmetry in the channel fixed by the HubbardStratonovich field. Working with the original fermionic degrees of freedom this bias can be avoided. However, the symmetrybroken phase is not easily accessible by a purely fermionic RG flow. Calculations starting from the symmetric phase run into divergences at a finite energy scale, that is before all degrees of freedom have been integrated out.
Most previous functional RG studies of interacting Fermi systems are based on a flow equation for the generator of the oneparticle irreducible (1PI) vertex functions [11]. This equation is equivalent to an infinite hierarchy of flow equations for the 1PI vertex functions. Truncating the hierarchy in a weak coupling expansion, the RG flow has been computed for various lowdimensional Fermi systems [8, 12, 13, 14, 15, 16]. Instead of the common choice of a momentum cutoff, it was found to be advantageous in certain situations to use as flow parameter the physical temperature [18], the strength of the interaction [19], or a cutoff in frequency space [20, 21].
Spontaneous symmetry breaking within the functional RG framework has frequently been studied by introducing bosonic fields for the order parameter, leading to a coupled flow of bosons and fermions [22, 23, 24]. A simple way to access spontaneous symmetry breaking in a purely fermionic setting is to combine the fRG with a meanfield calculation for the lowenergy degrees of freedom [17]. However, order parameter fluctuations are thereby discarded. To continue the fermionic functional RG flow into the symmetrybroken phase, an improved truncation of the hierarchy with an efficient treatment of selfenergy corrections, first formulated by Katanin [25], is very promising. Employing this improvement, the BardeenCooperSchrieffer (BCS) meanfield model can be solved exactly already by a loworder truncation [26], if a small external field is used to trigger the symmetrybreaking. The vaunted divergence is then regularized, but there is still a flow regime with strong coupling. The method was extended to the breaking of a discrete symmetry and finite temperatures [27]. Employing the interaction flow [19], the external field can be compensated by a counterterm, thus enabling the study of firstorder phase transitions which were previously out of the method’s reach [28]. Furthermore, this approach yields results at zero external field without having to pass through a strong coupling regime in the case of discretesymmetry breaking. Strong attractors for the flows of various counterterm strengths were found.
In this work, we apply the fermionic fRG to the attractive Hubbard model, using both the externalfield as well as the counterterm approach. The repulsive Hubbard model [29, 30, 31] was originally introduced to study ferromagnetism of itinerant electrons, but has since been employed to elucidate many aspects of strong correlations in solidstate physics. The attractive version exhibits swave superconductivity under certain conditions, including the case we study here, namely for a square lattice at zero temperature and weak or moderate coupling [32]. It is a good testbed because the correct order parameter is affected by fluctuations even in the weak coupling limit [33], while the model is fairly wellstudied, allowing comparisons of the results obtained using the new method. Lattice fermions with attractive interactions can nowadays be realized by cold atoms in optical lattices, and swave superfluidity in such systems has been observed very recently [34].
Our goals in this work are threefold. First, we would like to study the role of anomalous effective interactions with three ingoing particles and one outgoing particle (or vice versa), which arise in systems with both Cooper and forward scattering. We find that these interactions have a limited impact. Second, we would like to determine whether the application of the Katanintruncated fRG to the Hubbard model works as well as for meanfield models. For the momentumshell approach with a small external field, we obtain reasonable flows resulting in order parameter values comparable to the literature. It is possible to choose the strength of the external field two orders of magnitude below the final order parameter value without running into artificial divergences. For the interactionflow approach, we find remnants of the strongattractor behavior known from the meanfield models. Third, we compute the superconducting order parameter in the ground state of the attractive Hubbard model as a function of coupling and filling.
This paper is structured as follows. In section 2, the attractive Hubbard model is introduced formally. The formalism permitting the study of symmetry breaking in general and the anomalous effective interactions in particular is detailed in section 3. It is employed in the context of a generalized randomphase approximation (RPA) of the model’s behavior in section 4. Section 5 presents the result of the fRG studies, containing momentumshell and interaction flows as well as a study of the order parameter strength and a comparison to earlier calculations.
2 Attractive Hubbard model
We discuss the attractive Hubbard model on a twodimensional square lattice specified by the Hamiltonian
(1) 
denotes the number of lattice points, is the onsite attraction, and takes all values in the first Brillouin zone of the square lattice. The dispersion relation
(2) 
contains a probability amplitude for nearestneighbor hopping, a probability amplitude for nextnearestneighbor hopping, and a chemical potential . An external field
(3) 
is added to to trigger the instability towards superconductivity. This field is either kept small or turned off gradually in the fRG flow to approximate the situation with spontaneous symmetry breaking.
3 Nambu formalism
3.1 Functional integral and field substitution
We study the partition function in a functional integral representation
(4) 
within the Matsubara formalism, employing the grand canonical action
(5)  
where and are the Grassmann fields matching and , respectively. Furthermore, and volume as well as temperature factors are implicitly contained in the summation symbol , a convention we also adopt in the following. We employ the substitution
(6) 
to obtain
(7)  
3.2 Propagator and effective interaction
In (7), the external field appears on the offdiagonal of a matrix specifying the action’s quadratic part. According to the Dyson equation, inverting this matrix after subtracting the selfenergy leads to the propagator or oneparticle Green’s function
(8)  
(9) 
where and is the normal selfenergy. denotes the superconducting order parameter and contains the anomalous selfenergy as well as the external field. To simplify the diagrammatics, we antisymmetrize the interaction part
(10) 
where , if used as an index, is short for and and are exchanged in comparison with (7). If labels the rows and the columns of a matrix to the Nambu base , the antisymmetrized bare interaction reads
(11) 
By employing the Nambu field substitution (6), a formalism containing four potentially different normal (particlenumber conserving) and twelve potentially different anomalous (particlenumber nonconserving) effective interactions arises. Here, as we concentrate on conventional wave pairing, we require time reversal invariance to hold. This implies the following behavior of an effective interaction component under flipping of all its Nambu indices: if an odd number of the Nambu indices are , the sign of the effective interaction’s real part is flipped; otherwise, the sign of the imaginary part is flipped. Exploiting this symmetry, we assign symbols to the various effective interactions:
(12) 
The are collectively referred to as effective interactions as they have three incoming and one outgoing leg or vice versa if drawn according to the field representation of (5). They are generated in the symmetrybroken phase by diagrams of the type in figure 1.
is a anomalous effective interaction arising only in the symmetrybroken phase. , , and are normal effective interactions. Only and have nonzero bare values.
4 Meanfield and generalized RPA
Here we employ a meanfield approximation of the selfenergy and a generalized randomphase approximation (RPA) of the effective interactions. This represents a useful first step toward the fRG treatment, as it allows us to assess the basic behavior and the frequency and momentumdependence of normal and anomalous vertices. In particular, we find that anomalous 3+1 vertices have a sizable feedback on certain effective interactions.
4.1 Gap equation
The system’s selfenergy can be calculated approximatively by resumming a subset of perturbation theory diagrams as in figure 2.
This subset contains all diagrams corresponding to bubbles, sunrise diagrams, and trees of either or both. It is equivalent to selfconsistent HartreeFock theory. If we were to restrict the Hubbard Hamiltonian to only include the Cooper channel and forward scattering processes with zero momentum transfer, the resulting diagrammatic selfconsistency equation would be exact. Analogous statements can be made about the oneloop fRG treatment for such a restricted Hamiltonian[35]. The selfenergy equation can be expanded into two distinct parts,
The Matsubara sums are performed analytically. The result is
(13)  
(14) 
In this approximation, both and are constant in momentum space if the external field is. The Hartree contribution can be absorbed into the chemical potential and is subsequently ignored.
4.2 BetheSalpeter equation
A subset of all diagrams of the perturbation expansion for the effective interaction can be resummed into a BetheSalpeter equation as shown in figure 3.
This corresponds to a generalization of the randomphase approximation in the sense that particleparticle ladders in the representation are included in addition to the usual particlehole bubbles. Symbolically, this BetheSalpeter equation reads
(15)  
where, as a parameter, is short for . By the introduction of the loop matrix
(16) 
and by exploiting momentum and energy conservation while setting , (15) can be rewritten as
(17) 
where the Nambu indices of and are organized into pairs as in (11). only depends on in this approximation, as can be seen by solving for , which yields
(18) 
4.3 Numerical setup
The Matsubara sums in (18) can be performed analytically. The integration over the Brillouin zone and the matrix inversion on the righthand side of (18) is accomplished numerically. In the following, a system with , , , and is assumed. This choice of parameters corresponds to an approximately quarterfilled system far below its meanfield critical temperature. is chosen to be real. We mostly employ an external field . This keeps the numerics wellbehaved while having no appreciable influence on except near .
4.4 effective interactions
The are generated by diagrams as in figure 1. In the reduced meanfield models studied previously these vertices remain zero due to the restricted momentumstructure of the bare interaction where only particle pairs with total momentum equal to zero interact. However, in the general case with scattering of particle pairs with arbitrary total momentum they have to be considered. Due to the appearance of an anomalous propagator figure 1, these diagrams and hence the are nonzero only in the symmetrybroken phase. To obtain an overview of the momentum and frequency dependence of the effective interactions, we study the graphs in figure 4.
The righthand graphs show that the effective interactions exhibit no divergence in the channel. The at larger frequencies are suppressed by decreasing the external field. From the figure’s lefthand diagram, we learn that the ’s moduli are large around and around for .
4.5 Large interaction vertices and Goldstone boson
The instability toward superconductivity gives rise to a growth of certain combinations of the interaction vertices. From previous fRG studies on reduced meanfield models [26, 27] and from the solution of the BetheSalpeter equation, the meanfield picture of the leading components for zero total incoming momentum and frequency has been obtained as follows. There is one strongly growing combination of effective interactions that occurs in the flow equation for the anomalous selfenergy. It typically reaches its maximum at the critical scale or temperature and does not diverge below. In the BCS meanfield model, this combination is given by for real [26]. Another combination of the effective interactions, , diverges also below . This divergence is associated with the massless phase mode dubbed Goldstone boson present if a continuous symmetry, in our case a global symmetry, is broken. For cases where no continuous symmetry is broken as for the halffilled chargedensitywave meanfield model with a twofold degenerate ground state [27], the cutoffdependence of the effective interaction is comparable to the cutoff dependence of the combination of normal and anomalous effective interactions of a BCS meanfield model.
For the Hubbard interaction studied in this paper, the decay of these large vertices away from zerototalmomentum or frequency becomes important, as other channels not present in the meanfield models might pick up large contributions from the symmetrybreaking channel. In this section, we investigate and with regard to their momentum and energy dependence as well as to the impact of the effective interactions.
4.5.1 Momentum dependence
We analyze the momentum dependence of and at . We fit
(19) 
to data obtained from numerical calculations for small up to intermediate , see figure 5.
For the nondivergent part , the fit reproduces the numericallycalculated behavior. For small momenta, the agreement is better, while for larger momenta, deviations become larger, but vanish along the axes. This is due to an anisotropy in the numerical data which is not captured by the approximation (19).
For the divergent part , we note that as obtained from the fitting procedure is proportional to . This is confirmed by calculations for and . For intermediate momenta along the diagonal (not shown), the fit is not optimal as deviations of up to occur. The situation is better for large momenta, where deviations of about are common. The peak at is wellcaptured by the approximation (19), which permits us to confirm the dispersion relation for the lowlying Goldstone mode in the next section.
We consider the impact of the effective interactions. By setting loops combining an anomalous and a normal propagator to zero in (15), the effective interactions are set to zero. The modulus
(20) 
of the relative change of a given effective interaction combination upon switching on the effective interactions is plotted for and in figure 6.
We note that the magnitude of this change for the combination largely follows the magnitude of the effective interactions in the Brillouin zone. It varies between and . The combination corresponding to the system’s phase mode does not change to numerical accuracy (note the scale on the axis).
4.5.2 Energy dependence
The energy dependence of and is to be included in the approximate formula (19) by adding to the denominator, i.e. by considering the approximation function
(21) 
We restrict ourselves to the real part of the effective interactions for brevity. Contemplating the lefthand graph of figure 7,
we become aware of a shoulder close to zero frequency, which is not captured by (21). This implies a small nonlinearity of the dispersion relation of the amplitude mode.
The ansatz (21) reproduces the numericallyobtained results for upon properly adjusting the parameters, as the righthand part of figure 7 illustrates. The agreement is good for small Matsubara frequencies. In the inset, we see that the largefrequency behavior is not captured optimally, which is due to our choice of , having taken it from the fit in figure 5. However, the exponent for is confirmed by the fit. This together with the momentumdependence fit in section 4.5.1 confirms the linearity of the dispersion relation of the Goldstone mode connected to .
To understand the influence of the () on the energy dependence of and , we plot the relative change at as in (20), see figure 8.
We glean that for , the only appreciable change occurs at and is quite small in the limit . For , the converse is true: Only at is there no change, while even for small , a sizable change is observed which is only suppressed at large . This implies that it would not be prudent to neglect the , while the smallness of their influence on the part of suggests that their quantitative impact on the superconducting gap obtained by an fRG calculation is likely small. The latter supposition is motivated by the observation that drives the flow of the gap, as already shown for the BCS meanfield model in [26].
5 Functional renormalization group
In this section, we describe the fRG treatment of the attractive Hubbard model on the twodimensional square lattice at zero temperature. Our main goals are to study the flows produced as well as to obtain reliable quantitative results for the order parameter. Compared to the partial resummation in the last section, the fRG takes into account all oneloop diagrams. The consequential increase in numerical complexity is balanced by adopting a cruder discretization of momentum and frequency space. By taking into account all oneloop diagrams, the pairing channel is renormalized by other processes with nonzero total momentum such as spin and charge fluctuations. We will find that this causes a reduction of the superconducting order parameter with respect to the meanfield result. Furthermore, the fRG treatment allows for a dependence of the selfenergy around the Fermi surface.
5.1 Momentumshell and interactionflow procedures
We introduce a cutoff into our action by substituting
(22) 
We call
(23) 
momentumshell cutoff function [1] and
(24) 
interactionflow cutoff function [19]. The former name is reflected in the expansion of the system by an infinitesimal shell in momentum space if . The latter name can be understood by scaling the fields in (5) to , causing the noninteracting part of the action to be cutoffindependent and the interacting, quartic part to be proportional to . For the momentumshell flow, setting provides a solvable, noninteracting starting point. For , the original, fully interacting system is recovered. For the interaction flow, suppresses the interaction according to the rescaling of the fields outlined above and provides a solvable starting point. recovers the original, fully interacting system.
Employing the momentumshell cutoff, we introduce a small external field to select the symmetrybroken phase as the endpoint of the fRG calculation. This means our results in this case apply to a situation with spontaneouslybroken symmetry only in the limit . For the interactionflow cutoff, we add zero in the form to the quadratic part of the action, including in the noninteracting part and in the interacting part of the action. This approach has been shown to yield symmetrybroken results at zero external field [28] for .
5.2 fRG equations in the Katanin truncation
The flow equations which we employ are shown in figure 9,
where the singlescale propagator is used. The name can be understood by noting that because of the differentiation acting on , lives only on the energy scale for the momentumshell flow. These flow equations can be obtained by truncating according to Katanin [25] the infinite hierarchy of oneparticle irreducible (1PI) flow equations of [11]. Solving them yields the exact order parameter and effective interaction in the thermodynamic limit for any model with only twoparticle interactions which is meanfield exact [27, 36, 26].
Symbolically, the flow equation for reads
(25)  
where , , and . Note that the last two terms correspond to the second diagram pictured above, as the antisymmetrization is denoted explicitly in the symbolic equation. Prefactors arising from Fourier transformations are absorbed into the summation symbols.
The symbolic form of the flow equation for the selfenergy reads
(26) 
Here, the singlescale propagator is rewritten into a form better suited for numerical evaluation [27].
5.3 Numerical results
In the following, we simplify (25) and (26) and subsequently present results stemming from a numerical quadrature of the simplified flow equations. While we focus on the momentum and frequency structure of the effective interaction in section 4, in this section we are mainly interested in determining the order parameter beyond meanfield theory, taking into account all oneloop diagrams in the fRG. However, this comes at the price of an increase in numerical complexity, which we balance by adopting a cruder discretization of the effective interaction. We first drop the frequencydependence of . Consequently, no energy dependence of is generated in the flow according to (26). We can therefore do all Matsubara sums analytically. We furthermore drop the normal part of the selfenergy, thus neglecting Fermi surface shifts. All quantities can be considered real in our approximation. A parametrization taking into account partially the momentum dependence of the effective interaction is the socalled “patching” [9, 37]. The Brillouin zone is split into patches as shown in figure 10.
The functional value of the effective interaction at the midangle of the patch and on the Fermi surface is assumed to be also valid for all other momenta in the patch. Thereby, only the momentum dependence along the Fermi surface is taken into account. This is motivated by the finding that the dependence perpendicular to the Fermi surface is irrelevant in the RG sense for the momentumshell flow without symmetrybreaking ([5, 6]). We specify the momentum arguments of as triples of patch indices, so refers to any with in patch , in patch , and in patch . Note that to obtain momenta in the representation from these patch numbers, signs must be adjusted according to (6).
5.3.1 Momentumshell flows
We replace the step function from (23) by
(27) 
This simplifies the numerical calculation significantly, although one could suspect that the sharp cutoff might eliminate an integration. However, in the Katanin truncation, terms proportional to the derivative of the order parameter arise in addition to the terms proportional to the derivative of the cutoff function. Note that our results do not depend on in a reasonable range, which adjusts the sharpness of the cutoff.
In figure 11,
we plot flows of the order parameter in the first quadrant of the Brillouin zone. We note that the flow saturates for the two external field strengths employed. However, if we employ too small an external field, the effective interaction flows diverge and the method breaks down. Fortunately, the external field can be chosen at least a hundred times smaller than the final value of the superconducting gap. In calculations for meanfield models [27], the modification of the order parameter due to the external field is less than in this case.
From the resummation studies in section 4, we have learned that the Goldstone phase mode dominates for zero momentum transfer in the Nambu formalism, which is equivalent to zero total momentum if we transform back to the original fields according to (6). In figure 12, we plot several flows from this channel.
They dominate in magnitude all other effective interaction flows, including those not shown in the figure, and exhibit a spread of under variation of the momentum combination. They are quite similar to the flows found in [26]. Note that the maximal final value modulus increases stronger than proportionally to . This is a precursor of the divergence to arise for even smaller .
Complementing the phase mode is the effective interaction combination driving the order parameter. Flows of are plotted in figure 13.
The spread of the final values increases with decreasing , while the average remains approximately constant. The graphs resemble those found when a discrete symmetry is broken [27]. However, the monotonicity below the critical scale which had been established analytically for the meanfield model is violated by one of the flows in the righthand figure. Overall, all momentumshell flows exhibited resemble the flows for the meanfield models.
5.3.2 Interaction flows
In this section, we employ the interaction flow cutoff (24) with an initial anomalous selfenergy and a counterterm which we include in the bare propagator. In figure 14,
we show examples of the flow of the superconducting order parameter for various counterterms. We note that if we could treat the flow equations without truncation and further approximations, all counterterms should produce the same final selfenergy and vertices. Indeed, the countertermfRG for a halffilled chargedensitywave meanfield model exhibited a single strong attractor for the flow of the selfenergy [28].
For the more general model studied here, the truncation error becomes noticeable, and for certain values of the pairing counterterm, the interaction vertices still diverge. The flows are terminated if the maximal effective interaction exceeds four times the bandwidth, i.e., if the flow leaves the weakcoupling regime, or if reaches . We see in the left diagram that the flows no longer converge onto a single strong attractor as for the meanfield case [28]. Studying the final 5 of the flow shown in the right plot of figure 14 reveals that there exists a set of small counterterms for which the flows cross and leave the weakcoupling regime, and a minimal counterterm for which the effective interactions remain within the above limit up to . The flows for counterterms whose strength is just beyond this minimal counterterm terminate closer to each other than do the flows for larger counterterms. The minimal counterterm is approximately twice as strong as the final order parameter value obtained in the minimal counterterm flow. The behavior described above is interpreted as being analogous to the strongattractor behavior in the meanfield case. The overestimation of the final order parameter value for larger counterterms is attributed to the overestimation of the order parameter in the flow due to the neglect of the order parameter’s energy dependence and momentum dependence perpendicular to the Fermi surface. Consequently, we choose the minimal order parameter obtainable by terminating the flow at in the weakcoupling regime as our approximation for the physical order parameter. Another way to argue for the so chosen counterterm is that the continuous symmetry breaking necessitates a divergent vertex exactly at the end of the flow in order to render the Goldstone mode correctly. This infinite value is suppressed by our approximation, but the flow for the chosen counterterm yields the maximal final value for the Goldstone mode within the above limit. Finally, note that the order parameter obtained by these rules is approximately smaller than the value obtained with the momentumshell method in section 5.3.1, which is enhanced by a finite, albeit small, external field.
5.3.3 Order parameters
figure 15
shows the magnitude of the superconducting order parameter versus the chemical potential as calculated using the interaction flow method including a counterterm as described in the previous section. For small bare coupling , is maximal if the van Hove points lie on the bare Fermi surface, which is the case for . For larger bare coupling, the maximum gradually shifts to smaller . We compare the results from the fRG to results from meanfield theory in figure 16.
We note that the suppression of the meanfield values by the fluctuations taken into account by the fRG varies between a factor 1.6 and 2.0. This is quite similar to values found by secondorder perturbation theory [33]. Table 1

U [t] fRG fRG and MF 1.0 0.02 0.02 1.5 0.11 0.07 2.0 0.22 0.15 2.5 0.35 0.26 3.0 0.49 0.38
compares our results to results from a combination of fRG and meanfield theory [38] for finite . While the values are of the same order of magnitude, the results obtained by employing only the fRG are larger by between .
6 Conclusions
We have applied the fermionic oneparticle irreducible (1PI) version of the functional renormalization group (fRG) in the Katanin truncation to the symmetrybroken phase of a twodimensional, nonmeanfield model exhibiting superconductivity at zero temperature. The attractive Hubbard model has been chosen as the prototype system because it is sufficiently wellstudied to provide a good testbed for new manybody methods. In this case, particlehole as well as particleparticle loops are important for the symmetrybreaking properties of the system. The fRG takes both into account on an equal footing.
Our analysis has proceeded in three steps. First, the usual Nambu formalism has been extended to cover anomalous effective interactions featuring an odd number of incoming legs, whose behavior and importance was a focus of our interest. We have found that for the case of superconductivity, the effective interaction had to be extended into a matrix containing certain redundancies due to symmetries. This extension followed from a substitution of the fields in which the model is originally formulated by Nambu fields which are chosen specifically to produce an action of canonical shape while encompassing singlet Cooper pairs in the quadratic part. effective interactions were found to arise naturally within this formalism.
In the second step, a study of the system has been performed by means of a resummation of a subset of all perturbation theory diagrams. The strength of the symmetry breaking was determined by a meanfield gap equation. A class of diagrams containing bubble chains, ladders, and all chains comprising both was resummed into a BetheSalpeter equation. In the context of this equation, the impact of the effective interactions was studied by comparing effective interaction values calculated with and without the effective interactions. This impact was found to be zero for the zerofrequency part of the phase mode, and to be dependent on the external field for nonzero frequency for both phase and amplitude mode.
In the third step, the fRG has been numerically applied to the attractive Hubbard model. Two distinct variants of the fRG were employed. First, a momentumshell cutoff was used together with an external field providing a finite initial order parameter. The Goldstone mode divergence is regularized by the field. Despite the approximate nature of the flow and the explosive character of the fRG differential equation, flows reaching down to zero scale were found employing external field strengths two orders of magnitude smaller than the final order parameter values. The relative error of the results for the order parameter due to such an external field was estimated to be at most , judging from earlier calculations for meanfield models. For even smaller external field values, a regime of diverging flows was found. Second, we have calculated interaction flows, balancing an external field with a counterterm increasing continuously with the flow. We have noted that vestiges of strongattractor behavior found previously for a meanfield model are also present for the attractive Hubbard model flows. They permit the determination of order parameter values which are comparable to results from the literature. While these values proved to be smaller than values calculated by meanfield theory alone, as expected, they are generally slightly larger than results obtained by combining symmetricphase Wickordered fRG calculations and meanfield theory. This deviation is likely due to small differences between the 1PI and Wickordered fRG schemes, since symmetricphase 1PI calculations already yield larger critical scales than the corresponding Wickordered calculations. In summary, the 1PI fRG in the Katanin truncation was found to yield fairly accurate quantitative results for weak and moderate interaction strengths (up to ). We reckon that the diverging flows present for very small external fields in our calculations can be pushed back by improving the resolution of the discretization in both frequency and momentum space. This would furthermore improve the accuracy of the results. The method may also be employed to study symmetry breaking in interactingfermion models for which a stronger momentum dependence of the order parameter is expected, for example to study the repulsive Hubbard model. The method can also be applied to situations with competing order. The interactionflow scheme with possible symmetrybreaking in two different channels has recently been tested in a onedimensional model with competing longrange correlations. There it allowed for a correct determination of a quantum critical point between a chargegapped phase and a compressible phase with dominant superconducting correlations [39]. Finally, it will also be interesting to analyze how nonzero temperature affects the renormalization of the meanfieldlike instability encountered in the flows.
Acknowledgments
We thank Andrey Katanin, Jutta Ortloff, Matthias Ossadnik, Julius Reiss,
Manfred Salmhofer, and Philipp Strack for pleasant and valuable
discussions, and
Pawel Jakubczyk for providing very useful comments on the
manuscript.
References
 [1] K. G. Wilson. Renormalization Group and Critical Phenomena. II. PhaseSpace Cell Analysis of Critical Behavior. Phys. Rev. B, 4:3184, 1971.
 [2] G. Benfatto and G. Gallavotti. Renormalizationgroup approach to the theory of the fermi surface. Phys. Rev. B, 42:9967, 1990.
 [3] J. Feldman and E. Trubowitz. Perturbation theory for many fermion systems. Helv. Phys. Acta, 63:156, 1990.
 [4] J. Feldman and E. Trubowitz. The flow of an electronphonon system to the superconducting state. Helv. Phys. Acta, 64:213, 1991.
 [5] R. Shankar. Renormalization group for interacting fermions in d 1. Physica A, 177:530, 1991.
 [6] R. Shankar. Renormalizationgroup approach to interacting fermions. Rev. Mod. Phys., 66:129, 1994.
 [7] C. J. Halboth and W. Metzner. Renormalizationgroup analysis of the twodimensional hubbard model. Phys. Rev. B, 61:7364, Mar 2000.
 [8] C. Honerkamp, M. Salmhofer, N. Furukawa, and T. M. Rice. Breakdown of the LandauFermi liquid in two dimensions due to umklapp scattering. Phys. Rev. B, 63:35109, 2001.
 [9] D. Zanchi and H. J. Schulz. Weakly correlated electrons on a square lattice: A renormalization group theory. Europhys. Lett., 44:235, 1998.
 [10] V.N. Popov, Functional integrals and collective excitations (Cambridge University Press, Cambridge, 1987).
 [11] M. Salmhofer and C. Honerkamp. Fermionic Renormalization Group Flows —Technique and Theory—. Prog. Theor. Phys., 105:1, 2001.
 [12] S. Andergassen, T. Enss, and V. Meden. Kondo physics in transport through a quantum dot with Luttingerliquid leads. Phys. Rev. B, 73:153308, 2006.
 [13] S. Andergassen, T. Enss, V. Meden, W. Metzner, U. Schollwock, and K. Schonhammer. Functional renormalization group for Luttinger liquids with impurities. Phys. Rev. B, 70:075102, 2004.
 [14] T. Enss, V. Meden, S. Andergassen, X. BarnabeTheriault, W . Metzner, and K. Schonhammer. Impurity and correlation effects on transport in onedimensional quantum wires. Phys. Rev. B, 71:155401, 2005.
 [15] H. C. Fu, C. Honerkamp, and D.H. Lee. Renormalization group study of the electronphonon interaction in highT cuprates. Europhys. Lett., 75:146, 2006.
 [16] D. Rohe and W. Metzner. Pseudogap at hot spots in the twodimensional hubbard model at weak coupling. Phys. Rev. B, 71:115116, 2005.
 [17] J. Reiss, D. Rohe, and W. Metzner. Renormalized meanfield analysis of antiferromagnetism and dwave superconductivity in the twodimensional Hubbard model. Phys. Rev. B, 75:075110, 2007.
 [18] C. Honerkamp and M. Salmhofer. Temperatureflow renormalization group and the competition between superconductivity and ferromagnetism. Phys. Rev. B, 64:184516, 2001.
 [19] C. Honerkamp, D. Rohe, S. Andergassen, and T. Enss. Interaction flow method for manyfermion systems. Phys. Rev. B, 70:235115, 2004.
 [20] S. Andergassen. Functional renormalization group analysis of Luttinger liquids with impurities. PhD thesis, Universität Stuttgart, 2006.
 [21] T. Enss. Renormalization, Conservation Laws and Transport in Correlated Electron Systems. PhD thesis, Universität Stuttgart, 2005.
 [22] T. Baier, E. Bick, and C. Wetterich. Temperature dependence of antiferromagnetic order in the Hubbard model. Phys. Rev. B, 70:125111, 2004.
 [23] F. Schütz and P. Kopietz. Functional renormalization group with vacuum expectation values and spontaneous symmetry breaking. J. Phys. A, 39:8205, 2006.
 [24] P. Strack, R. Gersch, and W. Metzner. in preparation.
 [25] A. A. Katanin. Fulfillment of Ward identities in the functional renormalization group approach. Phys. Rev. B, 70:115109, 2004.
 [26] M. Salmhofer, C. Honerkamp, W. Metzner, and O. Lauscher. Renormalization Group Flows into Phases with Broken Symmetry. Prog. Theor. Phys., 112:943, 2004.
 [27] R. Gersch, C. Honerkamp, D. Rohe, and W. Metzner. Fermionic renormalization group flow into phases with broken discrete symmetry: chargedensity wave meanfield model. Eur. Phys. J. B, 48:349, 2005.
 [28] R. Gersch, J. Reiss, and C. Honerkamp. Fermionic functional renormalizationgroup for firstorder phase transitions: a meanfield model. New J. Phys., 8:320, 2006.
 [29] M. C. Gutzwiller. Correlation of electrons in a narrow band. Phys. Rev., 137:A1726, Mar 1965.
 [30] J. Hubbard. Electron correlations in narrow energy bands. P. Roy. Soc. Lond., 276:238, 1963.
 [31] J. Kanamori. Electron Correlation and Ferromagnetism of Transition Metals. Prog. Theor. Phys., 30:275, 1963.
 [32] R. Micnas, J. Ranninger, and S. Robaszkiewicz. Superconductivity in narrowband systems with local nonretarded attractive interactions. Rev. Mod. Phys., 62:113, 1990.
 [33] A. MartínRodero and F. Flores. Solution for the Unegative Hubbard superconductor including secondorder correlation effects. Phys. Rev. B, 45:13008, 1992.
 [34] J. K. Chin, D. E. Miller, Y. Liu, C. Stan, W. Setiawan, C. Sanner, K. Xu, and W. Ketterle. Evidence for superfluidity of ultracold fermions in an optical lattice. Nature 443:961, 2006.
 [35] C. Honerkamp and M. Salmhofer. Eliashberg Equations Derived from the Functional Renormalization Group. Prog. Theor. Phys., 113:1145, 2005.
 [36] R. Gersch, C. Honerkamp, D. Rohe, and W. Metzner. Fermionic renormalization group flow at all scales: Breaking a discrete symmetry. In Ferdinando Mancini and Adolfo Avella, editors, LECTURES ON THE PHYSICS OF HIGHLY CORRELATED ELECTRON SYSTEMS X: Tenth Training Course in the Physics of Correlated Electron Systems and High Tc Superconductors, page 236, Melville, New York, 2006. AIP.
 [37] D. Zanchi and H. J. Schulz. Weakly correlated electrons on a square lattice: Renormalizationgroup theory. Phys. Rev. B, 61:13609, 2000.
 [38] J. Reiss. Renormalised MeanField Analysis of the 2D Hubbard Model. PhD thesis, Universität Stuttgart, 2006.
 [39] M. Ossadnik and C. Honerkamp. in preparation.