many years ago , superconductivity was discovered in heavy - fermion compounds @xcite .
it was suggested @xcite that the superconductivity was due to collective electronic fluctuations and not due to electron - phonon interactions .
transport properties in the superconducting state were analyzed @xcite to show that superconductivity was in the d - wave symmetry .
it was also suggested that the d - wave symmetry is promoted by antiferromagnetic fluctuations @xcite with long enough correlation lengths .
this promotes scattering of fermions near the fermi - surface predominantly through angles around @xmath2 , which is essential for superconducting instability in the `` d - wave '' channel for a suitable fermi - surface @xcite .
the idea of long enough afm correlation lengths as essential for this mechanism is supported by the fact that in heavy - fermions , superconductivity occurs generally in the regime near the afm quantum critical point where the correlation lengths are long but the competing afm phase has lower condensation energy .    at the same time
, random phase approximation on the hubbard model was used to calculate the spin - fluctuation spectra and to suggest that d - wave superconductivity is promoted by such fluctuations @xcite .
the properties of the hubbard model have proven controversial in more elaborate calculations ; there are calculations which suggest that the ratio of the transition temperature @xmath0 to the typical electronic kinetic energy parameters @xmath3 is more than @xmath4 @xcite to less than @xmath5 @xcite .
since heavy - fermion properties require kondo effect of the f - orbital local moments and their magnetic interactions using the wide - band electrons , a multi - orbital model is obviously required @xcite .
the hubbard model was proposed as a sufficient model for the cuprate compounds @xcite .
but the discovery in under - doped cuprates @xcite of the predicted time - reversal breaking order parameter @xcite on the basis of a multi - orbital model raises doubts on the validity of the hubbard model for the cuprates . for pnictides ,
generalization of the hubbard model to multi - orbital situations and inclusions of hund s rule couplings appears essential .
we have a more modest goal in this paper than calculating spin - fluctuations from microscopic theory and using it to calculate properties of the superconductor . in recent years inelastic neutron scattering in the heavy fermion compounds cecu@xmath6si@xmath6 @xcite and ceirin@xmath7 @xcite have provided details of the afm fluctuation spectrum in the normal state . the primary aim of this paper is to estimate the superconducting transition temperature using the parameters provided by the experiments in these compounds .
to do so , we solve the eliashberg equations for d - wave superconductivity using a phenomenological afm spectral function with which the experimental data is in good accord .
the use of the eliashberg equations for quantitative calculations may be open to question because the migdal expansion parameter , which is of @xmath4 for the electron - phonon problem is of @xmath8 for such compounds if one assumes that the scale of the afm fluctuations extends to the order of the electronic bandwidth .
however , when the afm correlation length @xmath9 is large compared to the lattice constant @xmath10 or @xmath11 , the scale of the afm fluctuations is reduced correspondingly to @xmath12 . but
in the limit of large correlation lengths , new questions arise @xcite which are not important in the electron - phonon problem .
the most prominent among them are the role of inelastic scattering in depressing @xmath0 on the one hand @xcite , and the fact that the bcs type coupling constant @xmath13 appears to diverge when the characteristic fluctuation frequency @xmath14 and the bcs prefactor appears to go to @xmath15 .
an answer to these questions and various considerations which determine @xmath0 from afm interactions is possible from the numerical solution of the eliashberg equations .
we find that it is reasonable to conclude from a comparison of the calculated @xmath0 with experiments that afm fluctuations are responsible for d - wave superconductivity in the heavy fermion compounds . very importantly , with similar parameters we calculate the _ measured _
coefficient of the anomalous @xmath1 contribution to the resistivity in these compounds .
a claim to quantitative accuracy on both these quantities can however be made only to factors of @xmath16 .
we note here that if one adopts that the dimensionless measure @xmath17 for how high is the electronic fluctuation induced superconducting , the heavy fermions may be said to do very well indeed .
for example , in many cases , including the compounds studies here , this ratio is @xmath4 , similar to that of the cuprates .    following the proposals that afm fluctuations may also promote superconductivity in the cuprate compounds @xcite ,
there have been many discussions of the mechanism and many calculations based on the eliashberg equations .
a partial list includes the following @xcite .
the most complete of these calculations appear to us to be those carried out by monthoux and lonzarich ( ml ) @xcite , both for 2 and 3 dimensional models .
we present below calculations for the 2 dimensional square lattice model with a phenomenological spin - fluctuation spectrum , whose results are no different from those of ml for the range of parameters examined that are common .
a difference in the calculations is that we vary the parameters in the two dimensional model so that `` nesting '' at the afm wave - vector quantitatively changes .
the amount of nesting does have a significant effect on the results .
more important is that now that the afm fluctuation spectra is available , we can use the experimental parameters to test the ideas quantitatively .
we also discuss how to put limits on the parameters used based on sum - rule for the fluctuation spectra and show that they are inter - related .
results for the range of physical parameters that we find relevant for the heavy fermions is not available in the published results of ml .
this has bearing also on general conditions to determine the extent to which afm fluctuations give significant @xmath0 for relevant parameters in other compounds .
this paper is organized as follows : we present in sec .
( ii ) the models for fermi - surface and for the spin - fluctuations which we have investigated using the linearized eliashberg equations .
we discuss there the change of effective coupling constants with the afm correlation length using sum - rules so that the results for numerical solutions of the eliashberg equations presented later are presaged .
we present the results of the calculations in sec .
( iii ) and discuss the important conclusions immediately after the description of the models .
we also present , in an appendix , the explicit derivation of the coefficient of the @xmath18 resistivity from the measured form of the afm critical fluctuation spectra .
this is used in the text to estimate independently the value of a coefficient @xmath13 , which is important for the calculation of @xmath0 .
we give the parameters that have been deduced by inelastic neutron scattering for the heavy fermion compounds cecu@xmath6si@xmath6 @xcite and ceirin@xmath7 @xcite and compare the measured @xmath0 with the calculations .
we should emphasize that such a comparison is meant to be only illustrative of the physical principles involved ; no detailed quantitative agreement is to be expected , especially given that the electronic structure of these compounds is far more complicated than assumed in the models studied . however , enough details can be provided so that one can conclude that the idea of afm fluctuations near the quantum critical point in these compounds as the source of d - wave superconductivity is well supported .
for example using measured properties , different levels of assumed nesting in the band - structure need a coupling constant @xmath13 between 1.5 and 3 to get the measured @xmath0 . in this range of @xmath13 and for the measured afm correlation length
@xmath0 is close to being linear in @xmath13 .
this range of values is compared with the value of @xmath19 needed to get the measured coefficient of the @xmath18 resistivity , which is relatively insensitive to nesting .
one can assess the results from the fact that in the range of @xmath13 deduced , @xmath20 is found to be approximately linear in @xmath13 .
in our calculations we will consider two types of fermi - surfaces , a free electron fermi - surface and the others given by the tight binding spectrum in a two dimensional square lattice with nearest neighbor and next nearest neighbor hopping @xmath3 and @xmath21 respectively : @xmath22 the fermi - surface with the tight binding spectrum are shown in fig .
[ fig:01 ] for four values of the next nearest hopping @xmath23 and the afm wave - vector .
the nesting in the model changes as @xmath21 increases .
we will show detailed result for three fermi - surfaces , the free electron fermi - surface , the fermi surface ( fs1 ) with tight binding spectrum with @xmath24 and the fermi surface ( fs2 ) with tight binding spectrum with @xmath25 .
of the four fermi - surfaces shown in fig .
[ fig:01 ] , the one with @xmath24 has the worst nesting and the one with @xmath25 has the best nesting .
[ fig:02 ] shows the circular fermi - surface , fs1 , fs2 and the corresponding afm wave - vectors .
we will discuss using the results of ml together with ours , that if properly normalized density of states and fluctuation spectra are used , two dimensional and three dimensional models give similar results for @xmath0 provided one adjusts the ratio of the region of fermi - surface nesting to the total fermi - surface .
this is in general is always lower in three than in two dimensions .
it is also important to note , as discovered long ago @xcite , @xcite for the case of s - wave superconductors that @xmath20 is a rather gross quantity which depends to a very good approximation on the average density of states near the chemical potential only and not on details such as the number of fermi - surface sheets and shapes . for d - wave superconductors
, we we show below , it is important to also include effects of nesting of the fermi - surface near the afm wave - vectors .      the dynamical spin - susceptibility for itinerant fermions with afm correlations may be usefully divided into two parts : the normal contribution of non - interacting fermions @xmath26 and the part @xmath27 affected by afm correlations .
the two together obey the total magnetic moment sum - rule .
the non - interacting susceptibility can play only an insignificant role in promoting superconductivity @xcite and should be ignored .
the two contributions to the susceptibility can be distinguished by their momentum dependence .
the characteristic momentum dependence of the non - interacting spin fluctuations is on a scale of @xmath28 while that of the afm spin - fluctuations is much shorter .
a partial sum - rule on @xmath27 in terms of the ordered moment in the afm phase can be used to relate the integrated fluctuations to the afm correlation length .
double counting by using the sum - rule on the total susceptibility for the fluctuations and yet having ( free ) fermions interacting with such spin - fluctuations is incorrect as it over - counts the total degrees of freedom .
such considerations have usually been blithely ignored in most of the previous phenomenological work on this problem .
the fermions interact with spin - fluctuations with a phenomenological action @xmath29 @xmath30 will be chosen to have dimensions of inverse of energy ( after subsuming a factor of @xmath31 in its definition ) .
so @xmath32 is a coupling function of dimension of energy .
@xmath32 for heavy fermions is the exchange energy between the conduction electrons and the @xmath33- local moments .
its meaning for d - band problems is more ambiguous , and may be best inferred from independent experiments , for example the resistivity above @xmath20 .
a suitable phenomenological form for the dynamical spin - fluctuations due to afm correlations , with which experimental results @xcite can be fitted , is @xmath34 .
\label{chi}\end{aligned}\ ] ] where @xmath35 is the damping rate of the fluctuations , @xmath36 is the antiferromagnetic vector .
the correlation length @xmath9 is related to the deviation from the quantum critical point ( qcp ) by variation in pressure , doping , magnetic field , etc .
as well as by temperature .
@xmath35 , @xmath37 and @xmath9 may all be determined from experiments .
the temperature dependence of @xmath9 has been studied by renormalization group ( rg ) @xcite and by the self - consistent renormalization ( scr ) methods @xcite . @xmath38 near afm qcp in the 3 dimensions and dynamical critical exponent @xmath39 .
scr derives using the same dynamical critical exponent that near the magnetic qcp , @xmath40 .    in our calculations in two dimensions
@xmath9 will be assumed to be of the form @xmath41 where @xmath42 is the asymptotic @xmath43 value of the correlation length .
one of our results is that the temperature dependence of @xmath9 is of insignificant consequence in determining @xmath0 .
the linearized elisahberg equations give that the kernel for cooper pair coupling in the d - wave channel in a square lattice is proportional to the projection of @xmath44 to @xmath45 . in spin - fluctuations theories
, @xmath46 has only a smooth momentum dependence .
so , cooper - pair coupling prominently depends depends only on ( i ) the momentum dependence of @xmath47 determined by the correlation length @xmath9 ( ii ) the integrated weight in the momentum dependent part and ( iii ) the energy scale of the momentum dependent fluctuations .
the first is qualitatively obvious from the fact that a q - independent spin - fluctuation contributes zero to the cooper channel in the d - wave channel .
it is not possible to make quantitative statements on these effects without detailed calculations because the results also depend on the nesting in the band - structure near the afm @xmath37 .
we will show that the three ingredients in @xmath47 are not mutually independent .    to gain physical insight ,
the effect of @xmath9 on the integrated spectral weight may be discussed before detailed calculations through the the partial sum - rule on @xmath48 , which determines the effective coupling constant for superconductivity : @xmath49   \right ] .
\label{tmsr}\end{aligned}\ ] ] with the assumed lorentzian form , it is necessary to introduce an upper cut - off @xmath50 in the frequencies @xmath51 up to which the fluctuations extend .
actually , spin fluctuation are actually quite suppressed for @xmath52 and we can simply use @xmath53 in calculations of eliashberg equations .
it is important to take into account that there are four equivalent afm - vector for the two dimensional problem in the paramagnetic regime of the model , however strongly fluctuating it may be .
this has been taken into account in the sum - rule by multiplying the measured im@xmath54 by 4 . for @xmath55 ,
the number of equivalent afm vectors is larger and a correspondingly larger multiplicative factor should be used .    in the regime of very long correlation lengths , @xmath56 , i.e. close to the quantum - critical point
, the sum rule simply gives @xmath57 @xmath58 may to a first approximation be estimated from the ordered moment @xmath59 in nearby afm phase but more properly from integration of the relevant momentum and frequency range of the measured fluctuations in absolute units using polarized neutrons .
[ fig:00 ] shows the @xmath60-dependence of @xmath61 for @xmath62=1.0 , 2.0 and 10.0 .
let us now consider the sum - rule in the opposite limit , that the correlation length is small compared to the lattice constant , i.e. the system is very far from the quantum critical point
. then @xmath63\ ] ]    as already shown by ml and further elaborated below , for a given band - structure , the results of the eliashberg calculations for @xmath64 may be parametrized in terms of a dimensionless `` bare '' coupling constant @xmath13 and a correlation length @xmath9 , @xmath65 @xmath66 may be determined in terms of @xmath58 and therefore ( approximately ) to the ordered moment through the sum - rule .
we may define an effective coupling constant @xmath67 to incorporate the effect of the correlation length . using that the sum - rule becomes the total moment sum - rule in the limit of infinite correlation length and the maximum possible ordered moment , i.e. that of the afm insulator ( ignoring the zero - point effects )
, one concludes that in the limit of very large correlation lengths @xmath68 where @xmath33 is the fraction of the maximum possible ordered moment . from fig .
( [ fig:00 ] and from the detailed calculations presented in the next section , one deduces that the limit for @xmath69 is reached for @xmath70 , below which there is an exponential fall off of @xmath64 . for smaller correlation - lengths , fig .
( [ fig:00 ] ) shows that the @xmath67 for @xmath0 decreases with decreasing correlation length .    in the work of ml ,
@xmath13 values from about 5 to about 50 are used in the calculations with varying correlation lengths . actually , one obtains for the considerations of the sum - rules above that for spin-(1/2 ) problems , even the coupling constant @xmath69 is only of @xmath8 , because @xmath71 is of @xmath8 and so is the upper limit on the ratio @xmath62 .
an independent estimate of @xmath67 may be obtained from the normal state properties , for example the coefficient of the temperature dependence of the resistivity of non - fermi - liquid form in the quantum critical region .
again only @xmath13 of @xmath8 will be found consistent .
it is also important to note that these bcs type coupling constants do not carry information on the retardation effects due to the frequency dependence of the interaction ; these as well as the effects of inelastic scattering which are particularly important for anisotropic superconductors are properly treated through the numerical solution of the eliashberg equations .
the difference from electron - phonon induced superconductivity where a single parameter @xmath13 need by introduced @xcite should also be noted .
it should also be pointed out that in some heavy fermions , the quantum - critical fluctuations do not have the functional form given by the simple rg or scr approximations as above , but displays `` local criticality '' @xcite as suggested for the cuprates @xcite . in this paper , we only consider fluctuations which are well specified by the form given above .
the temperature dependence of the resistivity near the quantum - critical points has been derived several times @xcite . here , we rederive it paying special attention to the coefficient in front of the anomalous temperature dependent part .
an expression of the resistivity @xmath72 in the antiferromagnetic quantum critical region suitable for heavy fermions may be derived with the following formula derived from the boltzmann equation . @xmath73 where the integration is taken over the fermi - surface , and @xmath74 the fermi velocity .
this assumes that the _ actual _ electronic structure near the chemical potential is sufficiently complicated that in the temperature region of interest , vertex corrections which lead to emphasis on large momentum scattering for resistivity are unimportant . in that case the scattering rate which determines the resistivity is the same as the single - particle scattering rate averaged over the fermi - surface .
this is true in a multi - sheeted fermi - surface and is suitable for heavy fermions .
this is similar to the case of transition metals where the resistivity from electron - phonon scattering at low temperatures is @xmath75 in contrast to the nearly free - electron metals where it is @xmath76 . for weakly anisotropic single band scattering , as in the cuprates ,
the resistivity for large afm correlation lengths is close to the fermi - liquid temperature dependence although near the hot spots the scattering rate is nearly @xmath77 @xcite .    equation ( [ aq1 ] ) can also be expressed as follows .
@xmath78 where @xmath79 means the average over the fermi - surface , @xmath80 the renormalized effective mass .    here
, @xmath81 can be derived from the imaginary part of the self - energy . @xmath82 where the self - energy due to the antiferromagnetic quantum fluctuations is given as follows .
@xmath83 the result for the resistivity in the limit @xmath84 is derived in an appendix a. by explicitly calculating the self - energy given by eq .
( [ aq4 ] ) is @xmath85
the superconducting transition temperature is given by the linearized version of the eliashberg equations for the normal self - energy @xmath86 and the anomalous or pairing self - energy @xmath87 .
@xmath88 i \omega_n & = - \int_{\rm fs } \frac { { \rm d}^{d } s_{\vec{p } } } { ( 2 \pi)^d v_{\vec{p } } } \pi t \sum_{\omega_m } i\ { \rm sgn}(\omega_m ) g^2 \chi ( \vec{k } - \vec{p},i \omega_n - i\omega_m),\\      \label{leliash2 }      w ( \theta_{\vec k } , i \omega_n ) & = - \int_{\rm fs } \frac { { \rm d}^{d } s_{\vec{p } } } { ( 2 \pi)^d v_{\vec{p } } } \pi t \sum_{\omega_m } \frac { w ( \theta_{\vec p } , i \omega_m ) } { \vert \omega_m z ( \theta_{\vec{p } } , i \omega_m ) \vert } g^2 \chi ( \vec{k } - \vec{p},i \omega_n - i\omega_m).\end{aligned}\ ] ] here @xmath89 are the matsubara frequencies ; @xmath32 is a momentum - independent coupling matrix element , which has already been defined , @xmath90 is an angle parameterizing the fermi surface , @xmath91 is the density of states at angle @xmath90 .
the @xmath92-integral is over the fermi surface , @xmath93 is the unrenormalized velocity .
our principal general results for @xmath0 on the basis of solution of the linearized eliashberg equations in terms of @xmath13 and the parameters in @xmath48 are given in this section .
the numerical evaluation is done by first simplifying the eliashberg equations ( [ leliash1 ] ) and ( [ leliash2 ] ) , as far as possible analytically . the final expressions for the numerical evaluation , both for the circular fermi - surface and
the tight - binding fermi - surfaces are given in appendix b.    figure .
[ fig:03a ] show the @xmath94-dependences of @xmath64 for the circular fermi - surface on the bare coupling constant @xmath13 in the large correlation length limit on the right . for small bare coupling @xmath13 ,
the latter does have the bcs form while for @xmath95 , the dependence is approximately linear .
consistent with earlier discussions @xcite , @xmath64 shows a very shallow peak at around @xmath96 .
@xmath64 shows a drop - off as the correlation length decreases , while it also shows moderate decreases as the correlation length increases .
ml pointed out that this moderate decrease is caused by the rapid diverges of @xmath97 as the correlation length increases .
we note also that the quantum - classical crossover correction to the correlation length proportional to the factor @xmath98 in eqs .
( [ corr2 ] ) has a negligible effect on @xmath0 .
this will not be considered in any further calculations .    the principal message from fig . ( [ fig:03a ] )
is that the infinite correlation length result for @xmath0 is well obeyed up to @xmath99 with a very sharp fall off thereafter which will be seen later to be exponential . for large @xmath100 ,
no bcs type approximation for @xmath0 is valid .
the limit of very large correlation length is equivalent to the effective frequency of fluctuations @xmath14 , as may be seen from eq.([chi ] ) .
if we use the mcmillan @xcite type approximation , in which @xmath101 , the inverse of the average squared frequency of fluctuations , we get a divergent coupling .
( [ fig:03a ] ) gives a finite limit to @xmath64 , which depends on the bare coupling constant @xmath13 .
one may understand this result from the calculations of allen and dynes @xcite , deduced for s - wave eliashberg equation , that in the limit of a diverging coupling constant @xmath102 , where @xmath103 may be taken approximately to be the square - root of @xmath104 . -    next we show in fig .
( [ fig:04a ] ) the @xmath94-dependence of @xmath64 for the four fermi surface shown in fig .
( [ fig:01 ] ) . for @xmath24 , the worst nesting case ,
@xmath64 is similar to that for the circular fermi - surface .
improving the nesting condition increases , @xmath64 .
we show in ( [ fig:05a ] ) @xmath20 as a function of @xmath13 for the worst nesting fermi - surface and the best nested fermi - surface of those in fig .
( [ fig:01 ] ) . a increase of @xmath16 in @xmath0 for the similar values of @xmath13 is discerned from the worst to the best nesting conditions .    .
[ fig:05a ]    ml have also presented detailed results for calculations on a 3 d electronic dispersion with the symmetry of a cubic lattice .
they remark that other parameters being the same 2 dispersion gives higher @xmath0 than a three dimensional dispersion .
based on our results for changes in @xmath0 in the 2 d problem , we conclude that this is because of the much better nesting that is obtainable in model 2 d systems compared to the 3d systems for a given @xmath37 which spans the fermi - surface in some ( usually symmetry ) direction . in fact , we can place our 2 d - results for the very weakly nested fermi - surface over the ml results for the 3d fermi - surface and find for other parameters the same that the systematics of the results for @xmath0 as well as its value is very similar .
in this section , we compare the estimates of @xmath0 from the calculations with the experimental result in cecu@xmath6si@xmath6 and ceirin@xmath7 . for convenience , we show the measured intensity @xcite proportional to the dynamic structure factor @xmath105 in fig .
( [ fig : s(q ) ] ) for @xmath106 near @xmath107 .    .
[ fig : s(q ) ]    although the magnetic fluctuation spectrum found through inelastic scattering in cecu@xmath6si@xmath6 is well represented by the form of eq .
( [ chi ] ) , the electronic structure is far more complicated than assumed here or in the 3 d calculations of ml .
we have seen that @xmath0 , especially in the limit of large magnetic correlation lengths depends only on gross parameters like @xmath13 and secondarily on the amount of nesting .
the comparison can only be very limited and can only give insight into the orders of magnitudes expected and to the physics involved .
[ [ cecu_2si_2 ] ] cecu@xmath6si@xmath6 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    @xcite + to fit the phenomenological susceptibility to these results , the parameters take the following values : + _ afm wave vector _ : @xmath108= ( 0.22 , 0.22 , 1.46 ) + _ ordered moment in the afm phase _ :
@xmath109/ce _ correlation length : _
@xmath110 25  ; + _ characteristic energy scales _ : + @xmath111 1.5 [ mev ] ; + _ density of states _ : from the measured uniform @xmath112 paramagnetic susceptibility in the normal state , one deduces the @xmath113 ( 1/7 ) [ mev]@xmath114unit - cell .
+ _ afm spin - fluctuations parameter _ : @xmath66 : + the experimental results for @xmath115 in ref .
( ) are parametrized in terms of three quantities @xmath116 and @xmath117 . the correlation length @xmath9 in ref .
( ) is the same quantity used by us . for clarity
we give here the relation of the other two parameters to the parameters used by us . the conversion from the quantity @xmath118 , which we will call @xmath119 to our @xmath120
is obtained by equating the integral over all @xmath121 of equation s4 in ref .
( ) to the integral of the same physical quantity given in eq .
( [ tmsr ] ) . in the limit of @xmath122 ,
one gets @xmath123 .
the experimental result is @xmath124 .
this then gives @xmath125 .
+ the quantity @xmath117 is related to @xmath126 by @xmath127 .
+ _ transition temperature _ : @xmath128 0.6[k ] + cecu@xmath6si@xmath6 has a very anisotropic fermi - surface with very little dispersion along the tetragonal axis .
the fermi - surface in the plane is very complicated but we assume that just as in s - wave superconductivity @xcite , @xmath0 depends only on the average density of states at the fermi - surface , supplemented by knowledge of nesting of the fermi - surface near @xmath129 . among other things , our results below may be taken to be test of this assumption .
[ [ ceirin_5 ] ] ceirin@xmath7 + + + + + + + + + + + + + + + + + + + + + +    @xcite + although long - range magnetic order competing with superconductivity in ceirin@xmath7 has not been accessed in this stoichiometric compound , there are strong experimental results indicating that the compound lies in the vicinity to an afm quantum - critical point .
the resistivity of this material exhibits a non - fermi liquid behavior similar to that observed in cecoin@xmath7 , which is known to lie in the vicinity of an afm quantum - critical point which has been accessed by doping the compound .
moreover , the nuclear spin relaxation rate of ceirin@xmath7 is also similar to that of cecoin@xmath7 . the dynamical susceptibility has been recently deduced by nmr experiments @xcite in agreement with this conclusion .
+ to fit the susceptibility to these experimental result , the parameters take the following values . + _ afm wave vector _ : @xmath108= ( 0.5 , 0.5 , 0.5 ) + the chosen @xmath108 and ordered moment are taken to be that of the related compound cecoin@xmath7 @xcite + _ ordered moment in the afm phase _ : @xmath130 + _ correlation length _
+ @xmath131 at @xmath132[k]@xcite _ characteristic energy scales of afm : _ @xmath133 1.5[mev ] ; + _ transition temperature _ : @xmath0 = 0.4[k ]     the experimental results show that both cecu@xmath6si@xmath6 and ceirin@xmath7 lie not far from the asymptotic large correlation length limit and that their @xmath64 are both about 0.03 . for a circular fermi - surface , and
using the measure value of @xmath100 in the former , we may refer then simply to fig .
( [ fig:03a ] ) and find that @xmath134 gives the right value of @xmath0 . for the best nested fermi - surface , however , a @xmath135 is sufficient as shown in fig .
( [ fig:05a ] ) .
we may now try to estimate @xmath13 to see if these values are reasonable .
we do this in two different ways . to utilize the neutron scattering results for this purpose , we need to know @xmath32 besides the directly measured properties listed above .
the renormalizations in the heavy fermion problem are such that near the critical point the afm interaction between magnetic moments is of the same order as the heavy fermion bandwidth .
then @xmath32 is of the order of the effective fermi - energy , i.e @xmath136 .
then we may use the experimental values of @xmath137 and @xmath138 deduced from experiments above in eq.([lambda ] ) to get @xmath139 .
the above manner of estimation has forced us to guess the value of @xmath32 .
we can estimate the value of @xmath13 much better and independently from the non - fermi - liquid resistivity proportional to @xmath18 observed in the quantum critical regime of cecu@xmath6si@xmath6 , whose coefficient is proportional to @xmath13 .
the resistivity @xmath140 in the quantum - critical region for @xmath84 is given by eq.([rho ] ) .
using the values of cecu@xmath6si@xmath6 mentioned above , the resistivity is estimated @xmath141 [ @xmath142 from eq .
( [ rho ] ) .
the non - fermi liquid resistivity observed in cecu@xmath6si@xmath6 takes the form : @xmath143 @xcite where @xmath144 @xcite . from the comparison of the coefficient of @xmath18-term in the resistivity between theoretical and experimental results ,
@xmath13 is estimated as @xmath145 .
this should be considered an important evidence for the rather obvious idea that fluctuations that determine the normal state scattering also determine @xmath0 , and of the consistency of the present calculations .
the extent to which the calculations correctly estimate @xmath0 may be judged from the fact that in the range of @xmath13 from the different estimates for it , @xmath146 .
+ we comment briefly on an estimation of the condensation energy due to superconductivity and its comparison with the increase in energy of afm fluctuations on entering superconductivity @xcite . the latter has been estimated to be almost a factor of 20 larger than the superconducting condensation energy . the suggestion has been offered that this factor of 20 may be the increase in kinetic energy . in bcs theory for electron - phonon interactions , the absolute magnitude of the change in kinetic and in potential energy are both of the same order as the condensation energy .
so , a good reason has to be found for this factor of 20 .
we do not have a solution to this enigma .
we have presented a solution to the linearized eliashberg equations using a phenomenological spin - fluctuation spectrum and simple fermi - surfaces to highlight the important parameters that determine @xmath0 for d - wave symmetry .
careful attention has been paid to the partial sum - rule on the q - dependent part of the spin - fluctuation spectra to estimate the effective coupling constant which depends on parameters such as the total partial spectral weight , the correlation length and the upper frequency cut - off of the q - dependent spin - fluctuations .
these parameters are not independent and we show their relationship in the simple model studied . with regard to the electronic structure ,
a knowledge of the average density of states at the fermi - surface is sufficient for determining @xmath0 in the s - wave channel @xcite .
but for d - wave superconductivity through exchange of well correlated spin - fluctuations , this must be supplemented by a knowledge of nesting .
the results for the general solutions are employed for two heavy fermion compounds using their measured spin - fluctuation spectra and other quantities such as specific heat and magnetic susceptibility .
correct estimates for @xmath0 to factors of @xmath16 are obtained .
confidence in these results is bolstered by getting the correct observed temperature dependence of the anomalous @xmath18 resistivity with a coefficient using the same parameters , again correct to factors of @xmath16 .
this puts a semi - quantitative backbone to the surmise made long ago that d - wave superconductivity in such heavy fermions is promoted by large amplitude spin - fluctuations with large correlation lengths such as occur near some afm quantum critical points .
+ _ acknowledgements : _ part of the work by sn was done while visiting university of california , riverside with the aid of the global coe program ( g10 ) from the japan society for the promotion of science .
the work of km is partially supported by a grand - in - aid for scientific research on innovative area `` heavy electrons '' ( no . 20102008 ) and a grant - in - aid for specially promoted research ( no .
20001004 ) from the ministry of education , culture , sports , science and technology , japan .
the work of cmv is partially supported by nsf under grant dmr-1206298 . + * appendix a : derivation of resistivity near the antiferromagnetic quantum critical point * + an expression for the resistivity under assumptions suitable for heavy fermions with a multi - sheeted fermi - surface and/or sufficient impurity scattering @xcite is given by eq .
( [ rho ] ) in terms of the self - energy function eq .
( [ aq4 ] ) . here
, we derive the relation ( [ rho ] ) explicitly . substituting @xmath147 in the spectral representation into eq .
( [ aq4 ] ) and carrying out the @xmath148-summation , one gets @xmath149 taking the analytic continuation of @xmath150 , the imaginary is given as @xmath151 ^ 2 + \left ( \frac { \varepsilon - \xi_{\vec{p } - \vec{q } } } { \gamma_{\rm afm } } \right)^2   }   \left (   \tanh \frac{\xi_{\vec{p - q}}}{2k_{\rm b}t } + \coth \frac{\varepsilon - \xi_{\vec{p - q } } } { 2k_{\rm b}t } \right ) .
\label{aq6}\end{aligned}\ ] ] we now consider the behavior at around the antiferromagnetic quantum critical point , i.e. , @xmath152 . in a low temperature region where the non - fermi liquid behavior appears , @xmath153 gives the dominant contribution for eq .
( [ aq6 ] ) .
moreover , using the following relation , @xmath154 eq .
( [ aq6 ] ) is transformed as @xmath155 next , we consider the @xmath156-integration in eq .
( [ aq8 ] ) .
because the denominator in eq .
( [ aq8 ] ) has @xmath157 term , @xmath158 gives the dominant contribution in the @xmath159-integration .
therefore , one gets @xmath160 since the integrated function in eq .
( [ aq9 ] ) rapidly decays as @xmath161 increases , we take @xmath162 as @xmath163 and obtain following result by easy calculation .
@xmath164      here , we estimate the average over the fermi - surface in eq .
( [ aq11 ] ) assuming that the fermi - surface is spherical .
@xmath168 the dominant contribution in eq .
( [ aq12 ] ) comes from `` hot '' line where the relation @xmath169 is satisfied as shown in fig .
[ fig:11 ]    assuming that the dispersion near the fermi - surface is given by linear dispersion , we obtain @xmath170 , where @xmath171 is the deviation from the `` hot '' line . for one `` hot '' spot , the integration is estimated as @xmath172 changing the integration variable as @xmath173 , eq .
( [ aq13 ] ) is transformed as @xmath174 now , we take the upper limit of the integration as @xmath163 because we consider the low temperature region , and eq .
( [ aq14 ] ) can be calculated as @xmath175      substituting eq .
( [ aq16 ] ) into eq .
( [ aq11 ] ) , we obtain @xmath178 in this calculation , the @xmath36-vector is given by @xmath179 .
the @xmath36-vector of the cecu@xmath6si@xmath6 is observed as ( 0.215 , 0.215 , 0.1458 ) giving @xmath180 .
therefore , @xmath181 is estimated as @xmath182 , and we obtain the result used for the estimation of @xmath13 : @xmath183    * appendix b : final expressions for evaluation of @xmath0 * + _ circular fermi - surface _
+ for a circular fermi - surface , it is possible to do the momentum integrals in the eliashberg equations ( [ leliash1 ] ) and ( [ leliash2 ] ) analytically so that only a diagonalization in discrete frequency space needs to be done numerically .
the final expressions used for numerical evaluation for the normal and the anomalous self - energy are :      where @xmath185,\\        x & = & \tan^{-1 } \left [ \frac { \vert \vec{q } \vert \sin \theta_{\vec{q } } + \vert \vec{p } \vert \sin \theta_{\vec{p } } } { \vert \vec{q } \vert \cos \theta_{\vec{q } } + \vert \vec{p } \vert \cos \theta_{\vec{p } } } \right ] .
\label{f4}\\\end{aligned}\ ] ]    _ tight - binding fermi - surfaces _
+ with tight binding approximation , only some simplifications in the momentum integrals in the eliashberg equations ( [ leliash1 ] ) and ( [ leliash2 ] ) can be done analytically .
the final expressions used in this paper for numerical evaluation are @xmath186    for both circular and tight binding fermi - surfaces , the best numerical strategy to evaluate @xmath0 is to cast eqs .
( [ f3 ] ) and ( [ f6 ] ) in the form of an eigenvalue equation for the eigenvector @xmath187 @xmath188 \left [ \frac{w(i \omega_m ) } { \vert \omega_m \vert }   \right ] = 0 .
\label{a9}\end{aligned}\ ] ]    it should be noted that the matrix of eq .
( [ a9 ] ) is not hermitian because @xmath189 includes the renormalization factor @xmath190 .
if the angle dependence of @xmath97 can be neglected , we can define @xmath191 in a form which does not include @xmath97 , and we obtain the eigenvalue equation with a hermitian matrix for the eigenvector @xmath192 . in s - wave superconductor case , such a situation , namely angle - independent self - energy appears .
however , in the d - wave case , @xmath193 strongly depends on @xmath90 . on including @xmath97 in the kernel @xmath191 ,
the latter is no longer symmetric for the frequency exchange , @xmath89 and @xmath194 .    at high temperatures the eigenvalues of eq .
( [ a9 ] ) are close to the negative odd integers . as the temperature decreases , the largest eigenvalue increases and crosses zero at transition temperature @xmath195 .
+    99 f. steglich , j. aarts , c. d. bredl , w. lieke , d. meschede , w. franz , and h. sch@xmath196fer , phys .
lett . * 43 * , 1892 ( 1979 ) .
h. r. ott , h. rudigier , z. fisk , and j. l. smith , phys .
* 50 * , 1595 ( 1983 ) .
g. r. stewart , z. fisk , j. o. willis , and j. l. smith , phys .
lett . * 52 * , 679 ( 1984 ) . c. m. varma , in _ proceedings of the nato advanced summer institute on the formation of local moments in metals _ ,
vancouver , canada , 1983 ( edited by w. buyers ) , plenum press , new york ( 1984 ) ; comments on solid state physics , vol .
11 , 221 ( 1985 ) .
s. schmitt - rink , k. miyake and c. m. varma , phys .
57 * , 2575 ( 1986 ) .
p. j. hirschfeld , p. w@xmath197lfle , and d. einzel , phys .
b * 37 * , 83 ( 1988 )
. k. miyake , s. schmitt - rink , and c. m. varma , phys .
b * 34 * , 6554 ( 1986 ) . c. m. varma , rep . prog
. phys . * 75 * , 052501 ( 2012 ) .
d. j. scalapino , e. loh , jr . and j. e. hirsch , phys . rev .
b * 34 * , 8190 ( 1986 ) .
t. a. maier , et al .
lett , * 95 * , 237001 ( 2005 ) ; d. snchal , et al .
lett . , * 94 * , 156404 ( 2005 ) ; m. capone and g. kotliar , phys .
b * 74 * , 054513 ( 2006 ) .
takeshi aimi , masatoshi imada , j. phys .
* 76 * , 113708 ( 2007 ) .
c.m . varma and y. yafet , phys .
13 , 2950 ( 1976 ) .
p. w. anderson , science , * 235 * , 1196 ( 1987 ) .
for a review , see p. bourges and y. sidis , c. r. physique * 20 * , ( 2010 ) . c. m. varma , phys .
b * 55 * , 14554 ( 1997 ) ; m. e. simon and c. m. varma , phys .
lett . * 89 * , 270003 ( 2002 ) .
a. schroeder et al . , phys .
80 ( 1998 ) 5623 c.m .
varma et al . ,
63 , 1996 ( 1989 ) ; v. aji and c.m .
varma , phys .
b 79 , 18450 ( 2009 ) .
o. stockert , j. arndt , e. faulhaver , c. geibel , h. s. jeevan , s. kirchner , m. loewenhaupt , k. schmalzl , w. schmidt , q. si , and f. steglich , nature phys . * 7 * , 119 ( 2011 ) .
o. stockert , j. arndt , a. schneidewind , h. schneider , h.s .
jeevan , c. geibel , f. steglich , and m. loewenhaupt , physica b , * 403 * , 973 ( 2008 ) .
j. arndt , o. stockert , h. s. jeevan , c. geibel , and f. steglich , phys .
lett . * 106 * , 24601 ( 2011 ) .
s. kambe , h. sakai , and y. tokunaga , and r. e. walstedt , phys .
b * 82 * , 144503 ( 2010 )
.    a. millis , s. sachdev and c.m .
varma , phys .
b * 37 * 4975 ( 1988 ) .
see , for example the following and references therein , p. monthoux , a.v .
balatsky and d. pines , phys .
b * 46 * , 14803 ( 1992 ) ; d.j .
scalapino , in _ handbook of high temperature superconductivity _ , edited by j. schrieffer and j. brooks ( springer , new york , 2007 ) .
m. r. norman , phys .
b * 37 * , 4987 ( 1988 ) ; a. j. millis , phys . rev .
b * 45 * , 13047 ( 1992 ) ; p. monthoux , a. v. balatsky , and d. pines , phys .
. lett . * 67 * , 3448 ( 1991 ) ; p. monthoux and d. pines , phys .
* 69 * , 961 ( 1992 ) ; p. monthoux and d. pines , phys .
b * 47 * , 6069 ( 1993 ) ; p. monthoux and d. pines , phys . rev .
b * 49 * , 4261 ( 1994 ) ; p. monthoux , j. phys : condens .
matter * 15 * , s1973 ( 2003 ) ; p. j. williams and j. p. carbotte , phys .
b * 39 * , 2180 ( 1989 ) ; c. t. rieck , d. fay , and l. tewordt , phys . rev .
b * 41 * , 7289 ( 1990 ) ; h. chi and j. p. carbotte , phys . rev .
b * 49 * , 6143 ( 1994 ) ; t. moriya , y. takahashi , and k. ueda , j. phys . soc . jpn . * 59 * , 2905 ( 1990 ) ; a. abanov , a. chubukov and m.r .
norman , phys .
b*78 * , 220507(r ) ( 2008 ) .
p. monthoux and g. g. lonzarich , phys .
b * 59 * , 14598 ( 1999 ) . p. monthoux and g. g. lonzarich , phys .
b * 63 * , 054529 ( 2001 ) . j. a. hertz , phys . rev .
b * 14 * , 1165 ( 1976 ) .
a. j. millis , phys .
b * 48 * , 7183 ( 1993 ) .
t. moriya and k. ueda , adv .
phys . * 49 * , 555 ( 2000 ) ; t. moriya and k. ueda , rep .
* 66 * , 1299 ( 2003 ) ; t. moriya and t. takimoto , j. phys .
64 * , 960 ( 1995 ) .
t. moriya , _ spin fluctuations in itinerant electron magnetism _ ( springer - verlag , 1985 ) .
h. shishido , r. settai , d. aoki , s. ikeda , h. nakawaki , n. nakamura , t. iizuka , y. inada , k. sugiyama , t. takeuchi , k. kindo , t. c. kobayashi , y. haga , h. harima , y. aoki , t. namiki , h. sato , and y. @xmath198nuki , j. phys .
jpn . * 64 * , 960 ( 1995 ) .
m. yashima , h. mukuda , y. kitaoka , h. shishido , r. settai , and y. @xmath198nuki , phys .
b * 79 * , 214528 ( 2009 ) . c. stock , c. broholm , j. hudis , h. j. kang , and c. prtrovic , phys . rev .
lett . * 100 * , 087001 ( 2008 ) ; note that the analysis in this reference is pointed out in ref .
p. gegenwart , c. langhammer , c. geibel , r. helfrich , m. lang , g. sparn , and f. steglich , phys .
lett . * 81 * , 1501 ( 1998 ) .
h. schneider , z. kletowski , f. oster , and d. wohlleben , solid state commun . * 48 * , 1093 ( 1983 ) .
a. schroeder et al . , phys .
80 ( 1998 ) 5623 c.m .
varma et al . , phys .
63 , 1996 ( 1989 ) ; v. aji and c.m .
varma , phys .
b 79 , 18450 ( 2009 ) .
w. l. mcmillan , phys .
167 , 331 ( 1968 ) .
p. b. allen and r. c. dynes , phys .
b * 12 * , 905 ( 1975 ) .
see , for example , h.v .
lohneysen , a. rosch , m. vojta and p. woelfle , rev .
79 , 1015 ( 2007 ) .
r. hlubina and t. m. rice , phys .
b * 51 * , 9253 ( 1995 ) .


anderson localization  @xcite of electronic states in a disordered system is a path breaking observation that has extended its realm well beyond the electronic properties of disordered solid materials , and has been found out to be ubiquitous in a wide variety of systems such as the photonic  @xcite , phononic  @xcite , polaronic  @xcite , or plasmonic lattices  @xcite to name a few .
artificial , tailor made geometries developed using the improved fabrication and lithographic techniques have been instrumental in testing the validity of the basic theory of exciton localization in presence of disorder .
examples include the direct observation of localization of matter waves in recent times  @xcite .
localization is a result of quantum interference that has the strongest manifestation in one dimension where all the single particle states exhibit an exponential decay of amplitude as one moves away from a particular lattice point . for dimensions beyond one ,
the states retain their characteristic decay in amplitude , with the addition of a possibility of a metal - insulator transition in three dimensions .
the results have been substantiated by calculations of the localization length  @xcite , density of states  @xcite and their multi - fractality  @xcite .
extensive work has also revealed the intricacies of the single parameter scaling hypothesis  @xcite , its validity  @xcite or , variance  @xcite , and even violation  @xcite in low dimensional systems within a tight binding approximation .
the studies have subsequently been consolidated by direct experimental measurements of conductance distribution in quasi - one dimensional gold wires  @xcite .    in the last couple of decades
variations in the canonical cases of anderson localization have surfaced , particularly , within a tight binding description .
resonant tunnelling of electronic states has been shown to arise even in a one dimensional disordered system as a consequence of certain special kinds of positional correlation  @xcite .
discrete energy levels corresponding to _ extended eigenfunctions _
@xcite , appear within a typical point spectrum exhibited by the randomly disordered lattices rendering the two - terminal transport completely ballistic at special values of the fermi level .
such results are common even to infinite quasicrystalline lattices in one dimension where , the self similar character of the chain really allows an infinite number of unscattered ( extended ) eigenfunctions despite the absence of any translational invariance  @xcite .    in this context ,
rather unusual , yet meaningful question could be , is it possible to engineer , in a controlled fashion , an _ absolutely continuous _ portion in the energy spectrum of non - translationally invariant systems with the help of some kind positional correlations ?
the existence of real continuous bands are reported recently in quasi one dimensional or two dimensional systems only with diagonal disorder  @xcite .
we thus have a partial answer to the question , and it remains to be seen whether such bands of extended eigenfunctions can be generated and controlled with off - diagonal disorder are well .
we address this issue with respect to quasiperiodically ordered geometries . in this communication
we show that , in quasi - one dimensional ladder networks with off - diagonal quasi - periodic order it is possible to engineer absolutely continuous bands of energy eigenvalues .
we specifically work with the fibonacci quasiperiodic order , though the observation is by no means , restricted to them .
ladder networks within the framework of tight binding model have already been exploited to gain insight into charge transport in dna  @xcite , to unravel interesting physics related to polymer structures and ployacene lattice topology  @xcite .
comparatively speaking , controlled generation of absolutely continuous portions in the energy spectrum of the ladder networks , disordered or quasi - periodic , still remains to be explored .
if it works then such an effort may open up new physics , or allow one to explore the possibility of devising novel low dimensional filter networks .
the central result of this paper is that , one can _ create _ absolutely continuous subbands in the electronic spectrum of such quasi 1-d systems , by attaching a single quantum dot ( qd ) from one side to a selected set of lattice points in every strand of the ladder .
we discuss the results explicitly in terms of fibonacci quasi - periodic two and three arm ladder network models using a tight binding hamiltonian and a real space renormalization group ( rsrg ) method .
we also investigate whether a crossover in the fundamental character of the eigenfunctions , going from a typical _ critical _ nature to a bloch - like _ extended _ character can be engineered .
as will be shown , this indeed turns out to be quite a possibility in the case of the quasi - periodic ladder network ( qpln ) , particularly when the number of strands ` @xmath0 ' in the transverse direction increases .
the results hold even for an off - diagonally disordered ladder where a re - entrant insulator - metal - insulator transition becomes feasible .
the system is described by a tight binding hamiltonian , and is depicted in fig .
[ lattice1 ] .
the on - site potential at every lattice point is kept constant , while the nearest neighbor hopping integral along the arms of the ladder assumes two values dictated by the growth rule of the desired quasi - periodic sequence .
the electron can tunnel between the arms of the ladder by hopping from one site to its corresponding site on the other arm through the hopping integral @xmath1 ( see fig .
[ lattice1 ] ) .
the sites marked by @xmath2 in each arm of the ladder is coupled to a single level qd  @xcite attached from one side .
the qd couples to every @xmath2-site via a hopping integral @xmath3 .
there is no positional correlation , short range or long range in the conventional sense  @xcite .
attachment of qd s from a side can quench the conductance of a linear chain  @xcite , but an array of qd s in a ladder can have reverse effect , as will be seen
.    we will show in a completely analytical way that , for a particular relationship between the nearest neighbor hopping integrals along the ladder strand and the qd - strand coupling @xmath3 , and for a special value of the qd potential , we get _ absolutely continuous subbands _ in the spectrum .
this immediately raises a non - trivial issue that , the nature of the electronic spectrum and the eigenstates over a _ continuous range of eigenvalues _ can become sensitive to the _ numerical values _ of the hamiltonian parameters  a fact that is not commonly observed in the conventional anderson localization problem .
the full spectrum turns out to be an admixture of conventional _ critical eigenstates _
, characteristics of a fibonacci quasiperiodic chain , and absolutely continuous parts , giving rise to the possibility of a _ phase transition _ , particularly when the system grows in the transverse direction .
our results are analytically exact .    in section  [ sec2 ]
we describe the hamiltonian and the basic working principle for the two and three strand quasi - periodic ladders constructed following a fibonacci sequence . along with this , the real space renormalization scheme to obtain the energy spectrum and to judge the character of the single particle states is also presented . in section  [ sec3 ]
we present a result for the two terminal transmission coefficient of a finite sized qpln , and in section  [ sec4 ] we draw our conclusion .
in wannier basis the tight binding hamiltonian of a general @xmath0-arm ladder network reads  @xcite , @xmath4\ ] ] [ hamilton ] where , @xmath5 is , in general , an @xmath6 on - site potential matrix at the @xmath7-th vertical rung .
we shall explicitly discuss the two - strand ladder at first .
the cases @xmath8 can be trivially constructed .
@xmath9 is the creation ( annihilation ) operator , represented by rows ( columns ) of appropriate dimension , and @xmath10 is the @xmath6 hopping matrix representing hopping along and in between strands of the ladder , @xmath0 being equal to @xmath11 , @xmath12 or any integer .
the fibonacci ladder network is grown by attaching identical one dimensional fibonacci chains in the transverse direction .
each individual chain grows along the horizontal direction following the binary fibonacci sequence of two letters @xmath13 and @xmath14 ( representing two bonds ) .
the growth rule is , @xmath15 and @xmath16 , and the sequence begins with @xmath13 .
the on - site potentials at the vertices in each strand are designated by three symbols viz .
, @xmath17 , @xmath18 and @xmath19 for sites flanked by two @xmath13-bonds , @xmath13 on left and @xmath14 on right , and @xmath13 on right and @xmath14 on left respectively .
the attached qd is shown as  @xmath20 " in fig .
[ lattice1 ] , and the potential at the side coupled qd is designated by @xmath21 . the nearest neighbor hopping integral @xmath22 or @xmath23 along
ladder - arms following the fibonacci sequence .
@xmath24 is the inter - arm hopping , and @xmath25 designates the tunnel hopping between the qd and the @xmath2-sites of the ladder network . in our work , we shall assume @xmath26 , and @xmath27 .
so , this is purely a _ transfer model _
@xcite in the standard language of one dimensional quasicrystals
.    we shall analyze the problem using a method described earlier  @xcite in which an ` @xmath0'-strand ladder ( @xmath28 in the first phase of our discussion ) is decoupled , making a change of basis , into a system of @xmath0 independent linear chains .
the complete spectrum of the quasi @xmath29-d system of @xmath0 strands is then obtained by convolving the spectra of the individual linear chains .
in addition , we perform an rsrg analysis , which does not need any decoupling and yields a density of states spectrum that is in conformity with our results obtained using the decoupling scheme .
the two terminal transmission spectrum is also computed using the rsrg decimation scheme to highlight the effect of the side coupled qd s .
the analysis will be carried over to the three strand ladder , and a generalization to the @xmath0 strand case will be discussed using arguments .        at first , we renormalize the potential at the @xmath2-sites by decimating the side coupled qd .
this leads to a _
renormalized _ value of the on - site potential at the @xmath2-vertex , given by , @xmath30 .    , @xmath31 and @xmath32 are colored differently to distinguish between the nearest neighbor bond - environments . the qd s @xmath20 are coupled to the @xmath2-sites in each case .
the coupling between @xmath20 and @xmath2-site is denoted by @xmath3 and the inter - strand coupling is denoted by @xmath1.,width=302 ]    the schrdinger equation for the two - arm ladder can now be written equivalently in the form of the difference equations involving @xmath33 matrices : @xmath34 where , @xmath35 and , @xmath36 and , @xmath37 is the @xmath33 unit matrix . written explicitly , the matrix equation reads ,    @xmath38 \left ( \begin{array}{c } \psi_{i , a } \\ \psi_{i , b }    \end{array } \right ) =   \left ( \begin{array}{cccc } t_{i , i+1}^a & 0 \\   0 & t_{i , i+1}^b \end{array } \right ) \left ( \begin{array}{c } \psi_{i+1,a } \\ \psi_{i+1,b }    \end{array } \right )   + \left ( \begin{array}{cccc } t_{i , i-1}^a & 0 \\   0 & t_{i , i-1}^b \end{array } \right ) \left ( \begin{array}{c } \psi_{i-1,a } \\ \psi_{i-1,b }    \end{array } \right ) \label{differenceeq}\ ] ]    in eq .
, the potential @xmath39 takes on values @xmath40 corresponding to a renormalized @xmath2-site in each arm , while it is equal to @xmath41 for both the @xmath31 and @xmath32-sites in each arm .
the nearest neighbor hopping integrals @xmath42 are @xmath43 or @xmath23 depending on the local environment of the corresponding site in the arm @xmath44 . a look at the fig .
[ lattice1 ] will make it obvious that there are three such equations corresponding to the three different kinds of vertices viz .
, @xmath2 , @xmath31 and @xmath32 respectively .
the potential matrix @xmath5    corresponding to the renormalized @xmath2-sites is , @xmath45 and that for the @xmath31 and @xmath32 sites are @xmath46 we now diagonalize the @xmath47 matrix using a matrix @xmath48 .
the matrix @xmath48 makes a change over to a new basis defined by , @xmath49 .
as shown elsewhere  @xcite , such a change of basis leads to a decoupling of the set of eq .   into two independent equations , one corresponding to the @xmath50-arm alone , and the other , to the arm @xmath51 . written explicitly ,
the equations read , @xmath52 \phi_{i , a }   =   t_{l } \phi_{i+1,a } + t_{l } \phi_{i-1,a } \nonumber \\ & \left [ e - ( \epsilon_{0}+\xi ) \right ] \phi_{i , a }   =    t_{s } \phi_{i+1,a } + t_{l } \phi_{i-1,a } \nonumber \\ & \left [ e - ( \epsilon_{0}+\xi ) \right ] \phi_{i , a }   =   t_{l } \phi_{i+1,a } + t_{s } \phi_{i-1,a }   \label{arm1}\end{aligned}\ ] ] for the sites @xmath2 , @xmath31 and @xmath32 sequentially in arm @xmath50 . similar set of equations for @xmath2 , @xmath31 and @xmath32 sites in arm
@xmath51 will be , @xmath53 \phi_{i , b }   =    t_{l } \phi_{i+1,b } + t_{l } \phi_{i-1,b } \nonumber \\ & \left [ e - ( \epsilon_{0}-\xi ) \right ] \phi_{i , b }   =    t_{s } \phi_{i+1,b } + t_{l } \phi_{i-1,b } \nonumber \\ & \left [ e - ( \epsilon_{0}-\xi ) \right ] \phi_{i , b }   =    t_{l } \phi_{i+1,b } + t_{s } \phi_{i-1,b }   \label{arm2}\end{aligned}\ ] ]    needless to say , each of the sets eq .
and eq .   separately represents an independent fibonacci chain .
the spectrum of each chain bears the usual singular continuous character .
the spectrum of the original ladder network can be obtained by convolution of the density of states of the individual decoupled fibonacci chains . however , an analysis of these decoupled sub - systems yields rich spectral insight that would otherwise be difficult to obtain .
we are now in a position to throw some light on this issue .
first , we need to appreciate that , in the new basis each decoupled fibonacci chain is described by a set of three _ transfer matrices _ , viz . , @xmath54 , @xmath55 and @xmath56  @xcite . for the @xmath50-arm decoupled fibonacci chain , the explicit expressions of these matrices are , @xmath57/t_{l } & -1 \\   1 & 0   \end{array }
\right ) \nonumber \\
\mbox{\boldmath $ m_{\beta , a}$ } & = & \left ( \def\arraystretch{1.5}\begin{array}{cccc } [ e-(\epsilon_{0}+\xi)]/t_{s } & -t_{l}/t_{s } \\   1 & 0   \end{array } \right ) \nonumber \\
\mbox{\boldmath $ m_{\gamma , a}$ } & = & \left ( \def\arraystretch{1.5}\begin{array}{cccc } [ e-(\epsilon_{0}+\xi)]/t_{l } & -t_{s}/t_{l } \\   1 & 0   \end{array } \right ) \label{arma}\end{aligned}\ ] ] respectively .
a similar set of transfer matrices for the decoupled @xmath51-arm reads , @xmath58/t_{l } & -1 \\   1 & 0   \end{array }
\right ) \nonumber \\
\mbox{\boldmath $ m_{\beta , b}$ } & = & \left ( \def\arraystretch{1.5}\begin{array}{cccc } [ e-(\epsilon_{0}-\xi)]/t_{s } & -t_{l}/t_{s } \\   1 & 0   \end{array } \right ) \nonumber \\   \mbox{\boldmath $ m_{\gamma , b}$ } & = & \left ( \def\arraystretch{1.5}\begin{array}{cccc } [ e-(\epsilon_{0}-\xi)]/t_{l } & -t_{s}/t_{l } \\   1 & 0   \end{array } \right ) \label{armb}\end{aligned}\ ] ] @xmath54 and @xmath59
follow an arrangement in the fibonacci sequence  @xcite .    speaking in terms of
the individual decoupled fibonacci chains @xmath50 and @xmath51 , we already know that each such chain offers a singular continuous spectrum .
if we set @xmath60 , then each chain describes a separate fibonacci lattice with @xmath61 , and @xmath27 . the total transfer matrix for any @xmath62-th generation of such a chain is related to its predecessors at the @xmath63th and @xmath64th generations by the matrix map @xmath65 , which will exhibit a _ six cycle _  @xcite at @xmath66 for the decoupled @xmath50 and @xmath51 chains respectively .
the existence of a six cycle of the matrix map implies that the above energy values are indeed eigenvalues of the infinite decoupled chains , and hence , the infinite ladder network as well .
thus , the first fruit of the decoupling scheme is that , with @xmath67 we obtain at least two eigenvalues which definitely belong to the spectrum of the infinite ladder .
extending this idea to an @xmath0-strand fibonacci ladder we can easily identify @xmath0 different energy eigenvalues of the full quasi @xmath29-d fibonacci network in the absence of @xmath3 ( i.e. detaching the qd s ) in a completely analytical way .
the task becomes non - trivial when one consider the qpln as a whole , without decoupling .
now let us look at the more serious effect when we include the qd by making @xmath3 non - zero .
let us choose eq .  .
we set @xmath69 .
it can easily be verified that , with this potential , the commutator @xmath70 = 0 $ ] for the decoupled @xmath50-arm , _ irrespective of the energy @xmath71 _ of the electron , if @xmath72 .
this implies that , with these conditions the individual @xmath2-sites and the @xmath73 dimers in the decoupled @xmath50-arm can be arranged in any desired pattern , for example , in a perfectly periodic pattern .
the spectrum offered by the decoupled @xmath50-arm fibonacci chain will thus , under this choice of the qd potential @xmath21 and the tunnel hopping @xmath3 , will turn out to be _
indistinguishable _ from that of a perfectly periodic sequence of the @xmath2-site ( the ` renormalized ' @xmath2-site , to be precise ) and the @xmath73 doublet .
however , this argument needs to go through an acid test which we discuss below .
+ qd chain , and its renormalized version .
( b ) the @xmath73 binary ordered chain .
the ldgf is found at the site on the horizontal line in each case .
the qd is renormalized to produce an effective one dimensional periodic chain of @xmath74 sites , as already discussed in the text.,width=264 ]    owing to the commutation of @xmath54 and @xmath75 , the decoupled @xmath50-arm can be thought to be composed of two infinite periodic chains , one with @xmath2-sites coupled to the qd s , and the other , a periodic array of @xmath73 clusters .
the two chains are shown in fig .  [ chains ] .
each such infinite sublattice will have its own energy bands .
it is important to verify whether , under the special condition imposed on the strand - dot coupling @xmath3 , the periodic @xmath2-lattice and the periodic @xmath73 lattice share the same eigenvalue spectrum .
if yes , then we conclusively prove the existence of an absolutely continuous part in the energy spectrum of the original ladder network .    to check this , we first of all construct an effective one dimensional chain composed of renormalized @xmath2-atoms , obtained by decimating the side coupled dot .
the resulting chain is illustrated in fig .
[ chains](a ) with encircled lattice points .
each such lattice point has the renormalized on - site potential equal to @xmath76 .
each atom in the @xmath73 chain has the potential @xmath77 .
we now evaluate the local density of states ( ldos ) at any site of the renormalized ( encircled ) @xmath2-chain ( let s call it the @xmath74 chain ) and any site in the periodic @xmath73 chain ( fig .
[ chains](b ) ) .
the local densities of states are given , for these two infinite periodic lattices , and for a given set values of @xmath78 , @xmath1 , @xmath43 and @xmath23 , by @xmath79^{-1/2}$ ] and @xmath80^{-1/2}$ ] , where , @xmath81 and , @xmath82    if we now set @xmath83 , then a simple algebra yields that , the difference @xmath84 is given by , @xmath85 ( \lambda^2 + t_{l}^2 - t_{s}^2 ) } { ( e - \epsilon_{0 } - \xi ) ^2 } \label{den}\ ] ] which clearly vanishes for any value of the energy of the electron the moment we tune @xmath72 .
thus the ldos s of the two basic constituent periodic chains become identical _ independent of energy _ , resulting in a complete overlap of the spectra of these individual chains under the above resonance condition . as each of the @xmath2 and the @xmath73 sublattices
provide a two subband continuum in their energy spectra ( an obvious effect of the binary periodic lattice ) , it is obvious that these two spectral continua will retain their positions in the full spectrum of the fibonacci qpln when the qd potential and the qd - strand tunnel coupling are tuned to the above special values .
it should be appreciated that once we fix @xmath86 , the set of eq .   corresponding to the decoupled @xmath51 strand ,
represents a mixed model linear fibonacci chain .
its contribution to the full energy spectrum is the usual fragmented , multi - fractal one .
as the full spectrum of the two - arm ladder will be a convolution of both , we expect it to have absolutely continuous parts decorated with a spiky envelope , and a mixture of _ extended _ bloch - like states and the _ critical _ eigenfunctions which are characteristic of a quasi - periodic fibonacci sequence of potentials .
this is precisely what we observe in fig .
[ 2stranddos](a ) .
at this point one must be careful to interpret fig .
[ 2stranddos](a ) . in the full energy spectrum , critical and extended eigenstates
can never coexist .
so , even in the full spectrum of the un - decoupled qpln the entire range shown in fig .
[ 2stranddos](a ) by grey shade corresponds to extended bloch - like states only . beyond it
, every wave function will be critical in nature .
the contributions coming from each decoupled strand are only shown separately in this figure .
the above arguments are confirmed by an exact real space renormalization group ( rsrg ) calculation of the local density of states , first of each individual decoupled strands and then of the full two - arm fibonacci ladder without decoupling it . the first task is already discussed in literature  @xcite , and need not be talked about again in here .
we instead , provide the details of the rsrg scheme of the full ladder network .
we renormalize the full two - strand ladder by decimating the @xmath31-vertices .
this generates _ new _ hopping integrals across the diagonals of the fundamental rectangular plaquettes of renormalized ladder denoted by @xmath87 and @xmath88 in the equations to follow . to implement the decimation we need to assign three different symbols to the inter - arm hopping @xmath1 , viz . ,
@xmath89 , @xmath90 and @xmath91 connecting the @xmath2 , @xmath31 and @xmath32 vertices of the two strands .
of course , we begin with @xmath92 , and @xmath93 .
the recursion relations arising out of this decimation procedure are given by , @xmath94 where @xmath95 a small imaginary part is added to the energy @xmath71 , and the ldos at an @xmath2 , @xmath31 or @xmath32 site is obtained by calculating the respective local green s function @xmath96 , where @xmath97 or @xmath32 .
@xmath98 represents the fixed point of the respective on - site potential as the hopping integrals flow to zero under the rsrg iterations .
the ldos is given by , @xmath99\ ] ]    the results are presented in fig .  [ 2stranddos ] . in panel
( a ) the spectra of the decoupled @xmath50- and @xmath51-arms are shown separately , while panel ( b ) illustrates the spectrum of the full two - arm ladder as obtained from the rsrg study .
the existence of the absolutely continuous portions in the energy spectrum is obvious .
the character of the eigenfunctions belonging to the spectrum is checked from an observation of the flow of the the hopping integrals @xmath43 , @xmath23 and the diagonal hopping amplitude under rsrg iterations . for any energy eigenvalue , chosen from any portion of the absolutely continuous sub - bands ( grey shaded portion in fig .  [ 2stranddos](a ) ) all the hopping integrals , as obtained from eq .
, remain non - zero for an arbitrarily large number of iterations confirming the extended character of the wave functions . beyond the extended regime , the hopping integrals flow to zero under iteration indicating that the eigenfunctions are not extended . in fig .  [ 2stranddos](b ) we see the full ldos spectrum from the coupled two - strand qpln obtained from the rsrg method ( using the set of eq .  ) .
the light blue shade in the spectrum is just to highlight the region of non - zero ldos .
the extended eigenstates occupy the precise positions as depicted by the grey shaded part in fig .  [
2stranddos](a ) .
the eigenstates form a mixed spectrum of critical and extended states , as one traverses the energy axis from the left .
-vertex for a two - strand ladder network .
the upper panel shows the ldos spectra of the decoupled individual linear chains , one chain satisfying resonance condition .
the lower panel shows the ldos spectrum of the composite ladder network as obtained from the rsrg scheme , under resonance condition .
the values of the parameters chosen are @xmath100 , @xmath101 , @xmath102 , @xmath103 , @xmath104 , and @xmath105 , measured in unit of @xmath43.,width=264 ]    an additional support in regard of the continuous subbands and the extended eigenstates is also obtained from an observation of the flow pattern of the entire parameter space of the decoupled strands that satisfies the condition @xmath106 and @xmath107 .
for such a choice , under successive rsrg step , we always get @xmath108 and @xmath109 for the decoupled @xmath50-strand .
this is a typical flow in the parameter space that corresponds to the extended eigenstates in quasiperiodic fibonacci chains  @xcite . in the present case
it confirms the extended bloch - like character of the wave function for the entire energy regime ( given by the grey color in fig .  [ 2stranddos](a ) ) .
the motivation behind extending the above analysis to a three - arm fibonacci ladder stems    .,width=302 ]    from the curiosity to know whether , as one extends the ladder network in the transverse direction , any possibility of a re - entrant transition ( or , crossover ) from a singular continuous to an absolutely continuous energy spectrum becomes apparent .
for disordered ladders comprising the same @xmath2 and @xmath73 clusters this would imply the possibility of an _ insulator - metal _ transition that can be triggered by an appropriate choice of the numerical values of the qd potential @xmath21 and the tunnel hopping integral @xmath3 .    to check such a possibility we decouple a three - arm fibonacci ladder ( fig .  [ lattice2 ] ) .
we follow the same prescription for decoupling as outlined in the case of a two - strand ladder .
diagonalization of the potential matrix @xmath5 now yields three eigenvalues , viz .
, @xmath41 and @xmath110 , and the effective potential at the @xmath2-site in each decoupled arm will be @xmath111 ( for arm @xmath50 , say ) , @xmath112 ( for arm @xmath51 , say ) and @xmath113 ( for arm @xmath114 ) .
we now end up having three sets of equations , one for each decoupled strand .
for example , @xmath115   \phi_{i , a } = t_{l } \phi_{i+1,a } + t_{l } \phi_{i-1,a } \nonumber \\ & ( e - \epsilon_{0 } ) \phi_{i , a } = t_{s } \phi_{i+1,a } + t_{l } \phi_{i-1,a } \nonumber \\ & ( e - \epsilon_{0 } ) \phi_{i , a } = t_{l } \phi_{i+1,a } +   t_{s } \phi_{i-1,a }   \label{3decouple1}\end{aligned}\ ] ] corresponding to the @xmath2 , @xmath31 and @xmath32 sites respectively in the decoupled @xmath50-strand .
two sets of similar equations are also obtained for the decoupled @xmath51- and @xmath114-strands , with the effective potential at the @xmath2-site being equal to @xmath116 and @xmath117 respectively .
the effective potentials at the @xmath31 and @xmath32 sites in the @xmath51 and @xmath114 strand are @xmath110 as there is no side coupling to these sites .    if , as before , we set @xmath118 in eq .  , then @xmath70=0 $ ] for the decoupled @xmath50-strand .
the contribution to the full spectrum coming from this separate chain will be two absolutely continuous sub - bands . with the potential of the qd set at @xmath118 ,
the two other decoupled strands represent two different fibonacci chains in the mixed model .
the spectrum resulting from each of them will remain fragmented cantor type .
the overall spectrum , as we now have known , will be a convolution of the three .
the absolutely continuous portions will still hold their positions in the main spectrum , and will be flanked by the spiky , fragmented parts characteristics of critical eigenstates of a fibonacci lattice ( fig  [ 3stranddos ] ) .
the eigenfunctions belonging to the continuous sub - bands will , needless to say , be of extended character .
.,width=264 ]    one needs to appreciate that each spectrum offered separately by the decoupled @xmath51- and @xmath114-strands will have different  band - centres " around which the singular continuous spectra of the individual decoupled fibonacci chain will spread out . with an @xmath0-strand ladder model , there will be @xmath119 such centers from the @xmath119 decoupled strands around which cantor set energy eigenvalues will be distributed , the remaining strand contributing the absolutely continuous sub - bands .
the centres will be densely populated around the outer edge of the continuous sub - bands as @xmath120 .
thus , with increasing number of strands there is a possibility of observing a smooth critical to extended state crossover in the spectrum of such quasi @xmath29-d aperiodic ladder models .
clearly , the argument prevails even if we have an @xmath0-strand off diagonally disordered ladder network .
the side coupled qd s with appropriately tuned potential one can then expect a re - entrant _ insulator - metal _ transition .    before leaving this discussion , it is pertinent to mention that , the commutaion of transfer matrices leading to a complete , energy independent delocalization of electronic states have been reported recently in certain special examples of quasi - one dimensional lattices  @xcite .
but , those systems do not open up any possibility of an insulator - metal transition , in the current spirit of the work . at the same time
, it is useful to note that , the commutivity of transfer matrices has been used in the past to unravel discrete unscattered eigenstates in an array of quantum wells  @xcite or a disordered distribution of pschl - teller potential  @xcite .
the difference with the present communication and these two latter references lies in the energy - independent character of the matrix - commutation that plays the key role .
we now present partial results ( to save space ) for the two terminal transmission coefficient of a finite two strand ladder without decoupling it .
no change of basis is made , and the full finite qpln is treated within the rsrg scheme .
the basic working formula is already elaborated in the existing literature  @xcite .
the ladder is clamped between two semi - infinite perfectly ordered leads ( viz .
, the source and the drain ) at any two extremities ( fig .  [ lattice3 ] ) .
the leads are characterized by uniform on - site potential @xmath121 and nearest neighbor hopping integral @xmath122 .
the portion of the ladder clamped between the leads is renormalized to bring it down to a _ diatomic molecule _ , and the transmission coefficient is easily obtained  @xcite .
the procedure is standard and to save space again , we do not provide the details .    it should be appreciated that , the precise details of the transmission coefficient definitely depends on which points in the ladder the leads are attached to .
only one result is shown here in which the two semi - infinite ordered leads are connected to the extreme points of the upper strand of a @xmath123-th generation fibonacci ladder having @xmath124 bonds in each arm .
-th generation two - strand fibonacci ladder network under resonance condition .
both the incoming and outgoing leads , viz .
, the source and the drain are connected in the upper strand . system parameters chosen are @xmath100 , @xmath101 , @xmath102 , @xmath103 , @xmath104 , and @xmath105 , and the lead parameters are @xmath121 and @xmath122.,width=264 ]    one needs to appreciate that , even in the absence of any attached qd , a two strand finite qpln can exhibit high transmission spikes with @xmath125 close to unity .
this normally results from the finite size of the system , as well as from the fact that a ladder network consists of loops that brings in a flavor of long range connection , and may well be responsible for quantum interference leading to high transmittivity at selected values of electron - energy .
also , the attachment to the leads , and the difference between the hopping integral in the lead and that in the bulk of the system give rise to strong scattering effects at the contacts , affecting the overall transmission of the qpln system .
therefore , the results obtained from a finite size analysis should be carefully compared with the predictions which come from the anatomy of an infinite system .    nevertheless , the existence of the continuous parts in the energy spectrum of the infinite system is already showing up in the enhanced value of the transmission coefficient at appropriate values of the electron energy as can be seen in fig .  [ transport ] .
the apparent drop to close to zero value at certain regions even within the continuous zone is due to the scale at which the transmission coefficient is shown here .
however , the transmission coefficient exhibits wide fluctuations owing to the quasiperiodic nature of network .
the effect of quantum interference is manifestly observed here .
the finite value of the transmission coefficient even beyond the absolutely continuous part of the spectrum is due to the finite size of the system , and gets lower and lower , and the spectrum gets more spiky with increasing system size .
also , one can anticipate eigenstates with their ` localization length ' much larger compared to the size of the system as considered in the present calculation .
we have shown that fibonacci quasi - periodic ladder networks , where quasi - periodic order is introduced in the distribution of bonds , keeping the potential at every vertex same , are capable of producing absolutely continuous portions in their energy spectra .
this requires the attachment of quantum dots from one side to a special set of sites in each arm of the ladder .
a suitable adjustment of the dot potential , together with the dot - backbone coupling can generate absolutely continuous energy sub - bands , decorated at the flanks by the usual critical states  a characteristic of quasi - periodic lattices .
apart from opening up the possibility of engineering the transport of excitons in such multi - strand systems , a very basic issue of anderson localization is reviewed .
it seems that , the existence of localized eigenstates in similar disordered ladder networks may strongly depend on the numerical values of the parameters in the hamiltonian .
the analysis is carried over to other quasi - periodically ordered systems , such as the copper mean chain , the aperiodic thue morse chain and even to the disordered ladder networks with locally and non - locally coupled dots .
the conditions for the existence of a continuum in the energy spectrum becomes much more non - trivial in such systems .
these issues will be reported in a forthcoming article .
is grateful to dst , india for an inspire fellowship . a.c . acknowledges financial support from dst , india through a purse grant , and the abdus salam international centre for theoretical physics , trieste , italy for its hospitality and financial support .
60 p. w. anderson , phys .
* 109 * , 1492 ( 1958 ) .
b. kramer and a. mackinnon , rep .
* 56 * , 1469 ( 1993 ) .
e. abrahams , p. w. anderson , d. c. liciardello , and t. v. ramakrishnan , phys .
* 42 * , 673 ( 1979 ) .
e. yablonovitch , phys .
* 58 * , 2059 ( 1987 ) .
s. john , phys .
lett . * 58 * , 2486 ( 1987 ) .
f. r. montero de espinosa , e. jimnez , and m. torres , phys .
80 * , 1208 ( 1998 ) .
j. o. vasseur , p. a. deymier , g. frantziskonis , g. hong , b. djafari - rouhani , and l. dobrzynski , j. phys . : condens .
matter * 10 * , 6051 ( 1998 ) .
i. o. barinov , a. p. alodzhants , and s. m. arakelyan , quantum electron .
* 39 * , 685 ( 2009 ) .
m. grochol and c. piermarocchi , phys .
b * 78 * , 035323 ( 2008 ) .
a. tao , p. sinsermsuksaul , and p. yang , nature nanotechnol .
* 2 * , 435 ( 2007 ) .
a. christ , y. ekinci , h. h. solak , n. a. gippus , s. g. tikhodeev , and o. j. f. martin , phys .
b * 76 * , 201405 ( 2007 )
. b. damski _
lett . * 91 * , 080403 ( 2003 ) .
et al . _ , nature ( london ) * 453 * , 891 ( 2008 ) .
et al . _ , nature ( london ) * 453 * , 895 ( 2008 ) .
r. a. rmer and h. schulz - baldes , europhys
. lett . * 68 * , 247 ( 2004 ) .
a. eilmes , r. a. rmer , and m. schreiber , physica b * 296 * , 46 ( 2001 ) .
a. rodrguez , j. phys . a : math
. gen . * 39 * , 14303 ( 2006 ) .
a. rodriguez , l. j. vasquez , and r. a. rmer , phys .
b * 78 * , 195107 ( 2008 ) . a. rodriguez , l. j. vasquez , k. slevin , and r. a. rmer , phys . rev .
b * 84 * , 134209 ( 2011 ) .
s. d. pinski , w. schirmacher , and r. a. rmer , europhys . lett . *
97 * , 16007 ( 2012 ) .
lev i. deych , a. a. lisyanski , and b. l. alshuler , phys .
lett . * 84 * , 2678 ( 2000 ) .
m. titov and h. schomerus , phys .
lett . * 95 * , 126603 ( 2005 ) .
j. w. kantelhardt and a. bunde , phys .
b * 66 * , 035118 ( 2002 ) .
v. e. kravtsov and v. i. yudson , phys .
b * 82 * , 195120 ( 2010 ) . p. mohanty and r. a. webb , phys .
lett . * 88 * , 146601 ( 2002 ) .
d. h. dunlap , h .-
wu , and p. phillips , phys .
lett . * 65 * , 88 ( 1990 ) .
f. a. b. f. de moura and m. lyra , phys .
81 * , 3735 ( 1998 ) . s. sil , s. k. maiti , and a. chakrabarti , phys .
b * 78 * , 113103 ( 2008 ) .
a. rodriguez , a. chakrabarti , and r. a. rmer , phys .
b * 86 * , 085119 ( 2012 ) .
a. chakrabarti , s. n. karmakar , and r. k. moitra , phys .
b * 50 * , 13276 ( 1994 ) .
e. maci , phys .
b * 74 * , 245105 ( 2006 ) ; phys . rev .
b * 80 * , 125102 ( 2009 ) .
m. bravi , r. farchioni , g. grosso , and g. p. parravicini , phys .
b * 87 * , 035105 ( 2013 ) .
r. farchioni , g. grosso , and g. p. parravicini , eur .
j. b * 84 * , 227 ( 2011 ) .
b. kubala and j. knig , phys .
b * 67 * , 205303 ( 2003 ) .
r. farchioni , g. grosso , and g. p. parravicini , phys .
b * 85 * , 165115 ( 2012 )
. m. kohmoto , b. sutherland , and c. tang , phys .
b * 35 * , 1020 ( 1987 ) . j. a. ashraff and r. b. stinchcombe , phys . rev .
b * 37 * , 5723 ( 1988 ) .
a. chakrabarti , s. n. karmakar , and r. k. moitra , phys .
b * 39 * , 9730 ( 1989 ) .
a. chakrabarti , s. n. karmakar , and r. k. moitra , phys . lett .
a * 168 * , 301 ( 1992 ) .
b. pal , s. k. maiti , and a. chakrabarti , europhys . lett . * 102 * , 17004 ( 2013 ) .
i. gmez , f. dminguez - adame , and e. diez , physica b * 324 * , 235 ( 2002 ) .
a. rodrguez and j. m. cerver , phys .
b * 74 * , 104201 ( 2006 )
. a. d. stone , j. d. joannopoulos , and d. j. chadi , phys .
b * 24 * , 5583 ( 1981 ) .
a. chakrabarti , phys .
b * 74 * , 205315 ( 2006 ) .
b. pal and a. chakrabarti , eur .
j. b * 85 * , 307 ( 2012 ) .