## Effect of pairing fluctuations on the spin resonance in Fe-based superconductors

Alberto Hinojosa, 1 Andrey V. Chubukov, 1 and Peter W¨ olfle 2

1

Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA 2 Institute for Condensed Matter Theory and Institute for Nanotechnology, Karlsruhe Institute of Technology, D-76128 Karlsruhe, Germany

The spin resonance observed in the inelastic neutron scattering data on Fe-based superconductors has played a prominent role in the quest for determining the symmetry of the order parameter in these compounds. Most theoretical studies of the resonance employ an RPA-type approach in the particle-hole channel and associate the resonance in the spin susceptibility χ S ( q , ω ) at momentum Q = ( π, π ) with the spin-exciton of an s + -superconductor, pulled below 2∆ by residual attraction associated with the sign change of the gap between Fermi points connected by Q . Here we explore the effect of fluctuations in the particle-particle channel on the spin resonance. Particle-particle and particle-hole channels are coupled in a superconductor and to what extent the spin resonance can be viewed as a particle-hole exciton needs to be addressed. In the case of purely repulsive interactions we find that the particle-particle channel at total momentum Q (the π channel) contributes little to the resonance. However, if the interband density-density interaction is attractive and the π -resonance is possible on its own, along with spin-exciton, we find a much stronger shift of the resonance frequency from the position of the would-be spin-exciton resonance. We also show that the expected doublepeak structure in this situation does not appear because of the strong coupling between particle-hole and particle-particle channels, and Im χ S ( Q , ω ) displays only a single peak.

## I. INTRODUCTION

The spin resonance, observed by inelastic neutron scattering (INS) experiments first in the cuprates 1 and then in heavy-fermion 2 and Fe-based superconductors (FeSCs) , has been the subject of intense theoretical and 3 experimental studies over the past decade using both metallic 4-6,10 and near-localized strong-coupling scenarios 8 (see Ref. 9 for a review). The theoretical interpretations of the resonance can be broadly split into two classes. The first class of theories assumes that the spin resonance is a magnon, overdamped in the normal state due to the strong decay into low-energy particle-hole pairs, but emerging prominently in the superconducting state due to reduction of scattering at low energies 10,11 . In this line of reasoning the resonance energy Ω res is the magnon energy and as such it is uncorrelated with the superconducting gap ∆. However, the decay of magnons into particle-hole pairs is only suppressed at energies below 2∆, hence the magnons become sharp in a superconductor only if their energy is below 2∆. The symmetry of the superconducting state does not play a crucial role here. It is only relevant that the superconducting gap is finite at the Fermi surface (FS) points connected by the gap momentum Q .

Theories from the second class assume that the resonance does not exist in the normal state and emerges in the superconducting state as a feedback effect from superconductivity, like e.g., the Anderson-Bogolyubov mode in the case of a charge neutral single condensate component, a Leggett mode in the case of several gap components, a wave-like excitation of a spin-triplet order parameter, or a pair vibration mode in the case of a gap parameter possessing internal structure 21 . These 'feedback' theories can be further split into three subclasses. In the first the resonance is viewed as a spin- exciton, i.e., the pole in the dynamical spin susceptibility χ ( q , Ω) dressed by multiple particle-hole bubbles 4-6 . Such χ ( q , Ω) can be obtained by using a computational scheme based on the random-phase approximation (RPA). It was argued that, if the superconducting gap changes sign between FS points connected by Q , the residual attraction pulls the resonance frequency to Ω res &lt; 2∆, where the decay into particle-hole pairs is reduced below T c and vanishes at T = 0. As a result, at T = 0, χ ′′ ( q , Ω) has a δ -functional peak at Ω res . In this respect, if the resonance is an exciton, its existence necessary implies that the superconducting gap changes sign either between patches of the FS connected by Q or between different Fermi pockets again connected by Q . The role of the resonance in allowing to determine the structure of the gap in a number of different superconductors has been highlighted in . 7 Theories of the second subclass explore the fact that in a superconductor the particle-hole and particle-particle channels are mixed and argue that the strongest resonance is in the particleparticle channel and the measurements of the spin susceptibility just reflect the 'leakage' of this resonance into the particle-hole channel. The corresponding resonance has been labeled as the π -exciton 12 , where the π boson is a particle-particle excitation with total momentum Q (a 'pair density-wave' in modern nomenclature 13-15 ). Finally, theories of the third class explore the possibility that the resonance emerges due to coupling between fluctuations in particle-hole and particle-particle channel. Within RPA, such resonance is due to non-diagonal terms in the generalized RPA which includes both particle-hole and particle-particle bubbles. It was called a plasmon 16 to stress the analogy with collective excitations in an electronic liquid.

The interplay between the 'damped spin-wave' scenarios, the spin-exciton π -resonance, and the plasmon sce-

narios for 2D highT c cuprates has been studied in detail in the past decade. The outcome is that near and above optimal doping the resonance is best described as a spinexciton, with relatively weak corrections from coupling to the particle-particle channel 5,17,18 , while in the underdoped regime, where superconductivity emerges from a pre-existing pseudogap phase, both spin-exciton and spin-wave scenarios have been argued to account for the neutron resonance . 9 The situation is more complicated in 3D heavy-fermion systems 11,19 because there the excitonic resonance has a finite width even at T = 0 if the locus of FS points separated by Q intersects a line of gap nodes.

In this work we discuss the interplay between spinexciton, π -resonance, and plasmon scenarios in FeSCs. Previous studies of the resonance in FeSCs 6 focused only on the response in the particle-hole (spin-exciton) channel and neglected the coupling between particle-hole and particle-particle channels. Our goal is to analyze the effect of such coupling.

As we said, the resonance peak has been observed below T c in several families of FeSCs 3 . The spin response above T c in FeSCs is rather featureless away from small doping, which implies that the magnetic excitations in the normal state are highly overdamped and don't behave as damped spin waves. The full analysis of the spin resonance in FeSCs is rather involved as these systems are multi-band materials with four or five FSs, on which the superconducting gap has different amplitudes and phases. Still, the basic conditions for spin-excitonic resonance are the same as in the cuprates and heavy fermion materials: namely, the resonance emerges at momentum Q if the superconducting gap changes sign between FS points connected by Q . This condition holds if the superconducting gap has s + -symmetry, as most researchers believe, and changes sign between at least some hole and electron pockets. An alternative scenario 20 , which we will not discuss in this paper is that the superconducting state has a conventional, sign-preserving s ++ symmetry, and the observed neutron peak is not a resonance but rather a hump at frequencies slightly above 2∆.

In FeSCs that contain both hole and electron pockets, the resonance has been observed at momenta around Q = ( π, π ), which is roughly the distance between hole and electron pockets in the actual (folded) Brillouin zone. To account for the resonance and, at the same time, avoid unnecessary complications, we consider a minimal threeband model (one hole pocket and two electron pockets), and neglect the angular dependence of the interactions along the FSs. Including this dependence and additional hole pockets will complicate the analysis but we don't expect it to lead to any qualitative changes to our results.

We find that for repulsive density-density and pairhopping interactions, the resonance peak is, to a good accuracy, a spin-exciton. The π -resonance does not develop on its own, and the coupling between the resonant spin-exciton channel and the non-resonant π -channel only slightly shifts the energy of the excitonic resonance.

We also considered the case (less justified microscopically) when the interaction in the π channel is attractive, such that both spin-exciton and π resonance develop on their own at frequencies below 2∆. One could expect in this situation that the full dynamical spin susceptibility has two peaks. We found, however, that this happens only if we make the coupling between particle-hole and particle-particle channels artificially small. When we restored the original coupling, we found, in general, only one peak below 2∆. The peak is a mixture of a spinexciton and π -resonance and at least in some range of system parameters its energy is smaller than that of a spin-exciton and a π -resonance. This implies that, when both channels are attractive, the coupling between the two plays a substantial role in determining the position of the true resonance which can, at least partly, be viewed as a plasmon. A somewhat similar result has been obtained earlier for the cuprates 5,16,18 . For some system parameters we did find two peaks in Im χ S ( Q , ω ), but for one of them Im χ S has wrong sign. We verified that this indicates that for such parameters the system is unstable either against condensation of π excitations (i.e., against superconductivity at momentum Q ), or against the development of SDW order in co-existence with superconductivity.

The paper is organized as follows: In the next section we consider the model. In Sec. III we obtain the dynamical spin susceptibility within the generalized RPA scheme, which includes the coupling between particlehole and particle-particle channels. In Sec. IV we analyze the profile of χ S ( Q , ω ) first for purely repulsive densitydensity and pair-hopping interactions, and then for the case when we allow the density-density interaction to become attractive. We summarize our conclusions in Sec. V.

## II. THE MODEL

The FeSCs are multiband metals with two or three hole FS pockets centered around the Γ point (0 0) and , two elliptical electron pockets centered at ( π, π ) in the folded BZ with two iron atoms per unit cell. The electron pockets are elliptical and related by symmetry, while the hole pockets are C 4 -symmetric, but generally differ in size. Since we are only interested in studying the role of the particle-particle channel in the spin response function, for which the non-equivalence between hole pockets is not essential, we consider the case of two hole pockets and assume that they are circular and identical, and also neglect the ellipticity of electron pockets. Under these assumptions our model reduces effectively to only one hole pocket ( c fermions) and one electron pocket ( f fermions). The fact that there are actually two hole and two electron pockets only adds up combinatoric factors.

The free part of the Hamiltonian is

$$H _ { 0 } = \sum _ { { \mathbf k }, \sigma } \left ( \xi _ { { \mathbf k } } ^ { c } c _ { { \mathbf k } \sigma } ^ { \dagger } c _ { { \mathbf k } \sigma } + \xi _ { { \mathbf k } + { \mathbf Q } } ^ { f } f _ { { \mathbf k } + { \mathbf Q } \sigma } ^ { \dagger } f _ { { \mathbf k } + { \mathbf Q } \sigma } \right ), \quad ( 1 ) \quad \text{parti} \\ \text{prod}$$

where

$$\xi _ { \mathbf k } ^ { c } = \mu _ { c } - \frac { k _ { x } ^ { 2 } + k _ { y } ^ { 2 } } { 2 m _ { c } }, \quad \quad \quad ( 2 )$$

$$\xi _ { \mathbf k + \mathbf Q } ^ { f } = - \mu _ { f } + \frac { k _ { x } ^ { 2 } + k _ { y } ^ { 2 } } { 2 m _ { f } } \, \quad \quad ( 3 ) \, \quad$$

We do not study here how superconductivity develops from interactions, as that work has been done elsewhere 22 . Instead, we simply assume that the system reaches a superconducting state with s + -symmetry before it becomes unstable towards magnetism and take the superconducting gaps as inputs. In this state, the free part of the Hamiltonian is

$$H _ { 0 } ^ { S C } = \sum _ { { \mathbf k }, \sigma } \left ( E _ { { \mathbf k } } ^ { c } c _ { { \mathbf k } \sigma } ^ { \dagger } c _ { { \mathbf k } \sigma } + E _ { { \mathbf k } + { \mathbf Q } } ^ { f } f _ { { \mathbf k } + { \mathbf Q } \sigma } ^ { \dagger } f _ { { \mathbf k } + { \mathbf Q } \sigma } \right ), \quad ( 4 ) \quad \text{To m}$$

where the dispersions are E c k = √ ( ξ c k ) 2 +(∆ ) c 2 and E f k + Q = √ ( ξ f k + Q ) 2 +(∆ ) , and ∆ f 2 c = -∆ f ≡ ∆.

Now we consider interactions that contribute to the spin susceptibility. They consist of a density-density interband interaction u 1 and a correlated interband hopping u 3 . Intraband repulsion only affects the chemical potentials but does not otherwise contribute to the spin susceptibility. Interband exchange in principle contributes in the π channel but in renormalization group analysis it flows to small values 26 . In general, u 1 and u 3 depend on the angle in momentum space via coherence factors associated with the transformation from the orbital to the band basis 23 . However, this complication is not essential for our purposes and we take both interactions to be momentum independent. The interaction Hamiltonian is

/negationslash

/negationslash

$$H _ { \text{int} } & = u _ { 1 } \sum _ { [ 1 2 3 4 ], \sigma \neq \sigma ^ { \prime } } c _ { \text{p} _ { 1 } \sigma } ^ { \dagger } f _ { \text{p} _ { 2 } \sigma ^ { \prime } } ^ { \dagger } f _ { \text{p} _ { 3 } \sigma ^ { \prime } } c _ { \text{p} _ { 4 } \sigma } \\ & + u _ { 3 } \sum _ { [ 1 2 3 4 ], \sigma \neq \sigma ^ { \prime } } \left ( c _ { \text{p} _ { 1 } \sigma } ^ { \dagger } c _ { \text{p} _ { 2 } \sigma ^ { \prime } } ^ { \dagger } f _ { \text{p} _ { 3 } \sigma ^ { \prime } } f _ { \text{p} _ { 4 } \sigma } + \text{h.c.} \right ), \quad ( 5 ) \quad \text{//where} \\ & \quad \text{//visual} \\ & \quad \text{//kw.k.}$$

where the sum over momenta obeys momentum conservation as usual ( p 1 + p 2 = p 3 + p 4 ).

Because the interactions in the band basis are linear combinations of Hubbard and Hund interactions in the orbital basis, weighted with orbital coherence factors, the sign of u 1 and u 3 depends on the interplay between intraorbital and inter-orbital interactions 20,24 . The interaction u 3 contributes to the superconducting channel, and for an s + -gap structure must be repulsive. The sign of u 1 is a priori unknown. In most microscopic studies it comes out positive (repulsive), but in principle it can also be negative (attractive). We do not assume a particular sign of u 1 and consider first a case where u 1 is positive and then when it is negative. For the first case we show that a resonance can only originate from the particle-hole channel. For negative u 1 the π channel can produce collective modes as well, and we show that in general the resonant mode is a mix between spin exciton and a π -resonance.

## III. SUSCEPTIBILITIES AND RPA

We focus on susceptibilities at antiferromagnetic momentum Q which separates the centers of hole and electron pockets. Following similar work done on the cuprates , we define spin and 5 π operators as

$$\underset { \text{ take }, \text{ the } } ^ { \prime \text{ be-} } S ^ { z } ( \mathbf Q ) = \frac { 1 } { \sqrt { N } } \sum _ { \mathbf k } \left [ c _ { \mathbf k \alpha } ^ { \dagger } \sigma _ { \alpha \beta } ^ { z } f _ { \mathbf k + \mathbf Q \beta } + f _ { \mathbf k + \mathbf Q \alpha } ^ { \dagger } \sigma _ { \alpha \beta } ^ { z } c _ { \mathbf k \beta } \right ], \ ( 6 )$$

$$\intertext { a } _ { \substack { \pi ( \mathbf Q ) } } = \frac { 1 } { \sqrt { N } } \sum _ { \mathbf k } ^ { \mu \mu } \left [ c _ { \mathbf k \alpha } \sigma _ { \alpha \beta } ^ { x } f _ { \mathbf Q - \mathbf k \beta } \right ].$$

To make a closer connection to , the operator 5 π can be equivalently defined as π = 1 √ N ∑ [ k g k a k α σ x αβ a Q k -β ] , with | g k | = 1 2 and the sign of / g k is chosen so that it is positive near the hole FS and negative near the electron FS ( g k = -g Q k -).

For notational convenience we split S z into two operators such that S z = S c + S f , where

$$S _ { c } ( \mathbf Q ) = \frac { 1 } { \sqrt { N } } \sum _ { \mathbf k } c _ { \mathbf k \alpha } ^ { \dagger } \sigma _ { \alpha \beta } ^ { z } f _ { \mathbf k + \mathbf Q \beta }, \quad \quad ( 8 )$$

$$S _ { f } ( { \mathbf Q } ) = \frac { 1 } { \sqrt { N } } \sum _ { \mathbf k } ^ { \, \, \cdot \, } f _ { \mathbf k + { \mathbf Q } \alpha } ^ { \dagger } \sigma _ { \alpha \beta } ^ { z } c _ { \mathbf k \beta }.$$

We now define the susceptibilities χ ab (Ω m ) in terms of Matsubara frequencies as

$$\varsigma _ { a b } ( \Omega _ { m } ) = \int _ { 0 } ^ { 1 / T } \mathrm d \tau ^ { \prime } e ^ { i \Omega _ { m } \tau ^ { \prime } } \left \langle T _ { \tau } A _ { a } ( \tau ^ { \prime } ) A _ { b } ^ { \dagger } ( 0 ) \right \rangle, \quad ( 1 0 )$$

where A a = ( S , S c f , π, π † ) a . The actual spin susceptibility is given by χ S = χ 11 + χ 12 + χ 21 + χ 22 .

The bare susceptibilities χ 0 ab can be calculated in the usual way in terms of Green's functions and are given by bubbles made out of c and f -fermions, with different Pauli matrices in the vertices. At T = 0 and after performing analytic continuation to real frequency space, the (retarded) susceptibilities have the following form:

$$\begin{array} { c c c c } e & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & | & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ( & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

where χ 0 ab , A ab , and B ab are symmetric matrices. The expressions for A ab and B ab are presented in Table I in

Table I: Coefficients of bare susceptibilities. Note: k ′ = k + Q .

|   a |   b | A ab ( k )                    | B ab ( k )                    |
|-----|-----|-------------------------------|-------------------------------|
|   1 |   1 | ( v c k u f k ′ ) 2           | ( u c k v f k ′ ) 2           |
|   1 |   2 | - u c k v c k u f k ′ v f k ′ | - u c k v c k u f k ′ v f k ′ |
|   1 |   3 | u c k v c k ( u f k ′ ) 2     | - u c k v c k ( v f k ′ ) 2   |
|   1 |   4 | ( v c k ) 2 u f k ′ v f k ′   | - ( u c k ) 2 u f k ′ v f k ′ |
|   2 |   2 | ( u c k v f k ′ ) 2           | ( v c k u f k ′ ) 2           |
|   2 |   3 | - ( u c k ) 2 u f k ′ v f k ′ | ( v c k ) 2 u f k ′ v f k ′   |
|   2 |   4 | - u c k v c k ( v f k ′ ) 2   | u c k v c k ( u f k ′ ) 2     |
|   3 |   3 | ( u c k u f k ′ ) 2           | ( v c k v f k ′ ) 2           |
|   3 |   4 | u c k v c k u f k ′ v f k ′   | u c k v c k u f k ′ v f k ′   |
|   4 |   4 | ( v c k v f k ′ ) 2           | ( u c k u f k ′ ) 2           |

terms of coherence factors which are given by:

$$u _ { \mathbf k } ^ { c } = \sqrt { \frac { 1 } { 2 } \left ( 1 + \frac { \xi _ { \mathbf k } ^ { c } } { E _ { \mathbf k } ^ { c } } \right ) },$$

$$v _ { \mathbf k } ^ { c } = \sqrt { \frac { 1 } { 2 } \left ( 1 - \frac { \xi _ { \mathbf k } ^ { c } } { E _ { \mathbf k } ^ { c } } \right ) } \, \text{sgn} \, \Delta ^ { c }, \quad \ \ ( 1 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

and similar expressions for u f k ′ and v f k ′ .

To obtain the full susceptibilities χ ab we used the generalized RPA approach. Within this approach

$$\chi _ { a b } = ( 1 - \chi ^ { 0 } V ) _ { a c } ^ { - 1 } \chi _ { c b } ^ { 0 }, \text{ \quad \ \ } ( 1 4 )$$

where the sum over repeated indices is implied and V is given by

$$V = \frac { 1 } { 2 } \begin{pmatrix} u _ { 1 } & u _ { 3 } & 0 & 0 \\ u _ { 3 } & u _ { 1 } & 0 & 0 \\ 0 & 0 & - u _ { 1 } & 0 \\ 0 & 0 & 0 & - u _ { 1 } \end{pmatrix}. \quad \quad ( 1 5 ) \quad \text{fig} \\ \text{cpt} \quad \text{$\epsilon$}.$$

The solution for the full spin susceptibility can be written in a simpler form once we note that the matrix χ 0 ab has additional symmetry. Indeed, the functions A ab ( k ) and B ab ( k ) can be separated into parts that are even or odd with respect to ξ c k and ξ f k . If the momentum sums are evaluated only near the FSs, where the integration region can be chosen to be symmetric with respect to positive and negative values of ξ c k and ξ f k , then the odd parts cancel out and χ 0 ab acquires the following form:

$$\chi ^ { 0 } = \begin{pmatrix} a & b & c & - d \\ b & a & d & - c \\ c & d & e & f \\ - d & - c & f & e \end{pmatrix} \quad \quad ( 1 6 )$$

If we rotate the basis of operators from ( S , S c f , π, π † ) to ( S c + S , π f -π , S † c -S , π f + π † ) we find that both V and χ 0 become block diagonal, so the first 2 × 2 system decouples from the second one. In this case, the solution for the subset ( S c + S , π f -π † ) takes the simple form

$$\chi _ { S } = & \frac { \chi _ { S } ^ { 0 } + \delta \chi _ { S } ^ { 0 } } { 1 - u _ { S } ( \chi _ { S } ^ { 0 } + \delta \chi _ { S } ^ { 0 } ) }, \quad \quad ( 1 7 )$$

$$\chi _ { \pi } = & \frac { \chi _ { \pi } ^ { 0 } + \stackrel { \lambda } { \delta } \chi _ { \pi } ^ { 0 } } { 1 - u _ { \pi } ( \chi _ { \pi } ^ { 0 } + \delta \chi _ { \pi } ^ { 0 } ) },$$

$$\delta \chi _ { S } ^ { 0 } = & \frac { u _ { \pi } } { 1 - u _ { \pi } \chi _ { \pi } ^ { 0 } } \left ( \chi _ { S \pi } ^ { 0 } \right ) ^ { 2 }, \quad \quad \left ( 1 9 \right )$$

$$\delta \chi _ { \pi } ^ { 0 } = & \frac { \mu ^ { - n \wedge \pi } } { u _ { S } } \left ( \chi _ { S \pi } ^ { 0 } \right ) ^ { 2 }. \quad \quad \left ( 2 0 \right )$$

and u S = ( u 1 + u 3 ) / 4 and u π = -u / 1 4. The bare susceptibilities in this basis are given by

$$\chi ^ { 0 } _ { S } = \chi ^ { 0 } _ { 1 1 } + \chi ^ { 0 } _ { 1 2 } + \chi ^ { 0 } _ { 2 1 } + \chi ^ { 0 } _ { 2 2 } = 2 ( \chi ^ { 0 } _ { 1 1 } + \chi ^ { 0 } _ { 1 2 } ), \quad ( 2 1 )$$

$$\chi ^ { 0 } _ { \pi } = \chi ^ { 0 } _ { 3 3 } - \chi ^ { 0 } _ { 3 4 } - \chi ^ { 0 } _ { 4 3 } + \chi ^ { 0 } _ { 4 4 } = 2 ( \chi ^ { 0 } _ { 3 3 } - \chi ^ { 0 } _ { 3 4 } ), \quad ( 2 2 )$$

$$\chi ^ { 0 } _ { S \pi } = \chi ^ { 0 } _ { 1 3 } - \chi ^ { 0 } _ { 1 4 } + \chi ^ { 0 } _ { 2 3 } - \chi ^ { 0 } _ { 2 4 } = 2 ( \chi ^ { 0 } _ { 1 3 } - \chi ^ { 0 } _ { 1 4 } ). \quad ( 2 3 )$$

The first two bare susceptibilities contain contributions of products of two normal Green's functions and of two anomalous Green's functions, while χ 0 Sπ is composed of one normal and one anomalous Green's function (see Fig. 1 for some contributions of χ 0 Sπ to the spin response function).

Figure 1: Some mixed-channel contributions to the spin susceptibility and Raman scattering. The solid lines represent f fermions and the dashed lines c fermions, corresponding to quasiparticles from the electron and hole pockets, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0021v1-figure-2dbb7ea9-f9ea-44b7-9e84-a7020ae53498.jpg)

We obtained expressions for the real parts of the bare susceptibilities by replacing the momentum sums by integrals and evaluating them within an energy range from -Λ to Λ about the FSs. They are valid in the range 0 &lt; ω &lt; 2∆ in the limit when the broadening γ → 0.

$$R e \chi _ { S } ^ { 0 } ( \omega ) = & L + \frac { \omega ^ { 2 } } { 4 \Delta ^ { 2 } } R e \chi _ { \pi } ^ { 0 } ( \omega ), \quad \quad \quad ( 2 4 )$$

$$\text{Re} \chi _ { \pi } ^ { 0 } ( \omega ) = & \frac { m } { \pi } \frac { \stackrel { \text{$\mapsto \ell$} } { 4 \Delta ^ { 2 } } } { \omega \sqrt { 4 \Delta ^ { 2 } - \omega ^ { 2 } } } \arctan \left ( \frac { \omega } { \sqrt { 4 \Delta ^ { 2 } - \omega ^ { 2 } } } \right ), \\ \dots$$

$$\tilde { R e } \chi _ { S \pi } ^ { 0 } ( \omega ) = & \frac { \omega } { 2 \Delta } \tilde { R e } \chi _ { \pi } ^ { 0 } ( \omega ), & \cdot & \cdot \\ \pi ^ { \dagger } )$$

where L = m π log(2Λ / ∆) and we have neglected terms of order (∆ / Λ) . 2

where

Figure 2: Real part of the bare spin susceptibilities ( T = 0, units of ∆ -1 ) in the case of perfectly nested circular pockets ( m ≡ m c = m f = 1 100∆ ). In all numerical calculations µ = 10∆ and we include a finite broadening γ = ∆ 200 when / evaluating momentum integrals.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0021v1-figure-0d425560-917d-440e-a481-661f92058b3b.jpg)

In Fig. 2 we present the results of numerical calculations of the bare susceptibilities in the case of perfectly nested FSs ( µ c = µ f ≡ µ ). The susceptibility χ 0 S would diverge at ω = 0 in the absence of superconductivity, but becomes finite at a finite ∆. Conversely, χ 0 π at ω = 0 would be zero in the absence of superconductivity, but becomes non-zero because of ∆. Note that all three susceptibilities monotonically increase with frequency in the domain 0 ≤ ω ≤ 2∆ and diverge at ω = 2∆. The imaginary parts of the three bare susceptibilities (not shown in the plot) are infinitesimally small and undergo a discontinuous jump at ω = 2∆. If we make the electron pocket elliptical, the divergence in the real part is replaced by a local maximum.

We see that χ 0 S ( ω ) and χ 0 π ( ω ) have finite value at ω = 0. We recall that, in the absence of superconductivity, χ S (0) would diverge and χ π (0) would vanish at perfect nesting. At ω = 2∆ both bare susceptibilities diverge. The cross-susceptibility χ 0 Sπ ( ω ) vanishes at ω = 0 simply because it is composed from one normal and one anomalous Green's function but rapidly increases with ω and becomes comparable to χ S ( ω ) and χ π ( ω ) at ω ≤ 2∆. The cross-susceptibility between particle-hole and particle-particle channels has been recently analyzed for an s + -superconductor in the context of Raman scattering 25 . There, it was computed in the charge channel and was found to be very small due to near-cancellation between contributions from Fermi surfaces with plus and minus values of the superconducting gap. In our case we found that the contributions from hole and electron FSs add up rather than cancel. The difference is that in Raman scattering the side vertices in the susceptibility bubble have the spin structure given by δ α,β , while in our case the spin structure is, say, σ z α,β . For an s + - gap, Raman bubbles from hole and electron pocket have the same vertex structure but differ in the sign of the anomalous Green's function, hence the two contributions to cross-susceptibility have opposite signs, resulting in a cancellation that is complete in the case of perfect nesting and near-complete in the case of one circular and one elliptical pocket. This cancellation does not occur in our case because the σ z structure of the side vertices additionally flips the sign of one of the two diagrams, and the contributions to χ 0 Sπ ( ω ) from hole and electron FSs add constructively.

We next note that the terms δχ 0 S and δχ 0 π are precisely what is neglected when the particle-particle channel is not included in the calculation of the spin susceptibility. Setting these terms to zero reduces the expressions for the full χ S and χ π to the usual RPA results

$$\chi _ { S } = \frac { \chi _ { S } ^ { 0 } } { 1 - u _ { S } \chi _ { S } ^ { 0 } }, \quad \chi _ { \pi } = \frac { \chi _ { \pi } ^ { 0 } } { 1 - u _ { \pi } \chi _ { \pi } ^ { 0 } }. \quad ( 2 7 )$$

The effect of coupling the two channels can be seen more clearly by substituting (19) into (17), which yields

$$\text{lcula} _ { \text{effectly} } \quad \chi _ { S } = & \frac { \chi _ { S } ^ { 0 } ( 1 - u _ { \pi } \chi _ { \pi } ^ { 0 } ) + u _ { \pi } ( \chi _ { S \pi } ^ { 0 } ) ^ { 2 } } { ( 1 - u _ { S } \chi _ { S } ^ { 0 } ) ( 1 - u _ { \pi } \chi _ { \pi } ^ { 0 } ) - u _ { S } u _ { \pi } ( \chi _ { S \pi } ^ { 0 } ) ^ { 2 } } \\ \text{-} \text{-} \text{-} 0 \\ \text{-} \text{-} 0$$

The positions of resonance peaks are given by the zeroes of the denominator in this equation and we can see that the particle-hole and particle-particle channels are coupled through the mixed-channel susceptibility χ 0 Sπ ( ω ).

## IV. THE RESULTS

## A. Purely repulsive interaction, u 1 &gt; 0 , u 3 &gt; 0 .

For repulsive interactions u S &gt; u | π | &gt; 0 and u π &lt; 0. In the absence of χ 0 Sπ the resonance in χ S is present for any u S because the bare susceptibility χ 0 S is positive and diverges at ω = 2∆, hence the equation 1 -u S Re χ 0 S ( ω ) = 0 has a solution for 0 &lt; u S &lt; (Re χ 0 S (0)) -1 . In contrast, the fact that Re χ 0 π &gt; 0 means that no resonance originates from this channel. When χ 0 Sπ is included, we found in our numerical calculations that the effect of the particle-particle channel is that the peak in the imaginary part of the full susceptibility is shifted to a higher frequency (since Re { u S u π ( χ 0 Sπ ) 2 } &lt; 0). We show representative behavior of real and imaginary parts of the full spin susceptibility in Fig. 3

This result is also obtained when we consider an elliptical electron pocket, except that there is a minimum value for u S below which no resonance is observed. This is due to the fact that the bare susceptibilities have a local maximum instead of a divergence at ω = 2∆.

Figure 3: Imaginary part (main plot) and real part (insert) of the full spin susceptibility ( T = 0, units of ∆ -1 ) in the case of repulsive interactions. In this plot u S = 26∆ and u π = -13∆. The solid, black line indicates the full calculation. For comparison, the dashed, blue line indicates a calculation that neglects particle-particle contributions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0021v1-figure-b3ca8def-004b-4ce4-83e8-8e3e4108d95e.jpg)

## B. Partially attractive interaction, u 1 &lt; ,u 0 3 &gt; 0 .

We now consider an alternative case where the densitydensity interaction u 1 &lt; 0, hence u π = -u / 1 4 &gt; 0. We still assume u 3 &gt; u | 1 | such that u S &gt; 0. For positive u π , the π channel can acquire a resonance on its own since the equation 1 -u π Re χ 0 π ( ω ) = 0 necessarily has a solution at a frequency between 0 and 2∆ if u π &lt; (Re χ 0 π (0)) -1 . We assume that this inequality holds together with u S &lt; (Re χ 0 S (0)) -1 . If any of these two conditions are not satisfied, the system becomes unstable either towards π superconductivity with total momentum of a Cooper pair Q or towards magnetic order. In both cases, the analysis of the spin susceptibility has to be modified to include the new condensates. If the spinexciton and π channels were not coupled (i.e., if χ 0 Sπ ( ω ) was absent), we would find resonances in the spin and π channels at frequencies ω S and ω π , respectively, set by 1 -u S Re χ 0 S ( ω S ) = 0 and 1 -u π Re χ 0 π ( ω π ) = 0. This suggests the possibility that there may be two resonance peaks in the full spin susceptibility χ S ( Q,ω ) once we restore the coupling χ 0 Sπ ( ω ).

However, we found that for all values of u S and u π for which the pure s + -state is stable, there is only a single peak in the spin susceptibility at a frequency lower than both ω S and ω π . We show representative behavior in Fig. 4. The existence of a single peak is due to the fact that χ 0 Sπ is small at small frequencies, hence it does not prevent the increase of the real part of the spin susceptibility with increasing ω (see insert in Fig. 4) and only shifts the position of the lower pole ( ω S or ω π ) to a smaller value ω res . At the same time, at higher frequencies, χ 0 Sπ is no longer small relative to the other bare susceptibilities χ 0 S and χ 0 π . As a result, the denominator in χ S ( Q , ω ) in (28) passes through zero at ω = ω res and then remains negative all the way up to ω = 2∆ and does not cross zero for the second time.

To better understand this, we artificially add a factor /epsilon1 to χ 0 Sπ and consider how the solutions evolve as we progressively increase /epsilon1 between 0 and 1. At small /epsilon1 , the two solutions obviously survive and just further split from each other - the peak that was at a higher frequency shifts to a higher frequency and the other peak shifts to to a lower frequency. As /epsilon1 increases, the peak at a higher frequency rapidly moves towards 2∆. If we keep Im χ 0 ab strictly zero, this peak survives up to /epsilon1 = 1 with exponentially vanishing amplitude. If, however, we keep a small but finite fermionic damping in the computations of χ 0 ab , we find that the functions χ 0 ab ( ω ) increase but do not diverge at 2∆. In this situation, the higher frequency peak in χ S ( Q , ω ) vanishes already at some /epsilon1 &lt; 1.

We also considered the evolution of the two-peak solution with /epsilon1 in a different way: we postulated that the two peaks should be at ω res, 1 and ω res, 2 , both below 2∆ and solved the set of equations for u S and u π which would correspond to such a solution. At small /epsilon1 we indeed found some real u S and u π which satisfy 'boundary conditions' u S χ 0 S &lt; 1 and u χ π 0 π &lt; 1. However this holds only up to some /epsilon1 cr . At higher /epsilon1 the solutions for u s and u π become complex, which implies that the two-peak solution is no longer possible. At even higher /epsilon1 real solutions for u s and u π reappear, but they do not satisfy the boundary conditions. We searched for a range of ω res, 1 and ω res, 2 and for all values that we tested we found /epsilon1 cr &lt; 1, i.e., again there is only a single peak for the actual case of /epsilon1 = 1.

Another way to see that there is only one peak in the full χ S is to substitute the expressions for the real parts of the susceptibilities, Eqs. (24)-(26), into the denominator of Eq. (28) and express the real part of the term D = (1 -u S χ 0 S )(1 -u χ π 0 π ) -u S u π ( χ 0 Sπ ) 2 via χ 0 π = χ 0 π ( ω ). We obtain

$$\stackrel { \sigma \psi \psi } { S _ { \pi } } ( \omega ) \\ \text{and } \pi \\ \text{set by } \\ \text{This } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D } \quad \mathfrak { D }$$

Because (1 -u S L ) and (1 -u χ π 0 π (0)) are required to be positive for the stability of the paramagnetic state, at zero frequency, D is surely positive. At finite ω , the first term in (29) is positive, while the second one is negative and its magnitude monotonically increases with increasing ω . As a result, the denominator crosses zero only once, at some ω &lt; 2∆.

The single resonance peak is a mixture of a spinexciton and π -resonance and for the representative case shown in Fig. 4 its energy is smaller that that of spinexciton and a π -resonance. This implies that, when both channels are attractive, the coupling between the two plays substantial role in determining the position of the true resonance. From this perspective, the resonance at

u 1 &lt; 0 can, at least partly, be viewed as a plasmon. A somewhat similar result has been earlier obtained in the analysis of the resonance in the cuprates in the parameter range where π -resonance is allowed 5,16,18 .

Figure 4: Imaginary part (main plot) and real part (insert) of the full spin susceptibility ( T = 0, units of ∆ -1 ) in the case of u S = 30 5∆ and . u π = 130∆. In the absence of coupling we would observe resonance peaks in the spin and π channels at frequencies ω S ≈ 0 57∆ and . ω π ≈ 1 0∆, respectively, . indicated on the plot.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0021v1-figure-bc51402d-01bd-4cb7-8710-1c3dd72f53c8.jpg)

## V. CONCLUSIONS

We have studied the spin resonance at antiferromagnetic momentum ( π, π ) in an s + -superconducting state

1 J. Rossat-Mignod, L.P. Regnault, C. Vettier, P. Bourges, P. Burlet, J. Bossy, J.Y. Henry, and G. Lapertot, Physica C 185-189 , 86 (1991); H. A. Mook, M. Yethitaj, G. Aeppli, T. E. Mason, and T. Armstrong, Phys. Rev. Lett. 70 , 3490 (1993); H. F. Fong, B. Keimer, P. W. Anderson, D. Reznik, F. Dogan, and I. A. Aksay, Phys. Rev. Lett. 75 , 316 (1995); H.F. Fong et al, Nature 398 , 588 (1999); Ph. Bourges et al, Science 288 1234 (2000); H. He et al, Science 295 , 1045 (2002); G. Yu, Y. Li, E. M. Motoyama, and M. Greven, Nat. Phys. 5 , 873 (2009).

2 N.K. Sato, N. Aso, K. Miyake, R. Shiina, P. Thalmeier, G. Varelogiannis, C. Geibel, F. Steglich, P. Fulde, and T. Komatsubara, Nature 410 , 340 (2001); C. Stock, C. Broholm, J. Hudis, H.J. Kang, and C. Petrovic, Phys. Rev. Lett. 100 , 087001 (2008).

3 A. D. Christianson, E. A. Goremychkin, R. Osborn, S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. Todorov, H. Claus, D. Y. Chung, M. G. Kanatzidis, R. I. Bewley, and T. Guidi, Nature (London) 456 , 930 (2008); M. D. Lumsden, A. D. Christianson, D. Parshall, M. B.

of FeSCs by including contribution from the particleparticle channel, which in the superconducting state gets mixed with the particle-hole channel. We have shown that for purely repulsive interactions the inclusion of this channel does not qualitatively change the spin resonance, which remains predominantly spin-exciton and only slightly shifts to higher frequencies. For attractive density-density interaction, when both spin-exciton resonance in the particle-hole channel and π -resonance in the particle-particle channel are allowed, we found that strong coupling between the two channels destroys the two-peak structure and only one peak survives, whose frequency is smaller than would be that of a spin-exciton and π -resonance in the absence of the coupling. We argued that strong coupling between the particle-hole and particle-particle channels is peculiar to the spin susceptibility, while for the charge susceptibility, which, e.g., is relevant for Raman scattering, the coupling is much smaller.

We acknowledge useful discussions with R. Fernandes, S. Maiti, and Y. Wang. The work by A. H. and A.V.C. is supported by the DOE grant DE-FG02-ER46900. P. W. thanks the Department of Physics at the University of Wisconsin-Madison for hospitality during a stay as a visiting professor. P. W. also acknowledges partial support through the DFG research unit 'Quantum phase transitions'.

Stone, S. E. Nagler, G. J. MacDougall, H. A. Mook, K. Lokshin, T. McGuire, A. S. Sefat, R. Jin, B. C. Sales, and D. Mandrus, Phys. Rev. Lett. 102 , 107005 (2009); S. Chi, A. Schneidewind, J. Zhao, L. W. Harriger, L. Li, Y. Luo, G. Cao, Z. Xu, M. Loewenhaupt, J. Hu, and P. Dai, Phys. Rev. Lett. 102 , 107006 (2009); S. Li, Y. Chen, S. Chang, J.W. Lynn, L. Li, Y. Luo, G. Cao, Z. Xu, and P. Dai, Phys. Rev. B 79 , 174527 (2009); D.S. Inosov, J.T. Park, P. Bourges, D.L. Sun, Y. Sidis, A. Schneidewind, K. Hradil, D. Haug, C. T. Lin, B. Keimer, V. Hinkov, Nature Phys. 6 , 178 (2010); J.-P. Castellan, S. Rosenkranz, E. A. Goremychkin, D. Y. Chung, I. S. Todorov, M. G. Kanatzidis, I. Eremin, J. Knolle, A. V. Chubukov, S. Maiti, M. R. Norman, F. Weber, H. Claus, T. Guidi, R. I. Bewley, and R. Osborn, arXiv:1106.0771; Y. Qiu, W. Bao, Y. Zhao, C. Broholm, V. Stanev, Z. Tesanovic, Y. C. Gasparovic, S. Chang, J. Hu, B. Qian, M. Fang, and Z. Mao, Phys. Rev. Lett. 103, 067008 (2009); D.N. Argyriou, A. Hiess, A. Akbari, I. Eremin, M.M. Korshunov, J. Hu, B. Qian, Z. Mao, Y. Qiu, C. Broholm, and W. Bao, Phys. Rev. B 81,

- 4 N. Bulut, and D. J. Scalapino, Phys. Rev. B 47 , 3419 (1993); I. Mazin and V. Yakovenko, Phys. Rev. Lett. 75 4134 (1995); D. Z. Liu, Y. Zha, and K. Levin, Phys. Rev. Lett. 75 4130 (1995); N. Bulut and D. Scalapino, Phys. Rev. B 53 , 5149 (1996); A. J. Millis and H. Monien, Phys. Rev. B, 54 , 16172 (1996); D. Manske, I. Eremin, and K. H. Bennemann, Phys. Rev. B 63 , 054517 (2001); M.R. Norman, Phys. Rev. B 61 , 14751 (2000); ibid 63 , 092509 (2001); A. Chubukov, B. Janko and O. Tchernyshov, Phys. Rev. B 63 , 180507(R) (2001); A. Abanov, A. Chubukov, and J. Schmalian, Journal of Electron Spectroscopy and Related Phenomena 117, 129-151 (2001); F. Onufrieva and P. Pfeuty, Phys. Rev. B 65, 054515 (2002); I. Eremin, D.K. Morr, A.V. Chubukov, K.H. Bennemann, and M.R. Norman, Phys. Rev. Lett. 94, 147001 (2005); Ar Abanov, A. Chubukov, and M. Norman, Phys. Rev. B 78, 220507 (2008).
- 5 O. Tchernyshyov, M. R. Norman, and A. V. Chubukov, Phys. Rev. B 63 , 144507 (2001).
- 6 M.M. Korshunov, and I. Eremin, Phys. Rev. B 78 , 140509(R) (2008); T.A. Maier, and D.J. Scalapino, Phys. Rev. B 78 , 020514(R) (2008); T. A. Maier, S. Graser, D. J. Scalapino, and P. Hirschfeld, Phys. Rev. B 79 , 134520 (2009); J. Zhang, R. Sknepnek, and J. Schmalian, Phys. Rev. B 82 , 134527 (2010); S. Maiti, J. Knolle, I. Eremin, and A.V. Chubukov, Phys. Rev. B 84 , 144524 (2011).
- 7 D. J. Scalapino, Rev. Mod. Phys. 84 , 1383 (2012).
- 8 J. Brinckmann, and P. A. Lee, Phys. Rev. Lett. 82 , 2915 (1999); K. Seo, C. Fang, B.A. Bernevig, and J. Hu, Phys. Rev. B 79 , 235207 (2009).
- 9 for a review see M. Eschrig, Adv. Phys. 55 , 47 (2006).
- 10 D. Morr and D. Pines, Phys. Rev. Lett. 81 , 1086 (1998).
- 11 A.V. Chubukov and L.P. Gorkov, Phys. Rev. Lett. 101, 147004 (2008).
- 12 E. Demler, and S.-C. Zhang, Phys. Rev. Lett. 75, 4126 (1995); E. Demler, H. Kohno, and S.-C. Zhang, Phys. Rev. B 58 , 5719 (1998).
- 13 E. Fradkin and S. A. Kivelson, Nature Physics 8 , 865-866 (2012).
- 14 P. A. Lee, arXiv:1401.0519.
- 15 D. Agterberg and M. Kashuap, arXiv:1406.4959.
- 16 W.C. Lee et al, Phys. Rev. B 77, 214518 (2008) W.C. Lee and A.H. MacDonald, Phys. Rev. B 78, 174506 (2008).
- 17 I. Fomin, P. Schmitteckert, and P. W¨lfle, Phys. Rev. Lett. o 69 , 214 (1992).
- 18 Z. Hao and A. Chubukov, Phys. Rev. B 79 , 224513 (2009).
- 19 I. Eremin, G. Zwicknagl, P. Thalmeier, and P. Fulde, Phys. Rev. Lett. 101, 187001 (2008).
- 20 S. Onari and H. Kontani, Phys. Rev. B 84 , 144518 (2011).
- 21 D. Vollhardt and P. W¨lfle, 'The Superfluid Phases of Heo lium 3', Taylor and Francis, London, (1990).
- 22 see e.g., I.I. Mazin and J. Schmalian, Physica C, 469 , 614 (2009); S. Graser, T. A. Maier, P. J. Hirshfeld, D. J. Scalapino, New J. Phys. 11 , 025016 (2009); K. Kuroki, H. Usui, S. Onari, R. Arita, and H. Aoki, Phys. Rev. B 79, 224511 (2009); J-P Paglione and R.L. Greene, Nature Phys. 6 , 645 (2010); D.N. Basov and A.V. Chubukov, Nature Physics 7 , 241 (2011); P.J. Hirschfeld, M.M. Korshunov, and I.I. Mazin, Reports on Progress in Physics 74, 124508 (2011); H.H. Wen and S. Li, Annu. Rev. Condens. Matter Phys., 2 , 121 (2011); A.V. Chubukov, Annual Review of Condensed Matter Physics 3, 57 (2012).
- 23 A.V. Chubukov, M.G. Vavilov, and A.B. Vorontsov, Phys. Rev. B 80, 140515(R) (2009); A. F. Kemper, T. A. Maier, S. Graser, H.-P. Cheng, P. J. Hirschfeld, and D. J. Scalapino New J. Phys. 12 073030 (2010).
- 24 S. Raghu, Xiao-Liang Qi, Chao-Xing Liu, D. Scalapino, and Shou-Cheng Zhang, Phys.Rev.B 77 220503(R) (2008); A. V. Chubukov, D. Efremov, and I. Eremin, Phys. Rev. B 78 , 134512 (2008).
- 25 M. Khodas, A.V. Chubukov, and G. Blumberg, Phys. Rev. B 89 , 245134 (2014).
- 26 A. V. Chubukov, D. Efremov, and I. Eremin, Phys. Rev. B 78 , 134512 (2008).
Alberto Hinojosa Andrey V. Chubukov Peter Wolfle
Effect of pairing fluctuations on the spin resonance in Fe-based superconductors
The spin resonance observed in the inelastic neutron scattering data on
Fe-based superconductors has played a prominent role in the quest for
determining the symmetry of the order parameter in these compounds. Most
theoretical studies of the resonance employ an RPA-type approach in the
particle-hole channel and associate the resonance in the spin susceptibility
$\chi_S (\mathbf{q}, \omega)$ at momentum ${\bf Q} = (\pi,\pi)$ with the
spin-exciton of an $s^{+-}$ superconductor, pulled below $2\Delta$ by residual
attraction associated with the sign change of the gap between Fermi points
connected by ${\bf Q}$. Here we explore the effect of fluctuations in the
particle-particle channel on the spin resonance. Particle-particle and
particle-hole channels are coupled in a superconductor and to what extent the
spin resonance can be viewed as a particle-hole exciton needs to be addressed.
In the case of purely repulsive interactions we find that the particle-particle
channel at total momentum ${\bf Q}$ (the $\pi $channel) contributes little to
the resonance. However, if the interband density-density interaction is
attractive and the $\pi-$resonance is possible on its own, along with the
spin-exciton, we find a much stronger shift of the resonance frequency from the
position of the would-be spin-exciton resonance. We also show that the expected
double-peak structure in this situation does not appear because of the strong
coupling between particle-hole and particle-particle channels, and $\mathrm{Im}
\chi_S ({\bf Q}, \omega)$ displays only a single peak.