recent discovery of the fe - based superconducting compounds provided a great impetus to the research of superconductivity ( sc ) . since the first report on the superconducting transition at 7k with the doped laofep by kamihara et al .
@xcite , various substitutions ( mainly , p by as and la by ce , gd , sm , pr ) for this mother compound increase the superconducting temperature @xmath3 over 50k with sm(of)feas @xcite .
intensive investigations by many experimentalists and theorists have already revealed main metallic and superconducting properties of this group of materials .    from band calculations@xcite , it is agreed on that the 3d electrons of fe atoms are the main contributors to the conduction bands crossing the fermi surface ( fs ) . besides the degree of degeneracy
, the key feature of conduction bands is that it consist of hole band(s ) around @xmath4 point and electron band(s ) around @xmath5 point [ in the notation of the folded brillouin - zone ( bz ) scheme@xcite ] .
regarding the pairing symmetry , there are already many experiments : ( 1 ) knight shift below @xmath3 shows a clear drop indicating a spin singlet pairing @xcite ; ( 2 ) tunneling spectroscopy @xcite showed the zero - bias conductance peak ( zbcp ) - signature of a sign changing gap , but the interpretation of the shape of density of states ( dos ) @xmath6 is diverse ; ( 3 ) nuclear spin - lattice relaxation rate @xmath2 @xcite unanimously showed no coherent peak and @xmath7 dependence below @xmath3 , hence strongly suggesting a d - wave type gap ; and ( 4 ) specific - heat coefficient @xmath8 below @xmath3 ( ref .
@xcite )  although the measurement is not yet reaching low enough temperature  appears @xmath9 linear indicating the gap with lines of node .
all these experiments appear to be consistent with a d - wave type gap .
however , recent penetration depth measurements with prfeaso , sm(of)feas and nd(of)feas ( ref.@xcite ) strongly suggest a fully opened gap at low temperatures indicating a s - wave type pairing symmetry .    regarding the paring glues , the phonon interaction appears unlikely mainly because the electron - phonon coupling is estimated to be very weak ( @xmath10)@xcite . on the other hand , this series of materials , without doping , commonly has a spin density wave ( sdw ) transition at around @xmath11 k. when the superconductivity appears with doping , the sdw correlation is expected to remain , albeit the long range order disappears
recent neutron - scattering experiments with la(of)feas and ce(of)feas ( ref .
@xcite ) directly measured the antiferromagnetic ( afm)-type correlation of the fe d - electron spin moment .
the overall phase diagram with doping for ce(of)feas reveals a close correlation with an antiferromagnetism and superconductivity , suggesting the important role of magnetic fluctuations as a pairing glue .
it also shows that the generic phase diagram of these compounds shares the universal features with the high-@xmath3 cuprates , pu-115 superconductor @xcite , and various heavy fermion superconductors ; namely , the sc occurs in the neighborhood of the magnetic long range order when this magnetic order is suppressed . in particular
, the magnetic order is an afm type .
this universal phase diagram is very tantalizing because it appears to cover a wide class of unconventional sc materials with a range of @xmath3 from a few mk to 100 k and suggests that the afm fluctuation is a common thread and its characteristic energy roughly scales with the sc @xmath3 @xcite .    for the fe - based sc materials ,
several theoretical models were already proposed and most of them started with the orbital basis of the fe 3d - electrons including hubbard @xmath12 interaction(s ) and hund coupling(s ) @xmath13 @xcite . some of these studies @xcite found that the @xmath0s - wave gap as a dominant instability .
a d - wave gap also often appears as a second instability @xcite . in this paper , we took a phenomenological approach to investigate possible pairing states in the fe - pnitide superconductors .
the non - interacting part of hamiltonian is constructed by choosing a minimal set of topologically distinct two bands and the interaction part of hamiltonian is assumed from the experimental input @xcite , simulating a short range afm correlation . by solving the coupled bcs gap equations , we found the two sc gap solutions : a sign - changing s - wave gap and a double d - wave gap .
for the both sc states , we carried out the calculations of the sc properties such as the dos , temperature dependencies of the spin - lattice relaxation rate @xmath2 , knight shift , and penetration depth .
we particularly take into account of the interband coherence factors , unique to the two band model , in these calculations .
the results are discussed in comparison with the currently available experimental data .
( red ) and @xmath14 ( green ) bands and two gap solutions considered in the paper : ( a ) @xmath0s - wave gap and ( b ) double d - wave gap .
the width of the fs represents the local dos for each band within @xmath15 energy .
[ fig1],width=302 ]
we propose a minimal phenomenological two band model for the fe - based superconductors . for the non - interacting part of hamiltonian
, we observe that several band structure calculations @xcite of the fe - pnictide compounds reached the consensus that the fs of the doped compounds consists of two hole pockets and two electron pockets . to keep the essential physics , but avoiding unnecessary complexity , we choose only two topologically distinct bands : one hole band around @xmath4 point ( 0,0 ) and one electron band around @xmath5 point ( @xmath16 ) .
the main phenomenological assumption of our model is the interacting part of hamiltonian .
this pairing interaction @xmath17 is chosen to simulate a short range afm spin fluctuations peaking at the ordering wave vector @xmath18 .
this assumption is directly motivated by the experimental observations of the afm correlation in the fe - pnictide compounds by the neutron - scattering measurements @xcite .
the elastic neutron - scattering experiments for la(f@xmath19o@xmath20)feas and ce(f@xmath19o@xmath20)feas by dai and co - workers @xcite showed the long - range afm order of the fe3d - electron spins for the doping range of @xmath21 .
when this long - range afm order disappears beyond the doping around x=0.05 , sc appears up to the doping range x @xmath22 0.2 ( this is only the limit of the measured data ) .
this overall phase diagram appears to be generic for the fe - pnictides , and clearly shows that the afm correlation is the dominant magnetic correlation in this group of materials .
more importantly , this afm correlation is expected to continue to exist in the doping range where the sc phase occurs albeit becoming a short - rang one .
this speculation is supported by more recent inelastic neutron - scattering experiment in the doped ( bak)fe@xmath23as@xmath23 compound ( @xmath3=38k ) by christianson et al .
@xcite in this experiment , a clear magnetic - resonance peak is observed at the expected position of the afm correlation , i.e. , at ( @xmath24 ) momentum . with this series of experiments
it is clear that the afm correlation is the dominant magnetic correlation in the fe - pnictides among other competing magnetic correlations such as weak ferromagnetism @xcite , checkerboard afm @xcite , and afm stripe phase @xcite , which were theoretically proposed .
finally , the coupling matrix element is assumed to be a constant for simplicity .
the hamiltonian is written as    @xmath25    where @xmath26 and @xmath27 are the electron creation operators on the hole and the electron bands , respectively .
@xmath28 are the dispersions of the hole band and electron bands , respectively , defined as @xmath29 and @xmath30 . in this paper , we choose the band parameters as ( 0.30,0.24,-0.6 ) for hole band and ( 1.14,0.74,1.70 ) for electron band with the notation ( @xmath31 ) @xcite .
the pairing interaction @xmath32 is phenomenologically defined below .
it is all repulsive in momentum space and it represents a short range afm spin fluctuations as explained above .    @xmath33    where @xmath34 and @xmath35 are the two- dimensional momenta on the two dimensional bz and the parameter
@xmath36 controls the magnetic correlation length as @xmath37 ( @xmath38 is the unit - cell distance ) .
this interaction mediates the strongest repulsion when two momenta @xmath34 and @xmath35 are spanned by the ordering wave vector @xmath39 .
this condition is better fulfilled when the two momenta @xmath34 and @xmath35 reside each other on different bands in the model band structure ( see fig.1 ) . as a result ,
the sign - changing s - wave gap can form on each band as already suggested by several papers @xcite . however , this opposite sign gap on the hole and electron bands is not limited with the @xmath0s - wave state ( [ fig.1 ( a ) ] .
another possibility , which conforms to the lattice symmetry , is that each band develops a d - wave gap but with @xmath40 phase shift between two bands [ fig.1 ( b ) ] .
we call this type of gap as double d - wave gap .
we need to mention that our model did not include the screened coulomb interaction ( neither did the other theoretical investigations @xcite ) , which certainly exists in the fe - pnitide superconductors as well as in all metals in general .
the screened coulomb interaction is traditionally treated as `` coulomb pseudo potential '' @xmath41 in the conventional phonon- driven sc . however , the reliable estimate for its strength is practically impossible because even a small difference in @xmath41 would cause a large change in @xmath3 . in the fe - pnictides ,
if we are to determine @xmath3 theoretically , a quantitative estimate of @xmath41 is necessary .
we did not include it in our model interaction , first , because we do not know how to reliably estimate it in these compounds and , second , because the primary purpose of the present paper is not the prediction of the precise @xmath3 .
still we could investigate its generic effects on the different pairing symmetries such as @xmath0 s - wave and double d - wave gaps ; for example , how large value of @xmath41 is necessary to kill the @xmath0 s - wave pairing for a given strength of the afm interaction .
we think that this kind of analysis will dilute the focus of the present paper and therefore should be a separate investigation .
we briefly remark , however , on the general effects of the screened coulomb interaction .
the screened coulomb interaction becomes a short - range interaction in real space and therefore weakly momentum dependent in momentum space .
this type of interaction is almost harmless for the d - wave type pairing but extremely detrimental for the s - wave type pairing .
the pairing solutions in this paper should be considered with this point in mind .
now we solve the hamiltonian eq.(1 ) using the bcs approximation and the two band electrons need two sc order parameters ( ops )    @xmath42    after decoupling the interaction terms of eq.(1 ) using the above ops , the self - consistent mean field conditions lead to the following two coupled gap equations .
@xmath43 , \\
\nonumber % \delta_e ( k )   & = &   \\
\nonumber   - & \sum_{k^ { ' } } & [ v_{eh } ( k , k^ { ' } ) \delta_h ( k^ { ' } )   \chi_h ( k^ { ' } ) + v_{ee } ( k , k^{'})\delta_e ( k^ { ' } )   \chi_e ( k^ { ' } ) ] .
\\\end{aligned}\ ] ]    where @xmath44 , @xmath45 , etc are the same interaction defined in eq.(2 ) but the subscripts are written to clarify the meaning of @xmath44 = @xmath46 , @xmath45 = @xmath47 , etc . , and @xmath48 and @xmath49 specify the momentum @xmath50 located on the hole and electron bands , respectively . the pair susceptibilities are defined as @xmath51    where @xmath52 , and @xmath53 are the quasiparticle excitations and the dos of the hole and electron bands , respectively , and @xmath54 is the cutoff energy of the pairing potential @xmath55 .    when we solve the above gap equations eq(5 ) and eq.(6 ) , we numerically restricted the momenta @xmath56 and @xmath57 around the fss of the hole and electron bands within @xmath54 energy range . therefore , the fs shapes and the local dos @xmath58 of the realistic bands are faithfully taken into account in our gap solutions .
also no restriction on the functional forms of the gaps @xmath59 was imposed except the general symmetry depicted in fig.1 , so that the @xmath50-dependence of the gap functions @xmath60 will follow the characteristics of the bands and pairing interaction .
as explained in sec .
ii , the main pairing process with the afm spin - fluctuation mediated interaction @xmath55 is the interband pair hopping between the hole and the electron bands , in which a pair of electrons @xmath61 on the hole band scatters to a pair of electrons @xmath62 on the electron band and vice versa .
this process is particularly dominant when the size of the fs of each band is much smaller than the size of @xmath39 vector . considering only this interband pair process [ keeping only @xmath63 and @xmath64 terms in eqs.(5 ) and ( 6 ) ] , we observe the fact that the pair potential @xmath65 for the hole band electrons is provided by the pairs of electrons in the electron band and vice versa .
the physical consequence of it is that the relative sizes of the gaps and doss on each band are reversed ; namely , if @xmath66 , then @xmath67 holds in general .
this relation holds both for the @xmath0s - wave and for the double d - wave gap solutions and affects all superconducting properties such as tunneling dos , knight shift , @xmath2 , and penetration depth .    for all numerical calculations in this paper ,
we choose the parameters @xmath68 ( @xmath69 ) , @xmath70 , and @xmath71 ( average interaction @xmath72 ) .
our choice of band parameters produces @xmath73 and @xmath74 so that @xmath75 .
we think that these numbers represent the fe - based sc materials but should not be taken too seriously ; in particular the pairing strength @xmath71 is chosen freely for demonstration .
this solution for the fe - based sc is already proposed by several authors @xcite here we demonstrate that this solution is indeed realized with a simple phenomenological interaction , which mimics an afm spin fluctuations , on the minimal two band model representing the fe - based sc compounds .    as we described above , the reversed relation between the magnitude of the doss and the size of gaps holds more rigorously for the s - wave case and we suggest an approximate relation @xmath76 ( see the appendix for more detailed discussions ) .
this relation is a generic feature of the model .
therefore , given a substantial difference of dos between the hole and electron bands ( several band calculations @xcite indicate that this is true for the fe - based sc materials ) , at least two distinctively different sizes of the sc gaps should be observed in various experiments @xcite . in particular , because the band with a larger dos would dominate the physical properties but actually holds a smaller gap , this feature will modify various sc properties of the fe - based sc in unorthodox manner , such as @xmath77 value , temperature dependencies of various sc properties below @xmath3 , and the responses to impurities .
s - wave gap solutions @xmath78 and @xmath79 .
( b ) normalized dos of the hole band n@xmath80 ( red dotted line ) , electron band n@xmath81 ( blue dotted line ) , and the total n@xmath82 ( solid black squares ) .
[ fig2],width=302 ]    in fig.2 , the solution of the @xmath0s - wave gap and the corresponding dos are shown . as mentioned , @xmath83 for our bands .
accordingly the size of gaps of the hole band and the electron band is reversed as @xmath84 and @xmath85 .
the strongly momentum dependent pairing interaction and the realistic bands naturally induce an anisotropic modulation of the s - wave gaps with @xmath86 symmetry ; the modulation is stronger for the larger gap on the smaller dos band ( electron band around @xmath5 point ) .
compared to the case of the double d - wave solution , the average size of the @xmath0s - wave gap is larger by a factor of @xmath87 with the same pairing potential . therefore
, unless some other interactions are added , the ground state of our model is the @xmath0s - wave sc state .
this conclusion is already obtained by other authors @xcite with different models and approaches .
the separate and total doss plotted in fig.2(b ) show the main features of the @xmath0s - wave gap : two peak structure , the large dos with a small gap and the small dos with a large gap .
the overall shape of the total dos is not very much revealing compared to the current tunneling dos measurements @xcite .
however , it is too early to make a decisive conclusion with our calculations without including andreev scattering .
also the zbcp , the hallmark of a d - wave gap and observed in experiments with the fe - based superconductors @xcite , can equally be obtained with the @xmath0s - wave gap state .
we consider nuclear spin - lattice relaxation rate @xmath2 for the @xmath0s - wave gap .
several groups @xcite have reported that @xmath2 shows no coherence peak and the @xmath88 power law below @xmath3 , strongly suggesting an unconventional gap with lines of node such as a d - wave gap .
s - wave gap is known to have a constructive coherent factors for @xmath2 to induce the coherence peak over a temperature range below @xmath3 .
however , as mazin et al.@xcite envisaged , the sign - changing gaps between two bands provide a destructive coherent factor for the interband scattering which will largely cancel the intraband coherent factors . as a result the coherent peak of @xmath2 for the @xmath0s - wave gap sc will be substantially reduced .
the explicit formula that we used for the calculations is the following :    ( t ) of the @xmath0s - wave gap with @xmath89 = 1.5 .
( a ) separate term contributions of eq.(8 ) : total ( solid black square ) , hole band ( open red square ) , electron band ( open blue triangle ) , and interband term ( solid green square ) .
( b ) total @xmath2(t ) without ( solid black square ) and with ( open green pentagon ) damping .
( c ) the log - log plot of ( b ) .
the inset is a wide view [ fig3],title="fig:",width=340 ] ( t ) of the @xmath0s - wave gap with @xmath89 = 1.5 .
( a ) separate term contributions of eq.(8 ) : total ( solid black square ) , hole band ( open red square ) , electron band ( open blue triangle ) , and interband term ( solid green square ) .
( b ) total @xmath2(t ) without ( solid black square ) and with ( open green pentagon ) damping .
( c ) the log - log plot of ( b ) .
the inset is a wide view [ fig3],title="fig:",width=302 ]    @xmath90 \nonumber \\ &
+ & \bigl [   n_h ^2(0)\left\langle re     \frac{\delta_h(k)}{\sqrt{\omega^2-\delta_h ^2(k ) } }    \right\rangle_{k}^2    \nonumber \\    & + & 2 n_h ( 0 ) n_e ( 0 ) \left\langle re       \frac{\delta_h(k)}{\sqrt{\omega^2-\delta_h ^2(k ) } }      \right\rangle_{k } \left\langle re       \frac{\delta_e(k^{'})}{\sqrt{\omega^2-\delta_e ^2(k^ { ' } ) } }      \right\rangle_{k^ { ' } }   \nonumber \\    & + &   n_e ^2(0 ) \left\langle re       \frac{\delta_e(k)}{\sqrt{\omega^2-\delta_e ^2(k ) } }      \right\rangle_{k}^2 \bigr]\biggr \rbrace.\end{aligned}\ ] ]    for the temperature dependence of the gaps @xmath91 , we use a phenomenological formula , @xmath92 . by choosing the values of @xmath93
, we can partially take into account of the strong coupling superconductivity effect .
@xmath94 is not a sensitive parameter for final results ; we take @xmath95 in this paper .
fig.3(a ) shows the contributions to the @xmath2 relaxation rate from each terms of eq.(8 ) : the hole band , the electron band , and the interband terms .
it shows that the cancellation of the coherence factors is not perfect in general unless the conditions @xmath96 as well as @xmath97 are fulfilled .
nevertheless , due to the large cancellation by the interband coherence factor , the height of the coherence peak is very much reduced [ compare the total @xmath2 and the hole band only @xmath2 in fig.3(a ) ] .
small amount of impurities can easily wash out this reduced coherence peak as shown in fig.3(b ) ; the damping rate @xmath98 is enough to completely kill the coherence peak .
the subtle part is to fit the @xmath7 power law below @xmath3 .
it requires to tune @xmath99 ratio . in fig.3 , @xmath100 ( automatically , it makes @xmath101 which is quite a large value ) is used for the best fit .
fig.3(c ) shows that this pseudo-@xmath88 behavior is not extended to the very low - temperature region as in the d - wave case because this @xmath88 behavior in the @xmath0s - wave gap is not an intrinsic property of the lines of nodes . at low temperatures
, there appears the exponential drop inevitably due to the full gaps , and then it finally reaches the impurity - induced @xmath9-linear region because we added some amount of impurities to kill the coherent peak .
all these details put rather stringent conditions to confirm the @xmath0s - wave gap state with experiments .    s - wave gap .
the total ( solid black square ) , hole band ( open red circle ) , and electron band ( open blue triangle ) contributions are shown separately .
( a ) @xmath102 and ( b)@xmath89 = 0.5 .
[ fig4],width=302 ]     of @xmath0s - wave gap and its separate contributions from the hole and electron bands .
( a ) @xmath102 and ( b)@xmath89 = 0.5 .
[ fig5],width=302 ]    now we consider the knight shift which is the measure of uniform susceptibility in sc phase .
because it is a @xmath103 probe , there is no interband contribution and the total knight shift is just sum of the contributions from each band as follows :    @xmath104.\end{aligned}\ ] ]    in fig .
4 , the normalized knight shift ( uniform spin susceptibility ) is plotted and it shows the typical flat behavior of a s - wave gap at low temperatures .
the contributions from the hole and electron bands show separately the feature of the larger dos with small gap and the smaller dos with a larger gap .
4(a ) is the results with @xmath89 = 1.5 , the best gap-@xmath3 ratio to fit @xmath88 behavior of @xmath2 in fig.3 .
fig4(b ) showed the results with @xmath89 = 0.5 to demonstrate the convex shape of knight shift which was reported by matano et al.@xcite for pr(fo)feas .
impurities does not change much of this feature unlike in the case of d - wave gap .    finally , we consider the penetration depth .
the static response function to the electromagnetic fields is the following : @xcite    @xmath105    the @xmath106 limit of this kernel @xmath107 is directly proportional to the superfluidity density or 1/@xmath108 in the london limit . for our two band model ,
total kernel is the sum of @xmath109 and @xmath110 with the proper weighting factor proportional to the dos n@xmath111 of each band and there is no interband screening current .
@xmath112 is the non - local parameter and can be rewritten in more convenient form as @xmath113 .
@xmath114 and @xmath115 are the coherence length and the penetration depth at zero temperature , respectively .
apparently if @xmath116 is small compared to 1 , the non - local effect becomes negligible .
a typical value of @xmath117 for ybco was estimated about 0.01 , for example @xcite . for the fe - based superconductors , we believe that @xmath117 is not much larger than the values of the high-@xmath3 cuprates . also for a s - wave case , the non - local effect does not change much of the temperature dependence of @xmath118 except the overall magnitude . therefore , we take @xmath119 for the calculations of the penetration depth in the @xmath0s - wave case . however , this effect can induce an important modifications in the d - wave case , which will be discussed in sec.iiib .    fig.5 shows the normalized superfluidity density @xmath120 and separate contributions from the hole and electron bands .
the exponentially flat region appears at low temperatures due to the full gap opening , which is consistent with recent experiments.@xcite relatively narrower region of the flat part ( for t @xmath121 ) compared to the ordinary s - wave gap is another feature due to the smaller gap with the larger dos of the @xmath0s - wave gap sc .
a subtle part here is to fit the high- temperature region ( 0.3 t@xmath122 ) . with @xmath123 = 1.5 ( the same value used for the @xmath2 fit ) , this part becomes too convex [ fig.5(a ) ] in comparison to the experiments .
@xcite a smaller gap-@xmath3 ratio can make it concave as shown in fig.5(b ) ( with @xmath89 = 0.5 ) ; this concave feature was recently observed by martin et al .
@xcite     and @xmath79 .
( b ) normalized dos of the hole band n@xmath80 ( red dotted line ) , electron band n@xmath81 ( blue dotted line ) , and the total n@xmath82 ( black squares ) .
[ fig6],width=302 ]    in summary , the @xmath0s - wave gap state provides the most consistent descriptions for the penetration depth experiments .
however , it explains @xmath2 only for a limited temperature range , even with a fine tuning of @xmath124 ratio and impurities .
knight shift of any shape can be fit with two band parameter ( this is also true with the double d - wave gap ) .
the tunneling dos does not provide a decisive conclusion .
( t ) of the double d - wave gap .
( a ) each term contributions of eq.(8 ) : total ( solid black square ) , hole band ( open red square ) , electron band ( open blue triangle ) , and interband term ( open green triangle ) .
( b)log - log plot of total @xmath2(t ) for @xmath125 ( black square ) and 1.0 ( green circle ) .
[ fig7],width=302 ]      in fig.6 , the gap solution and the corresponding dos of the double d - wave gap are shown .
as mentioned , our model bands have @xmath83 , and consequently gap in the hole band @xmath126 is smaller than the one of the electron band @xmath127 .
the sizes of the maximum gaps are @xmath22 5 times smaller than the @xmath0s - wave gap solutions .
therefore , the double d - wave gap solution is not the best sc state for our phenomenological model with an antiferromagnetic pairing interaction .
this result is in agreement with other theoretical studies @xcite .
this conclusion may change with the correlation length of the afm fluctuations , the sizes of the fs of the hole and electron bands , etc .
but we numerically found that @xmath0s - wave gap solution is favored compared to the double d - wave gap solution for most cases .
as discussed in sec.ii , however , the screened coulomb interaction may change this conclusion .
the detailed studies about this issue will be discussed in a separate paper .    to complete the comparisons
, we calculated the same sc properties of the double d - wave gap state .
the separate and total doss shown in fig.6(b ) demonstrates the main features of the double d - wave gap : the large dos band with a small gap and the small dos band with a large gap .
this result shows a similar feature of the tunneling dos measurement by wang et al .
@xcite except the zbcp , which does not show up in our simple dos calculation but should appear when the tunneling conductivity is properly calculated with andreev scattering process .
we consider the nuclear spin - lattice relaxation rate @xmath2 . as in the case of @xmath0s - wave gap
, there are three contributions for the total @xmath2 relaxation rate : hole band , electron band , and interband scattering terms .
the formula is the same as eq.(8 ) but the last three terms should be dropped because the fs average of @xmath59 vanishes in this case .
we use the same form of temperature dependent gap function as @xmath128 . in fig.7
, @xmath129=1.75 is used for the best @xmath88 fit below @xmath3 . however , in the double d - wave gap , @xmath130 provide reasonably good fits , showing a more tolerance than the @xmath0s - wave gap state .     and
( b ) @xmath89 = 1.0 [ fig8],width=302 ]    fig.7(a ) shows the separate contributions from each channel together with the total contribution . as in the @xmath0s - wave case
, @xmath131 provides the largest contribution and @xmath132 provides the smallest contribution .
there is also the interband term @xmath133 .
in contrast to the @xmath0s - wave case , all three terms display a similar temperature dependence and no coherence peaks .
fig.7(b ) shows the same @xmath134 ( black squares ) in log - log plot .
the overall features of @xmath2 are the ones of the typical d - wave sc state : no coherence peak near @xmath3 and @xmath135 below @xmath3 , and consistent with the current nmr experiments @xcite . at very low temperatures ,
t - linear behavior starts to appear due to a small damping for the numerical calculations ( @xmath136 ) . for comparison
, we also show @xmath134 ( open green circles ) with @xmath137 = 1.0 , a smaller gap-@xmath3 ratio ; it exhibits a substantial convex part below @xmath3 and then starts displaying the universal @xmath88 behavior before entering the impurity dominating region .    in fig.8
, we show the result of the uniform spin susceptibility which is measured as knight shift .
fig.8(a ) shows the results with @xmath125 .
the hole band contribution is dominant as in @xmath2 and the electron band contribution show the steeper drop just below @xmath3 because of the larger gap-@xmath3 ratio @xmath138 .
the overall behavior of the total @xmath139 below @xmath3 shows a typical d - wave behavior such as @xmath9-linear at low temperatures .
fig.8(b ) shows the results with @xmath140 .
a smaller gap-@xmath3 ratio makes the knight shift convex as in the @xmath0s - wave case and observed by matano et al.@xcite for pr(fo)feas .
this result demonstrates that this convex behavior of knight shift is irrelevant to the gap symmetry but a generic feature of the two - gap ( or multigap ) sc .
but it reveals that the gap-@xmath3 ratio @xmath141 is much smaller than the standard bcs value , where @xmath142 refers to the gap of the band with largest dos .
now we calculate the penetration depth . as we discussed in sec.ii , most of experiments , up to now , report a flat temperature dependence of @xmath143 at low - temperature region and suggest a fully gapped sc state @xcite .
a naive double d - wave gap state has no chance to explain this flat behavior at low temperatures .
therefore , we consider a non - local effect of the electromagnetic response of the double d - wave gap superconductor as a possible cause to modify the typical temperature dependence . in order to include the effect of the non - local electrodynamics
, we use the fully @xmath144dependent kernel @xmath145 [ eq.(10 ) ] and put it into the integral formula for @xmath146 with the specular boundary condition ,    @xmath147    where @xmath148 are the normalized kernels as @xmath149 and @xmath150 is a dimensionless momentum .
the results with a diffusive boundary condition are qualitatively the same ; therefore , they will not be discussed . for the non - local parameter @xmath117
, we think it to be much smaller than 1 for the fe - based superconductors , but here we take it as a free parameter and see how large value is needed to fit experimental data .
fig.9 and fig.10 show the normalized total superfluidity density @xmath151 and separate contributions from the hole and electron bands for the double d - wave gap state .
fig.9 used @xmath152=1.75 and fig.10 used @xmath153=1.0 . in each figure , panel ( a ) is a local limit ( @xmath154 ) and the panel ( b ) is a non - local limit ( @xmath155 ) .
the local cases display the typical d - wave behavior at low temperatures , i.e. , the linear decrease in @xmath9 from @xmath156 .
the extreme non - local cases ( @xmath157 ) introduce a substantial round - off ( @xmath158 ) region at low temperatures , which is , however , not an exponentially flat behavior as the recent experiments claim .
further , even a rough fitting requires an unreasonably large non - local parameter @xmath117 .     of double d - wave gap and its separate contributions from the hole and electron bands with @xmath159 .
( a ) @xmath160 and ( b ) @xmath117 = 0.5 .
[ fig9],width=302 ]     of double d - wave gap and their separate contributions from the hole and electron bands with @xmath161 .
( a ) @xmath160 ; ( b ) @xmath117 = 0.5 .
[ fig10],width=302 ]    in summary , the double d - wave gap state can provide consistent descriptions for tunneling dos , @xmath2 , and knight shift .
however , there is an intrinsic difficulty to explain the flat behavior of the penetration depth at low temperatures .
also , in our model with an afm mediated pairing interaction only , the double d - wave gap solution is energetically less favored than the @xmath0s - wave gap solution .
we demonstrated that a minimal model with a phenomenological pairing interaction of the afm spin fluctuations can allow both the @xmath0s - wave gap and the double d - wave gap solutions with the realistic bands of the fe - based sc compounds . with the same parameters ,
the @xmath0s - wave gap solution is energetically more favorable by a factor of @xmath225 times , so that it has a better chance to be realized in the fe - based sc compounds .    in both cases
, we found that the approximate relation @xmath76 holds because it is a generic feature of the two gap sc when an interband pair scattering is the dominant pairing interaction .
this relation appears for all sc properties in subtle way , which modifies the value of @xmath162 and other sc properties in unorthodox way . numerically solving the coupled gap equations for the two bands
, we found the detailed structure of the gap functions @xmath163 , which showed an anisotropy ( @xmath22 20 @xmath164 ) of the @xmath0s - wave gaps .
we also calculated the key sc properties , for both gap states , such as tunneling dos , @xmath2 , knight shift , and penetration depth and discussed them in comparison with experiments .
when we calculated these quantities , we paid a special attention to the interband coherence factor which is a unique feature of multigap sc .
this interband coherence factor particularly produced an important modification to the @xmath2 relaxation rate of the @xmath0s - wave gap state .
the @xmath0s - wave gap state provides the most consistent descriptions for the penetration depth experiments : the flat low- temperature behavior.@xcite besides the low - temperature behavior , the high - temperature ( 0.3 @xmath165 ) behavior  due to a large difference of the gap sizes @xmath166 and @xmath167 and their corresponding dos @xmath168  can be either concave or convex .
however , @xmath2 experiments can only be fitted for a limited temperature range , even with a fine tuning of @xmath124 ratio and impurities .
the @xmath0s - wave gap state is not inconsistent with the knight shift and the tunneling dos data , but overall does not provide any decisive merit in comparison with the double d - wave gap .    the double d - wave gap state , although it is energetically less favorable in our model unless additional interactions are added , provides the best fit to the @xmath2 experiments .
however , it has a difficulty to explain the penetration depth experiments for the low temperature flat behavior .
it requires an unreasonably large non - local effect to fit the low temperature part ; it is still not exponentially flat but only @xmath158 .
if this low temperature part of @xmath146 is , indeed , confirmed to be exponentially flat , the double d - wave gap state should be ruled out .
tunneling dos and knight shift can be fit with the double d - wave gap state as much as with the @xmath0s - wave gap state
.    in conclusion , quantitative calculations , carried out in this paper , with the two most promising sc gap states can serve as guidelines for sorting out the possible pairing states of the fe - based sc in comparison with the current and future experiments .
for that , very low temperature measurements and systematic studies with the amount of impurities will provide decisive information to determine the correct gap symmetry .
_ note added  _ recently , we have known that similar studies of @xmath2 for the @xmath0s - wave state were carried out by two groups @xcite where only the interband scattering process was analyzed and by another recent paper @xcite where both the interband and intraband processes were considered as in our paper .
this work was supported by the kosef through the grants no .
krf-2007 - 070-c00044 , and no .
krf-2007 - 521-c00081 .
the bottom baseline is @xmath171 and the top baseline is @xmath170 .
[ fig1],width=340 ]    in the main text , we claimed the approximate relation @xmath172 as a generic feature of the two- band model with a dominant interband interaction .
this kind of relation will have direct and important implications to the experimental observations .
however it is pointed out by mazin @xcite that this is not a rigorous identity in general . in this appendix
, we clarify the degree of the validity of this relation .
in general , these equations do not yield a simple algebraic relation between @xmath184 and @xmath170 , but we can obtain the simple relations for the limiting cases .
first , for the extreme weak coupling limit , ie .
, when @xmath185 , the two logarithmic terms become asymptotically equal as @xmath186 , and we obtain the same relation as the @xmath178 case eq.(a.5 ) . on the other hand , for strong - coupling limit ,
, when @xmath187 ( which is certainly an unphysical limit ) , @xmath188/\delta_{h , e } ) \approx \omega_{afm}/\delta_{h , e}$ ] and we obtain the relation ,      having found the results of the two limiting cases , we can guess that the gap ratio @xmath184 should be in between these two limiting ratios .
for example , we can attempt an expansion with @xmath190}$ ] starting from the weak coupling limit @xcite , and we obtain , in the first order of @xmath191 ,      where @xmath193 is a dimensionless coupling constant . for practical use , we numerically solve eqs .
( a.1 ) and eq.(a.2 ) and plot the ratio @xmath194 as a function of @xmath195 in fig.a.1 .
the ratio @xmath184 becomes a universal curve when it is normalized by the distance between @xmath196 and @xmath197 .
the result indeed shows that when @xmath198 , it is in between two limiting ratios , @xmath196 and @xmath197 as we expected from the above analysis .    in reality ,
there exist two complications .
first , the intraband couplings @xmath199 and @xmath200 need to be included . a little analysis of eqs .
( 5 ) and ( 6 ) as well as of numerical results reveals that this effect always enhances the gap ratio toward the limit @xmath197 .
another complication arises from the fact that there are more than two bands in real materials @xcite . applying the results of the above analysis
, we can suggest the following approximate relations .
first , we classify the bands of the real materials into two groups : the hole bands around @xmath4 point and the electron bands around @xmath5 point , respectively . then in the strong - coupling limit ,              g. f. chen , z. li , d. wu , g. li , w.z .
hu , j. dong , p. zheng , j.l .
luo , n.l .
wang , phys .
lett . * 100 * , 247002 ( 2008 ) ; x. h. chen , t. wu , g. wu , r. liu , h. chen , and d. fang , nature ( london ) * 453 * , 761 ( 2008 ) ; z. ren , et al . , chin .
lett . , * 25 * , 2215 ( 2008 ) .
e. manousakis , jun ren , and e. kaxiras , arxiv:0806.3432 ( unpublished ) ; 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 , arxiv:0807.2009 ( unpublished ) .
k. matano , z.a .
ren , x.l .
dong , l.l .
sun , z.x .
zhao , and guo - qing zheng , europhys . lett . * 83 * 57001 ( 2008 ) ; h .- j .
grafe , d. paar , g. lang , n. j. curro , g. behr , j. werner , j. hamann - borrero , c. hess , n. leps , r. klingeler , and b. buechner , phys .
* 101 * , 047003 ( 2008 ) .
t. y. chen , z. tesanovic , r. h. liu , x. h. chen , and c. l. chien , nature ( london ) , * 453 * , 1224 ( 2008 ) ; y. wang , l. shan , l. fang , p. cheng , c. ren , and h. wen , arxiv:0806.1986 ( unpublished ) ; l. shan , y. wang , x. zhu , g. mu , l. fang , c. ren , and h. wen , europhys .
letters , * 83 * , 57004 ( 2008 ) .
h. mukuda , n. terasaki , h. kinouchi , m. yashima , y. kitaoka , s. suzuki , s. miyasaka , s. tajima , k. miyazawa , p.m. shirage , h. kito , h. eisaki , and a. iyo , arxiv:0806.3238 , j. phys .
( to be published ) ; y. nakai , k. ishida , y. kamihara , m. hirano , and h. hosono , arxiv:0804.4765 , j. phys .
( to be published ) .
l. malone , j.d .
fletcher , a. serafin , a. carrington , n.d .
zhigadlo , z. bukowski , s. katrych , and j. karpinski , arxiv:0806.3908 ( unpublished ) ; k. hashimoto , t. shibauchi , t. kato , k. ikada , r. okazaki , h. shishido , m. ishikado , h. kito , a. iyo , h. eisaki , s. shamoto , and y. matsuda , arxiv:0806.3149 ( unpublished ) ; c. martin , r. t. gordon , m. a. tanatar , m. d. vannette , m. e. tillman , e. d. mun , p. c. canfield , v. g. kogan , g. d. samolyuk , j. schmalian , and r. prozorov , arxiv:0807.0876 ( unpublished ) .      c. de la cruz et al . ,
nature ( london ) * 453 * , 899 ( 2008 ) ; j. zhao , q. huang , c. de la cruz , s. li , j. w. lynn , y. chen , m. a. green , g. f. chen , g. li , z. li , j. l. luo , n. l. wang , and p. dai , arxiv:0806.2528 ( unpublished ) .
q. han , y. chen , and z. d. wang , europhys .
lett . * 82 * , 37007 ( 2008 ) ; tao li , arxiv:0804.0536 .
( unpublished ) ; v. cvetkovic , and z. tesanovic , arxiv:0804.4678 ( unpublished ) ; g. baskaran , arxiv:0804.1341 ( unpublished ) ; p. a. lee , and xiao - gang wen , arxiv:0804.1739 ( unpublished ) ; x. dai , z. fang , y. zhou , and fu - chun zhang , arxiv:0803.3982 ( unpublished ) .


supersymmetry  ( susy ) is a well - motivated extension of the standard model  ( sm ) that provides an explanation for the stability of the hierarchy between the weak scale and the planck scale  @xcite .
however , experiment excludes susy from being an exact symmetry at low energies .
if supersymmetry is _ softly _ broken , it is possible to push the superpartner masses up enough to make the theory consistent with experiment while still preserving the electroweak - gravity hierarchy . for this to work
, the masses of the sm superpartners should not be significantly larger than the electroweak scale , putting them within reach of upcoming particle collider experiments such as the lhc .
a major drawback of this scenario of softly - broken low - energy supersymmetry is that the supersymmetry breaking operators generically introduce many new sources of quark and lepton flavor mixing , leading to flavor - changing neutral currents  ( fcnc ) in conflict with experiment  @xcite .
one way to cure this problem is to take no - scale boundary conditions in which the soft scalar masses and couplings all vanish at a boundary input scale , with the gaugino soft masses being the only significant source of supersymmetry breaking in the visible sector at this scale  @xcite .
such boundary conditions can arise from gaugino mediation  @xcite , or conformal running effects  @xcite .
low - scale scalar soft terms , which are necessary to obtain a viable phenomenology , are generated primarily from the gaugino masses in the course of renormalization group  ( rg ) running from the boundary input scale down to near the electroweak scale .
this dominant contribution is controlled by gauge interactions , implying that the regenerated scalar soft terms , and thus the squark and slepton masses , will be approximately flavor - universal at the low scale .
a slight generalization of this pure no - scale picture , which still leads to scalar soft terms that are roughly flavor universal at the low scale ( for the first and second generations ) , is to allow non - vanishing higgs soft scalar masses at the boundary scale . in the context of gaugino mediation
this can arise if the higgs multiplets propagate in the bulk  @xcite , while for conformal running , it can emerge from hidden - sector interactions related to the origin of the @xmath7 term  @xcite .
these _ higgs exempt no - scale _  ( hens ) models often have an important phenomenological advantage over pure no - scale constructions : in large regions of the parameter space hens models have a neutralino lightest superpartner particle  ( lsp )  @xcite .
this differs from the situation in strict no - scale models where typically the predominantly right - handed stau is the lightest sm superpartner particle  @xcite .
while such models can be viable if the true lsp is a gravitino  @xcite , a neutralino lsp with exactly conserved @xmath8-parity tends to make for a better thermal dark matter candidate .
hens models can push the right - handed stau mass above that of the lightest neutralino by the influence of the non - vanishing higgs soft masses on the rg running of the other soft parameters .
the primary contribution to the slepton soft terms comes from the hypercharge @xmath9 term @xmath10 that appears in the rg equations for all the mssm scalar soft masses  @xcite .
this contribution increases the right - handed slepton soft masses in the course of running down provided @xmath11 is negative . in the portion of the hens model parameter space that is also consistent with electroweak symmetry breaking , the higgs scalar soft masses
must typically be very tachyonic at the boundary scale to obtain a neutralino lsp , and the uncolored sfermions still end up being relatively light  @xcite
.    while these hens models are compelling , the tachyonic higgs soft masses at the boundary scale lead to concerningly large tachyonic low - scale higgs soft masses .
these , combined with the small ( but positive ) slepton soft masses , suggest that the true vacuum of the theory might be a charge - and - color - breaking  ( ccb ) minimum , or there may exist unbounded - from - below  ( ufb ) directions that are only stabilized far out in field space by higher - dimensional operators  @xcite .
the presence of such features need not exclude these regions of the parameter space provided our sm electroweak vacuum is metastable and long - lived relative to the age of the universe .    in the present work we study the vacuum structure of hens models , and
compute the lifetime of the sm vacuum state when it turns out to be metastable .
we concentrate on the stability of the sm vacuum at zero temperature with respect to vacuum tunneling .
thermal effects in the early universe can also potentially induce thermal transitions between different vacua , and modify the shape of the potential itself .
however , these thermal effects tend to stabilize the origin of the field space to which the sm vacuum is connected , favoring this vacuum over others that lie further out in the field space  @xcite .
thermal effects are especially effective in delaying the formation of vacua that break color due to the large thermal corrections from the strong gauge and top quark yukawa couplings  @xcite .
scalar field evolution during and after inflation may also populate non - standard vacua , although the precise result depends on the details of inflation and lies beyond the scope of this paper  @xcite . by focusing solely on the @xmath12 constraints we obtain conservative and unambiguous bounds on the hens parameter space that do not depend on the cosmological history .
motivated by gauge coupling unification , throughout this work we take the input scale to be @xmath13 and we assume a universal gaugino mass @xmath14 at this scale .
the set of independent free parameters after imposing consistent electroweak symmetry breaking is therefore term is also free , but we fix it to be positive .
this slightly increases the model prediction for @xmath15 , as favored by the experimentally measured value of this quantity  @xcite . ]
m_1/2 ,  ,  m_h_u^2(m_c ) ,  m_h_d^2(m_c ) .
[ inputpars ] investigations of no - scale models with non - universal gaugino masses can be found in refs .
similar input soft terms have also been studied in the context of refs .
@xcite    the outline of this paper is as follows . in section  [ hensvac ]
we investigate the vacuum structure of hens models .
we discuss the general features of vacuum tunneling and outline the relevant criteria we use to judge which of the possible ccb and ufb vacua are the most dangerous in section  [ bounce ] . in section  [ constraint ]
we map out the portions of parameter space that are phenomenologically consistent and give rise to a sm vacuum state that is sufficiently long - lived .
section  [ conc ] is reserved for our conclusions . in appendix  [ appa ]
we outline our numerical technique for estimating the vacuum tunneling rate .
we present in appendix  [ appb ] an exact tunneling solution for an inverted - parabola barrier .
with the hens model input parameters of eq .  , relatively small slepton soft masses as well as large tachyonic higgs soft masses obtain near the electroweak scale in much of the phenomenologically allowed parameter space  @xcite .
this spectrum of soft parameters frequently implies the existence of ufb directions and ccb minima that are much deeper than the standard electroweak vacuum .
indeed , we find that nearly the entire allowed parameter space in the hens model ( before imposing vacuum stability constraints ) has at least one ufb direction .
we investigate the existence and nature of such potentially dangerous vacuum features in the present section .    given the soft breaking spectrum that arises in the hens scenario , the results of ref .
@xcite suggest that the most dangerous vacuum feature will be a sleptonic ufb-3 direction .
this direction has @xmath16 , with @xmath17 , @xmath18 , @xmath19 , and @xmath20 all non - zero . turning on expectation values for these fields
, there exists a @xmath9- and @xmath21-flat direction that is only lifted by quadratic supersymmetry breaking operators . to obtain @xmath21-flatness ,
only the @xmath21-term of @xmath22 must be cancelled .
this can be arranged by taking || = |_l|=|_r| = , [ stauufb ] with the relative phase of @xmath17 and @xmath18 chosen appropriately .
@xmath9-flatness is then obtained by setting    @xmath21-flat field configuration provided |h_u^0| > -|| .
[ hulow ] the scalar potential along this direction in field space then becomes  @xcite v_ufb-3 = ( m_h_u^2+m_l_i^2)|h_u^0|^2 + ||(m_l_3 ^ 2+m_e_3 ^ 2+m_l_i^2)|h_u^0| - .
[ vufb3 ] when @xmath23 is negative , the potential becomes unbounded in the limit @xmath24 .
it will ultimately be stabilized by loop corrections or higher - dimensional operators ( that we have not included in eq .  ) at a location that is very deep and far out in field space relative to the electroweak vacuum .
this sleptonic ufb-3 direction is particularly dangerous in the hens models on account of the large and negative values of @xmath25 and the smaller values of @xmath26 and @xmath27 that emerge in the low - energy spectrum .
these properties imply that the barrier against tunneling from the electroweak vacuum near the origin out to the deeper ufb-3 direction , arising from the linear term in @xmath28 in eq .  , will not be especially large .
the barrier will be further weakened by larger values of @xmath29 which enhance the coupling @xmath30 .
other similar ufb-3 directions may be present in the theory , but they will generally have larger barriers due to the larger values of the squark soft masses or the smaller values of the first- and second - generation lepton yukawa couplings .
when @xmath28 does not satisfy the bound given in eq .
, the lowest - energy @xmath21-flat direction in the potential has @xmath31 , and is given by  @xcite @xmath32 this potential is no longer @xmath9-flat , and is stabilized at a finite value of @xmath28 . if this point occurs with @xmath28 less than the bound of eq .  , it is a constrained local ccb minimum . on the other hand ,
when this constrained local extremum has @xmath28 larger than the bound of eq .
, it represents a saddle point that is unstable under flowing to a non - zero value of @xmath33 .    in practice
, we find that for smaller values of @xmath28 the potential can be reduced further by relaxing the @xmath21-flatness constraint of eq .  .
the effect of dropping the @xmath21-flatness constraint is that the minimal potential for a given value of @xmath28 is deformed slightly away from the precise ufb-3 form of eqs .
( [ vufb3],[vufb3small ] ) , but that the same qualitative features remain .
in particular , there usually remains a local extremum at @xmath34 and @xmath35 .
we shall designate this local extremum as ccb-4 .
if this extremum occurs at smaller values of @xmath28 , on the order of the electroweak scale , it can be a local ccb minimum . when the ccb-4 extremum occurs with a value of @xmath28 much larger than the electroweak scale , it is generally a saddle point that flows in the @xmath36 direction to a genuine ufb-3 direction of the form given in eq .  .
even when it is only a saddle , the ccb-4 point plays an important role in determining the tunneling rate from the electroweak minimum to the ufb-3 direction , as we will discuss below .    a stable ccb-4 minimum with stops
can also arise for smaller values of @xmath29 and a correspondingly larger @xmath37 yukawa coupling . in this case
we have @xmath38 while @xmath39 and @xmath40 are all non - zero .
the relevant potential is @xmath41 this potential is generally stable against excursions in the @xmath42 direction on account of the larger squark soft masses that arise in the hens model . with @xmath43 at the low scale , this potential often has a ccb-4 minimum with both @xmath44 and @xmath45 non - zero .
however , as we show below , the barrier to tunneling to this minimum from the standard electroweak minimum to this ccb-4 minimum is usually safely large , again on account of the larger values of the squark soft masses as well as the less negative values of @xmath46 . for similar reasons
, we expect that the rate for tunneling to the other potential ccb minima discussed in ref .
@xcite will typically be less constraining than the rate to tunnel to a stau ufb-3 direction .
the existence of vacua deeper than the standard electroweak minimum in hens models implies there is a danger of tunneling into one of these phenomenologically unacceptable states . at the very least , the lifetime for this tunneling must be greater than the age of the universe .
the vacuum - to - vacuum transition rate associated with tunneling can be calculated using path integral methods  @xcite . in the semiclassical approximation ,
the lifetime of the vacuum is found to be @xmath47 } \label{liftim}\ ] ] where @xmath48 is a dimension - four prefactor to be discussed below , @xmath49 denotes the _ bounce solution _ for the @xmath50-@xmath51 field , and @xmath52 is the euclidean action , @xmath53=\int d^4x_e\left(\left|\nabla\bar{\phi}_i\right|^2 + u(\bar{\phi}_1, ... ,\bar{\phi}_i)\right)=t[\bar{\phi}_i]+v[\bar{\phi}_i ] .
\label{bouact}\ ] ] the _ bounce solution _ for the fields @xmath49 is the extremum of the euclidean action that is @xmath54-symmetric and obeys the the following boundary conditions : @xmath55 where @xmath56 is the euclidean distance .
these boundary conditions correspond to a field that begins in a metastable vacuum and tunnels through a barrier separating it from a deeper vacuum state , emerging with zero kinetic energy .
we focus on @xmath54 symmetric solutions because these are expected to have the least action , and therefore dominate the tunneling probability  @xcite .
even with the simplification of an @xmath57 symmetry , the equation of motion for the bounce can not in general be solved analytically .
we compute the bounce solution numerically using the _ improved action method _  @xcite .
the details of this method are outlined in appendix  [ appa ] .
even more difficult to compute is the non - exponential pre - factor @xmath48 in eq .
( [ liftim ] )  @xcite . on general grounds ,
we expect it to be on the order of the mass scale setting the size of the potential barrier . in low - energy supersymmetry , a good estimate for this number is @xmath58 , which we take to be the case throughout the rest of this article .
the precise value of this pre - factor is unlikely to affect our qualitative conclusions as the multiplicative uncertainty in its value is much more slowly - varying than the exponentiated large values of the action that lead to acceptable lifetimes .
our choice for the pre - factor is also conservative , in that choosing a larger number here would only exclude more points . with this pre - factor , it is found that the lifetime of the sm vacuum will be greater than the age of the universe , @xmath59 , provided @xmath60 > 400 $ ] .
the value of the bounce action for tunneling between a pair of local vacua depends on the relative depth of the minima , the number of distinct minima , the height of the barriers between them , and the relative size of the field values within them . in general
, the bounce solution represents a configuration of fields that simultaneously minimizes these opposing contributions , with the kinetic term favoring slowly varying fields , and the potential term preferring to reach the deepest minimum as quickly as possible .
some intuition for this balance can be obtained by examining special cases in which the bounce solution can be derived analytically .
one such example is a potential consisting of piecewise linear segments  @xcite .
motivated by the ufb-3 potential in eq .  , we present another exact tunneling solution in appendix  [ appb ] for an inverted parabolic potential .
as expected , our analytic tunneling solution shows that the bounce action increases with the size of the barrier .
this solution also implies that the depth of the minimum to which one is tunneling to ceases to matter once it becomes very deep .
thus , we can safely compute the rate to tunnel into a ufb-3 direction without knowing where it is ultimately stabilized provided the corresponding minimum is very deep relative to the height of the barrier .
on the other hand , the relative depth of the minima is relevant to the tunneling rate when they are nearly degenerate , as can be seen from the analytic _ thin - wall _ approximate solution that can be safely applied in this case  @xcite .
one important feature not captured by our one - dimensional analytic tunneling solution is that the bounce action tends to also increase when there are more independent fields involved in the tunneling process , as each one of them contributes non - negatively to the kinetic portion of the bounce action . for these reasons ,
numerical solutions are used to get the details correct .
the possibility of tunneling from the sm vacuum to a phenomenologically unacceptable vacuum places strong constraints on the parameter space of the hens model . to investigate these constraints ,
we have calculated the bounce action for tunneling to a ccb-4 minimum or a ufb-3 direction for a series of representative hens parameter sets .
our strategy is to fix @xmath14 and @xmath29 , and scan over the input values of @xmath61 and @xmath62 that lead to an acceptable low energy spectrum .
we focus on the values @xmath63 for @xmath64 and @xmath65 for @xmath66 . in each of these scans we have consider points that meet the phenomenological constraints laid out in ref .
@xcite , including the current lep lower bounds on the lightest higgs boson .
we also keep the points with a charged slepton as the lightest mssm superpartner found in ref .
@xcite , noting that they would require a more complicated cosmology and likely r - parity violation or a gravitino lsp  @xcite . by keeping such points ,
our results are also applicable to minimal gaugino mediation subject to our assumptions about the gaugino masses and the input ( compactification ) scale .     and
@xmath62 for @xmath67 and @xmath68 .
all points shown are consistent with collider phenomenology .
the points enclosed below by the dash - dot line have a neutralino lsp .
the solid line separates the region with a ccb-4 minimum or saddle point ( left ) from that which only has a ufb-3 ( right ) .
@xmath69 is cosmologically safe .
[ m12500b30],scaledwidth=70.0% ]    we first consider the cases with @xmath68 , where the tau yukawa coupling is enhanced .
this has the effect of opening up a ccb-4 minimum or ccb-4-like saddle point that flows to a ufb-3 direction . in fig .
[ m12500b30 ] we show ranges of the bounce action for tunneling out of the sm vacuum to a ccb-4 minimum or a ufb-3 direction for @xmath67 as a function of the input values of @xmath61 and @xmath62 .
all points to the left of the solid line in this figure have either a ccb-4 minimum or saddle point .
the points to the right have either a stable sm vacuum or a ufb-3 direction with no ccb-4 stationary point .
to the left and the right of the solid line we show the bounce action for tunneling along the direction leading to the ccb-4 stationary point or to the ufb-3 direction , whichever is smaller .
we emphasize the ccb-4 stationary point here , even when it is only a saddle point , because we find that the most dangerous lowest - action tunneling path from the electroweak minimum is typically one that passes near this point with @xmath70 .
the dot - dashed line in fig .
[ m12500b30 ] indicates the upper border of the portion of parameter space in which the lsp is the lightest neutralino .
( see ref .
@xcite for more details . ) . in fig .
[ m121000b30 ] we show the same quantities as in fig .  [ m12500b30 ] for @xmath71 .     and @xmath62 for @xmath72 , @xmath64 and @xmath73 .
all points shown are consistent with collider phenomenology .
the points enclosed below by the dash - dot line have a neutralino lsp .
the solid line separates the region with a ccb-4 minimum or saddle point ( left ) from that which only has a ufb-3 direction .
@xmath69 is cosmologically safe .
[ m121000b30 ] , scaledwidth=70.0% ]    both fig .
[ m12500b30 ] and fig .
[ m121000b30 ] show that for @xmath68 the lifetime for tunneling to a stau ufb-3 direction or ccb-4 minimum is shorter than the age of the universe , corresponding to @xmath74 , over a large portion of the otherwise acceptable parameter space . as expected , the newly disallowed regions are those with large and negative values of @xmath25 and @xmath46 at the input scale @xmath75 . of these two soft masses , only @xmath25 appears in the potential relevant for the ufb-3 direction or the ccb-4 minimum , and it has the stronger effect . indeed ,
the isocontours of @xmath52 coincide roughly with lines of constant @xmath25 for smaller values of @xmath46 . for smaller values of @xmath25
the sm vacuum becomes sufficiently long - lived to describe our universe .
the bounce action for these points is much larger simply because the effective width of the barrier is larger .
note that this stable region includes the minimal gaugino mediation point , @xmath76 at @xmath75 .    comparing the plots for @xmath67 and @xmath77
, we see that larger @xmath14 tends to yield a slightly more stable sm vacuum .
this arises simply because increasing the input gaugino masses also increases the low - scale slepton soft masses through rg running . on the other hand
, larger gaugino masses also permit more negative values of @xmath25 , so there remain significant parameter regions in which the tunneling rate is too fast . from fig .
[ m12500b30 ] , where @xmath67 , we see that nearly the entire region in which the lightest superpartner is a neutralino is ruled out by our vacuum stability considerations . for @xmath71 , there is a small region in which the lightest superpartner is a neutralino and the electroweak vacuum is sufficiently long - lived .
the bounce action for tunneling to a stau ufb-3 direction as a function of the hens model parameters @xmath25 and @xmath46 .
the other hens parameters have been fixed to be @xmath78 , @xmath66 , and @xmath73 .
all points shown are consistent with collider phenomenology .
the points between the two dash - dotted lines have a neutralino lsp .
@xmath69 is cosmologically safe.,scaledwidth=70.0% ]    we turn next to the cases with @xmath66 , for which the tau yukawa coupling is less enhanced . in figs .
[ m12500b10 ] and [ m12300b10 ] we show ranges of the bounce action for tunneling from the sm vacuum to the ufb-3 direction for @xmath79 and @xmath80 as a function of the input values of @xmath25 and @xmath46 at @xmath75 .
no ccb-4 local extremum is found for any of the points scanned over since the deviation from @xmath21-flatness is a sizeable effect in this case .
the dot - dashed lines in these plots enclose the portion of parameter space in which the lsp is the lightest neutralino .
( see ref .
@xcite for more details . ) from these plots we see that the lifetime of the standard electroweak vacuum against tunneling to a ufb-3 direction is safely large for both @xmath81 .
this is the result of the smaller value of the tau yukawa coupling @xmath82 , which gives rise to a larger barrier against tunneling , eq .  .
as before , with all else equal @xmath52 increases with @xmath14 and so higher values of @xmath14 remain safe .     the bounce action for tunneling to a stau ufb-3 direction as a function of the hens model parameters @xmath25 and @xmath46 .
the other hens parameters have been fixed to be @xmath83 , @xmath66 , and @xmath73 .
all points shown are consistent with collider phenomenology .
the points between the two dash - dot lines have a neutralino lsp .
@xmath69 is cosmologically safe.,scaledwidth=70.0% ]    at smaller values of @xmath29 the top quark yukawa coupling grows larger , and we should check that the tunneling rate to ccb-4 minima associated with the stops is adequately small . as was previously discussed , this direction will only occur if @xmath46 is large and negative . in fig .
[ m12500b10ccb ] we show contours of the bounce action for tunneling to the stop ccb-4 minimum , as well as the regions in which the sm minimum is the true minimum .
only a very few points at the largest and most negative values of @xmath46 are excluded , while in the great majority of the parameter space the sm vacuum is the true minimum .
we find a similar result for @xmath84 .     the bounce action for tunneling to a stop ccb-4 minimum as a function of the hens model parameters @xmath25 and @xmath46 .
the other hens parameters are @xmath78 , @xmath66 and @xmath73 .
all these points are consistent with collider phenomenology .
the points between the two dash - dot lines have a neutralino lsp .
@xmath69 is cosmologically safe.,scaledwidth=70.0% ]
in the present work we have examined the constraints placed on the hens model from vacuum stability . due to
the large tachyonic higgs soft masses that can emerge in this model , there often arise local vacuum states deeper than the standard electroweak minimum .
many points that are consistent with collider phenomenological constraints ( described in ref .
@xcite ) are ruled out because they lead to an overly short - lived sm vacuum .
the most dangerous vacuum feature is a ufb-3 direction involving the stau fields .
we have also found a new ccb-4 saddle point that facilitates tunneling to the ufb-3 direction . as a result ,
vacuum tunneling rates tend to be too fast for larger values of @xmath29 . at lower values of @xmath29 ,
tunneling to a ccb-4 direction involving stop fields rules out a very small portion of the parameter space , with the rest of the parameter space being safely long - lived .
we conclude that the hens models with a neutralino lsp and larger values of @xmath68 are mostly ruled out subject to our assumptions about the input scale and gaugino mass universality .
on the other hand , minimal gaugino mediation and the hens models without a neutralino lsp may still have a sufficiently long - lived electroweak vacuum state at larger values of @xmath29 . for lower values of @xmath29 , such as @xmath85 , the constraints from vacuum tunneling are much weaker , and most of the parameter space that is consistent with collider lower - energy bounds remains viable .
we thank f. adams , p. kumar , j. mason , and d. poland for helpful discussions . this work is supported in part by the u.s .
department of energy under doe grant de - fg02 - 95er40899 .
j.e . would also like to thank the cern th department for its hospitality while this work was in progress .
the bounce action is a stationary point of the euclidean action given in eq .
subject to the boundary conditions of eqs .
( [ bcond1],[bcond2 ] ) .
the corresponding equations of motion for the @xmath57 symmetric solution are @xmath86 where @xmath50 runs over the independent fields .
these equations are a set of non - linear coupled differential equations with an _ a priori _ unknown starting point .
these conditions together make it a very difficult problem to solve and require numerical techniques  @xcite .
the technique we use in the present work is called the _ improved action method _  @xcite . in this method ,
additional terms are added to the action that are identically zero for the bounce solution .
the advantage of adding these terms is that they make the bounce solution a minimum of this modified action and not just an extremum .
the term that does this is found by making the change of variable @xmath87 in eq .
( [ bouact ] ) .
because the bounce is the extremum of the action , the first derivative of the scaled action with respect to @xmath88 will be zero for @xmath89 .
this gives the following condition : @xmath90 + 2\sigma^4u[\bar{\phi}_i]\right)|_{\sigma=1 } = 0 \label{scainv}\ ] ] this relation illustrates that the potential term must be negative
. the kinetic term can not be negative because it is the integral of a sum of squares . since the potential term scales as @xmath91 and the kinetic term scales as @xmath92 , eq .
defines a maximum .
thus , we have determined the maximal direction of the saddle point . by adding to the action
the absolute value of this quantity to a positive power , the saddle point of the action can be turned into a minimum . in this case
the improved action is @xmath93=t[\bar{\phi}_i]+u[\bar{\phi}_i]+\lambda \left|t[\bar{\phi}_i]+2u[\bar{\phi}_i]\right|^n\label{modact}\ ] ] where @xmath94 and @xmath95 are positive constants .    to solve for the bounce with this improved action
, we take an initial profile for the vevs @xmath96 with the kinetic and potential terms @xmath97=2\pi^2\delta^4\sum_{m=1}^{l-1 } ( \rho_{m+1}-\rho_m)\rho_m^3\left[\sum_{i}^n \frac{(\phi_i^{m+1}-\phi_i^m)^2}{2(\rho_{m+1}-\rho_m)^2\delta^2}\right ] \\
u[\phi_i]=2\pi^2\delta^4\sum_{m=1}^{l-1 } ( \rho_{m+1}-\rho_m)\rho_m^3u(\phi_{1}^m, ... ,\phi_{n}^m).\end{aligned}\ ] ] @xmath98 is a parameter determined by eq .  .
inspired by the thin - wall approximation , we take the following initial guess for the bounce solution @xmath99 the coefficients @xmath100 and @xmath101 can be solved for by applying the boundary conditions @xmath102 and @xmath103 , where @xmath104 are the field values in the sm minimum and @xmath105 are the values in the vacuum to which the tunneling connects .
since @xmath105 is _ a priori _ unknown , this leaves @xmath105 , @xmath106 , and @xmath107 as free parameters .
these parameters are determined by first substituting the field profile in eq .   into the modified action for each @xmath96 .
the @xmath105 will not in general be the field configuration of the minima , but rather some field points inside the well . in the case of a ufb-3 direction
there is no minimum , and @xmath105 will be some point on the runaway downslope with a potential energy less than that at sm minimum .
the exact value of @xmath105 as well as @xmath106 and @xmath107 are determined using a minimization routine that finds the coefficients that minimize the modified action . in ref .
@xcite the authors used the thick - wall approximation as an initial guess . in our case , the guess given in eq .
works better numerically because it is adaptable to both thick- and thin - wall potential profiles .
once the initial profile , coefficients and all , is determined , we randomly vary each lattice site .
the variations are stopped when further iterations do not reduce the modified action . to ensure that we arrive at the bounce solution
, we choose a value of @xmath94 that ensures @xmath108/2/v[\phi_i])^{1/2}<1.001 $ ] .
the smallest value of @xmath94 able to maintain this condition and used to optimize the code was close to @xmath109 with @xmath110 .
we exhibit here an exact bounce solution for a single field tunneling through an inverted parabolic potential . to the best of our knowledge
, this simple solution has not been presented elsewhere in the literature , although a related solution for a linearized potential was obtained in the pioneering work of ref .
@xcite .
consider the run - away potential v ( ) = a|| - b||^2 , where @xmath111 is a real scalar field .
the origin is a metastable minimum for @xmath112 .
focussing on @xmath113 , the equation of motion for an @xmath57-symmetric bounce is + + ( 2b - a ) = 0 , [ del ] where @xmath114 is the euclidean distance .
the boundary conditions for the bounce are ( = 0 ) & = & 0[bc1 ] + ( r ) & = & 0[bc2 ] + ( = r ) & = & 0[bc3 ] for potentials in which the false vacuum is flat , @xmath115 , as we assumed in eq .  . here
, the potential has a singular first derivative at the false vacuum at the origin . as a result
, the bounce reaches the origin at finite @xmath8 and remains there asymptotically  @xcite .
we have checked that deforming the potential into a concave parabola very near the origin to resolve this singularity leads to @xmath115 , but otherwise has only a small ( and smooth ) effect on the solution presented here away from the origin .
is a linear second - order differential equation that can be solved analytically .
the general solution before imposing boundary conditions is ( ) = + , where @xmath116 and @xmath117 are bessel functions of the first kind , x = , and @xmath118 and @xmath119 are constant coefficients .
the values of @xmath118 and @xmath119 , as well as @xmath8 , are fixed by the boundary conditions , eqs .
( [ bc1]-[bc3 ] ) .
the first of these , eq .  , gives 0&= & .()|_x=0 + & = & -_x0 ( ) , where we have made use of the recursion properties of bessel functions  @xcite .
for this limit to be non - singular , we must have @xmath120 . on the other hand
, @xmath121 vanishes at the origin so no constraint is imposed on @xmath118 by this condition .
applying eq .
, we find 0 & = & c_1.|_x = x_+ + & = & -c_1j_2(x_+),where we have defined @xmath122 . from this , we conclude that @xmath123 must be a zero of @xmath124
. in particular , for the minimal action bounce solution , it should be the first non - trivial zero : x_+ 5.1356  finally , let us apply the condition of eq .   to fix @xmath118 , c_1 = - .
s_e & = & 2 ^ 2_0^rd^3 [ se ] + & = & x_+^4 + & & ( 215).the value of the field @xmath111 at the escape point where it emerges from tunneling is given by ( 0 ) = - _ x0 . for comparison ,
the maximum of the potential occurs at @xmath125 and the width of the potential barrier is @xmath126 .
thus , the escape point is well beyond the peak , but only by a factor of a few .
the exact solution presented here can also be extended to potentials that are piecewise segments of parabolas , in analogy to the solution for piecewise - linear potentials presented in ref .
( when the local potential is a concave parabola , the general local solution consists of modified bessel functions of the first kind . )
this offers the possibility of obtaining closed - form expressions for bounce solutions to a wide variety of potentials by approximating them with segments of parabolas . in practice
, we find that the matching conditions between adjacent parabolic segments often lead to complicated implicit transcendental equations .
the generalization to thermal transitions over piecewise - parabolic barriers is also straightforward .      for general reviews ,
see : + s.  p.  martin , hep - ph/9709356 ; d.  j.  h.  chung , l.  l.  everett , g.  l.  kane , s.  f.  king , j.  d.  lykken and l.  t.  wang , phys .
rept .   * 407 * , 1 ( 2005 ) [ hep - ph/0312378 ] .
f.  gabbiani , e.  gabrielli , a.  masiero and l.  silvestrini , nucl .
b * 477 * , 321 ( 1996 ) [ arxiv : hep - ph/9604387 ] . m.  misiak , s.  pokorski and j.  rosiek , adv .
high energy phys .   * 15 * , 795 ( 1998 ) [ arxiv : hep - ph/9703442 ] . for a review ,
see a.  b.  lahanas and d.  v.  nanopoulos , phys .
rept .   *
145 * , 1 ( 1987 ) ; a.  b.  lahanas ,  no scale supergravity : a viable scenario for understanding the susy cern - th-7092 - 93 _ lectures given at international school of subnuclear physics : 31th course : from supersymmetry to the origin of space - time , erice , italy , 4- 12 jul 1993 _      a.  e.  nelson and m.  j.  strassler , jhep * 0009 * , 030 ( 2000 ) [ arxiv : hep - ph/0006251 ] ; a.  e.  nelson and m.  j.  strassler ,  exact results for supersymmetric renormalization and the supersymmetric jhep * 0207 * , 021 ( 2002 ) [ arxiv : hep - ph/0104051 ] . m.  a.  luty and r.  sundrum , phys .
d * 65 * , 066004 ( 2002 ) [ arxiv : hep - th/0105137 ] .
m.  luty and r.  sundrum , phys .
d * 67 * , 045007 ( 2003 ) [ arxiv : hep - th/0111231 ] .
m.  ibe , k.  i.  izawa , y.  nakayama , y.  shinbara and t.  yanagida , phys .
d * 73 * , 015004 ( 2006 ) [ arxiv : hep - ph/0506023 ] .
m.  ibe , k.  i.  izawa , y.  nakayama , y.  shinbara and t.  yanagida , phys .
d * 73 * , 035012 ( 2006 ) [ arxiv : hep - ph/0509229 ] .
m.  schmaltz and r.  sundrum , arxiv : hep - th/0608051 .
m.  dine , p.  j.  fox , e.  gorbatov , y.  shadmi , y.  shirman and s.  d.  thomas , phys .
d * 70 * , 045023 ( 2004 ) [ arxiv : hep - ph/0405159 ] .
a.  g.  cohen , t.  s.  roy and m.  schmaltz , jhep * 0702 * , 027 ( 2007 ) [ arxiv : hep - ph/0612100 ] ; t.  s.  roy and m.  schmaltz , phys .
d * 77 * , 095008 ( 2008 ) [ arxiv:0708.3593 [ hep - ph ] ] .
h.  murayama , y.  nomura and d.  poland , phys .
d * 77 * , 015005 ( 2008 ) [ arxiv:0709.0775 [ hep - ph ] ] .
g.  perez , t.  s.  roy and m.  schmaltz , arxiv:0811.3206 [ hep - ph ] .
m.  schmaltz and w.  skiba , phys .
d * 62 * , 095005 ( 2000 ) [ arxiv : hep - ph/0001172 ] ; m.  schmaltz and w.  skiba , phys .  rev .
d * 62 * , 095004 ( 2000 ) [ arxiv : hep - ph/0004210 ]
. w.  buchmuller , j.  kersten and k.  schmidt - hoberg , jhep * 0602 * , 069 ( 2006 ) [ arxiv : hep - ph/0512152 ] ; w.  buchmuller , l.  covi , j.  kersten and k.  schmidt - hoberg , arxiv : hep - ph/0609142 .
s.  p.  martin and m.  t.  vaughn , phys .
d * 50 * , 2282 ( 1994 ) [ erratum - ibid .
d * 78 * , 039903 ( 2008 ) ] [ arxiv : hep - ph/9311340 ] .
j.  f.  gunion , h.  e.  haber and m.  sher , nucl .
b * 306 * , 1 ( 1988 ) .
m.  claudson , l.  j.  hall and i.  hinchliffe , nucl .
b * 228 * , 501 ( 1983 ) .
j.  a.  casas , a.  lleyda and c.  munoz , nucl .
b * 471 * , 3 ( 1996 ) [ arxiv : hep - ph/9507294 ] .
a.  riotto and e.  roulet , phys .
b * 377 * , 60 ( 1996 ) [ arxiv : hep - ph/9512401 ] .
a.  kusenko , p.  langacker and g.  segre , phys .
d * 54 * , 5824 ( 1996 ) [ arxiv : hep - ph/9602414 ] .
m.  s.  carena , m.  quiros and c.  e.  m.  wagner , phys .
b * 380 * , 81 ( 1996 ) [ arxiv : hep - ph/9603420 ] ; m.  carena , g.  nardini , m.  quiros and c.  e.  m.  wagner , arxiv:0809.3760 [ hep - ph ] .
m.  dine , l.  randall and s.  d.  thomas , phys .
lett .   * 75 * , 398 ( 1995 ) [ arxiv : hep - ph/9503303 ] .
j.  r.  ellis , k.  a.  olive and y.  santoso , phys .
b * 539 * , 107 ( 2002 ) [ arxiv : hep - ph/0204192 ] ; j.  r.  ellis , t.  falk , k.  a.  olive and y.  santoso , nucl .
b * 652 * , 259 ( 2003 ) [ arxiv : hep - ph/0210205 ] ; h.  baer , a.  mustafayev , s.  profumo , a.  belyaev and x.  tata , jhep * 0507 * , 065 ( 2005 ) [ arxiv : hep - ph/0504001 ] .
s.  r.  coleman , phys .
d * 15 * , 2929 ( 1977 ) [ erratum - ibid .
d * 16 * , 1248 ( 1977 ) ] . c.  g.  .
callan and s.  r.  coleman , phys .
d * 16 * , 1762 ( 1977 ) .
s.  r.  coleman , v.  glaser and a.  martin , commun .
phys .   * 58 * , 211 ( 1978 ) .