# Screening of the Raman response in multiband superconductors – application to iron-pnictides

###### Abstract

We performed model calculations of Raman responses for multiband 2D superconductors. The multiband effects of screening in the symmetry channel were investigated analytically and numerically for a band structure model mimicing ARPES data on iron-pnictide materials. An acceptable agreement between our model calculations and recent experimental data is demonstrated by modification of the band structure parameters.

###### pacs:

74.25.nd, 74.70.Xa, 74.20.Mn^{†}

^{†}thanks:

## I Introduction

A large variety of compounds belonging to the iron-based superconductors have been found since the discovery of superconductivity in in the beginning of 2008 kamiharadisc (). They can be classified into four different families, divided by stoichiometric differences, exhibiting superconducting (SC) behavior induced by hole-doping, electron-doping as well as solely the application of high pressure (for a review see reviewishida ()). All compounds are layered materials with a common layer of an iron-pnictide in between different spacing layers or simply layers of . So the location of the SC inducing physics is the tetragonal -containing layer leading to highly two dimensional (2D) electronic structure singhcylfsdft1 (). Angle resolved photo-emission spectroscopy (ARPES) finds several bands crossing the Fermi surface arpesnature (); arpespnictidenature (); arpespnictidegap () (FS) which is a novelty since, with the exception of MgB, superconductors have been most often treated as single band materials. The mechanism that leads to superconductivity including the pairing symmetry of the SC order parameter is still under debate (a theoretical review is provided by Ref. mazintheoreticalreview ()). Gap symmetries both with nodes (-wave or -wave) and without nodes (-wave) are still under consideration mazintheoreticalreview ().

Electronic Raman scattering (for a review see reviewraman ()) could clarify some of those unresolved issues. In the light scattering experiment different parts of the Brillouin zone (BZ) can be emphasized with different weights through the selection of incident and outgoing light polarizations with respect of the crystallographic directions. Thus by manipulating the directions of the polarization vectors of the incident and scattered photons in a Raman experiment, it is possible to determine the symmetry of the pair-breaking excitations related to the symmetry of the SC order parameter reviewraman (). Furthermore, the low energy tail of the Raman responses can be used to distinguish between a fully gaped FS or a SC gap with nodes allowing the excitation of quasi particles at arbitrarily small energies.

The theory of electronic Raman spectroscopy in superconductors kleindierker () has been very successful in applications to high-T cuprates reviewraman (); strohmcardona () and now allows to make predictions for the Raman responses of the iron-pnictides. Recently model calculations of expected Raman responses for a superconductor with multi-band order parameter and including the long-range Coulomb interactions responsible for the multi-band screening effects have been performed by Boyd et. al. boydramanpnictide (). These calculations used multiple circular FS sheets as a band structure and for the majority of the work the Raman vertices are obtained by an expansion into FS harmonics for cylindrical FS. Taking the screening of the charged background into account but leaving out vertex corrections capturing effects of final state interactions a variety of different gap distributions on the FS sheets with and without nodes were used. Thus many different forms were obtained for the Raman responses which show the possible scenarios expected from the data. By comparing the singularities in different symmetry channels of such calculated Raman responses with the peak structure of actual Raman data the reconstruction of the momentum and frequency dependence of the superconducting order parameter, critical to understanding of the superconducting mechanism in these materials, is possible. In another study chubukovgapsym (); kleinviewpoint () it has been established by analyses of the vertex corrections that for an extended -wave gap symmetry, the Raman response has a true collective mode resonance peak below the fundamental gap . Since the collective mode is predicted only for this particular pairing symmetry, conscientious polarization resolved low-frequency Raman scattering study of the iron-pnictide materials may provide a way to unambiguously distinguish between various suggested gap symmetries.

Very recently the first electronic Raman scattering experiments on , a member of the 122-family, have been performed hackldata () showing a change of the Raman response in the SC state only in the and channel. A clear peak of -singularity shape in the latter one and a broader peak in the former one at a slightly higher energy are observed and the low energy shape is interpreted as evidence for accidental gap nodes.

In this paper model calculations of Raman responses for a multiband free electron model are performed. In chapter II the general theory of Raman scattering in the SC state is presented and explicit forms of the Raman vertices and an analytical expression for the Raman susceptibility are calculated for a 2D tetragonal symmetry and a constant SC gap. Then the screening correction is taken into account adding an extra term in the channel. Multiband effects in this symmetry channel for two and more than two bands with an example of a four band response are presented in chapter III. The general effects found in chapter III are used in chapter IV to explain the Raman responses for a band structure model inspired by ARPES data. First a constant gap and then a -dependent gap of extended -wave symmetry is used for which the single band structure is investigated with an only angular dependent model gap. Finally the Raman response with the extended s-wave gap is compared to experimental data and it is shown that the band structure model can be modified to come to acceptable agreement. It is important to note that final state interactions leading to vertex corrections are neglected in this work. The consequences of including vertex corrections that could give rise to resonance modes or notably modify the Raman response are discussed in Refs. chubukovgapsym (); blumbergmgb2 (). However, this work focuses on multiband screening effects and features of vertex corrections are left for future publications.

The fact that in the applied assumptions analytical calculations are possible allows us to find the origin of the obtained general mutliband effects in the screened channel. With that the shape of the numerically calculated Raman responses for the particular band structure with extended s-wave symmetry can be explained and further applications to other multiband SC might follow. Furthermore the nature of the main features of the single band Raman responses with the extended s-wave gap are identified analytically with a simplified model gap. Thus our work adds further details to the previous study by Boyd et. al. boydramanpnictide ().

## Ii Theory of electronic Raman scattering in the superconducting state

The Raman process is inelastic scattering of light in which the energy shift between the incident photon (, ) and the scattered photon (, ) is measured. In the case of electronic Raman scattering this two photon process creates charge density fluctuations in the sample and trough a virtual interband transition (in the non resonant case) an electron is excited from a state below the Fermi energy () to a state above. In a clean system the energy difference between those two states corresponds to the energy loss of the photon. It can be shown that the cross section of the scattered light is proportional to the imaginary part of the Raman susceptibility strohmcardona () which is defined in linear response theory in the following way

(1) |

Here , ( is Boltzmann’s constant), is the partition function and is the Hamiltonian of the unperturbated system. The perturbation Hamiltonian is proportional to the effective density operator strohmcardona ()

(2) |

which contains fermionic annihilation and creation operators in the band. Furthermore is proportional to the so called Raman vertex which depends on the polarization direction of the incident and scattered photon and thus on the experimental setup (light polarizers). This quantity will be specified and discussed in detail in the next paragraph. With the wavelength of the incident and scattered photon being in or near the visible spectrum the corresponding wave vectors and therefore also their difference q is much smaller than ( is the SC gap and is the component of the Fermi velocity along q). This makes the limit a reasonable simplification that leads to the explicit form of the Raman susceptibility in the SC state

(3) |

The quantity is the Tsuneto function tsuneto ()

(4) | ||||

with the SC gap , the SC quasi particle energy and a small broadening factor . This function represents the creation of quasi particles out of the SC condensate. The first term is proportional to a SC coherence factor and is unity at the Fermi momentum. The further away from the quasi particles are created the smaller this factor gets so the major contribution to the Raman response originates from the vicinity of the . The second term is a statistical factor that takes into account already thermally excited quasi particles. Since the following calculations are made under the assumption that this factor is simply unity. The third term takes care of the necessary energy conservation and excludes a creation of quasi particles inside the SC gap.

Now the Raman vertex will be specified by an analytical expression. In the nonresonant case, where the electronic excitation does not include a real intermediate state, the Raman vertex can be calculated with an effective mass approximation strohmcardona () in the following form

(5) |

with denoting the electrons mass and and being the light polarization vectors of the scattered and incident photon, respectively. In this approximation the Raman vertex is proportional to the second derivative (curvature) of the band energy dispersion in the particular -direction selected by the two light polarization vectors which is proportional to the inverse effective mass of the band in - theory. Taking into account that a periodic crystal transforming according to a special point symmetry group is investigated group theoretical arguments can be used to further simplify the expression in Eq. (5). Any operator of the form (=initial, =final) can be decomposed according to the point symmetry group of the crystal which will be taken to be in this case. This leads to the decomposition of into projected operators corresponding to the irreducible representation of the point group reviewraman ().

(6) | ||||

which results, together with Eq. (5), in the Raman vertices for the four symmetry channels

(7) | ||||

The different symmetry channels can be investigated by choosing certain combinations of the light polarization vectors through the setup of light polarizers in the experiment. Using an effective mass approximation (Eq. (5)) to obtain the Raman vertices of Eq. (LABEL:eq:parabolverteces) is a different approach than the one in Boyd et al. boydramanpnictide () where an expansion into FS harmonics is done leading to angular dependent vertices. In this work a 2D free electron model for the band structure will be used which results in a constant vertex for the channel since the second derivative of the energy band dispersion is proportional to the constant inverse effective mass of the band. For a parabolic band all other symmetry channels are vanishing which is a consequence of the effective mass approximation in combination with a free electron model. This leads to the focus of the following investigations being on the fully symmetric channel. A tight binding model, which is the next possible step, will not lead to vanishing vertices for the non fully symmetric channels (except , the symmetry of the antisymmetric tensor strohmcardona (), which is usually negligible).

Using a 2D free electron model for the energy dispersion of the band, which is determined by the effective mass and the chemical potential , the Raman susceptibility can be calculated analytically. First the Tsuneto function of Eq. (4) is inserted into Eq. (3), the limit of is taken and the sum is replaced by a 2D -integration. Evaluating only the imaginary part of the integral and inserting the result into the Kramers-Krönig relation leads to the analytical expression

(8) |

with the definition

(9) |

and with and being the 2D density of states (DOS) at which is a constant proportional to (detailed calculation can be found in masterthesis ()). Since the cross section of the Raman process or Raman efficiency is proportional to the one band response is obtained by taking the imaginary part of Eq. (8). For () is vanishing and the square root is entirely real leading to a vanishing Raman response. This is due to the fact that no quasi particles can be created without breaking a Cooper pair which has a binding energy of . At () the denominator hits zero and the function diverges as a square root singularity, a consequence of the diverging SC quasi particle DOS at this energy. Above () the square root is entirely imaginary and in only remains which is a constant of leading to for the single band Raman response.

So far it was not taken into account that in the solid a charged background of electrons is present that reacts through long range Coulomb interaction on the charge density fluctuations imposed by the Raman process. The electron background rearranges in a way to screen the perturbative electric fields and thus works against the Raman process. This screening effect can be included in the theory within a random phase approximation like sum strohmcardona () that leads in the limit of small q (compared to the inverse coherence length and the Fermi wave vector ) to an additional term in the Raman susceptibility. This term includes a k sum over the 1st BZ of that vanishes if is not fully symmetric hence exhibiting sign changes in between different parts of the 1st BZ leading to those parts canceling out each other. Thus only the fully symmetric channel gets screened and the screened Raman susceptibility is

(10) |

with

(11) |

and

(12) |

The second term in Eq. (10) is called screening term and consists of a normalization by a sum over terms shown in Eq. (12) and the square of a sum over terms shown in Eq. (11). Considering only one free electron band the Raman response vanishes.

(13) | ||||

This means that the screening of a free electron background is perfect and all charge density fluctuations are completely compensated by the rearrangement of the free electrons.

## Iii Multiband effects in the screened channel with constant SC gaps

### iii.1 Two bands in symmetry

In contrast to the rather trivial one band case interesting features appear as soon as two free electron bands are considered. Taking a closer look at the second term in Eq. (10) reveals that through the square in the numerator cross terms of the complex functions containing real and imaginary parts corresponding to distinct bands contribute to the two band Raman response. Having two bands furthermore adds complexity by the possibility of choosing different effective masses (magnitude and sign) and SC gaps. The general form of the screened Raman response shown in Eq. (13) extended for two bands can be transformed into a simpler form

(14) | ||||

For the sake of generality the effective masses (determining and ) and the SC gap (present in ) are taken to be different for both bands performing the transformation in Eq. (14). We note here that for a two band model the vertex corrections due to final state interactions can notably modify the Raman response (compare Eq. (14) to Eq. (6) in Ref. blumbergmgb2 ()). The first thing that becomes obvious in the form of the last line of Eq. (14) is that for equal effective masses (equal ’s) the two band Raman response is vanishing independent of the choice for the two SC gaps. So the screening is again perfect.

Further investigations of the different two band features require a case by case study of all possible combinations of the two tunable quantities, namely the effective mass and the SC gap. First the case of two equal SC gaps () is considered. In order to see the effect of screening in this case one has to take a closer look at the screening term which is

(15) |

with

(16) |

In Eq. (15) all common constants are absorbed in . The absolute values of the effective masses in Eq. (16) have their origin in the fact that the DOS is always a positive quantity in contrast to the Raman vertices for which the sign of the effective mass remains. In the numerator of the inverse effective screening mass the absolute value of the effective masses cancels out leaving only a sign dependence. This feature has its origin in the two dimensionality of the band structure for which the DOS is a constant proportional to and will not remain for 1D or 3D systems where the DOS includes an energy dependence. For two bands with equal signs this leads to a finite screening term and thus a partial screening (if the effective masses have different absolute values) while for unequal signs the numerator of in Eq. (16) is zero. Thus for a combination of a hole-like and an electron-like band with equal SC gaps the screening term vanishes and the two band Raman response is entirely unscreened no matter how different the absolute values of the effective masses of the two bands are. This case is displayed in the left column of Fig. 1 where the unscreened Raman response (black) and the screened one (red) show no difference demonstrating the vanishing screening term. Additionally it has to be noted that the square root singularity previously found in the unscreened one band Raman response remains in the screened two band Raman response for equal SC gaps.

Now the case of unequal SC gaps is investigated. Reinserting the explicit form of (see Eq. (9)) into the simplified form of the screened two band Raman response (last line in Eq. (14)) results in

(17) |

with the new short representations

(18) | ||||

The form in Eq. (17) shows that the singularities at , present in the unscreened case, are both removed in the screened case since the square root terms and hit zero separately always leaving one of them finite. This case is illustrated in the right column of Fig. 1 with different combinations of effective masses for the screened Raman response (red, orange, yellow). For equal effective masses (orange) a broad dome, peaked in the middle of twice the two SC gap values, appears. When both bands have different effective masses (red, yellow) this shape changes into a sharper peak located in the vicinity of twice the SC gap corresponding to the band with the smaller effective mass.

### iii.2 More than two bands in symmetry

Having seen the interesting features originating from the interplay of two bands through their cross terms in the screening term the next logical step is to investigate what happens by including more bands into the Raman response. Writing down a general form for more than two bands (analog to the upper two lines of Eq. (14)) it becomes clear that the screening vanishes only in very special cases. Only if all SC gaps are equal and an even number of bands consisting of the same amount of hole-like and electron-like bands is considered the Raman response will be unscreened independent of the absolute values of the effective masses. In all other cases the strength of the screening will depend on the combination of effective masses of the contributing bands. A transformation analog to the one done in Eq. (14) including free electron bands leads to

(19) |

Here , and run from to . This simplified form shows that the screened band Raman response consists of a sum over all possible combinations of two band terms (with ) each normalized by all bands. Furthermore it is obvious that terms including two bands with equal Raman vertices vanish and that the dominating terms in the sum are the ones including two bands with opposite signs since for them the Raman vertices in the bracket add up.

Now the issue of the removal of the singularities previously seen for the two band response with unequal SC gaps will be studied for more than two bands. In order to simplify this investigation the example of a four band Raman response, instead of the general band form in Eq. (19), will be considered first. For four bands there are six possible combinations of two band terms that add up in the sum of the total Raman response. Each of those six terms can be written as

(20) |

with

(21) | ||||

If the two bands and in Eq. (20) have different SC gaps and will hit zero at different energies. Thus the term will remain finite when and vice versa leading to a finite denominator of Eq. (20) for all energies. For equal SC gaps in Eq. (20) and are equal and will vanish at the same energy resulting in a vanishing denominator and hence a singularity at that energy. This means that as soon as there are two bands with equal SC gaps contributing to the four band Raman response there is one term in the sum that is singular at twice this SC gap. An illustration of this is seen in Fig. 2 where different cases for the SC gaps are considered. The unscreened Raman response (black) shows the location of the singularities and the screened one (red) displays which ones get removed through screening. Whenever a gap is present more than once the screened four band Raman response contains a singularity because of the above given explanation. The generalization of the four to the band case can easily be done by adding further terms to the denominator of Eq. (20). The argumentation will remain the same as above and the discovered result that only singularities corresponding to uniquely present gap values are removed holds true.

## Iv Application to a multiband free electron band structure based on ARPES data for an iron-pnictide superconductor

### iv.1 Band structure model

Here the general investigations of free electron Raman responses in the SC state will be applied to a particular band structure model inspired by ARPES measurements. In the band structure was obtained arpespnictidenature (); arpespnictidegap () consisting of two hole like, almost parabolic bands around the -point and a propeller-like structure around the -point. This propeller has an electron-like parabolic band in its center which is surrounded by four hole-like blades. In order to mimic this structure within a free electron model the two bands around the -point and the one around the -point are approximated by a 2D parabola which has the form

(22) |

with . The blades are approximated by 2D ellipsoids, having two different effective masses in the - and -directions, represented by

(23) | ||||

Here , with and , which is extracted from arpespnictidenature (). The effective masses and chemical potentials used for the different bands are given in TABLE 1 privatcom (). Eq. (22) and Eq. (23) with the parameters from TABLE 1 lead to the band structure shown in the left of Fig. 3 and its corresponding FS displayed in the top middle of Fig. 3.

band | inner -par. | outer -par. | blade | |

[] | -2.1 | -6 | 0.8 | -0.61, -2.46 |

[] | 35 | 45 | -17 | 10 |

Combining Eq. (LABEL:eq:parabolverteces) with Eq. (22) and Eq. (23) results in explicit Raman vertices for the used band structure model. For the parabolic bands we obtain

(24) | ||||

and the Raman vertices for the elliptical bands are

(25) | ||||

For the used band structure all bands contribute to the fully symmetric channel while the channel only has a contribution from the elliptical blades and the other symmetry channels completely vanish. This is a consequence of the free electron model in combination with the effective mass approximation used inhere.

At this point it is interesting what the new aspects are that come into play by including elliptical bands. The effect that is a direct consequence of the form of the Raman vertex and the 2D DOS, both different for parabolic and elliptical bands, is the vanishing screening for two bands of opposite signs and equal constant SC gaps. Since the Raman vertex of an elliptical band in symmetry is not just proportional to the inverse effective mass but to the sum of the two different effective masses in the particular -direction and the 2D DOS is proportional to the square root of the product of the two effective masses the inverse effective screening mass differs from Eq. (16). For two elliptical bands it changes to

(26) |

where the plus describes two bands of the same sign and the minus two bands of opposite signs and and for band are denoted by and for a shorter notation. Eq. (26) shows that for two elliptical bands with equal signs the screening is not vanishing (just like it is the case for two parabolic bands) and that its strength depends on the absolute values of the effective masses. Furthermore it can be seen that for two bands of opposite signs and different ratios between and the screening term remains finite while for the same ratio of the effective masses both terms in the numerator of Eq. (26) cancel each other and the screening is vanishing.

### iv.2 Screened Raman response with a constant gap

In order to be able to calculate the Raman response of the introduced band structure shown in Fig. 3 the distribution of the SC gap size along the different bands has to be specified. As a first step the SC gaps are assumed to be constant on each band which makes it possible to directly apply the general investigations of chapter III. ARPES measurements find two different SC gaps arpespnictidegap () of which the smaller one is located on the outer -parabola and the bigger one is present on all the other bands. The magnitude found for the smaller gap is roughly half the size of the bigger one. Thus the distribution of the SC gaps used here consists of a smaller gap on the outer -parabola and a bigger one on the inner -parabola, the -pocket and the blades both related by .

The screened (red) and unscreened (black) Raman response for these assumptions presented in Fig. 4 is reminiscent of the one shown in the lower right of Fig. 2. In both plots the singularity at twice the smaller SC gap is removed since it is uniquely present and the one at twice the bigger SC gap remains because three bands exhibit this same gap value. As previously mentioned the Raman response just has a contribution by the blades leading to a simple square root singularity structure like in the left hand side of Fig. 1 with the singularity located at twice the bigger SC gap. The other two symmetry channels have vanishing Raman vertices and hence show no response at all in this model.

### iv.3 Raman response with an extended s-wave gap

A next step in the calculation of Raman responses for the chosen band structure is to drop the assumption of constant SC gaps and to use one general -dependent SC gap that determines the gap magnitude on each band dependent on its location in the 1st BZ. Experimental results are at this point still controversial even for compounds of the same stoichiometric family (122 in this case) with evidence for gap nodes in some experiments nodesin122a () and measurements of finite gaps having no nodes in other experiments nonodesin122 (); arpespnictidegap (). However the candidate favored by the majority of publications, for example chubukovgapsym (); korshunovspmgap (); mazinspmgap (), is the extended s-wave gap which is represented by the following analytical expression chubukovgapsym ()

(27) |

Here runs from to . This gap exhibits a diamond shaped node in the 1st BZ, shown as a red line in the top middle panel of Fig. 3, which does not coincide with any of the bands at the FS. Thus all bands are fully gaped but there is a variation of the gap magnitude on either one of them. This angular variation on each band at the FS (after a transformation into polar coordinates) is displayed in the bottom right panel of Fig. 3. The extended s-wave gap has its maximum at the -point and its minimum of the same magnitude but negative sign () is located at the -point. Measured values of the maximum gap arpespnictidegap (); Dingarpesgap (); muonspinrotgap () range in between and supporting the choice of for the following calculations.

Combining Eq. (3) and Eq. (4) at and noting that the Raman vertex is a constant in the applied model leads to

(28) |

for the Raman susceptibility of the band. For the numerical calculations the sum is replaced by a 2D -integration in polar coordinates and a constant broadening parameter of representing a clean system is used. The obtained unscreened single band Raman responses are displayed in Fig. 5 in the corresponding coloring to the 1st BZ FS in Fig. 3. The plots show three main features that are seen best in the magenta response of the blade since this is the band exhibiting the largest gap variation at the FS. There is a discontinuous falloff at about corresponding to which denotes the top of the band in the SC state above which no quasi particles can be excited. This feature is seen at higher energies for the other bands (since they have larger ’s) and thus is not present in this plot. At twice the maximum SC gap there is a -singularity and at twice the minimum SC gap a threshold can be seen. These two features are also obtained for the inner -parabola (yellow) but can only be found at a much smaller energy resolution for the other two bands since those have a much smaller gap variation (see lower right of Fig. 3). The nature of the -singularity and the threshold are found in analytical investigations of the imaginary part of Eq. (28) using a simplified gap of the form with running from to . Noting that the integration must only run over the real part of the integrand leads to differences in the integration boundaries at twice the minimum and maximum SC gap resulting in the two distinct features (for a detailed investigation see masterthesis ()).

Knowing the single band Raman response functions is now used to calculate all possible combinations of the screened two band Raman responses and the screened Raman response of all bands. According to Eq. (19) the former ones allow to explain the origin of the features seen in the latter one since it consist of a sum of the renormalized two band terms. In the upper left of Fig. 6 the three combinations of two hole-like bands are shown and in the upper right the corresponding plot for the combinations of the three hole-like bands with the electron-like -pocket is illustrated. The total screened Raman response, containing all bands in the 1st BZ, is displayed in the bottom panel of Fig. 6. Its main feature is a sharp peak at about whose origin lies in the three electron-hole two band responses. They all show a peak slightly below twice the maximum SC gap of the -pocket and those add up to the observed peak in the total Raman response. The three hole-hole two band terms generally show a weaker response since the screening is more effective for them leading to a rather featureless broad continuum down to the energy of the threshold of the blade. Since all bands have different crossing lines with the FS the maximum SC gaps are different on all of them and thus all -singularities are removed corresponding to the argumentation in chapter III B. The Raman response again only consists of the four blades which results in a Raman response entirely equal to the magenta single band Raman response of the blade seen in Fig. 5 but multiplied by . All other symmetry channels are again vanishing.

### iv.4 Comparison with experiment

The recent measurements of Raman responses hackldata () on only partially resemble the obtained Raman responses with an extended s-wave gap. While the negligible contribution of the SC state to the and response in the experimental data is consistent with the vanishing response in those channels for the used band structure model the and responses are quite different. The structure in the response of the data points towards a power law behavior for small energies and thus supports a SC gap with nodes and not the threshold observed in this work. Even stronger disagreement is seen in the channel where no sharp peak appears in the data but a rather broad dome, peaked at higher energies than the peak in the response, is found. However the -singularity shape of the peak of the here presented Raman response is also present in the data.

One way to tune our band structure model in order to get closer to the measured data is to shift the blades further towards the node line of the gap. To be consistent with the measured peak positions the effective masses of all bands, except the outer -parabola, and the chemical potential of the -pocket are changed and the maximum gap value is set to . All this and a larger broadening parameter leads to the FS topology (FS) and the Raman responses illustrated in Fig. 7. The threshold in the channel (green) is now at small energies. The total shape is close to the experimental one but to obtain the proper intensity the shown Raman response is multiplied by 5 since the spectral weight of the channel is quite small in the calculated response. This is a consequence of a lack of contribution by other bands than the blades to the channel and might also be due to some resonance effects in the experiment that are not taken into account in the calculations within effective mass approximation. The Raman response (blue) in the right hand side of Fig. 7 is now also quite close to the measured one except for the still quite obvious peak at about . In general the agreement between the Raman responses of the data on and the modified free electron band structure is quite good which shows that the obtained shape of the data does not necessarily have to be explained with gap nodes. Furthermore it is not clear that the measured compound is in the clean limit still leaving the possibility of rather strong impurity scattering that could also explain the low energy tail of the data. With the introduction of possible scattering centers into the layer (where the important physics for SC is assumed to take place) by the replacement of by this scenario is quite reasonable. Certainly more data on other compounds is needed to clarify these issues.

## V Conclusion

In this work we performed model calculations of Raman responses for a multiband superconductor in a 2D free electron model. For a constant SC gap parameter analytical calculations were used to investigate multiband effects in the screening correction term in the fully symmetric channel. It was shown that for two parabolic bands with equal SC gaps and opposite signs of effective masses the screening term vanishes which also was proven to hold true for two elliptical bands with equal ellipticity. Moreover the removal of the singularities through screening for two bands with different SC gap sizes was presented. This feature remained in the Raman response of multiple bands for all contributing bands exhibiting a uniquely present SC gap size and the example of a four band response was discussed. Additionally is was shown through analytical calculation that multiband responses consist of a sum of two band terms normalized by all contributing bands.

The general effects found for a constant SC gap parameter were then used to explain the SC Raman response of a model free electron band structure inspired by ARPES measurements on iron-pnictide materials. Using a SC gap of extended s-wave symmetry the single band response features, a threshold at twice the minimum SC gap value and a -singularity at twice the maximum SC gap value, were clarified. Furthermore we identified that the main features of the four band Raman response of the model band structure with the extended s-wave gap were originating from the two band terms consisting of bands with opposite signs in the effective masses. A comparison of this four band Raman response with experimental data was done and by modifying the band structure parameters an acceptable agreement was demonstrated.

## Acknowledgements

The authors would like to thank S. Borisenko, D. Evtushinsky and A. Koitzsch for ARPES related discussions and for the band structure parameters. The work of CS was supported in part by the Deutscher Akademischer Austauschdienst (DAAD).

## References

- (1) Y. Kamihara, T. Watanabe, M. Hirano and H. Hosono, J. Ame. Chem. Soc. 130, 3296 (2008).
- (2) K. Ishida, Y. Nakai and Hideo Hosono, J. Phys. Soc. Japan 78, 062001 (2009)
- (3) D. J. Singh, Phys. Rev. B 78, 094511 (2008)
- (4) D. H. Lu, M. Yi, S.-K. Mo, A. S. Erickson, J. Analytis, J.-H. Chu, D. J. Singh, Z. Hussain, T. H. Geballe, I. R. Fisher and Z.-X. Shen, Nature 455, 81 (2008)
- (5) V. B. Zabolotnyy, D. S. Inosov, D. V. Evtushinsky, A. Koitzsch, A. A. Kordyuk, G. L. Sun, J. T. Park, D. Haug, V. Hinkov, A. V. Boris, C. T. Lin, M. Knupfer, A. N. Yaresko, B. Buchner, A. Varykhalov, R. Follath, S. V. Borisenko, Nature 457, 569 (2009)
- (6) D. V. Evtushinsky, D. S. Inosov, V. B. Zabolotnyy, A. Koitzsch, M. Knupfer, B. Büchner, M. S. Viazovska, G. L. Sun, V. Hinkov, A. V. Boris, C. T. Lin, B. Keimer, A. Varykhalov, A. A. Kordyuk and S. V. Borisenko, Phys. Rev. B 79, 054517 (2009)
- (7) I.-I. Mazin and J. Schmalian, arXiv:0901.4790v3 (unpublished)
- (8) T. P. Devereaux and R. Hackl, Rev. Mod. Phys. 79, 175 (2007)
- (9) M. V. Klein and S. B. Dierker, Phys. Rev. B 29, 4976 (1984)
- (10) T. Strohm and M. Cardona, Phys. Rev. B 55, 12725 (1997)
- (11) G. R. Boyd, T. P. Devereaux, P. J. Hirschfeld, V. Mishra and D. J. Scalapino, Phys. Rev. B 79, 174521 (2009)
- (12) A. V. Chubukov, I. Eremin and M. M. Korshunov, Phys. Rev. B 79, 220501(R) (2009)
- (13) M. V. Klein, Physics 2, 46 (2009)
- (14) B. Muschler, W. Prestel, R. Hackl, T. P. Devereaux, J. G. Analytis, Jiun-Haw Chu and I. R. Fisher, Phys. Rev. B 80, 180510(R) (2009)
- (15) G. Blumberg, A. Mialitsin, B. S. Dennis, M.V. Klein, N. D. Zhigadlo, and J. Karpinski, Phys. Rev. Let. 99, 227002 (2007)
- (16) T. Tsuneto, Phys. Rev. 118, 1029 (1960)
- (17) C. Sauer, Master thesis, Rutgers The State University of New Jersey (2009), online at: http://girsh.rutgers.edu/Papers/Thesis/MasterThesisChristophSauer.pdf
- (18) The values are obtained through private communication with A. Koitzsch at Institute for Solid State Research, IFW-Dresden, PO Box 270116, 01171 Dresden, Germany
- (19) R. T. Gordon, C. Martin, H. Kim, N. Ni, M. A. Tanatar, J. Schmalian, I. I. Mazin, S. L. Bud’ko, P. C. Canfield and R. Prozorov, Phys. Rev. B 79, 100506(R) (2009)
- (20) K. Terashima, Y. Sekiba, J. H. Bowen, K. Nakayama, T. Kawahara, T. Sato, P. Richard, Y.-M. Xu, L. J. Li, G. H. Cao, Z.-A. Xu, H. Ding and T. Takahashi, Proc. Nat. Acad. Sci. 106, 7330 (2009)
- (21) M. M. Korshunov and I. Eremin, Phys. Rev. B 78, 140509(R) (2008)
- (22) I. I. Mazin, D. J. Singh, M. D. Johannes and M. H. Du, Phys. Rev. Let. 101, 057003 (2008)
- (23) H. Ding, P. Richard, K. Nakayama, K. Sugawara, T. Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. Takahashi, Z. Wang, X. Dai, Z. Fang, G. F. Chen, J. L. Luo and N. L. Wang, Europhys. Let. 83, 47001 (2008)
- (24) T. J. Williams, A. A. Aczel, E. Baggio-Saitovitch, S. L. Budko, P. C. Canfield, J. P. Carlo, T. Goko, J. Munevar, N. Ni, Y. J. Uemura, W. Yu and G. M. Luke, Phys. Rev. B 80, 094501 (2009)