in the early universe , a scalar field often oscillates coherently after inflation  @xcite .
coherent oscillation of a scalar field may have various cosmological effects such as generation of curvature perturbation , dark matter , baryon asymmetry and unwanted relics .    in a simple model of a real scalar field
@xmath4 with a canonical kinetic term and a scalar potential of @xmath5 , it is known that the scalar field stops at the initial position @xmath6 until the hubble parameter becomes equal to the effective mass of the scalar .
after that , @xmath4 begins a coherent oscillation for @xmath7 with @xmath8 being the critical exponent , which depends on the background equation of state ( eq .  ( [ cond1 ] ) ) . in this case , the mass scale of @xmath4 becomes larger and larger compared with the hubble scale and such a coherent scalar dissipates its energy due to the hubble expansion .
the amplitude of the coherent oscillation @xmath9 decreases as  @xcite @xmath10 with @xmath11 being the scale factor of the universe .
such a scalar field can act as a curvaton  @xcite , which generates the primordial density perturbation of the universe  @xcite . in particular , curvaton models with @xmath12
were studied in some literature  @xcite .    on the other hand , for @xmath13 ,
the scalar field never oscillates and it just follows a so - called scaling solution  @xcite ) scaling solution .
we do not regard an oscillating solution as a scaling solution even if its amplitude satisfies a simple power law like ( [ scale1 ] ) . ]
@xmath14 this is an attractor and it was pointed out that such a scalar field can not generate curvature perturbation since the initial fluctuation of @xmath6 is erased during the scaling regime  @xcite .    in this paper , we revisit the scalar dynamics for @xmath15 . in ref .
@xcite it is stated that in this critical case , the scalar field is not attracted to a scaling solution but it is something like between a scaling and oscillation .
nevertheless , they say that for @xmath15 both the oscillating solution ( [ scale1 ] ) and scaling solution ( [ scale2 ] ) show the same behavior with the cosmic expansion and hence the curvature perturbation also vanishes .
however , we will see that the actual dynamics of the scalar field for @xmath15 is far from both the oscillating solution ( [ scale1 ] ) and scaling solution ( [ scale2 ] ) .
even averaged features of the scalar dynamics can not be described by these solutions .
instead , a peculiar feature of the scalar dynamics shows up , which we call a `` pseudo scaling solution '' .
we will show it both in a numerical and analytical way and see that the curvature perturbation does not vanish in this case .
the pseudo scaling solution can have a variety of physical implications .
as we noted above , the curvature perturbation can be generated by a curvaton even if it follows the pseudo scaling solution , and hence it may offer the seed of the large scale structure of our universe .
another example is a symmetry breaking field such as the peccei - quinn field  @xcite . in the supersymmetric ( susy ) axion model ,
it is possible that the peccei - quinn field follows the pseudo scaling solution in the early universe . as we will see later
, we can naturally solve the isocurvature perturbation as well as the domain wall problem in such a case .
the pseudo scaling behavior also has some implications on the baryon number asymmetry produced by the affleck - dine mechanism  @xcite .
therefore , it is worth studying detailed features of the pseudo scaling solution .
the organization of this paper is as follows . in sec .
[ sec : mono_pot ] , we discuss the dynamics of the scalar field with a monomial potential in a general background equation of state . here
the  pseudo scaling " scalar dynamics is extensively studied .
also , some comments on the hubble mass term are given there .
cosmological implications of the pseudo scaling solution are explained in sec .
[ sec : cos ] .
we discuss the curvature perturbation and the domain wall problem there .
the last section is devoted to conclusions and discussion .
[ sec : mono_pot ]      let us consider the lagrangian @xmath17 where the scalar potential is given by the monomial one @xmath18 where @xmath12 is an integer .
, @xmath4 should be interpreted as @xmath19 . ]
the equation of motion of @xmath4 is given by @xmath20 with @xmath21 being the hubble parameter .
we take @xmath22 and assume @xmath23 . in the matter ( radiation)-dominated ( md ( rd ) ) universe , @xmath24 . a more general equation of state may be obtained in the inflaton oscillation dominated universe ( see app .
[ sec : app ] ) .
let us define @xmath25 with @xmath11 being the scale factor of the universe and @xmath26 the initial time .
we have taken @xmath27 .
then the equation of motion becomes  @xcite @xmath28 where the prime denotes derivative with respect to @xmath29 and @xmath30 the equation ( [ eqofm1 ] ) represents the motion of the rescaled field @xmath31 in the effective potential @xmath32 in these variables , @xmath31 has a minimum at @xmath33 where @xmath34 thus if @xmath35 the sign of the friction term is positive , and hence the oscillation of @xmath31 damps and it relaxes to the minimum of the potential @xmath33 .
let us suppose the condition ( [ cond1 ] ) is satisfied so that there is a nontrivial minimum of the effective potential and @xmath31 relaxes there . , and @xmath36 for @xmath37 .
thus , @xmath38 is ensured and there is a nontrivial minimum of the potential if @xmath39 .
it does not always hold if there is a hubble mass term in the potential ( see eq .  ) . ]
this represents the scaling solution : in terms of the original field @xmath4 , @xmath40 in this case , independently of the initial condition , @xmath4 finally follows the scaling solution without oscillation .
note also that the scaling solution implies that the effective mass scale of the scalar decreases in proportion to the hubble parameter : @xmath41 .    on the other hand ,
if the inequality ( [ cond1 ] ) is inverted , the friction term of ( [ eqofm1 ] ) has a negative sign , and hence @xmath31 never relaxes to its potential minimum .
this corresponds to an oscillating solution .
this behavior is also understood as follows .
assuming an oscillating solution for @xmath4 with potential @xmath42 , we obtain @xmath43 by making use of the virial theorem with @xmath9 being the oscillation amplitude of @xmath4 .
thus the effective mass of @xmath4 scales as @xmath44 .
if this can continue to exceed the hubble parameter @xmath45 , the oscillation survives .
this condition is written as @xmath46 this is just an inversion of the condition ( [ cond1 ] ) .
we have seen that a scalar field approaches to a scaling solution for @xmath47 and oscillating solution for @xmath48 independently of the initial condition .
however , the case of @xmath15 is more nontrivial . in this case
, the friction term of ( [ eqofm1 ] ) vanishes : @xmath49 in terms of @xmath31 , it does not universally approach to some fixed solutions and continues to oscillate without damping .
the resulting dynamics significantly depends on the initial condition of @xmath4 . if the initial condition is chosen so that @xmath31 oscillates around @xmath33 without crossing @xmath50 , @xmath4 also does not cross @xmath51 .
the resulting dynamics of @xmath4 looks like neither scaling nor oscillating solution ; we call this a `` pseudo scaling solution '' .
hereafter we consider the case of @xmath15 .    before solving the equation of motion
, we must specify the initial condition .
the natural initial condition is set just after inflation as @xmath52 and @xmath53 , where @xmath54 and @xmath55 are the hubble parameter and the cosmic time at the end of the inflation , respectively .
we assume @xmath56 without loss of generality . here
@xmath6 should satisfy @xmath57 since otherwise @xmath4 rolls down the potential during inflation .
this condition roughly means @xmath58 .
therefore , in terms of @xmath31 , it is initially located close to the origin and rolls down the effective potential ( [ veff ] ) toward the minimum . in particular
, we consider the case @xmath59 by noting that the initial velocity of @xmath31 is given by @xmath60 we obtain the conserved energy density of @xmath31 : @xmath61 here we have used the approximation @xmath62 which is valid as long as @xmath63 .
since the energy density is negative , @xmath31 does not overshoot the origin @xmath50 , which means that @xmath4 does not oscillate around @xmath51 .
the rescaled field @xmath31 just oscillates around @xmath33 without damping .. although it holds during inflation , the relation is modified after inflation depending on the effective equation of state of the inflaton oscillation dominated universe as @xmath64 for @xmath65 .
the approximation @xmath66 as an initial condition is valid for inflation models in which the inflaton mass scale is much larger than the hubble scale , such as new inflation .
we have numerically checked that @xmath67 is indeed negative for such a case . on the other hand , for chaotic inflation , our approximation is not justified and @xmath67 may be positive .
however , the dynamics is almost identical to the case of negative @xmath67 after @xmath4 is replaced with @xmath19 even if @xmath4 once overshoots @xmath51 as long as @xmath68 is satisfied . ]     for @xmath15 for the potential @xmath0 for @xmath69 in the left panel and @xmath70 in the right panel .
we have normalized @xmath4 by the planck mass @xmath71 .
the time is normalized by @xmath72 , where @xmath73 is a natural mass scale of this system .
the red solid line represents a numerical solution , which corresponds to a `` pseudo scaling '' solution , while the green dashed line shows a hypothetical scaling solution for representative purpose .
, title="fig : " ]   for @xmath15 for the potential @xmath0 for @xmath69 in the left panel and @xmath70 in the right panel .
we have normalized @xmath4 by the planck mass @xmath71 .
the time is normalized by @xmath72 , where @xmath73 is a natural mass scale of this system .
the red solid line represents a numerical solution , which corresponds to a `` pseudo scaling '' solution , while the green dashed line shows a hypothetical scaling solution for representative purpose .
, title="fig : " ]    the dynamics looks so simple in terms of the rescaled field @xmath31 .
this `` pseudo scaling '' solution is more nontrivial in terms of the original field @xmath4 .
let us see fig.[fig : ps ] where we have numerically solved the equation of motion of @xmath4 with initial conditions @xmath74 for @xmath69 in the left panel and @xmath70 in the right panel .
we have normalized @xmath4 by the planck mass @xmath71 .
the time is normalized by @xmath72 , where @xmath73 is a natural mass scale of this system .
the red solid line corresponds to a numerical solution , which corresponds to a pseudo scaling solution , while the green dashed line shows a scaling solution @xmath75 for representative purpose . for comparison
, we also show time evolution of @xmath76 for cases with @xmath77 ( scaling ) , @xmath78 ( pseudo scaling ) and @xmath79 ( oscillating ) for @xmath80 in fig .
[ fig : ps2 ] .
it is clear from this figure that the pseudo scaling solution resembles neither the scaling solution nor the oscillating solution even qualitatively .
to understand this behavior , let us solve the equation ( [ eqofm2 ] ) more carefully .
our goal is to understand main features of the pseudo scaling solution and estimate how the periodicity depends on the initial condition @xmath81 . ) , but the following discussion is enough for our purpose . ] here it is important to note that the dynamics of @xmath31 is dominated by the region @xmath82 since @xmath83 in the case of our interest .
is satisfied . ]
it means that we can almost neglect the second term of eq . in order to understand main features of the dynamics . in the region @xmath82 , the equation of motion
is approximated as @xmath84 this has a general solution consisting of the growing and decaying mode , @xmath85 where @xmath86 and @xmath87 are determined by the initial conditions . in terms of @xmath4 , this solution corresponds to @xmath88 let us divide the time into the first phase @xmath89 , the second phase @xmath90 and the third phase @xmath91 . here
, we define @xmath92 as the time at which the motion of @xmath4 starts .
the time when the motion of @xmath4 stops is denoted as @xmath93 . some time after @xmath94 the field @xmath4 again starts to move , and we define that time as @xmath95 ( see fig.[fig : ps ] ) .
the @xmath96 term in the effective potential becomes important only for some short periods at around @xmath97 and @xmath98 .
[ [ first - phase ] ] first phase + + + + + + + + + + +    in the time interval @xmath99 , @xmath86 and @xmath87 are determined by the initial conditions : @xmath100 this solution breaks down at @xmath101 where @xmath31 becomes close to @xmath33 .
after that @xmath31 passes through @xmath33 and goes back to @xmath102 within order one unit time in terms of @xmath29 . we can estimate @xmath103 and @xmath92 as @xmath104 thus the time @xmath92 is roughly the epoch at which the hubble parameter becomes equal to the effective mass of @xmath4 .
it is obvious that @xmath4 remains almost constant during this period .
[ [ second - phase ] ] second phase + + + + + + + + + + + +    next , let us consider the interval @xmath90 .
just what we have to do is to determine the constants @xmath86 and @xmath87 in ( [ sol ] ) with appropriate initial conditions .
we set the initial condition at @xmath105 ( and @xmath106 ) , with @xmath107 taken somewhere between @xmath108 and @xmath33 , say , @xmath109 .
the precise position of @xmath107 is not important since the dynamics is dominated by the region @xmath82 as we noted before .
since @xmath67 is conserved during the whole process , we can express @xmath110 as @xmath111 where @xmath112 the sign of @xmath113 is flipped because @xmath31 changes the direction of motion due to the higher order term in @xmath114 at @xmath115 .
then , defining @xmath116 and @xmath117 as the constants @xmath86 and @xmath87 in eq .   at @xmath118
, we obtain @xmath119 hence @xmath120 .
then the end of the second period , @xmath121 , is the time at which the decaying mode @xmath122 becomes equal to the growing mode @xmath123 .
thus we obtain @xmath124 as @xmath125 where in the last equality we have approximated @xmath126 and @xmath127 . here
we have included an @xmath128 constant @xmath129 which comes from the uncertainty of @xmath107 as well as effects of the @xmath96 term .
we have numerically checked that the @xmath108 dependence of @xmath93 is precisely given by eq .   for @xmath63 , or @xmath130 .
it should be noticed that @xmath93 is orders of magnitude longer than the inverse mass scale of the scalar field . from ( [ solphi ] )
, we can see that @xmath4 scales as @xmath131 during this time interval , so we have @xmath132 and eventually at @xmath94 it becomes @xmath133 note that during this regime , @xmath4 satisfies @xmath134 , so the potential term becomes irrelevant and the effective mass scale @xmath135 becomes smaller and smaller than the hubble parameter .
the energy density of @xmath4 is dominated by the kinetic term and scales as @xmath136 during this regime .
[ [ third - phase ] ] third phase + + + + + + + + + + +    lastly , let us consider the third period @xmath91 .
this is just a continuation of the second period , where the growing mode dominates over the decaying mode .
thus the solution is @xmath137 this is valid until it becomes close to @xmath33 .
thus the end of the third period , @xmath138 is estimated as here is exactly the same as that in eq .   since the dynamics of @xmath31 is periodic . ]
@xmath139 therefore it is seen that the relative time interval @xmath140 is equal to @xmath141 . during this time interval , @xmath142 remains almost constant as is clear from ( [ solphi ] ) .
[ [ section ] ]    after @xmath143 , we have exactly the same dynamics as the second phase and the following third phase since there is no damping in the equation of motion of @xmath31 .
this one `` pseudo - oscillation '' consisting of one set of the second and third period takes much longer time than the inverse mass scale of the scalar field and even than the hubble time scale . in terms of @xmath4 ,
the field value reduces by an amount of @xmath144 during this one pseudo oscillation .
the scaling of @xmath4 consists of @xmath145 regime and @xmath146 regime , both are completely different from the scaling solution @xmath147 .
the scaling solution is obtained only after averaging the pseudo oscillation , but the period of the pseudo oscillation is much longer than the hubble time scale , hence using the scaling solution would lead to qualitatively wrong conclusions on the physical implications such as the curvature perturbation .
finally , we make some comments on qualitative differences between the scaling solution for @xmath148 and the pseudo scaling solution for @xmath149 . in the case of the scaling solution , the kinetic term , the hubble friction term and the potential term in the equation of motion are all comparable .
therefore , if the potential is not monomial , the motion of @xmath4 soon deviates from the scaling solution once other terms than the original one dominate the potential . on the other hand , in the case of the pseudo scaling solution , the kinetic term and the hubble friction term dominate over the potential term .
this can be seen by the fact that @xmath150 and @xmath151 are the solutions of the equation @xmath152 for @xmath153 .
, it means that the quadratic term dominates the effective potential . ] thus , even after the potential is dominated by other terms than the one with @xmath149 , the motion of @xmath4 still follows the pseudo scaling solution for some period until the potential term becomes comparable to the kinetic and friction terms .
changes , e.g. due to the inflaton decay , the pseudo scaling behavior terminates soon.]we will see this in more detail in sec .
[ sec : cos ] .      here
we briefly comment on the case where the hubble mass term exists in the potential : @xmath154 the equation of motion of @xmath4 is given by @xmath155 the equation of motion of the rescaled field @xmath31 is the same as ( [ eqofm1 ] ) except that @xmath156 is now replaced with @xmath157 in particular , for @xmath15 , the equation of motion becomes @xmath158 as long as @xmath159 , the existence of the hubble mass term does not qualitatively affect the pseudo scaling dynamics studied in the previous subsection .
if @xmath160 , the crucial difference from the case of @xmath159 is that the initial position of @xmath4 can not be taken freely . in particular , for positive @xmath161 , we have @xmath162 and there is no dynamics .
let us consider the case of negative hubble mass : @xmath163 . in this case
, @xmath4 initially sits at the minimum of the potential @xmath164 where @xmath165 this is the same order of @xmath33 during inflation , and hence @xmath31 begins oscillation around @xmath33 with an oscillation time scale of order one in terms of @xmath166 .
note that the temporal minimum of the potential @xmath167 exhibits the same scaling as the scaling solution @xmath75 . thus @xmath76 follows the temporal minimum @xmath168 with oscillation around it . in this case , this oscillation time scale is close to the hubble time scale .
this is qualitatively different from the case of @xmath159 in the previous subsection .
let us see the dynamics more closely .
in order to evaluate the initial amplitude @xmath169 around the temporal minimum @xmath167 , we must specify the origin of the hubble mass term of @xmath4 .
it arises from the planck - suppressed interaction as @xmath170 where @xmath171 and @xmath172 are constants of order unity , @xmath173 denotes the inflaton filed and @xmath174 is the inflaton potential . during a slow - roll inflation , the potential energy of the inflaton dominates ,
hence we have @xmath175 . after inflation
, the inflaton begins a coherent oscillation and the kinetic energy and potential energy are equally distributed , roughly speaking . in the case of harmonic oscillation of the inflaton
, we have @xmath176 . therefore , unless @xmath177 , the temporal minimum @xmath167 just after inflation is different from that during inflation .
this shift of the potential minimum occurs in the time scale of inflaton mass , rather than the hubble scale , so @xmath4 can not track the temporal minimum  @xcite .
thus the @xmath4 oscillation around the temporal minimum @xmath167 is induced just after inflation , and the amplitude is of the order of @xmath178 .
, @xmath4 can not track the temporal minimum  @xcite . ]
another term is coupling to the ricci scalar @xmath179 : @xmath180 if this is a dominant origin of the hubble mass term , @xmath181 during inflation and @xmath182 after inflation for the harmonic inflaton oscillation . oscillates between positive and negative values . ]    anyway , the oscillation amplitude around the temporal minimum @xmath169 decreases as the universe expands according to @xmath183 because of the number density conservation . by using this , we obtain the following scaling : @xmath184 the oscillation amplitude relatively decreases compared with the position of the temporal minimum for @xmath13 , while the ratio is constant for @xmath15 .
thus @xmath4 safely relaxes to its temporal minimum for @xmath13 . for @xmath15 ,
@xmath4 oscillates around @xmath167 and never cross @xmath51 thereafter if @xmath185 initially .
this fact is important to study the formation of topological defects , such as domain walls in the axion model@xcite .
so far we have considered the dynamics of scalar field with @xmath15 having only the non - renormalizable potential @xmath0 and the hubble mass term .
more realistic models may have other terms in the potential . in this section
we consider cosmological implications of the scalar dynamics with @xmath15 based on the potential @xmath186 the presence of the hubble mass term like @xmath187 does not change the result much as long as @xmath188 . in this type of the potential , the dynamics of @xmath4
is eventually dominated by the mass term at @xmath189 , even if @xmath4 obeys a pseudo scaling law at first .
this is easily seen by rewriting the equation of motion as @xmath190 since @xmath191 , the effect of the mass term becomes important at @xmath189 .
note that even if the potential is dominated by the mass term , the scalar field does not always feel the mass term if @xmath4 obeys a pseudo scaling law as we explained at the end of sec .
[ sec : ps ] . , where @xmath192 is defined below eq .
]    fig .
[ fig : mass ] shows numerical results for the time evolution of @xmath76 for @xmath78 for @xmath193 with a mass term . in the left ( right )
panel we have assumed positive ( negative ) mass squared .
the red solid ( green dashed ) line corresponds to @xmath194 @xmath195 .
other parameters and initial conditions are the same as those in fig .
[ fig : ps ] .
these results clearly show that the pseudo scaling law is maintained even after the mass term dominates the potential and the oscillation by the mass term starts at at @xmath196 .
first consider the abundance and the curvature perturbation of @xmath4 .
since @xmath4 is light during inflation , it obtains a quantum fluctuation of @xmath197 and @xmath4 may act as a curvaton - like field .
however , it is stated in ref .  @xcite that the curvature perturbation vanishes for both @xmath148 and @xmath15 because of the scaling behavior of the curvaton field .
if @xmath4 obeys a scaling solution ( [ scaling ] ) , the field value of @xmath76 at some uniform density slice does not depend on the initial position @xmath6 : this is easily checked as @xmath198 hence the @xmath6 dependence vanishes . here
@xmath199 is the time at which the effective mass becomes equal to the hubble scale ( see eq .
( [ t1 ] ) ) .
thus curvature perturbation is not generated by @xmath4 with a scaling law .
this is true for @xmath13 in which the scaling solution is an attractor , but for @xmath15 we have seen that the actual dynamics is far from the scaling law . for @xmath200 ,
@xmath4 scales as @xmath201 the initial @xmath6 dependence remains and hence the perturbation @xmath202 is not always erased .    for concreteness ,
let us suppose that @xmath4 decays into radiation at @xmath203 with a constant decay rate @xmath204 after the dynamics is dominated by the quadratic mass term for @xmath200 and oscillates around the origin .
we also assume that the inflaton energy density scales as @xmath205 until the curvaton enters in the quadratic oscillation regime , and decays into radiation before the curvaton decays , i.e. , we assume @xmath206 with @xmath207 being the inflaton decay rate . in this subsection
we mainly consider the case of positive mass squared @xmath208 and give only results for the case of negative mass squared @xmath209 at the end . for details on the case of negative
mass squared , we refer readers to app .  [
sec : hill ] .
the energy density of the curvaton at the onset of its oscillation , @xmath210 , is evaluated as @xmath211 where we have defined @xmath212 .
it represents the boundary at which the mass term begins to dominate the potential .
therefore , the fraction of the @xmath4 energy density at @xmath213 is given by @xmath214 this expression assumes that @xmath4 is subdominant at the decay .
otherwise , it is equal to @xmath215 .
now , let us evaluate the curvature perturbation generated in this scenario . from the expression ( [ rhom ] ) , we can see that the initial @xmath6 dependence explicitly remains .
thus the curvature perturbation of the order of @xmath216 is generated if @xmath4 eventually dominates the universe .      in the case of positive mass squared @xmath240 , to the second order in @xmath241
, we obtain the curvature perturbation as @xmath242\left(\frac{\delta\phi_i(\vec x)}{\phi_i } \right)^2,\end{aligned}\ ] ] where we have neglected the inflaton fluctuation .
the non - linearity parameter of the local type is given by , @xmath243 \frac{n_\phi^4}{(n_i^2+n_\phi^2)^2}.\ ] ] here we also show the results in the case of negative mass squared @xmath209 . if the mass term dominates the dynamics during the regime @xmath244 , the curvature perturbation is given by @xmath245,\end{aligned}\ ] ] where @xmath246 is a numerical factor .
the non - linearity parameter is given by @xmath247 for more details , see app .
[ sec : hill ] .
this is one example of the property of the curvature perturbation from the pseudo scaling curvaton .
although these expressions depend on how the background evolves and when the curvaton decays , it generally holds that the curvature perturbation is generated for @xmath15 in which @xmath4 obeys a pseudo scaling law .
next , we consider implications on the domain wall formation .
consider the case of negative mass squared : @xmath248 .
the potential has a @xmath249 symmetry @xmath250 if @xmath1 is even .
if @xmath4 obtains a large initial value during inflation and @xmath251 , it oscillates around the symmetric point @xmath252 after inflation and large fluctuations of @xmath4 particles ( and other particles which interact with @xmath4 ) are induced due to the efficient particle production .
such large fluctuations may cause the non - thermal restoration of the symmetry and @xmath4 may be trapped at the origin ( @xmath252 )  @xcite .
the trapping effect becomes diluted as the universe expands , and @xmath4 finally rolls down to the minimum of the potential .
eventually , the domain wall is formed and the universe is expected to be overclosed .
this is the rough picture of the domain wall problem we consider here .
a physically well - motivated example is the peccei - quinn ( pq ) field , which often leads to the axionic domain wall problem  @xcite due to the nonthermal restoration of the pq symmetry  @xcite .    on the other hand ,
for @xmath253 , the scalar field @xmath4 obeys the pseudo scaling law . as we saw earlier , the energy density of @xmath31 is negative for a reasonable initial condition ( see eq .  ) .
hence @xmath4 never crosses the origin during the pseudo scaling and no particle production occurs .
, no efficient particle production is expected to occur since the period of the oscillation is much longer than the cosmic time in the case of the pseudo scaling solution .
] moreover , even after the dynamics is dominated by the mass term , @xmath4 never crosses the origin because the mass term also prevents @xmath31 from crossing the origin ( see eq .  ) .
the field @xmath4 just oscillates around @xmath254 and the domain wall problem does not occur .
therefore , the pseudo scaling solution provides a natural solution to the domain wall problem .
here we emphasize that the pseudo scaling solution for @xmath253 is more powerful in solving the domain wall problem than the scaling solution for @xmath255 .
to see this point , let us consider the following potential : @xmath256 where we assume @xmath257 and @xmath258 .
first consider the scaling solution @xmath255 . in this case , if the initial value of @xmath4 is large enough , the field @xmath4 obeys the scaling solution after inflation .
if @xmath259 , the dynamics is dominated by the @xmath260 term and @xmath4 starts to oscillate around the origin when @xmath261 .
thus , for sufficiently large @xmath262 , the domain wall problem occurs even if @xmath4 obeys the scaling solution at first .    on the other hand , in the case of the pseudo
scaling solution @xmath253 , the dynamics becomes dominated by the @xmath260 term when @xmath263 where @xmath264 is defined as , and hence the dynamics is dominated by the @xmath260 term before the mass term becomes effective for @xmath265 . ]
@xmath266 we can deduce this as follows . in terms of @xmath31 ,
the effective potential is given by @xmath267 the amplitude of @xmath31 is of the order of @xmath33 , and hence the @xmath268 term becomes effective once the condition @xmath269 is satisfied .
this condition is equivalent to eq .  .
therefore , for example , if the @xmath260 term dominates the dynamics for the regime @xmath270 , the field value of @xmath4 at @xmath271 is given by @xmath272 if the initial amplitude satisfies @xmath273 , @xmath4 starts to oscillate with the amplitude much smaller than @xmath274 .
thus , the domain wall problem is less likely to occur even in this case compared to the scaling solution .
so far we have focused only on the domain wall problem , but similar implications can be obtained for other types of topological defects .
we have shown that a scalar field with potential @xmath0 with @xmath15 exhibits a peculiar behavior which we call a pseudo scaling .
it resembles neither oscillating nor scaling solution .
we have also shown that the primordial curvature perturbation is generated by this pseudo - scaling curvaton .
the critical exponent @xmath8 depends on the background expansion law : @xmath275 @xmath276 in the md ( rd ) universe .
even a renormalizable term can lead to a pseudo - scaling dynamics , i.e. , @xmath277 is obtained for @xmath278 , where the hubble parameter is given by @xmath279 .
this observation has implications on some realistic models and cosmological scenarios . in supersymmetric ( susy ) theories , for example , it is often the case that a scalar field has a non - renormalizable potential .
an example is the affleck - dine ( ad ) field  @xcite which consists of the flat direction in the mssm and have the superpotential of the form  @xcite @xmath280 where @xmath4 denotes the ad field and the scalar potential is then given by @xmath281 where @xmath282 and @xmath283 represent the effect of susy breaking .
if there is no hubble mass term or its coefficient is small , and the last term dominates the potential , the ad field may enter the pseudo - scaling regime .
the @xmath284 flat direction , which is lifted up by the @xmath285 superpotential , matches the condition for the pseudo scaling if the background is matter - dominated as in the inflaton oscillation dominated era .
in such a case , the lepton or baryon asymmetry generated by the ad mechanism will be significantly suppressed compared with the naive expectation , e.g. , the case of negative hubble mass term .
another example is the susy axion model , in which the superpotential is given by  @xcite @xmath286 where @xmath4 and @xmath287 are the pq field having pq charge @xmath288 and @xmath289 respectively .
the scalar potential is @xmath290 where we have taken @xmath291 for simplicity .
if @xmath292 , @xmath4 can have a large vev which may be suitable as a solution to the strong cp problem  @xcite .    in the case where the pq symmetry is broken before inflation ,
the axion field usually obtains a large fluctuation and may cause the isocurvature perturbation problem .
one way to avoid the isocurvature problem is to consider the situation where the pq field takes a very large value during inflation  @xcite .
in such a case , the isocurvature perturbation gets suppressed due to the large expectation value of the pq field . however , since the initial field value is large , particles are expected to be effectively produced due to the oscillation of the pq field after inflation . as a result , the pq symmetry may be non - thermally restored and hence the domain wall problem may occur after the qcd phase transition  @xcite . in this paper , we have seen that if the pq scalar obeys the pseudo scaling law after inflation , no such a disastrous situation is realized .
thus , for example , if @xmath293 and the pq scalar starts to move during the inflaton oscillation dominated era ( with the mass term ) , the domain wall problem as well as the isocurvature perturbation problem can be naturally solved .
this work was supported by the grant - in - aid for scientific research on scientific research a ( no.26247042 [ kn ] ) , young scientists b ( no.26800121 [ kn ] ) and innovative areas ( no.26104009 [ kn ] , no.15h05888 [ kn ] ) .
this work was supported by world premier international research center initiative ( wpi initiative ) , mext , japan .
the work of y.e . and m.t .
was supported in part by jsps research fellowships for young scientists and the program for leading graduate schools , mext , japan .
in this appendix we derive the expansion law @xmath294 or @xmath22 in the inflaton oscillation dominated universe .
we assume that the inflaton potential around the potential minimum is approximated as @xmath295 where @xmath173 denotes the inflaton field measured from its potential minimum .
the equation of motion of the inflaton field is given by @xmath296 from this we obtain @xmath297 where @xmath298 is the inflaton energy density .
on the other hand , we have @xmath299 taking time average over the inflaton oscillation time period and assuming that the hubble parameter is much smaller than the inflaton mass scale , we obtain the virial theorem @xmath300 where @xmath301 indicates time average over the inflaton oscillation period . by using this relation ,
we obtain @xmath302 thus the energy density scales as @xmath303 on the other hand , the freedman equation @xmath304 implies @xmath305 .
therefore , we obtain the expansion law @xmath306 for @xmath307 @xmath308 , the usual md ( rd)-like universe @xmath193 @xmath309 is realized .
more unconventional expansion laws are also obtained , for example , @xmath310 for @xmath311 and @xmath312 for @xmath313 . in terms of @xmath314 , the critical exponent @xmath315 is given by @xmath316
in this appendix , we give a formal solution of the equation of motion  .
note that the energy density of @xmath31 @xmath317 is conserved .
thus , the kinetic term of @xmath31 is given by @xmath318 by integrating this , we obtain the formal solution of @xmath31 as @xmath319 for @xmath320 region and @xmath321 for @xmath322 region . here
@xmath323 and @xmath324 are some initial values of @xmath29 and @xmath31 , respectively . in particular , the period of the oscillation of @xmath31 is given by @xmath325 where @xmath326 are the solutions of the equation @xmath327 with @xmath328 .. for @xmath329 , we must replace @xmath330 by zero and also double the right hand side of the equation when @xmath1 is even . ]
note that @xmath326 as well as @xmath331 depends on the initial condition .
however , it is difficult to solve this equation analytically .
moreover , the approximation given in sec .
[ sec : ps ] is sufficient for our purpose in this paper .
in this section we briefly review the curvature perturbation generated by the curvaton with a hilltop potential of the form @xmath332 as shown in sec .  [
sec : cos ] , if the initial position @xmath6 is set @xmath333 , where @xmath334 is the minimum position of the potential and @xmath335 , @xmath4 obeys a pseudo scaling law until @xmath21 becomes comparable to @xmath282 . for concreteness ,
let us suppose that this breaking of pseudo scaling law happens at @xmath91 .
the field value at this stage is given by ( [ phit2 ] ) : @xmath336 and this can be much smaller than @xmath192 .
the scalar dynamics thereafter is the same as that with a scalar potential ( [ tachyon ] ) with initial position @xmath337 close to the origin . in other words ,
the pseudo scaling scalar dynamics naturally provides us with a hilltop initial condition . in the following ,
we treat @xmath338 as if it is the `` initial '' condition of the dynamics , although the true initial position is @xmath6 . without loss of generality . ]    as in sec .
[ sec : abu ] , let us assume @xmath339 .
the abundance of the curvaton at @xmath213 is evaluated as @xmath340 this expression assumes that @xmath4 is subdominant at the decay . otherwise , it is equal to 1 . here
@xmath341 and @xmath342 denote the hubble parameter and the curvaton field value at the onset of its oscillation , respectively .
what is nontrivial in the hilltop scenario is that @xmath341 strongly depends on the initial condition @xmath337 in the hilltop limit @xmath343@xcite .
let us express @xmath341 and @xmath342 in terms of @xmath337 . before the onset of oscillation
, we can neglect the second term of ( [ tachyon ] ) .
the formal solution of the equation of motion with a tachyonic mass term @xmath344 and @xmath22 is written as @xmath345 where @xmath346 is the modified bessel function of the first kind .
in general , the solution is a superposition of the modified bessel function of the first and second kind .
however , without loss of generality we can consider the region @xmath347 and @xmath348 , and hence the bessel function of the first kind is sufficient to describe the dynamics .
the solution is conveniently approximated as @xmath349   & ~~~{\rm for}~~~mt \ll 1,\\       \displaystyle \phi _ * \frac{2^\nu\gamma(\nu+1)}{\sqrt{2\pi } } \frac{e^{mt}}{(mt)^{3p/2}}&~~~{\rm for}~~~mt \gg 1 .
\end{cases }      \label{phi_sol}\ ] ] first note that the velocity of the scalar field is given by , we have a solution for a general potential : @xmath350 . ]
@xmath351 since the oscillation epoch is significantly delayed in the hilltop limit @xmath343 , i.e. , @xmath352 , the condition of the start of the oscillation is given by @xmath353 on the other hand , from eq .
( [ phi_sol ] ) , we obtain @xmath354 for reasonable choice of @xmath355 , we have @xmath356 .
this justifies our assumption @xmath357 . then following relations
are derived , @xmath358 using these equations and noting that @xmath359 , we obtain @xmath360 where we have taken only leading terms in @xmath361 .
hence the curvature perturbation ( [ zeta ] ) is given by @xmath362,\end{aligned}\ ] ] to the leading order in @xmath361 . if @xmath337 would be a true initial condition set during inflation , then @xmath363 and this expression gives a correct estimate for the curvature perturbation from the hilltop curvaton . in this case , the non - linearity parameter of @xmath364 for @xmath365 is a generic prediction@xcite .    in our pseudo - scaling curvaton case ,
the true initial position is @xmath6 and it is @xmath202 that obtains a quantum fluctuation of @xmath366 . from the relation ( [ phi_ini ] ) , we obtain @xmath245 .
\end{aligned}\ ] ] the non - linearity parameter is given by @xmath247 therefore , if the inflaton contribution to the curvature perturbation is negligible , we have @xmath367 for @xmath365 .
thus in this case the non - linearity parameter may become sizable .
99 a.  h.  guth , phys .
d * 23 * , 347 - 356 ( 1981 ) ; a.  a.  starobinsky , phys .
b * 91 * ( 1980 ) 99 ; k.  sato , mon .  not .
soc .   * 195 * , 467 - 479 ( 1981 ) ; a.  d.  linde , phys .
b * 108 * ( 1982 ) 389 ; a.  albrecht and p.  j.  steinhardt , phys .
* 48 * , 1220 ( 1982 ) .
m.  s.  turner , phys .
d * 28 * , 1243 ( 1983 ) . s.  mollerach , phys .
d * 42 * , 313 ( 1990 ) ; a.  d.  linde and v.  f.  mukhanov , phys .
d * 56 * , 535 ( 1997 ) [ astro - ph/9610219 ] ; k.  enqvist and m.  s.  sloth , nucl .
b * 626 * , 395 ( 2002 ) [ hep - ph/0109214 ] ; d.  h.  lyth and d.  wands , phys .
b * 524 * , 5 ( 2002 ) [ hep - ph/0110002 ] ; t.  moroi and t.  takahashi , phys .
b * 522 * , 215 ( 2001 ) [ erratum - ibid .
b * 539 * , 303 ( 2002 ) ] [ hep - ph/0110096 ] .
k.  dimopoulos , g.  lazarides , d.  lyth and r.  ruiz de austri , phys .  rev .
d * 68 * , 123515 ( 2003 ) [ hep - ph/0308015 ] .
k.  enqvist , s.  nurmi , g.  rigopoulos , o.  taanila and t.  takahashi , jcap * 0911 * , 003 ( 2009 ) [ arxiv:0906.3126 [ astro-ph.co ] ] ; k.  enqvist , s.  nurmi , o.  taanila and t.  takahashi , jcap * 1004 * , 009 ( 2010 ) [ arxiv:0912.4657 [ astro-ph.co ] ] .
c.  t.  byrnes , k.  enqvist , s.  nurmi and t.  takahashi , jcap * 1111 * , 011 ( 2011 ) [ arxiv:1108.2708 [ astro-ph.co ] ] .
k.  mukaida , k.  nakayama and m.  takimoto , phys .
d * 89 * , no . 12 , 123515 ( 2014 ) [ arxiv:1402.1856 [ astro-ph.co ] ] .
a.  r.  liddle and r.  j.  scherrer , phys .
d * 59 * , 023509 ( 1999 ) [ astro - ph/9809272 ] .
r.  d.  peccei and h.  r.  quinn , phys .
lett .   * 38 * , 1440 ( 1977 ) .
i.  affleck and m.  dine , nucl .
b * 249 * , 361 ( 1985 ) .
m.  dine , l.  randall and s.  d.  thomas , nucl .
b * 458 * , 291 ( 1996 ) [ hep - ph/9507453 ] .
k.  harigaya , m.  ibe , m.  kawasaki and t.  t.  yanagida , arxiv:1507.00119 [ hep - ph ]
. k.  nakayama , f.  takahashi and t.  t.  yanagida , phys .
d * 84 * , 123523 ( 2011 ) [ arxiv:1109.2073 [ hep - ph ] ] .
a.  a.  starobinsky , jetp lett .   *
42 * , 152 ( 1985 ) [ pisma zh .
.  fiz .   * 42 * , 124 ( 1985 ) ] ; m.  sasaki and e.  d.  stewart , prog .
phys .   * 95 * , 71 ( 1996 ) [ astro - ph/9507001 ] ; m.  sasaki and t.  tanaka , prog .
phys .   * 99 * , 763 ( 1998 ) [ gr - qc/9801017 ] ; d.  h.  lyth , k.  a.  malik and m.  sasaki , jcap * 0505 * , 004 ( 2005 ) [ astro - ph/0411220 ] ; m.  sasaki , j.  valiviita and d.  wands , phys .
d * 74 * , 103003 ( 2006 ) [ astro - ph/0607627 ] .
l.  kofman , a.  d.  linde and a.  a.  starobinsky , phys .
lett .   * 76 * , 1011 ( 1996 ) [ hep - th/9510119 ] .
i.  i.  tkachev , phys .
b * 376 * , 35 ( 1996 ) [ hep - th/9510146 ] . for a review ,
see m.  kawasaki and k.  nakayama , ann .
nucl .  part .
sci .   * 63 * , 69 ( 2013 ) [ arxiv:1301.1123 [ hep - ph ] ] .
s.  kasuya , m.  kawasaki and t.  yanagida , phys .
b * 409 * , 94 ( 1997 ) [ hep - ph/9608405 ] ; phys .
b * 415 * , 117 ( 1997 ) [ hep - ph/9709202 ] ; m.  kawasaki , t.  t.  yanagida and k.  yoshino , jcap * 1311 * , 030 ( 2013 ) [ arxiv:1305.5338 [ hep - ph ] ] .
t.  moroi , k.  mukaida , k.  nakayama and m.  takimoto , jhep * 1306 * , 040 ( 2013 ) [ arxiv:1304.6597 [ hep - ph ] ] ; jhep * 1411 * , 151 ( 2014 ) [ arxiv:1407.7465 [ hep - ph ] ] .
t.  gherghetta , c.  f.  kolda and s.  p.  martin , nucl .
b * 468 * , 37 ( 1996 ) [ hep - ph/9510370 ] .
h.  murayama , h.  suzuki and t.  yanagida , phys .
b * 291 * , 418 ( 1992 ) . for a review ,
see j.  e.  kim , phys .
rept .   * 150 * , 1 ( 1987 ) ; j.  e.  kim and g.  carosi , rev .
* 82 * , 557 ( 2010 ) [ arxiv:0807.3125 [ hep - ph ] ] .
a.  d.  linde and d.  h.  lyth , phys .
b * 246 * , 353 ( 1990 ) ; a.  d.  linde , phys .
b * 259 * , 38 ( 1991 ) .
m.  kawasaki , k.  nakayama and f.  takahashi , jcap * 0901 * , 026 ( 2009 ) [ arxiv:0810.1585 [ hep - ph ] ] .
m.  kawasaki , t.  kobayashi and f.  takahashi , phys .
d * 84 * , 123506 ( 2011 ) [ arxiv:1107.6011 [ astro-ph.co ] ] .


two - band superconductivity was first proposed in the late 1950 s and further studied in the 1960 s and 1970 s as a possible explanation for understanding superconductivity in s- and d- band metals .
however , it was not until the discovery of superconductivity at relatively high temperature in mgb@xmath5 in 2001 , that the most promising example of two - band superconductivity was found .
after nearly four years of effort and the development of an extensive literature based on this material , it is now evident that this compound is an electron - phonon mediated superconductor , but with two distinct bands and hence two energy gaps .
indeed , this is firmly established from detailed bandstructure calculations that provide the electron - phonon spectral densities which , when in turn are used in a two - band eliashberg formalism , give rise to predictions which are in accord with experiment . in fact
, the picture that emerges is one of excellent agreement with the data , with the match being as good as is found between theory and experiment in conventional one - band superconductors@xcite .
a detailed review and summary of the level of agreement found for the thermodynamic properties , bcs ratios , penetration depth and temperature - dependent energy gaps is found in ref .
with such an extensive comparison between rigorous eliashberg theory and experiment completed , it is now appropriate to move on to the transport properties where less theoretical work has been done .    in this paper
, we consider the critical current which is of fundamental interest in the discussion of transport . while most technological applications focus on the bulk critical current density , which requires understanding of flux lattice pinning , etc .
, here we are presenting the critical current in thin films .
this geometry is more relevant to understanding issues associated with the fundamental superconducting state , rather than focussing on material details that lead to vortex pinning , etc . , and hence can be expected to provide further insight into the nature of multiband superconductivity in these materials .
we note in passing , at this point , that some authors use the term depairing current for the critical current although , technically , the depairing current is the current at which the first cooper pair is broken and the gap decreases , and the critical current is the maximum current possible , which is greater than the depairing current . in simple one - band s - wave bcs theory the values are almost the same and hence there is the tendency to use the terms interchangeably . here , these two quantities can be quite different and so we will use the term critical current to refer to the global maximum obtainable in the current density .
the full temperature - dependence of the critical current @xmath2 was first studied in conventional one - band s - wave bcs superconductors by rogers@xcite and bardeen@xcite .
further work following on this was done by parameter@xcite for the clean limit , and by maki@xcite for the dirty limit , where the latter author introduced a description of the superfluid current density using the green s function formalism .
a study of the effect of of a current on the quasiparticle states density of states of a superconductor was done by fulde@xcite .
the most general work was that of kupriyanov and lukichev@xcite , who calculated the full temperature dependence of the critical current in bcs theory for arbitrary impurity scattering , which returned the usual ginzburg - landau results for the critical current near @xmath6@xcite .
further considerations led to an initial attempt to include strong inelastic scattering@xcite followed by a full eliashberg calculation and study of strong coupling effects for conventional electron - phonon superconductors done by ourselves@xcite .
thus , the formalism and basic one - band results have been well - established in both the weak - coupling bcs limit and strong - coupling eliashberg regime , at all temperatures and for arbitrary impurity scattering .
more recently there has been a renewed theoretical effort on the critical current as applied to exotic gaps , which are anisotropic , such as d - wave@xcite and f - wave@xcite .
the two - band case is , in some sense , a highly anisotropic system and hence the study of the critical current in this instance also complements these recent works .    with regard to mgb@xmath5 ,
specifically , there has been some recent experimental work by kunchur@xcite , where the critical current has been measured in these systems . on the theoretical side , koshelev and
golubov@xcite have gone beyond ginzburg - landau to use the usadel equations to study this system and find interesting effects such as a kink in the temperature - dependence of the c - axis critical current that reflects the underlying two - band nature .
this type of feature is observed in other properties such as the penetration depth which also shows an inflection point in the temperature dependence , however , these effects may disappear with increased interband coupling .
finally , a more recent work has been done on the topic of the non - linear term in the superfluid current in mgb@xmath5 and its effect on the nonlinear microwave response , which is relevant for device applications@xcite .    in our work here , we use the general green s function formulation presented by maki@xcite where the cooper pairs are given a finite center of mass momentum @xmath1 .
this @xmath1 leads to a boost in the quasiparticle energy ( sometimes called a doppler shift ) .
the superfluid velocity @xmath7 is taken as uniform in the case of discussing currents in thin films or wires , unlike the doppler shift effect in the context of superfluid currents in the vortex state , where @xmath8 varies with spatial position .
this latter topic has been the subject of a large number of papers in the recent literature relating to d - wave superconductors@xcite .
a very fine physical description of the effect of this doppler shift on the quasiparticle density of states and the resulting effect on the superfluid density is presented in a paper by xu et al.@xcite , where they consider the non - linear meissner effect in d - wave superconductors from a starting point that discusses the basics of the superfluid current density .
thus , with a body of literature developing on the topic of critical currents and superfluid currents , it is of interest to elucidate the features in the critical current due to two - band superconductivity , where unusual effects can occur due to a transfer from one band to the other of the dominant contribution to superfluid current and hence the critical current .
our paper is structured in the following manner . in section  ii
, we introduce the basic theoretical equations which we have evaluated .
we then examine , in section  iii , results for renormalized bcs theory near @xmath6 and at @xmath9 and apply these results to mgb@xmath5 . in section
iv , we present our eliashberg results for finite temperature and with impurities , evaluating the critical current for mgb@xmath5 and other model parameters .
we also compare our results to the data for mgb@xmath5 .
finally , we form our conclusions in section  v.
while various approaches have been used in the past for calculating the superfluid current , we will use the green s function method as presented by maki@xcite .
thus the superfluid current density @xmath10 is given in terms of finite temperature green s functions as : @xmath11\quad , \label{eq : js1}\ ] ] where @xmath12 is the electron charge , @xmath13 , the electron mass , @xmath14 is the electron momentum , and @xmath15 is the temperature .
the green s function is given in terms of the matsubara representation by @xmath16^{-1 } , \label{eq : greensfn}\ ] ] where @xmath17 , @xmath18 is the chemical potential , and the @xmath19 and @xmath20 are the pauli matrices for the particle - hole and spin spaces , respectively .
the @xmath21 is the applied superfluid momentum and it results in a shift of the quasiparticle energies , given to first order in @xmath1 by @xmath22 . here ,
@xmath23 and @xmath24 are the imaginary axis renormalized superconducting order parameter and renormalized matsubara frequencies , respectively , and are defined below .
the integral over energy @xmath25 for the combined eqns .
( [ eq : js1 ] ) and ( [ eq : greensfn ] ) can be done and the result , generalized to two bands , is given in terms of a sum of two partial currents @xmath26 and @xmath27 , for the first and second band , respectively : @xmath28 with @xmath29 where @xmath30 , @xmath31 is the electron density and @xmath32 is the fermi velocity of the @xmath33th band ( @xmath34,2 ) . the nonlinear eliashberg equations for @xmath35 and @xmath36 have been generalized to two bands and with the inclusion of the effect of the @xmath1
, they are given as@xcite : @xmath37\nonumber\\ & \times&\int_{-1}^{1}\frac{dz}{2 } \frac{\tilde\delta_j(m)}{\sqrt{(\tilde\omega_j(m)-is_jz)^2+\tilde\delta_j^2(m)}}\nonumber\\ & + & \pi\sum_j(t^+_{ij}-t^-_{ij})\int_{-1}^{1}\frac{dz}{2 } \frac{\tilde\delta_j(n ) } { \sqrt{(\tilde\omega_j(n)-is_jz)^2+\tilde\delta_j^2(n)}}\nonumber\\ \label{eq : del}\end{aligned}\ ] ] and @xmath38 where @xmath39 and @xmath40 are the ordinary and paramagnetic impurity scattering rates , respectively , and the electron - phonon kernels @xmath41 with phonon energy @xmath42 enter through @xmath43 here , the @xmath44 indexes the @xmath44th matsubara frequency @xmath45 , with @xmath46 , where @xmath47 .
the @xmath48 are coulomb repulsions , with a high energy cutoff @xmath49 usually taken to be about six to ten times the maximum phonon frequency . as written , eqns .
( [ eq : js2])-([eq : z ] ) are for three dimensions with isotropic fermi surface . for general anisotropic fermi surface ,
a fermi surface integral would remain , which would need to be done numerically@xcite .
this goes beyond the scope of this paper , however , such details are not expected to change the qualitative features of our results but rather lead to quantitative changes only .    in the following sections , we present both eliashberg results and renormalized bcs ( rbcs ) results , the latter of which can be quite successful at capturing the essential features of the eliashberg calculations without the full numerical complications of the more sophisticated calculations even for two - band models@xcite . to develop the rbcs results we use the two - square - well approximation or @xmath50 model written for two - bands ( for details see ref .
@xcite ) , where now @xmath51\delta_j(t)}{\sqrt{(\omega_m - i\bar s_jz)^2+\delta_j^2}}. \label{eq : rbcsdel}\ ] ] in this expression , @xmath52 , @xmath53 and @xmath54 with @xmath55 taken to represent either the debye frequency or some other characteristic energy scale representing the phonons in the system . likewise ,
the current from eq .
( [ eq : js2 ] ) can be reduced to @xmath56 equations  ( [ eq : rbcsdel ] ) and ( [ eq : rbcsjs ] ) will be used in the next section to develop analytic results for comparison with full numerical calculation .
in order to evaluate the critical current near @xmath6 , we need to know the effect of the current on the gap near @xmath6 . in this limit
, both @xmath57 and @xmath1 will be small .
after long but straightforward algebra , to first order in @xmath58 where @xmath59 is the reduced temperature , the gap equation of eq .
( [ eq : rbcsdel ] ) reduces to the form @xmath60 with a second equation obtained from eq .
( [ eq : bcs1 ] ) by switching indices @xmath61 . here
@xmath62 and @xmath63 where @xmath64 and @xmath65 .
for @xmath66 as @xmath67 and with no superfluid momentum @xmath68 , we recover the equation for the critical temperature ( symmetric in @xmath61)@xcite : @xmath69 with @xmath70 . throughout this paper
, we will find it instructive to consider several simplifying limiting cases .
the completely decoupled case corresponds to @xmath71 and the assumption that @xmath72 . in this case
, the first channel determines the @xmath6 , which is given by @xmath73 . as an aside
, we note that the small separable @xmath74 anisotropy model studied extensively in the older literature@xcite maps on to this model if @xmath75 , @xmath76 and @xmath77 , with the anisotropy @xmath74 assumed small . in this case , eq .
( [ eq : bcs3 ] ) simplifies greatly since @xmath78 so that @xmath79 = ( 1+\lambda)/[\lambda(1+a^2)]$ ] which gives @xmath80)$ ] .
this is the well - known result that anisotropy increases the value of the critical temperature over its isotropic value ( @xmath81 ) .
returning to eq .
( [ eq : bcs1 ] ) which gives the gap just below @xmath6 , we need to specify @xmath82 .
it is given ( after more lengthy algebra ) as @xmath83 and an equivalent expression for @xmath84 is obtained from eq.([eq : bcs4 ] ) by switching indices @xmath61 .
solving for @xmath85 , we find @xmath86 with @xmath87 \label{eq : bcs6}\ ] ] and @xmath88 .
\label{eq : bcs7}\ ] ] likewise , the gap @xmath89 , is formed from the above three equations by exchanging @xmath61 . when @xmath90 , these equations properly reduce to those obtained previously by nicol and carbotte@xcite . in the limit of decoupled bands , @xmath91 and @xmath92 .
some care is required when treating the equivalent equations for @xmath93 and @xmath94 because the combination @xmath95 is replaced by @xmath96 , which is zero .
nevertheless , we can show @xmath97 , which is expected on physical grounds as the second band does not contribute near @xmath6 , where @xmath98 for @xmath99 , and @xmath100 $ ] , a well - known result .
we can also recover other well - known results for the small separable anisotropy model defined previously .
note first that when the fermi velocities are equal , @xmath101 in eq .
( [ eq : bcs7 ] ) reduces to @xmath102 , with @xmath103 given by eq .
( [ eq : bcs6 ] ) .
further , @xmath104 and @xmath105 which leads to ( eq .  ( [ eq : bcs5 ] ) ) @xmath106 \label{eq : bcs8}\ ] ] and an identical equation for the second gap with the first factor of @xmath107 replaced by @xmath108 . for @xmath90 , eq .  ( [ eq : bcs8 ] ) reduces to the known result that the average gap @xmath109 with @xmath110 and @xmath111 . for the isotropic case ( equivalent to one band ) @xmath81 and eq .
( [ eq : bcs8 ] ) reduces to @xmath112 as is well known .    the expression for the contribution to the current from the @xmath33th band
is given in renormalized bcs by eq .
( [ eq : rbcsjs ] ) .
expanding near @xmath113 for small @xmath1 and @xmath57 gives @xmath114 the total current is the sum over both bands and can be written as @xmath115 \frac{7\zeta(3)}{4(\pi t_c)^2 } eq_s , \label{eq : bcs13}\ ] ] where @xmath116 .
substituting into eq .
( [ eq : bcs13 ] ) the new expression from eq .
( [ eq : bcs5 ] ) ( and equivalent for @xmath61 ) for @xmath85 , one obtains two terms . the first proportional to @xmath1 and
the second to @xmath117 .
this gives the total current @xmath0 as a function of @xmath1 .
the critical current @xmath118 is obtained as the maximum of @xmath119 as a function of @xmath1 . after finding the extremum value of @xmath1
, we obtain for @xmath118 @xmath120^{3/2 } \biggl [ \frac{n_1}{m_1^*}\frac{1}{\chi_1^\prime}+\frac{n_2}{m_2^ * } \frac{1}{\chi_2^\prime}\biggr]^{-1/2}. \label{eq : bcs14}\end{aligned}\ ] ] we first note that the standard ginzburg - landau @xmath121 temperature dependence of the critical current remains unmodified in the two - band model . as a first check on our expression in eq .
( [ eq : bcs14 ] ) , we can recover the known one - band result when both bands are assumed to be identical by taking @xmath122 , @xmath123 , and the value of @xmath124 . here
@xmath125 is the total electron density per unit volume equal to the sum @xmath126 .
we obtain the standard ginzburg - landau result@xcite : @xmath127 note that all effective mass renormalization @xmath128 factors have cancelled in eq .
( [ eq : bcs15 ] ) . applying the separable gap anisotropy model to eq .
( [ eq : bcs14 ] ) leads to a modification of eq .
( [ eq : bcs15 ] ) by a multiplicative factor of @xmath129 .
this gap anisotropy will decrease @xmath118 over its value for @xmath130 .
decoupled bands give the usual expression of eq .
( [ eq : bcs15 ] ) since the second band makes no contribution near @xmath113 , but now @xmath44 , @xmath13 , and @xmath131 are those of the first band only . in the general case , we need to return to eq .
( [ eq : bcs14 ] ) which tells us how material parameters affect the slope of @xmath132 which varies linearly with @xmath58 near @xmath6 .
if we normalize the current , as we will in the next section , to a slope of negative one near @xmath6 , this will mean that the current at @xmath9 will vary with material parameters in contrast to the one - band bcs curve for which it is a universal number for the clean limit case that we are discussing here .
a well known result of ginzburg - landau theory is that near @xmath6 , @xmath133 , where @xmath134 is the thermodynamic critical field and @xmath135 is the london penetration depth , both for the case of zero current . in a previous paper@xcite ,
we have given expressions for @xmath134 and @xmath136 as @xmath99 .
using these with eq .
( [ eq : bcs14 ] ) , we arrive at the expression relating the two - band @xmath4 , @xmath3 and @xmath2 : @xmath137 where @xmath138 is the velocity of light and @xmath139 is a complicated material - dependent parameter .
it is given by @xmath140^{-1 } \label{eq : bcs18}\ ] ] where we have found it convenient to change from electron density @xmath31 , to density of states at the fermi surface @xmath141 .
this modifies the usual ginzburg - landau relation by a constant factor . for the decoupled case where
@xmath142 , @xmath143 , and @xmath144 , @xmath139 becomes equal to 1 .
there is no change in the ginzburg - landau relation and this makes sense since near @xmath145 we are dealing with the single band ( the first one ) . for the isotropic gap case , where @xmath122 for any @xmath146 , from eqs .
( [ eq : bcs6 ] ) and ( [ eq : bcs7 ] ) we obtain @xmath147 and @xmath148 , the average of the squares of the fermi velocities in the two bands , which gives @xmath149 for identical bands , this expression gives 1 , as it must .
this also holds if either the two fermi velocities are the same or the two density of states match , but not if both parameters for one band differ from those of the other band .
we see that , in general , @xmath139 can be either greater or less than one and the ginzburg - landau relation ceases to hold .    for identical bands but with gap anisotropy in the separable small anisotropy @xmath74 model , @xmath150 and @xmath151 , with @xmath104 and @xmath105 which gives @xmath152 to order @xmath74 .
this result is expected since , as is well known , the @xmath74 anisotropy correction to @xmath134 and @xmath153 are both of the form @xmath154 which gives a factor of @xmath129 in the ratio @xmath155 .
this agrees perfectly with the anisotropy factor we derived earlier for @xmath118 .
thus , in this case , the ginzburg - landau relation will hold .
we wish to comment on the use of the ginzburg - landau relation @xmath133 for comparison with data at all temperatures , which is done by substituting two - fluid model forms for @xmath134 and @xmath135 .
the ginzburg - landau result should only hold near @xmath6 and indeed we find that when we form the ratio of @xmath155 using the @xmath134 and @xmath135 from the full eliashberg procedure , the resulting @xmath118 curve rapidly deviates from the eliashberg calculated @xmath118 .
the results of using such a phenomenological procedure can give curves which are lower or higher than the true result , and will not even reproduce the qualitative feature of a kink in the @xmath118 that we have found in the two - band calculation .
we conclude that the ginzburg - landau procedure should only be used near @xmath6 as is appropriate to its derivation and no reliable use can be made of it phenomenologically for lower temperatures .
we now turn to low temperature and examine results of rbcs at @xmath9 . at zero temperature , the renormalized bcs gap equation ( [ eq : rbcsdel ] )
takes on the form @xmath156 where the barred quantities are renormalized by the factor of @xmath157 of eq .
( [ eq : rbcsz ] ) and @xmath158    ( 250,200 )    these equations need to be solved numerically .
results are shown in the four lefthand panels of fig .
the top frame is for @xmath159 , @xmath160 , @xmath161 ( nearly decoupled bands ) and @xmath123 .
because @xmath162 is so small , the larger gap @xmath163 exhibits the known behavior as a function of @xmath164 of a one - band s - wave superconductor@xcite with a small modification . in one - band ,
the gap would be unchanged until @xmath165 at which point it would begin a sharp drop toward zero which is reached at @xmath166 . here
, the same behavior occurs but there is a very slight drop from the @xmath90 value of @xmath57 at @xmath167 due to the coupling to the second band .
the lower gap behaves differently .
it retains its @xmath90 value until @xmath167 at which point it also begins a rapid drop but retains a small finite value until the point when @xmath168 .
this small tail is due to the coupling @xmath169 which guarantees that the second gap is nonzero , although it can be very small , as long as the first gap is finite . in the second frame @xmath169
has been increased to 0.1 .
this results in a considerable integration of the two gaps . at @xmath90 the lower gap ratio @xmath170
has increased from the 0.4 value of the top frame to 0.57 .
more significantly for the present paper , as @xmath1 increases both gaps retain their @xmath90 value until @xmath171 , which is the point at which the second gap starts to decrease , but this drop is not as fast as in the top frame because of the large interband coupling . at the same point as the second gap starts to decrease , the first is also reduced but by much less . beyond @xmath165 both gaps drop sharply as expected and reach zero together at @xmath172 , a value which is smaller than the one band value of 1.35 .
this shift in zero along the horizontal axis is also accompanied by a shift to the left , relative to the decoupled case , of the inflection points just described , but these shifts are small .
these are signatures of the two - band nature of our system . in the third frame of fig .
[ fig1 ] ( lefthand side ) , we have retained @xmath173 but now have increased the ratio @xmath174 to 2.0 instead of 1.0 .
this has the effect of reducing the region of @xmath1 in which the second gap remains constant at its @xmath90 value and is due to the increase in @xmath175 .
for the last frame , we have kept the parameters as in the first frame but have increased @xmath176 to a value of 0.7 .
this has the effect of increasing the @xmath90 value of the second gap and also increasing the range over which it stays constant .
additionally , the tail beyond @xmath167 is not as small as it is in the top frame due to less gap anisotropy .
( 250,200 )    in renormalized bcs , eq .
( [ eq : rbcsjs ] ) for the partial current can be evaluated analytically at zero temperature to obtain @xmath177^{3/2}\biggr],\nonumber\\ & & \qquad\qquad\qquad \qquad { \rm for }   \frac{\bar s_i}{\delta_{i0}}>1\end{aligned}\ ] ] results for the current are presented in the righthand side frames of fig .
[ fig1 ] for the four sets of parameters considered on the lefthand side . in all cases , the dotted lines give the results for the partial contribution to the current coming from the second band and the dashed for those from the first .
the sum of the two partial contributions is the solid curve .
a first feature to note is that for the three frames for which @xmath173 ( nearly decoupled ) , the dashed curves are essentially identical .
the straight line segment has slope 0.5 in our units and the zero is very nearly at @xmath178 as in the one - band case .
note that on the vertical axis , we have used the total electronic density of states @xmath179 with @xmath180 .
if we had used only @xmath181 then the slope of the straight line segment would have been 1 which corresponds to what is expected in the one - band case . a second point to note
is that increasing @xmath169 to 0.1 as in the second frame from the top has not had a large effect on the partial current from the first band except to reduce somewhat the value of @xmath1 at which the current goes to zero .
this is accompanied by a small shift to the left in the position of the peak and a small reduction in its height . by contrast , the effect of the parameter changes are much more significant on the dotted curve ( the current in the second band ) . in the top frame ,
the slope of the straight line portion of the dotted curve is a little larger than 0.5 because of the renormalization factors on @xmath175 relative to @xmath182 .
more significantly the partial current starts to deviate from linearity at @xmath183 at which point it drops rapidly towards zero but remains finite ( although small ) until the current in the first band has reached zero .
this also applies to the third frame although in that case the current in the second band has been increased in value because of the increase in @xmath184 by a factor of two . in all cases
the total current ( solid curve ) exhibits two peaks .
the peak arising at lower values of @xmath1 has a significant contribution from both bands while the other peak at higher values of @xmath1 is due almost entirely to the partial current of the first band when @xmath169 is very small . on the other hand ,
when there is more integration of the two systems ( second frame from top ) both bands can make a significant contribution . in the top three frames ,
the first peak is lower than the second while in the fourth frame , the opposite holds .
the height of the first peak can be made higher by increasing the value of @xmath184 as in the third frame compared with the first ( or equivalently , increasing @xmath185 as the density of states enters in combination with @xmath131 and so will have a similar effect , which is redundant to discuss here ) or by increasing @xmath176 as in the fourth , again compared with the first .
it can also be increased by increasing the offdiagonal coupling as in the second frame .
when this is done , however , the contribution of the second band to the total current in the second peak is also increased .
finally we note that the peak heights in the two middle frames are only very slightly different and the lower @xmath1 peak could be made to be the highest through a small tuning of the parameters . as the highest peak determines the value of @xmath2 , this transfer from one peak to the other by varying parameters is very interesting and will produce nonstandard effects at finite temperature .      before leaving our discussion of rbcs at @xmath9 ,
it is of interest to consider the specific case of mgb@xmath5 now well - established to be a two - band electron - phonon superconductor .
microscopic parameters have been calculated from extensions of bandstructure calculations to obtain the electron - phonon spectral densities .
the coulomb repulsion parameters that enter the gap equations ( [ eq : del ] ) and ( [ eq : z ] ) are also known .
the parameters are @xmath186 , @xmath187 , @xmath188 , @xmath189 , @xmath190 , @xmath191 , @xmath192 , @xmath193 , with @xmath194 mev . the ratio of the two density of states @xmath195 and of the fermi velocities @xmath196 . to switch to our notation , @xmath197 .
we start with rbcs results . in the top frame of fig .
[ fig2 ] , we show the dependence of @xmath163 and @xmath198 normalized to @xmath199 as a function of @xmath164 . on comparison with some of the results of fig
.  [ fig1 ] , we note that mgb@xmath5 shows considerable integration of the two bands and also both gaps become zero at @xmath200 which is considerably reduced from the one - band bcs value .    an interesting quantity that could be measured
is the quasiparticle density of states @xmath201 in the presence of a current , which was first examined for one - band s - wave superconductors by fulde@xcite .
such measurements have recently been reported for pb in ref .
@xcite . for each band separately @xmath202\end{aligned}\ ] ] with @xmath203 . in the bottom frame of fig .
[ fig2 ] , we give results for @xmath204 as a function of @xmath205 for various values of @xmath1 as indicated in the figure .
the solid curve is essentially the zero current limit and shows the expected superposition of two bcs quasiparticle density of states with square root singularities at the two gaps .
as the current is increased in the linear region only ( cross - referencing with the upper frame of fig .
[ fig3 ] which shows @xmath0 for mgb@xmath5 ) , the singular structure is smeared and the gap region fills as does the region between the two gaps .
gapless behavior appears at a value of @xmath1 near the first peak in the total current due to the density of states in the small band going gapless .
this is in contrast to what is found in a one - band s - wave superconductor@xcite where the gapless behavior occurs at the critical current .
( 250,200 )    the superfluid current @xmath0 versus @xmath1 for the rbcs calculation with mgb@xmath5 parameters is plotted in the upper frame of fig .
the line types are defined in the same way as for the @xmath0 curves in fig .
we note that this graph is very similar to that shown in the second frame ( righthand side ) of fig .
[ fig1 ] , where it has been discussed in relation to other cases . as in that graph , the peak in the total current seen at low @xmath1 is lower than the one at higher @xmath1 , so that this second peak determines the critical current .
it is of considerable interest to compare these renormalized bcs results with those from full eliashberg calculations based on the @xmath206 spectra given in ref .
very similar values of the parameters can also be found in ref .
such results are presented in the lower frame of fig .
there are several differences that are worthy of note .
the full eliashberg results are at @xmath207 and the rbcs are for @xmath9 .
this difference in temperature is small and is not expected to lead to any significant differences .
examination of the partial currents shows that the dashed curves are almost the same in both plots except that strong coupling effects have shifted the point of zero current to a slightly smaller value of @xmath1 .
this is accompanied by a shift to the left of the position of the peak , a reduction in its height and an increase in the rounding , all small effects .
these differences are seen to be slightly more prominent in the dotted curve .
this shows that rbcs can be used with confidence as a first approximation .
however , from this point on , we will consider only eliashberg results as we have available to us the electron - phonon spectral functions and we wish to consider the full effects of temperature and impurities .
( 250,200 )
at the end of the last section , we showed a comparison between rbcs and eliashberg theory ( fig .
[ fig3 ] ) . given that the electron - phonon spectral functions for mgb@xmath5 are available , it is appropriate to work with the full theory for @xmath2 as a function of temperature and impurity scattering .
we begin with a consideration of temperature effects . in the upper frame of fig .
[ fig4 ] , we show the current @xmath208 normalized by @xmath209 as a function of @xmath164 for nine values of temperature from @xmath207 to 0.9 in steps of 0.1
. the parameters used are those appropriate to mgb@xmath5 with the upper curve reproduced from the bottom frame of fig .
[ fig3 ] . as the temperature
is increased , the peak with larger value of @xmath1 remains the highest and therefore determines the critical current given as the solid curve in the lower frame .
here , @xmath2 has been normalized by @xmath210 so that the slope at @xmath145 is negative one by arrangement for the typical critical current plot of @xmath211 , which is designed to bring out the ginzburg - landau behavior near @xmath6 .
we see a smooth ( concave downward ) increase in @xmath118 as @xmath212 is reduced .
the value of the critical current at @xmath213 is @xmath214 .
this is lower than the one - band bcs value of 0.72 .
this reduction could reflect the two - band nature of mgb@xmath5 as well as strong coupling corrections@xcite due to our use of the full eliashberg equations ( [ eq : del ] ) and ( [ eq : z ] ) to calculate the current given by eq .
( [ eq : js2 ] ) .
however , reference to the lower frame of fig .  [ fig3 ] shows that the second band contributes only about 10% to the total critical current at @xmath9 and this remains true at @xmath6 .
consequently , this curve is very nearly the single band result and deviates from the classical bcs curve almost entirely because of retardation effects .
also shown is the data of kunchur@xcite , plotted as solid dots .
there is good agreement between the data and the solid curve which is quite remarkable as the calculation has no free parameters , in principle , as they are all given in the literature . including impurity scattering has an effect on the temperature dependence of @xmath118 , as we will discuss later in relation to fig .
unfortunately , information on the individual intraband scattering rates in these films for @xmath215 and @xmath216 bands , separately , is not available to us . including some impurity scattering in the @xmath215 band
would futher improve the agreement with experiment .
the other critical current curves in fig .
[ fig4 ] are for different ratios of @xmath174 as indicated in the figure caption .
in addition to allowing for some uncertainty in the reported @xmath131 values in the literature , this also mimics the possibility of currents associated with the c - axis where the @xmath216-band is known to have a much smaller @xmath131 than the @xmath215-band .
likewise , we have not included the details of the full fermi surface average in the expression for the current , however , it is known that the @xmath216-band is quasi - two - dimensional as opposed to the more three - dimensional @xmath215-band .
dahm and scalapino@xcite find significant geometrical corrections associated with such averaging which could also accentuate the difference between the two bands .
this may also be approximately captured by our use of different @xmath131 ratios .
all the curves are normalized to have a slope of negative one at @xmath145 .
as @xmath184 is increased , the normalized critical current at zero temperature increases and can become considerably larger than its one - band bcs value .
there is also a change in the temperature variation of the curves , the dashed curve shows a clear kink around @xmath217 while the long - dashed curve shows concave upward behavior at intermediate temperatures .
this can be understood better when the pattern of behavior shown in fig .
[ fig5 ] is considered .
( 250,200 )    the top frame gives the temperature evolution for the case @xmath218 and the bottom is for a ratio of 4 .
the notation is as for the top frame of fig .
[ fig4 ] and for the same nine reduced temperatures . in the top frame , it is the peak at lower @xmath1 which determines the critical current at @xmath213 while at higher @xmath212 , it is the peak with higher @xmath1 .
this crossover occurs between @xmath219 and @xmath220 and leads to the kink seen in the dashed curve of fig .
[ fig4 ] ( lower frame ) .
if we look at the individual contributions from each band ( not shown here ) to the total current , we find that at @xmath9 it is the second band which makes the dominant contribution because of its large value of @xmath131 but near @xmath6 this contribution has decayed sufficiently that it only makes a 30% contribution to the total .    in the lower frame of fig .
[ fig5 ] , the peak at smaller value of @xmath1 has greatly increased in magnitude over its value in the top frame .
it has also moved to a lower value of @xmath1 as expected from fig .  1 .
however , now there is no crossover from smaller to larger @xmath1 peak and hence no kink .
nevertheless , the dotted curve is very different from a one - band bcs variation .
we can understand better how this can arise by once more looking at the partial contribution to the critical current from each band .
at @xmath9 it is the second band which dominates while at @xmath6 it accounts only for 2/3 of the total with the remaining 1/3 contribution from band one .
this difference in admixture of the two bands with temperature is seen to be enough to cause an upward curvature in the dotted curve at intermediate temperatures as well as to increase the value of the normalized critical current at @xmath9 to a value considerably larger than the bcs value of 0.72 .
increasing @xmath184 further as in the long - dashed curve of the lower frame of fig .
[ fig4 ] increases the @xmath9 value even further as the peak at lower @xmath1 value becomes even more prominent .
details of the temperature variation of the critical current clearly depend on details of the microscopic parameters of the two bands involved , such as values of @xmath221 .
for example , it is possible to shift the kink in @xmath118 to higher values of reduced temperature and to make it sharper by considering other parameter sets .
( 250,200 )    to illustrate this , in fig .
[ fig6 ] we show results for a lorentzian @xmath222 spectrum defined in detail in our previous paper@xcite . in this model , @xmath159
and we take @xmath223 to decrease the gap anisotropy as compared with mgb@xmath5 .
we have also taken @xmath173 ( nearly decoupled bands ) .
the top frame gives the current @xmath224 for the case of @xmath218 , for the nine values of reduced temperature previously chosen . with this @xmath131 ratio ,
the peak with lower value of @xmath1 remains highest at small @xmath212 and the crossover of the maximum to the other peak occurs at a reduced temperature between 0.6 and 0.7 .
this manifests itself as a sharp kink in the solid curve for @xmath2 ( bottom frame ) at this same temperature .
also shown in this same frame are results for the critical current when @xmath225 ( dotted ) and @xmath226 ( dot - dashed ) .
for the larger value of @xmath184 , the normalized current at @xmath9 is increased over the @xmath218 case as expected and for @xmath227 it is reduced but remains above bcs . in this latter case , the two peaks in the current versus @xmath1 are closer together than they are in the top frame of fig .  [ fig6 ] and seem to almost merge before the crossover from lower to higher @xmath1 peak occurs , which means that the kink in the dotted curve can hardly be seen .
it is clear that a rich pattern of behavior exists for the current versus @xmath1 at different reduced @xmath212 in a two - band model .
this translates into quite distinct temperature dependence for @xmath118 as compared with the one - band bcs canonical behavior .
some addition insight into this complicated behavior can be obtained from a consideration of simplified but analytic bcs results in the limiting case of two decoupled , well - separated bands , i.e. , for small values of the gap anisotropy @xmath228 . under this assumption ,
@xmath144 , @xmath229 , @xmath230 , and @xmath231 so that the critical current near @xmath6 given by formula ( [ eq : bcs14 ] ) reduces to @xmath232 and is determined entirely by band one , as previously commented upon .
as the temperature is lowered , however , two possible circumstances can arise .
the critical current can remain determined by the same band ( number one ) in which case we can write approximately @xmath233 a second possibility is that a crossover to the peak at lower @xmath1 value occurs and we have approximately @xmath234 for the first case , @xmath235 which is the one - band bcs result where we have normalized the critical current in such a way that @xmath211 has a slope of negative one at @xmath145 .
this conforms with what we have done in the lower frames of figs .
[ fig4 ] and [ fig6 ] .
for the second case , by assumption , there is a crossover and @xmath236^{2/3}. \label{eq : bcse}\ ] ] this simplified approximate formula shows clearly that the anisotropy @xmath237 on its own reduces the value of @xmath238 below its universal bcs value of 0.72 while a large value of @xmath239 relative to @xmath240 or @xmath184 relative to @xmath241 increases it and that this reduced dimensionless quantity is no longer universal .
it can depend on the band structure parameters as well as on the gap anisotropy .
next we consider the effect of impurities on the current .
it is instructive to begin by expanding eq .
( [ eq : js2 ] ) for @xmath224 to lowest order in @xmath1 .
doing so leads to the well - known and physically expected result that for small superfluid velocity @xmath242 where @xmath243 is the superfluid density of the @xmath33th band given for a strong coupling superconductor by @xmath244^{3/2}}\biggr ] \label{eq : bcsg}\ ] ] in bcs at @xmath9 , this becomes @xmath245 and @xmath246 for band 1 and band 2 , respectively .
the slope of @xmath247 as a function of @xmath8 gives @xmath248 with eq .
( [ eq : bcsg ] ) valid for any temperature and any impurity content .
the impurities enter only through the eliashberg gap equations ( [ eq : del ] ) and ( [ eq : z ] ) .
as the temperature is increased , the superfluid density is reduced and so the slope of @xmath224 seen in figs .
[ fig4 ] to [ fig6 ] is correspondingly reduced .
impurities also reduce the superfluid density . in bcs , the known result for each of the two bands separately , to lowest order , is @xmath249 with @xmath250 and @xmath251 , the coherence length and impurity mean free path in the @xmath33th band , respectively .
this applies to intraband scattering which reduces the initial slope of the partial current coming from each band .
( 250,200 )    results are presented in fig .
[ fig7 ] , where we also consider the interband case @xmath252 . in the top frame ,
we show @xmath0 versus @xmath1 for the reduced temperature @xmath253 in the lorentzian model used in fig .
the solid curve applies to the pure case and is for comparison .
we show only three cases , one where @xmath254 ( dotted curve ) , @xmath255 ( dashed curve ) , and @xmath256 , where @xmath257 is the @xmath6 without impurities ( in the case of @xmath258 and @xmath259 , @xmath6 will be unaffected , but not so for interband impurities ) .
intraband impurities reduce only the partial contribution to the current coming from that band .
it reduces its value on the vertical axis and extends the x - axis to higher values of @xmath1 leaving the second band largely unaffected . of course
, since both bands contribute significantly to the total current in the region of the peak corresponding to the smaller value of @xmath1 , its magnitude is reduced even in the dotted curve . on the other hand , for the dashed curve
, there is little change in the region of the second peak which is dominated by the first band .
by contrast interband impurities modify both bands in a similar fashion .
the corresponding critical current as a function of @xmath260 is shown in the lower frame .
the solid line is the pure case of fig .
as expected , a finite value of @xmath258 has reduced the critical current below its value in the pure case at all temperatures , however , when normalized to a slope of magnitude one near @xmath6 , the effect pushes the curve up at low @xmath15 .
a kink remains and the crossover temperature is only slightly shifted toward higher values . for the dashed curve ,
the significant modifications arise only below the kink temperature .
by contrast the main effect of the interband impurities is to wash out the kink structure leaving the rest of the curve largely unaffected ( of course the @xmath6 would show a small shift downwards relative to @xmath257 , not shown here ) .
in conventional one - band s - wave superconductors , the energy gap @xmath57 remains unaffected by the flow of a supercurrent up to a critical value of the superfluid velocity @xmath7 ( ie .
, @xmath261 ) at which point pairbreaking becomes possible . in this region ,
the current @xmath0 is linear in @xmath8 with proportionality constant equal to the superfluid density . beyond this critical value of @xmath8 , the rapid reduction in the gap in addition to a quasiparticle backflow
current leads to @xmath262 at @xmath263 . for two nearly decoupled bands
, we would expect that the total current to exhibit a two - peak structure as a function of @xmath8 .
a first one near the critical value of the superfluid velocity for band two and the other near the larger critical value of @xmath8 corresponding to the first band alone .
the first peak has in general a significant contribution from both bands ( unless @xmath264 ) while the second peak could be due mainly to the first band .
the relative admixture depends on microscopic parameters such as density of states @xmath265 and fermi velocities @xmath32 , as well as gap anisotropy .
when @xmath184 is made large as compared to @xmath241 , the first peak moves to lower values of @xmath8 and can also become larger than the second at low temperature .
as the temperature is increased the superfluid density in each band is reduced but that in the second band is more strongly affected because of its smaller gap value . consequently , since the critical current is determined by the height of the highest peak , one can have a crossover from the first to second peak as @xmath15 is raised from @xmath266 to @xmath6 and this will reflect itself as a kink in the temperature dependence of the critical current . such a feature is also found in the dirty limit calculations presented in ref .
@xcite for the critical current along the c - axis .
when there is a significant interband coupling or the two gaps are close in value ( isotropic case ) the two peaks in the total current cease to be independent and a rich pattern of behavior emerges for the temperature and impurity dependence of the two peaks .
in particular , as the temperature is raised , the first peak can remain the determining peak for @xmath2 up to @xmath6 or it can become the smaller of the two at a crossover temperature which can be increased or decreased through judicious choices of parameters .
another possibility is that the two peaks merge into one or that the second peak is the dominant one at all temperatures .
application to the specific case mgb@xmath5 , which involves no adjustable parameters , gives good agreement with the in - plane data of kunchur@xcite .
if the resulting critical current @xmath211 is normalized to have slope of negative one as a function of @xmath58 near @xmath145 , the corresponding value of the critical current at @xmath9 will differ from its classic one - band bcs value of 0.72 .
strong coupling effects resulting from our use of the full eliashberg equations , with appropriate electron - phonon spectral densities describing the electron - phonon interaction , are known@xcite to always reduce the normalized value of @xmath267 while two - band effects can reduce or increase it depending on the detailed values of the microscopic parameters involved .
in particular , a large value of @xmath184 relative to @xmath241 can lead to values that can be even larger than one .
this is expected for the case of current orientated along the c - axis for the specific parameters representative of mgb@xmath5 . in this case , because the @xmath216 band ( band one in our notation ) is nearly two - dimensional , the fermi velocity of the @xmath215 band in the c - direction will be much greater than that for the @xmath216 band .
the superfluid density in each subband , which determines the initial slope of the linear in @xmath8 regime for the supercurrent , is reduced with the increase in intraband impurity scattering .
this fact can be used to manipulate the position along the @xmath8 axis and size of the two peaks in the total current and consequently the temperature dependence of the resulting critical current .
introducing impurities in the first band will reduce the critical current at all temperatures since this band always makes a contribution to the total current . on the other hand ,
if they are introduced only in the second band , this will leave the critical current largely unaltered around @xmath113 .
interband impurity scattering affects both bands although not necessarily in exactly the same way and leads to the smoothing out of the kink in @xmath268 versus @xmath212 when it exists .    in conclusion
, we predict a complex pattern of behavior for the total current as a function of @xmath8 , temperature and impurity content which can be used to restrict further the values of the microscopic parameters involved in two - band superconductors , such as mgb@xmath5 .
ejn acknowledges funding from nserc , the government of ontario ( prea ) , and the university of guelph .
jpc acknowledges support from nserc and the ciar .
this research was supported in part by the national science foundation under grant no .
phy99 - 07949 .