iron - based high-@xmath0 superconductors@xcite have a complex fermi surface that is determined by several bands , an effect resulting from the hybridization of the @xmath1 orbitals of iron.@xcite band structure calculations have shown that the bands that define the two hole pockets around the @xmath2 point have mostly @xmath3 and @xmath4 character , while the two electron pockets around the m point have @xmath5 , @xmath4 , and a smaller amount of @xmath6 contributions.@xcite for this reason it is important to understand superconductivity in multiorbital systems in general terms .    among the first to address this complex problem several years ago were suhl _
et al._@xcite using a model consisting of two orbitals , one @xmath7 and one @xmath8 , that did not hybridize with each other .
thus , in this case each band was determined by one single orbital .
they showed that bcs pairing@xcite could occur in each band and , since in the most general case the electron - phonon interaction would have different strengths for electrons in the different bands , it was proposed that two different superconducting gaps could arise .
almost 50 years were needed to observe experimental evidence of this phenomenon . in 2001 , superconductivity with @xmath9  k was observed in mgb@xmath10.@xcite despite the high-@xmath11 , it became clear that the bcs mechanism@xcite was at play and for the first time two different superconducting gaps were observed.@xcite as shown in ref .
, the fermi surface ( fs ) is determined by two bands : the @xmath12 band formed by the @xmath13 orbitals of b , and the @xmath14 band constituted by a linear combination of the @xmath15 and @xmath16 b - orbitals .
although three orbitals determine the fs , it is interesting to notice that only two different bcs gaps are observed .
this occurs because two of the three orbitals hybridize with each other and determine one single band , which couples strongly to the lattice phonons .
this opens a large superconducting gap on the @xmath14 fs .
the other orbital , @xmath13 , does not hybridize and forms the @xmath12 band that couples weakly to the lattice phonons determining a second , smaller , superconducting gap at the fs of the @xmath12 band .
thus , the number of different gaps that can arise in a multiorbital system is related to the degree of hybridization among the orbitals .
also note that in this early effort _
interband hopping of pairs of electrons belonging to the same band was included _ but the possibility of interband pairing , i.e. , _ pairs formed by electrons belonging to two different bands _ , was not considered .    in this paper , the subject of superconductivity in multiorbital systems
is revisited , in particular to shed light on the possible symmetry of the pairing operator of the pnictides superconductors .
the motivation is that the pairing operators that have been discussed the most thus far@xcite assume that only intraband pairing should occur , namely the two electrons of the cooper pair belong to the same band.@xcite however , numerical simulations@xcite performed on a two - orbital model@xcite for the pnictides favor an interorbital pairing operator that , when transformed to the band representation , results not only in intraband pairing but it includes interband pairing as well.@xcite for the pnictides , interband hopping of pairs formed between electrons in _ the same _ band is often denoted as `` interband superconductivity''.@xcite the situation discussed in the present paper is _ different _ and involves cooper pairs where the two electrons come from two different bands , which we will call `` interband pairing '' . using symmetry arguments and mean - field approximations ,
the plausibility and physical meaning of such an interband pairing in multiorbital systems will be discussed .
interband pairing has previously been addressed in the context of quantum chromodynamics ( qcd ) and cold atoms,@xcite heavy fermions,@xcite cuprates@xcite , and bcs superconductivity.@xcite in the case of heavy fermions , it was argued that interband pairing could occur if two fermi surfaces arising from different bands are very close to one other,@xcite while in qcd and cold atoms it was presented as a possibility for the case of sufficiently strong attractive pairing interactions , or for weaker attractions among particles with very different masses.@xcite as it will be discussed for a simple model in sec .
[ sec : interband ] , three different regimes , shown schematically in fig .
[ regimes ] , can result from a purely interband pairing as a function of the strength of the pairing potential @xmath17 : ( 1 ) a normal regime where the ground state is not superconducting ( namely in purely interband pairing an infinitesimal attraction does not lead to superconductivity ) , ( 2 ) an exotic superconducting `` breached '' regime where gaps open at the normal fermi surfaces while new fermi surfaces defining regions containing unpaired electrons are created , and ( 3 ) a superconducting regime resembling bcs states , at large attractive coupling.@xcite     of an interband pairing attraction .
the label `` normal '' denotes a non - superconducting state .
the case `` breached '' is an exotic regime with superconductivity and gaps , coexisting with fermi surfaces ( or several nodes ) and electrons that do not pair .
`` bcs '' is the large attraction region , where the ground state resembles that of a bcs superconductor and all electrons participate in the pairing.,title="fig:",width=302 ] 0.3 cm    the paper is organized as follows . in sec .
[ sec : pairings ] , the general form of pairing operators in multiorbital systems will be presented , remarking how the symmetry is determined by the spatial and the orbital characteristics of the operator .
the interorbital pairing operator with @xmath18 symmetry obtained numerically in a two - orbital model for the pnictides is discussed , emphasizing that this operator presents a mixture of intra and interband pairing in the band representation . in sec .
[ sec : interband ] , a simple toy model with pure interband pairing attraction is introduced .
this simplified model is discussed in order to illustrate the effects of interband pairing on observables , such as the occupation number and the spectral functions @xmath19 .
the stability of the interband paired state is also discussed .
the occupation number , spectral functions , and the stability of the @xmath18 pairing state are the subject of sec .
[ sec : interorbital ] which is directly related to the physics of pnictides , while sec .
[ sec : conclusions ] is devoted to the conclusions .
in single - orbital models , the symmetry of a spin - singlet pairing state is completely determined by the properties of its spatial form factor .
more specifically , the pairing operator will have the form : @xmath20 where @xmath21 destroys an electron with momentum @xmath22 and spin projection @xmath14 and @xmath23 is the form factor that transforms according to one of the irreducible representations of the crystal s symmetry group .
thus , @xmath23 determines the symmetry of the operator .
these form factors depend on the lattice geometry and generally they may be very complex . however , in materials with short pair coherence lengths , such as the high-@xmath0 cuprates , the assumption that the two particles that form the pair can be very close to one other is usually made .
the cu - oxide planes in the cuprates have the symmetry properties of the @xmath24 group and the case @xmath25 , which transforms according to the irreducible representation @xmath26 , provides the well - known @xmath8-wave symmetry pairing .    in multiorbital systems , on the other hand
, a spin singlet pairing operator will have both spatial and orbital degrees of freedom and it will be given by @xmath27 where @xmath28 destroys an electron with momentum @xmath22 , in orbital @xmath29 , and with spin projection @xmath14 , @xmath23 is the spatial form factor as indicated above , and @xmath30 is a matrix in the space spanned by the orbitals involved . the dimension of @xmath31 is equal to the number of orbitals that are considered to be of relevance . in this case ,
notice that the symmetry of the pairing operator would in general be determined by the product of the symmetry properties of @xmath32 and the symmetry of the orbital contribution @xmath30 .
only if @xmath30 is the identity matrix do the orbital contribution become trivial , because the identity matrix transforms according to @xmath33 .
thus , this is the only case where @xmath23 fully determines the symmetry of the pairing operator , as in the single - orbital example .
the minimum model for the pnictides considers the two orbitals @xmath5 and @xmath4 , which are strongly hybridized.@xcite all the possible pairing operators , up to nearest - neighbor distance , that are allowed by the lattice and orbital symmetries have been already calculated.@xcite numerical simulations performed on the two - orbital model suggest that the favored pairing operator at intermediate couplings , where the state is both magnetic and metallic ,  @xcite has symmetry @xmath18 and is given by eq .
( [ 2 ] ) with @xmath34 , which transforms according to @xmath33 , and @xmath35 which transforms according to @xmath18  @xcite ( where @xmath36 are pauli matrices ) .
thus , the non - trivial symmetry under rotations arises from the orbital portion of the operator .
this pairing operator has been studied at the mean - field level in ref .  .
in the orbital representation , the bogoliubov - de gennes hamiltonian matrix is given by : @xmath37 with    @xmath38    and @xmath39 where @xmath40 , with @xmath41 being the strength of the pairing interaction and @xmath42 the mean - field parameter obtained by minimizing the energy .
since the two orbitals are hybridized via @xmath43 , in the band representation the hamiltonian matrix becomes : @xmath44 where @xmath45 and @xmath46 are given by @xmath47 @xmath48 and @xmath49 and @xmath50 are the elements of the change of basis matrix @xmath51 given by @xmath52 with @xmath53 .
remember that @xmath45 and @xmath46 , are functions of the momentum @xmath22 , and @xmath54 .
thus , it is clear that in the band representation , in addition to the intraband pairing given by @xmath45 , there is also interband pairing given by @xmath46 . among pairing operators compatible with the symmetry of the model ,
the ones that do not lead to interband pairing not only can not mix orbitals , but have to contain @xmath55 , i.e. , the identity matrix.@xcite one example for such an operator would be the @xmath56 pairing.@xcite    the discussion above describes general properties of hybridized multi - orbital systems . if the orbitals are hybridized , but not related to one another by symmetry , there is no reason to expect that the coupling between the electrons in each orbital , and the interaction that produces the pairing , will have the same strength for all the orbitals and lead to a unit matrix in the orbital sector .
then , it is expected that interband pairing will arise in general and , thus , it is important to understand its consequences , providing the main motivation for the present manuscript .
to address qualitatively the issue of interband pairing , postponing the matters of stability to the next subsection ( sec .
[ sec : stable ] ) , let us consider the following two - bands simplified model with interband pairing : @xmath57 where @xmath58 label two bands that are not hybridized , @xmath14 is the spin projection , and for simplicity @xmath59 which gives parabolic bands that are degenerate at @xmath60 with energy @xmath61 , and with a chemical potential @xmath62 .
this can be considered as a crude representation of the two hole - pocket bands around the @xmath2 point in the pnictides , but more importantly presents a simple toy model where the effects of interband pairing can be studied . as before ,
the parameter @xmath63 is the product of an attractive potential @xmath41 between electrons in the two different bands and a mean - field parameter @xmath42 determined by minimizing the total energy .
the band dispersion without the interaction is presented in the inset of fig .
[ bands ] .    ) , for the indicated values of @xmath64 ( defined in the text ) as a function of the momentum @xmath65 .
the case shown is for @xmath66 , @xmath67 , and @xmath68 .
inset : non - interacting band dispersion for the same parameters . , title="fig:",width=302 ] 0.3 cm    the bogoliubov - de gennes matrix expressed in the basis expanded by @xmath69 has the form : @xmath70 this matrix can be diagonalized becoming @xmath71 in the basis expanded by @xmath72 and the change - of - basis matrix is given by @xmath73 with @xmath74 .
the four energy eigenvalues are @xmath75 , \label{13}\ ] ] where the first positive ( negative ) sign corresponds to the eigenvalues of the upper ( lower ) block , labeled @xmath76 and @xmath77 ( @xmath78 and @xmath79 ) in eq .
( [ 12a ] ) and in fig .
[ bands ] .
the second sign differentiates between the two solutions in each block .
figure  [ bands ] shows the eigenvalues for @xmath80 and @xmath81 .
when @xmath80 , then @xmath82 and @xmath83 , the two bands define two circular fermi surfaces with radius @xmath84 and @xmath85 where they cross the chemical potential ( @xmath62 ) .
this is illustrated schematically in fig .
[ diag ] .    .
@xmath84 and @xmath85 indicate the fermi momentum of the two bands when @xmath80 , while @xmath86 and @xmath87 indicate the position of the fermi momenta for the case of a finite but small pairing potential .
the shaded rings indicate the regions with width @xmath88 and @xmath89 in momentum space where electrons can pair . the white region in between the rings contains unpaired electrons in band 1.,title="fig:",width=302 ] 0.3 cm    notice that number operators can be defined in the two bases that are being considered here , i.e. @xmath90 and @xmath91 , which of course become equivalent when @xmath64=0 .
thus , in the basis @xmath90 the number operator is , @xmath92 while in basis @xmath91 the number operator is ( @xmath93=a , b ) @xmath94 the total electronic occupation of the system is given by @xmath95 then , for @xmath80 we find that @xmath96 for @xmath97 , @xmath98 for @xmath99 since in this region @xmath100 while @xmath101 , and finally @xmath102 for @xmath103 since both orbitals are totally filled with electrons .
these results are represented by the dashed lines in fig .
[ nkab12 ] .     for ( a ) the whole system ; ( b ) band b ,
( c ) band a ; ( d ) orbital 1 ; and ( e ) orbital 2 for the indicated values of the pairing potential @xmath64 , and @xmath66 , @xmath67 , and @xmath68.,title="fig:",width=302 ] 0.3 cm      in the nontrivial case of @xmath64 different from zero , the bands @xmath104 and @xmath105 result from the hybridization of the @xmath106 and @xmath107 bands due to @xmath64 .
it is interesting to observe that an internal gap opens at the crossing of bands @xmath76 with @xmath77 @xmath108 the chemical potential and between @xmath78 and @xmath79 @xmath109 the chemical potential , as indicated with circles in fig .
[ bands ] where results for @xmath81 are displayed .
these are very important differences with respect to conventional bcs calculations were all the action is restricted to the original fermi surfaces : for multiorbital models , gaps can open in other portions of the band structure as well .
while the @xmath104 bands are separated by a gap , the bands @xmath105 still cross the chemical potential determining two fermi surfaces at @xmath110 and at @xmath111 , even in this pairing state ( again , this is different from the one - orbital pairing standard bcs ideas ) .
the new fermi surfaces are shown schematically in fig .
[ diag ] : the interior fs has expanded and the exterior one has contracted . in fact , it will be shown below that the effect of @xmath64 is to try to equalize the two fermi surfaces , as this pairing attraction grows in magnitude . calculating @xmath112 we observe that @xmath113 for @xmath114 and @xmath115 , where @xmath116 is the element of the change - of - basis matrix in eq .
( [ 12b ] ) .
this agrees with the bcs expression for @xmath112 but it jumps discontinuously to @xmath98 for @xmath117 . such a result is shown by the solid lines in fig .
[ nkab12](a ) .
the jumps indicate the existence of the two fermi surfaces , which are here present even in the paired state .
thus , some electrons in the region in between the two fermi surfaces may behave like normal unpaired electrons .
a better understanding of the electronic behavior can be achieved by studying the electronic density in the two bases @xmath91 and @xmath90 .
we find that @xmath118 , as shown in fig .
[ nkab12](c ) , which is the standard bcs behavior in agreement with the fact that bands @xmath104 are separated by a gap .
thus , all the electrons in this band participate in the pairing and they do not have a fs . on the other hand , it can be shown that @xmath119 for @xmath114 and @xmath115 , but for @xmath117 there is a discontinuous change of behavior to @xmath120 [ fig .
[ nkab12](b ) ] .
thus , the fermi surfaces are determined by electrons in the band @xmath78 .
however , note that in this region @xmath121 is an increasing function of @xmath122 while @xmath123 is a decreasing function of @xmath122 which satisfies @xmath124 for all @xmath122 .
this behavior can be better understood by calculating the ( photoemission ) spectral functions @xmath19 , which allow us to obtain @xmath125 .     for the non - interacting case , @xmath80 .
( b ) same as ( a ) but for @xmath81 .
note that there are still fermi surfaces in this intermediate coupling case , namely gaps open far from the chemical potential but not at the fermi level .
( c ) same as ( a ) but for @xmath126 .
the presence of a gap for all the momenta shown is now clear.,title="fig:",scaledwidth=40.0% ]   for the non - interacting case , @xmath80 .
( b ) same as ( a ) but for @xmath81 .
note that there are still fermi surfaces in this intermediate coupling case , namely gaps open far from the chemical potential but not at the fermi level .
( c ) same as ( a ) but for @xmath126 .
the presence of a gap for all the momenta shown is now clear.,title="fig:",scaledwidth=40.0% ]   for the non - interacting case , @xmath80 .
( b ) same as ( a ) but for @xmath81 .
note that there are still fermi surfaces in this intermediate coupling case , namely gaps open far from the chemical potential but not at the fermi level .
( c ) same as ( a ) but for @xmath126 .
the presence of a gap for all the momenta shown is now clear.,title="fig:",scaledwidth=40.0% ] 0.3 cm    in the non - interacting case , shown in fig .
[ akw0](a ) , the spectral function shows two peaks , corresponding to the two bands @xmath106 and @xmath107 , for each value of @xmath122 .
this is the expectation for free electrons in a non - interacting multiorbital system .
the two non - interacting fs s are located where each band passes across the chemical potential .
when the pairing interaction becomes finite , the spectral function develops four peaks for each value of @xmath122 as it can be observed in fig .
[ akw0](b ) for @xmath81 .
this is also the expected result since the bcs interaction generates a `` shadow '' or bogoliubov band for each band present in the non - interacting system .
for example , close to @xmath127 the bands @xmath76 and @xmath78 have almost all the spectral weight , given by @xmath128 which is very close to 1 in this region , and follow a dispersion similar to the non - interacting bands @xmath106 and @xmath107 , while the bands @xmath79 and @xmath77 appear with very small spectral weight given by @xmath129 .
the latter are the bogoliubov or `` shadow '' bands .
these shadow bands appear above the chemical potential for large values of @xmath122 , as expected .
but what happens in the intermediate region @xmath117 ?
we see that the band @xmath130 crosses with @xmath131 at @xmath86 determining a fs .
most of the spectral weight below the chemical potential belongs to the band @xmath131 which contains the unpaired electrons but the band @xmath132 has appreciable shadow spectral weight indicating that there are also some paired electronic population as indicated in figs .  [ nkab12](b ) and ( c ) . as @xmath133 increases
, spectral weight is transferred continuously from @xmath131 to @xmath132 , behavior associated with the internal gap opened by the pairing interaction , so that when @xmath133 approaches @xmath87 most of the spectral weight below the chemical potential is in @xmath132 , paired electrons , and @xmath131 , unpaired electrons , has shadow spectral weight as seen in figs .
[ nkab12](b ) and ( c ) . at @xmath134 the second fs
is determined by the crossing of @xmath135 at the chemical potential indicated by the sudden jump in @xmath136 .
thus , the unpaired electrons in this region coexist with paired electrons and the spectral functions do not resemble the non - interacting ones .    it is also illuminating to analyze what happens with @xmath112 in the basis @xmath90 . while @xmath137 in the regions where there are no unpaired electrons , we find discontinuities associated with fermi surfaces in the two distributions that are given by @xmath138 and @xmath139 for @xmath117 [ figs .
[ nkab12](d , e ) ] .
this indicates that the pairing interaction @xmath64 has been able to promote some electrons from above to below @xmath85 in orbital 2 .
also , electrons have been transferred from their original location in orbital 1 , in the neighborhood of @xmath85 and above @xmath84 , to both orbitals 1 and 2 around @xmath84 .
these are the electrons in 1 and 2 that have become paired ( the pairing is indicated by the shadowed circular regions in fig .
[ diag ] ) .
but the interaction was not strong enough to provide pairing partners to all the extra electrons originally in orbital 1 and , thus , they have been left unpaired in between the two paired regions , as indicated in fig .
[ diag ] .
thus , the interband pairing attraction creates pairs of electrons belonging to different orbitals within an interval @xmath140 around each of the two original fermi surfaces .
the width of the pairing region increases with @xmath64 .
the pairing partners are obtained by promoting electrons with momentum @xmath141 and @xmath142 in both orbitals and by moving electrons from the more populated to the less populated orbital .
this creates the conditions to pair electrons near both fermi surfaces .
the electrons in band 1 that could not find promoted partners remain unpaired . whether this state is stable or not depends , of course , on the balance between kinetic and pairing energies , which will be discussed in sec .
[ sec : stable ] .      as the interaction @xmath64 increases further ,
the number of unpaired electrons is reduced .
this means that @xmath86 and @xmath87 become closer to each other making the size of the intermediate region with unpaired electrons in fig .
[ diag ] smaller .
eventually , the two momenta become the same @xmath143 for @xmath144 , and for @xmath145 the region with unpaired electrons vanishes and a full gap opens in the system whose physics now resembles bcs , except for the fact that the pairs are constituted by electrons from different orbitals .
the electronic population of the system in such a case , e.g. at @xmath126 , is presented in fig .
[ nk1 ] and the corresponding spectral functions are shown in fig .  [
akw0](c ) .
it can be shown that now @xmath146 for all values of @xmath133 and all the particles around the two noninteracting fermi surfaces now participate in the pairing .
the spectral functions show four peaks , i.e. bogoliubov bands for all values of @xmath133 , as it can be observed in fig .
[ akw0](c ) .    , and for @xmath66 , @xmath67 , and @xmath68 .
the case @xmath126 illustrates the `` strong '' pairing attraction regime where gaps open in the original fermi surface.,title="fig:",width=302 ] 0.3 cm    thus , notice that while band 1 contained more electrons than band 2 for @xmath80 [ see figs .
[ nk1](b ) and ( c ) ] , the interband pairing mechanism transfers electrons from one band to the other so that both bands have the same number of electrons in the superconducting state .
consequently , the smaller fs expands and the larger one shrinks .
then , when the pairing becomes strong enough , the two fermi surfaces become equalized and no unpaired electrons remain . whether this situation can be achieved will depend on the strength of the interaction and the energy balance , as discussed in the next subsection .
if the two non - interacting fermi surfaces are very close to each other in momentum space , even relatively weak pairing interaction could effectively be strong enough to make the interband pairing resemble bcs pairing as in the case of large @xmath64 in the present example .
as discussed in the introduction , the possibility of interband pairing has been previously discussed in the context of qcd and cold atomic matter .
similar effects on the fs as found in the present study were described , although the physics was different because in the qcd context each band contained different kinds of particles and the pairing was thus not able to promote particles from the majority to the minority band.@xcite the issue of stability was explored in the qcd framework , and it was found that a purely interband paired state could be stabilized for pairing attractions @xmath108 a certain cut - off value , which could become very small for a large difference between the masses of the two paired species.@xcite    in the case of our model , however , we have found ( see below ) that the purely interband - paired state only becomes stable when the attraction is sufficiently strong that no unpaired particles are left , i.e. when the two shaded regions overlap and the unpaired region in fig .
[ diag ] vanishes .
this means that , although the pairs would be formed by electrons in different orbitals , the physics would be analogous to bcs . a gap will be opened in the full fermi surface of the simple model studied here .    in order to study the issue of stability ,
let us assume that the interaction term responsible for the interband attraction is given by @xmath147 where @xmath148 and @xmath149 is the number of sites . performing the standard mean - field approximation : @xmath150 and @xmath151 and making the substitution @xmath152 ( and an analogous substitution for the product of annihilation operators ) ,
the mean - field results are obtained .
as usual , the fluctuations around the average given by @xmath153 are assumed to be small .
defining @xmath154 we obtain the following mean - field hamiltonian : @xmath155 equation  ( [ eq:12 ] ) can be recovered by defining @xmath156 , and disregarding the constant last term of eq .
( [ mf ] ) .
we can calculate the total energy @xmath157 for eq .
( [ mf ] ) as a function of @xmath40 .
if for a given @xmath158 the energy has a minimum , this indicates that the interband - paired state is stable . note that having a term linear in @xmath42 in eq.([mf ] ) is not sufficient to conclude the appearance of a superconducting state at small @xmath42 , since the sign of the coefficient of the linear term can change sign with @xmath64 .
a similar situation occurs for the magnetic state of undoped pnictides : a finite hubbard @xmath51 must be reached to stabilize the `` striped '' state @xcite .
returning to superconductivity , there are two regions of interest : ( i ) @xmath159 , which corresponds to the case in which two fermi surfaces are present in the paired state ; ( ii ) @xmath145 , which corresponds to the fully gapped case . in fig .
[ fig : stab](a ) , @xmath160 vs. @xmath64 is shown for different values of @xmath41 .
it can be observed that a second minimum develops for @xmath161 and it becomes stable for @xmath162 .
the minimum always occurs for @xmath145 which means that it corresponds to the case in which there are no unpaired electrons in the system .
the results shown in the figure are robust in the sense that changes in the values of @xmath163 , @xmath164 or the chemical potential were not found to stabilize the state with unpaired electrons .
thus , in this respect an attraction that is only interband can only lead to a stable superconducting state in the strong attraction region .    ) per site vs. @xmath40 for different values of @xmath41 .
the case @xmath144 is indicated with dashed lines : it separates the regions corresponding to nodal and nodeless states .
( b ) mean - field energy ( eq .  [ mf ] ) per site with intraband pairing with strength @xmath165 vs. @xmath40 , for different values of @xmath41.[fig : stab],title="fig:",scaledwidth=40.0% ] ) per site vs. @xmath40 for different values of @xmath41 .
the case @xmath144 is indicated with dashed lines : it separates the regions corresponding to nodal and nodeless states .
( b ) mean - field energy ( eq .  [ mf ] ) per site with intraband pairing with strength @xmath165 vs. @xmath40 , for different values of @xmath41.[fig : stab],title="fig:",scaledwidth=40.0% ]      the results of the previous paragraphs may seem negative with respect to the relevance of the `` intermediate '' state with simultaneous coexistence of pairing and fermi surfaces .
however , as pointed out in sec .
[ sec : pairings ] , most of the pairing operators allowed by the lattice and orbital symmetry in the pnictides are characterized by a _
mixture _ of both intra and interband pairing .
thus , it is important to consider such a situation in our simple model as well . in the case of the @xmath18 pairing operator , eqs .
( [ 6 ] ) , ( [ 7 ] ) , and ( [ 8 ] ) indicate that the pairing is purely interband only for @xmath166 or @xmath12 and @xmath167 or @xmath12 , because @xmath168 along these lines . thus , let us now consider our simple model in a brillouin zone ( bz ) defined by @xmath169 and with the addition of _ intraband _ pairing with intensity @xmath165 that is vanishing for @xmath166 or @xmath12 and @xmath167 or @xmath12 , i.e. , the pairing is purely interband only along those directions .
the energy bands now behave as in fig .
[ bands ] only along @xmath170 and @xmath171 , while the dispersion along any other direction is shown in fig .
[ bandasmix ] for @xmath80 and @xmath81 .    )
( along the main diagonal ) for the indicated values of the pairing potential @xmath64 , with the addition of an intraband pairing with strength @xmath165 , as described in the text.,title="fig:",width=302 ] 0.3 cm    performing the mean - field approximation similarly as explained above , we have found that the superconducting state now becomes stable on _ both _ sides of the original critical value @xmath144 .
figure  [ fig : stab](b ) shows that the pairing state is stabilized for @xmath172 but the value of @xmath42 where the minimum is located is such that @xmath173 and , thus , two nodes will be present along the @xmath174 and @xmath175 axes .
increasing the value of @xmath41 , the minimum eventually occurs for @xmath145 . for these larger values of @xmath41
, there would consequently be no nodes .
then , in this section it has been shown using a simple model that the interorbital paired state can become stable if the attraction @xmath41 is sufficiently strong . in this case
, it is the nodeless case that is stable even for purely interband attraction .
in addition , the very interesting novel phase with coexisting nodes and unpaired electrons in the majority band also requires intraband pairing to be stable , with strength similar to that of the interband , at least in parts of the bz .
in the previous section , a simple model was presented , both with exclusively interband pairing and with both inter and intraband pairing , and it was found that the intraband pairing stabilizes the state with a mixture of superconductivity and metallicity . as mentioned in the introduction , it is expected that in the most general cases the pairing operators allowed by the symmetry of the lattice and of the orbitals will , in the band representation , have both intra and interorbital pairing .
thus , now we will present and discuss , at the mean field level , the occupation number and the spectral functions for the pairing operator obtained from the numerical study of the two - orbital model for the superconducting state of the pnictides introduced in sec .
[ sec : pairings ] .
for the two - orbital model eq .
( [ eq : bcs_ours ] ) with pairing interaction @xmath81 ( continuous lines ) and @xmath80 ( dashed lines ) .
parameters : @xmath176 , @xmath177 , @xmath178 , @xmath179 .
( b ) same as ( a ) but for the orbital @xmath5 .
( c ) same as ( a ) but for the orbital @xmath4 .
( d ) same as ( a ) but for band 1 .
( e ) same as ( a ) but for band 2.[fig : nkxy05],scaledwidth=40.0% ]       for the two - orbital model eq .
( [ eq : bcs_ours ] ) without pairing interaction for the orbital @xmath5 ( orange / light ) and the orbital @xmath4 ( red / dark ) .
the ` floor ' indicates the fs in red .
( b ) same as ( a ) but for band 1 ( orange / light ) and band 2 ( red / dark)[nkxysurv0],title="fig:",scaledwidth=40.0% ]   for the two - orbital model eq .
( [ eq : bcs_ours ] ) without pairing interaction for the orbital @xmath5 ( orange / light ) and the orbital @xmath4 ( red / dark ) .
the ` floor ' indicates the fs in red .
( b ) same as ( a ) but for band 1 ( orange / light ) and band 2 ( red / dark)[nkxysurv0],title="fig:",scaledwidth=40.0% ]    let us start with the non - interacting case in which @xmath80 .
the spectral functions along high - symmetry directions in momentum space are presented in fig .
[ fig : akb2g_00 ] and , as expected , they reproduce the non - interacting band dispersion .
we also show the total occupation number @xmath112 along the same directions [ dashed lines in fig .
[ fig : nkxy05](a ) ] , as well as the occupation number for each orbital @xmath180 [ dashed lines in panel ( b ) ] and @xmath181 [ dashed lines in panel ( c ) ] and the orbital occupation in the quadrant of the first bz defined by @xmath182 in fig .  [ nkxysurv0](a ) .
it is clear that the electron - like and hole - like fermi surfaces are determined by an admixture of the two orbitals .    on the other hand , in the band representation
, band 1 determines the electron pockets while band 2 forms the hole pockets as it can be seen from the behavior of @xmath183 and @xmath184 ( indicated by the dashed lines in panels ( d ) and ( e ) of fig .
[ fig : nkxy05 ] ) and by the light ( orange ) and dark ( red ) surfaces in fig .
[ nkxysurv0](b ) where the fs is also indicated .
it is clear that the electronic occupation of band 1 is smaller than the electronic population of band 2 , so that unpaired electrons would be expected to belong predominantly to band 1 , as in the simple model of the previous sec .
[ sec : interband ] .
it is interesting to notice that in the orbital representation , on the other hand , the electrons are equally distributed among the @xmath185 and @xmath186 @xmath1 orbitals .
let us discuss what occurs when the pairing interaction becomes nonzero . to simplify the discussion , define the following points in momentum space : @xmath187 , @xmath188 , @xmath189 , and @xmath190 . as remarked in ref .  , the @xmath18 pairing operator always has nodes along the @xmath191 direction because the spatial form factor @xmath192 vanishes along that line .
but , as soon as @xmath64 is finite , a gap opens along the @xmath193 direction [ notice that along this direction the pairing is purely intraband since in eq .  ( [ 6 ] ) , @xmath194 and thus @xmath195 in eq .
( [ 8 ] ) ] .
along @xmath196 , @xmath197 , @xmath198 , and @xmath199 nodes associated to the different number of electrons in band 1 and band 2 remain [ notice that along these directions the pairing is purely interband since @xmath168 because eq .
( [ eq : bcs_ours ] ) and eq .
( [ 6 ] ) become identical to one other ] . when the pairing interaction @xmath64 becomes strong enough to make @xmath146 , as described in the simplified model presented in the sec .
[ sec : interband_strong ] , these nodes vanish . along any other direction in the bz
a mixture of intra and interorbital pairing will be present.@xcite    let us first consider a relatively small pairing @xmath81 . in fig .
[ fig : nkxy05](a ) , @xmath112 is presented along high symmetry directions . along @xmath200
there is no pairing and , thus , @xmath112 is unchanged from the non - interacting case shown in the figure with dashed lines . along @xmath198 ,
where only interband pairing occurs , no effects are observed at the electron pocket fs but a rounding in @xmath112 indicating pairing is observed at the hole pocket fs .
however , @xmath112 shows discontinuities at two points indicating the existence of nodes . along the diagonal direction @xmath201 , where all the pairing is intraband
, it is found that @xmath112 exhibits standard bcs behavior at both hole fermi surfaces indicating the opening of gaps .
finally , along @xmath196 it can be observed a rounding of @xmath112 at the hole fermi surfaces , indicating pairing , and a sharp jump at the electron fs .
we can further analyze the pairing in the orbital representation .
the occupation number for the orbitals @xmath185 and @xmath186 is shown in figs .
[ fig : nkxy05](b - c ) . along @xmath200 , where the pairing is zero ,
we observe how the fs for the electron pocket at @xmath202 @xmath203 is totally determined by electrons in the orbital @xmath185 @xmath204 and how the population of each orbital varies smoothly between the two fermi surfaces , always satisfying @xmath205 . from @xmath206 to @xmath207 , @xmath208 for @xmath209 indicating that the electrons at the hole pocket fs are paired in a wide region around it , but the pairing region is very narrow around the electron pocket because the pairing is reduced by the small value of @xmath23 . along the diagonal ,
i.e. from @xmath207 to @xmath2 , standard intraband pairing occurs at both hole fermi surfaces and , thus , @xmath208 while from @xmath2 to @xmath206 @xmath208 for @xmath210 indicating pairing around the hole fs and almost no pairing occurs at the electron pocket fs .
the behavior @xmath211 for all @xmath133 for which pairing occurs and @xmath212 unchanged from the non - interacting value for all other @xmath133 , was observed for all values of @xmath64 .
for this reason , figures for @xmath212 in the orbital representation will not be shown for the additional values of @xmath64 discussed below .
the population of the different bands is presented in fig.[fig : nkxy05 ]  ( d - e ) .
the fermi surfaces along @xmath200 are clearly determined only by band 1 , which is the band that forms the electron pockets .
it can also be observed that there is almost negligible pairing at the electron fs along @xmath198 , but there is clear interband pairing at the hole fs .
this is an indication that , due to the spatial variation of the pairing interaction , the attraction is much stronger at the hole pockets than at the electrons pockets .
also notice that at the hole pocket fs the pairing is interband and thus @xmath146 , but this does not happen along the diagonal direction @xmath201 where intraband pairing occurs , and there are more paired electrons belonging to band 2 than to band 1 . along the direction @xmath213
again we observed a stronger pairing effect at the hole fs than at the electron one .      in this section ,
we discuss the form of the spectral functions @xmath214 , which can be measured in angle resolved photoemission spectroscopy ( arpes ) experiments , for weak to strong interorbital pairing .
the spectral function for @xmath81 is depicted in fig .
[ fig : akb2g_05 ] .
as discussed above , the interorbital @xmath18 operator leads to _ intra_-band coupling along the @xmath193 line , where one consequently clearly sees the hole pockets to be gapped . along @xmath196 ,
the pairing is purely _
inter_-band , and one finds a fermi surface on both the hole and electron pockets , indicating that @xmath81 corresponds to the `` breached '' phase in fig .
[ regimes ] .
however , the spectral weight determining the hole - pocket node is weak and might thus be missed in the analysis of experiments .
the node on the electron pocket fs , on the other hand , is robust and it should be observed , if present . the same occurs along @xmath198 ; the signal for the node at the electron pocket fs should be robust while the one at the hole pocket fs will be weak .
as the pairing interaction increases , the nodes resulting from the interband pairing should get closer to each other , as discussed in sec .
[ sec : interband ] .
figures  [ fig : nk36](a - c ) show @xmath112 for @xmath215 along high - symmetry directions in the band representation .
one can see that there are still unpaired electrons , and @xmath215 consequently falls into the `` breached '' region schematically represented in fig .
[ regimes ] .     for the two - orbital model eq .
( [ eq : bcs_ours ] ) with pairing interaction @xmath215 ( continuous lines ) and @xmath80 ( dashed lines ) .
parameters : @xmath176 , @xmath177 , @xmath178 , @xmath179 .
( b ) same as ( a ) but for band 1 .
( c ) same as ( b ) but for band 2 .
( d ) total occupation number @xmath112 for the two - orbital model eq .  ( [ eq : bcs_ours ] ) with pairing interaction @xmath216 ( continuous lines ) and @xmath80 ( dashed lines ) .
( e ) same as ( d ) but for band 1 .
( f ) same as ( e ) but for band 2.[fig : nk36],title="fig:",scaledwidth=40.0% ]   for the two - orbital model eq .  ( [ eq : bcs_ours ] ) with pairing interaction @xmath215 ( continuous lines ) and @xmath80 ( dashed lines ) .
parameters : @xmath176 , @xmath177 , @xmath178 , @xmath179 .
( b ) same as ( a ) but for band 1 .
( c ) same as ( b ) but for band 2 .
( d ) total occupation number @xmath112 for the two - orbital model eq .  (
[ eq : bcs_ours ] ) with pairing interaction @xmath216 ( continuous lines ) and @xmath80 ( dashed lines ) .
( e ) same as ( d ) but for band 1 .
( f ) same as ( e ) but for band 2.[fig : nk36],title="fig:",scaledwidth=40.0% ]    it is interesting to notice that while @xmath217 , on the other hand @xmath218 indicating that the reconstruction around the hole pockets is much larger than around the electron pockets .
this can also be observed in fig .
[ fig : nk36](b ) , where we observe unpaired electrons along @xmath198 , but not along @xmath196 .
the effect occurs in part due to the smaller value of @xmath23 at the electron pockets but also because , due to the band dispersions , the price in kinetic energy for interband pairing is much larger at the electron pockets than at the hole pockets .
the spectral density in fig .  [ fig : akb2g_3 ] further illustrates that the pairing interaction is more effective along @xmath196 than along @xmath198 . only shadow spectral weight crosses the chemical potential along @xmath196 , while strong spectral weight crosses the chemical potential twice along @xmath198 and leads to two nodes .      finally , let us consider a much stronger value of the pairing , such as @xmath216 , in the bcs region of fig .
[ regimes ] , for which the only nodes observed are those along @xmath191 , due to the vanishing of the pairing operator . figures  [ fig : nk36](d - f ) , show that @xmath112 is discontinuous only along @xmath191 , while it is smooth along all the other directions indicating pairing .
note that @xmath146 along @xmath196 and @xmath198 where interband pairing occurs , while along @xmath193 , @xmath212 is smooth but different for each band because the pairing is intraband .
the behavior of the spectral functions displayed in fig .
[ fig : akb2g_6 ] shows that spectral weight only crosses the chemical potential along the @xmath191 direction .
in this situation , in the folded bz , nodes should occur only at the points where the two electron pockets cross with each other , as indicated in ref .  .
in the rest of the bz an anisotropic gap will be observed .
finally , let us discuss the important issue of the stability of the @xmath18 pairing state .
it will be assumed , following the notation in appendix a of ref .  ,
that eq .
( [ eq : bcs_ours ] ) has arisen from an interorbital attractive potential of the form @xmath219 and that @xmath220 where @xmath221     for different values of the attraction @xmath222 as a function of @xmath223.[fig : ebg],scaledwidth=40.0% ]    @xmath64 in eq .
( [ 5 ] ) is then given by @xmath223 and the mean - field energy @xmath157 can be calculated for a given @xmath222 as a function of @xmath64 .
the results are presented in fig .
[ fig : ebg ] . in this figure
, a minimum for @xmath158 can be seen for @xmath224 .
this indicates that if the pairing attraction overcomes a finite threshold level , the mean field results presented here with @xmath225 will be stable .
thus , _ if the attraction exists , both multinodal ( breached ) and states with nodes only along x - y ( strong coupling ) are possible . _
summarizing , in this manuscript the possibility of intra and interband pairing in multiorbital systems has been discussed .
while interband pairing has previously been studied in the context of qcd , cold atoms,@xcite heavy fermions,@xcite cuprates@xcite , and bcs superconductivity,@xcite most of the pairing operators proposed for the pnictides are based on the premise that the pairing has to be purely intraband , i.e. , both electrons in the cooper pairs belonging to the same band , compatible with the assumption that in these superconductors all the action must occur at the fermi surfaces .
however , symmetry considerations show that in models for the pnictides interorbital pairing is allowed with the two members of the cooper pair belonging to different bands.@xcite this is not surprising since the bands that determine the electron and hole fermi surfaces consist of hybridized orbitals . these interorbital pairing operators give rise to not only intra but also interband pairing when the band representation is used .
in addition , numerical calculations in a minimal two orbital model for the pnictides favor one of these non - trivial pairing operators.@xcite as a consequence , a clear discussion of the role of interband pairing is necessary .
the explicit calculations shown here of the electronic occupation in the orbital and the band representations , as well as the calculation of the spectral functions showing the distribution of bogoliubov bands , may offer guidance in the interpretation of arpes experiments .
in particular , most arpes measurements determine the fs in the normal state and study the opening of the superconducting gap by monitoring @xmath226 as they lower the temperature.@xcite notice that this approach would miss the nodes associated with the `` breached '' phase since in the superconducting state a gap would be observed in @xmath226 while the node would be detected in @xmath227 .
thus , experimentalists should investigate the possibility of nodes at points in momentum space that do not belong to the normal state fs and they must keep in mind that some of the nodes may be determined by shadow bands with very small spectral weight .    the most recent experimental results with the polarization dependence of the arpes spectra for bafe@xmath228co@xmath229as@xmath10 provided the allowed contribution of each of the five @xmath1 orbitals to the electron and hole fermi surfaces.@xcite while discrepancies with proposed four and five orbital models are remarked in that publication , it is interesting to observe that the minimal two - orbital model addressed here does not contradict the arpes findings if we disregard the additional @xmath230 hole pocket fs that they present , which is reasonable because it has @xmath231 character , an orbital not included in the minimal model considered here .
in fact , along their @xmath193 direction , which corresponds to our @xmath196 , our first hole fs is purely @xmath5 as it is their @xmath232 hole fs ; our second hole fs , which arises upon folding our extended fs along @xmath191 is purely @xmath4 , as it is their @xmath233 hole fs , and our electron fs is purely @xmath4 as it is their @xmath234 electron fs .
our second electron fs ( obtained upon folding ) has a purely @xmath5 character , while their electron fs @xmath235 appears to be mostly @xmath4 and @xmath6 , but with some amounts of @xmath5 as well .
along the diagonal direction @xmath196 , which corresponds to our @xmath193 , the two hole fermi surfaces are a symmetric admixture of @xmath5 and @xmath4 , exactly as in the two - orbital model.@xcite thus , these similarities between the experimental results and the band composition of the simple two - orbital model offer encouragement towards exploring whether the main physics of the pnictides can be captured with such a minimum number of degrees of freedom .    we have also mentioned in the text that the symmetry of the lattice and of the orbitals introduce constraints on the possible pairing operators . in general , a purely intraband pairing , such as the proposed @xmath236 state,@xcite would occur only if the coupling of the electrons with the source of the attraction is _ identical _ for all orbitals . considering the different spatial orientations of the orbitals ,
it is not obvious that this should be the case , since , as discussed in the introduction , phonons couple differently to electrons in the @xmath13 and @xmath15 boron orbitals in the case of mgb@xmath10 .
thus , if indeed the coupling results to be the same for all @xmath1 orbitals we could use this fact to elucidate the pairing mechanism ; but , if this is not the case , we would expect some degree of interband pairing , at least in some regions of the brillouin zone and , thus , it is important to study the experimental and theoretical consequences of such a possibility .
it is interesting to observe that in the undoped limit there are also three regimes : normal , coexisting magnetism and metallicity , and magnetism without a fermi surface .
for details see r. yu , k. t. trinh , a. moreo , m. daghofer , j. a. riera , s. haas , and e. dagotto , phys .
b * 79 * , 104510 ( 2009 ) .
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 , epl * 83 * , 47001 ( 2008 ) .
l. wray , d. qian , d. hsieh , y. xia , l. li , j.g .
checkelsky , a. pasupathy , k.k .
gomes , c. v. parker , a.v .
fedorov , g.f .
chen , j.l .
luo , a. yazdani , n.p .
ong , n.l .
wang , m.z .
hasan , phys .
b * 78 * , 184508 ( 2008 ) .
v. b. zabolotnyy , d. s. inosov , d. v. evtushinsky , a. koitzsch , a. a. kordyuk , j. t. park , d. haug , v. hinkov , a. v. boris , d. l. sun , g. l. sun , c. t. lin , b. keimer , m. knupfer , b. buechner , a. varykhalov , r. follath , and s. v. borisenko , nature * 457 * , 569 ( 2009 ) .
y. zhang , b. zhou , f. chen , j. wei , m. xu , l.  x. yang , c. fang , w.  f. tsai , g.  h. cao , z.  a. xu , m. arita , c. hong , k. shimada , h. namatame , m. taniguchi , j.  p. hu , and d.  l. feng , arxiv:0904.4022 , 2009 .