since its discovery fifty years ago @xcite anderson localization has remained a very active field of research ( see for instance the reviews @xcite ) . according to the scaling theory @xcite , there is no delocalized phase in dimensions @xmath20 , whereas there exists a localization / delocalization at finite disorder in dimension @xmath21 .
to get some insight into this type of transition , it is natural to consider anderson localization on the cayley tree ( or bethe lattice ) which is expected to represent some mean - field limit .
the tight - binding anderson model on the cayley tree has been thus studied by various authors . in @xcite , the recursion equation for
the self - energy has been studied to establish the mobility edge as the limit of stability of localized states . in @xcite , the exponent @xmath19 governing the divergence of the localization length
was shown to be @xmath22 . in @xcite ,
the supersymmetric formalism has been used to predict the critical behavior of some disorder - averaged observables ( see also @xcite where similar results have been obtained for the case where anderson tight - binding model is replaced by a non - linear @xmath23-model ) . in @xcite , the problem was reformulated in terms of recursions on riccati variables to obtain weak - disorder expansions .
other studies have focused on random - scattering models on the cayley tree @xcite .
more recently , the interest in anderson localization on the cayley tree has been revived by the question of many - body localization @xcite , because the geometry of the fock space of many - body states was argued to be similar to a cayley tree @xcite .    in the present paper
, we consider the anderson tight binding model on the cayley tree already studied in @xcite and we study numerically the statistical properties of landauer transmission and eigenstates as a function of the disorder strength @xmath0 and of the number @xmath1 of generations .
we find that several probability distributions propagate as traveling waves with a fixed shape , as the number @xmath1 of tree generations grows .
the fact that traveling waves appear in disordered models defined on trees has been discovered by derrida and spohn @xcite on the specific example of the directed polymer in a random medium and was then found in various models @xcite . for the case of anderson localization on trees , the analysis of @xcite concerning
the distribution of the self - energy in the localized phase is actually a traveling wave analysis ( see appendix a for more details ) , although it is not explicitly mentioned in these terms in @xcite .
the fact that traveling waves occur has been explicitly seen in the numerical study of the transmission distribution for a random - scattering model on the cayley tree ( see fig . 8b of @xcite ) and has been found within the supersymmetric formalism @xcite . in the field of traveling waves and front propagation
( see the reviews @xcite ) , there exists an essential separation between two classes : in pulled fronts , the velocity is determined by the form of the tail of the front and thus by the appropriate linearized equation in the tail region , whereas in pushed fronts , the velocity is determined by the bulk properties and thus by the non - linear dynamics in the bulk region .
it turns out that for disordered systems defined on trees , the traveling waves that appear usually correspond to cases where it is the tail of the probability distribution that determines the velocity @xcite .
in particular for anderson localization on the cayley tree , the traveling wave propagation is also determined by the tails @xcite . from this traveling wave point of view
, the localization / delocalization anderson transition thus represents a traveling / non - traveling critical point .
such a traveling / non - traveling phase transition for a branching random walk in the presence of a moving absorbing wall has been studied recently in @xcite and we find here very similar behaviors in the critical region .
the paper is organized as follows . in section [ trans ] , we describe the statistical properties of landauer transmission . in section [ eigen ] , we discuss the statistical properties of eigenstates , as measured by the entropy and inverse participation ratios ( i.p.r . ) .
we summarize our conclusions in section [ conclusion ] . in appendix
a , we translate the analysis of ref .
@xcite concerning the distribution of the self - energy in the localized phase into a traveling - wave tail analysis for the landauer transmission discussed in section [ trans ] . in appendix b
, we explain how such a similar tail analysis can be performed in the delocalized phase for auxiliary variables that are closely related to inverse participation ratios .
finally , in appendix c , we recall the results of ref .
@xcite concerning the traveling / non - traveling phase transition for a branching random walk in the presence of a moving absorbing wall , since these results are used as a comparison in the text to understand the finite - size scaling properties in the critical region .
in quantum coherent problems , the most appropriate characterisation of transport properties consists in defining a scattering problem where the disordered sample is linked to incoming wires and outgoing wires and in studying the reflexion and transmission coefficients .
this scattering theory definition of transport , first introduced by landauer @xcite , has been much used for one - dimensional systems @xcite and has been generalized to higher dimensionalities and multi - probes measurements ( see the review @xcite ) . in dimension @xmath24 , the transfer matrix formulation of the schrdinger equation yields that the probability distribution of the landauer transmission becomes asymptotically log - normal @xcite , i.e. one has @xmath25 where @xmath26 represents the localization length and where @xmath27 is a sample - dependent random variable of order @xmath28 distributed with a gaussian law .
although it is often assumed that this log - normal distribution persists in the localized phase in dimension @xmath29 , recent numerical studies @xcite are in favor of the following scaling form for the logarithm of the conductance @xmath30 with exponents of order @xmath31 and @xmath32 @xcite , whereas @xmath33 from eq .
[ trans1d ] . for the cayley tree that we consider in this paper
, we will find below that the fluctuation exponent vanishes @xmath34 , and we will discuss the probability distribution of the variable @xmath27 . but let us first recall the appropriate scattering framework for the cayley tree @xcite .
we consider the anderson tight - binding model @xmath35 where the hopping between nearest neighbors @xmath36 is a constant @xmath37 and where the on - site energies @xmath38 are independent random variables drawn from the flat distribution @xmath39 the parameter @xmath0 thus represents the disorder strength .
starting at generation @xmath40 and ending at generation @xmath41 ( on the figure @xmath42 ) is attached to one incoming wire and to @xmath43 outgoing wires . in section [ trans ] , we discuss the properties of the total transmission @xmath44 where @xmath45 is the reflexion amplitude of the incoming wire , and @xmath46 the transmission amplitudes of the outgoing wires.,height=226 ]    we consider the scattering geometry introduced in @xcite and shown on fig .
[ figscattering ] : the finite tree of branching ratio @xmath47 is attached to one incoming wire at its root ( generation @xmath40 ) and to @xmath43 outgoing wires at generation @xmath41 .
one is interested into the eigenstate @xmath48 that satisfies the schrdinger equation @xmath49 inside the disorder sample and in the wires where one requires the plane - wave forms @xmath50 these boundary conditions define the reflexion amplitude @xmath45 of the incoming wire and the transmission amplitudes @xmath46 of the @xmath51 outgoing wires . to satisfy the schrdinger equation of eq .
[ schodinger ] within the wires with the forms of eq .
[ psiqoutside ] , one has the following relation between the energy @xmath52 and the wave vector @xmath53 @xmath54 to simplify the discussion , we will focus in this paper on the case of zero - energy @xmath55 and wave - vector @xmath56 @xmath57 because the zero - energy @xmath55 corresponds to the center of the band where the delocalizion first appears when the strength @xmath0 of the disorder is decreased from the strong disorder localized phase .    inside the cayley tree @xmath58 , the schrdinger equation of eq .
[ schodinger ] involves one ancestor denoted by @xmath59 and @xmath47 descendants denoted by @xmath60 @xmath61 whereas for the last generation @xmath41 there is only one ancestor and one descendant ( outgoing wire ) @xmath62    as explained in @xcite , it is convenient to introduce the riccati variables @xmath63 that represent the ratio of the wave function of two neighboring sites . on the outgoing wires ,
these riccati variables are fixed by eq .
[ psiqoutside ] to be @xmath64 the schrdinger equation of eq .
[ schodingerbord ] gives the first recursion @xmath65 whereas the schrdinger equation of eq .
[ schodingerinside ] gives the recursion inside the tree for @xmath58 @xmath66 on the other hand , the value of the riccati variable for the origin of the tree @xmath40 is fixed by the incoming wire of eq .
[ psiqoutside ] @xmath67 and thus the reflexion coefficient @xmath45 can be obtained via @xmath68 from the @xmath69 obtained via the recursion of eq .
[ riccatiinside ] . from the conservation of energy ,
the total transmission @xmath70 is related to the reflexion coefficient @xmath71 @xmath72    as explained in @xcite , the criterion for the localization / delocalization phases is then the following :    \(a ) if the riccati variable @xmath69 at the root of the tree converges towards a real random variable as @xmath10 , the reflexion is total @xmath73 and the transmission vanishes @xmath74 .
\(b ) if the riccati variable @xmath69 keeps a finite negative imaginary part as @xmath10 , the reflexion is only partial @xmath75 and the transmission remains finite @xmath76    we refer to @xcite for the results of a weak disorder expansion within this framework , and for a numerical monte - carlo approach to determine the mobility edge in the plane @xmath77 . here
we will instead study the statistical properties of the transmission @xmath70 at zero energy @xmath55 as a function of the disorder strength @xmath0 ( eq [ flat ] ) and of the number @xmath1 of generations .
but before discussing the disordered case , let us first describe the finite - size properties of the pure case .      in the pure case where all on - site energies vanish @xmath78 , all branches are equivalent and there is no dependence on @xmath79 . the recursions for riccati variables of eqs [ riccatioutgoing ] , [ riccatifirst ] and [ riccatiinside ] give @xmath80 i.e. one obtains the simple alternation between odd and even generations @xmath81 the reflexion coefficient @xmath45 thus vanishes exactly for any even size @xmath82 @xmath83 and the transmission coefficient in each branch @xmath84 reads @xmath85 ( note that for pure trees of uneven size , the reflexion would not vanish .
this is why in the disordered case , we will only consider trees of even sizes @xmath41 where in the corresponding pure case , the reflexion vanish exactly for any finite tree ) .
it is instructive to write now the corresponding wave function as a function of the generation n , @xmath86 this exponential decay of the wavefunction in this delocalized case is very peculiar to the tree geometry : it is imposed by the energy conservation and by the exponential growth of the number of sites with the generation @xmath87 .
if one wishes to study numerically real trees , one is limited to rather small number of generations of order @xmath88 ( see for instance the study @xcite based on exact diagonalization ) because the number of sites and thus of random energies grows exponentially in @xmath1 . from the point of view of convergence towards stable probability distributions via recursion relations , it is thus better to use the so - called pool method that will allow us to study much larger number of generations ( see below eq . [ pool5 ] ) .
the idea of the pool method is the following : at each generation , one keeps the same number @xmath89 of random variables to represent probability distributions . within our present framework , the probability distribution @xmath90 of the riccati variables @xmath91 at generation @xmath87
will be represented by a pool of @xmath89 complex values @xmath92 .
it is convenient from now on to change the notation @xmath93 with respect to fig .
[ figscattering ] so that @xmath40 now corresponds to the contacts with the outgoing wires ( see fig .
[ figscattering ] ) .
the numerical results presented below have been obtained with the following procedure :    \(i ) the random variables of the initial pool are given by ( eq . [ riccatifirst ] ) @xmath94 where @xmath95 are independent random variables drawn with eq .
[ flat ] .
\(ii ) from the random variables @xmath96 of the pool at generation @xmath97 , the pool at generation @xmath87 is constructed as follows . for each @xmath98 ,
one generates a new random energy @xmath99 with the law of eq .
[ flat ] and one draws @xmath47 random indices @xmath100 among the pool of generation @xmath101 to construct the variable @xmath102 as follows(eq . [ riccatiinside ] ) @xmath103    \(iii ) from these pools of riccati variables , one may compute via eqs [ reflexioncoef ] and [ deftotaltrans ] a pool of total transmission @xmath104 for trees of @xmath87 generations using @xmath105    for the anderson model on the bethe lattice , the pool method has been already used , in particular in @xcite with a pool @xmath106 with a number @xmath107 of generations , in @xcite with pools up to @xmath108 , and in @xcite with pools up to @xmath109 with a number @xmath110 of generations .
the pool method is also very much used for disordered systems on hierarchical lattices ( see for instance@xcite ) ,    in the remaining of this section , we present the numerical results obtained with a pool of size @xmath111 we have also results for a pool of size @xmath112 with a number of generations @xmath113 to see how the results change with the pool size .
however , the number of generations @xmath114 for this bigger pool @xmath115 has turned out to be less precise in the critical region .
all figures shown below thus corresponds to data obtained with the pool of size @xmath116 .    as is usual with the pool method @xcite , the location of the critical point depends on the pool , i.e. on the discrete sampling with @xmath89 values of probability distributions .
it is expected to converge towards the thermodynamic critical point only in the limit @xmath117 ( see the discussion of section [ fspool ] in appendix b ) .
nevertheless , for each given pool , the critical behaviors with respect to this pool - dependent critical point usually allows a good measure of critical exponents @xcite .
for instance , for the pool of size @xmath116 , the critical value @xmath118 of the disorder strength where the localization - delocalization occurs for the landauer transmission is of order @xmath119 for the pool of size @xmath115 , we find that it is higher and of order @xmath120 .
this rather important shift of the pseudo - critical point with the pool size which has already been seen in @xcite , seems to be due to the very slow logarithmic convergence of traveling wave velocity in the presence of cut - off ( see the discussion of section [ fspool ] in appendix b ) .
so we stress that here , in contrast to @xcite , our goal is not to determine the true thermodynamic mobility edge @xmath121 , but instead to understand the critical behaviors of the finite - pool results with respect to the pool - dependent critical point of eq .
[ critipoole5 ] .
we will first discuss the behavior of the typical transmission @xmath122 defined by @xmath123 as a function of the number @xmath1 of generations and disorder strength @xmath0 , before we turn to the distribution around this typical value .       in the localized phase @xmath125 : ( a ) linear decay of @xmath126 as a function of the number @xmath1 of generations ( see eq .
[ transloc ] ) .
( b ) behavior of the slope @xmath127 as a function of the disorder strength @xmath0 : it vanishes linearly @xmath128 with @xmath129 ( see eq .
[ critipoole5 ] ) and @xmath6 ( see eq .
[ xiloc ] ) , title="fig:",height=226 ]   in the localized phase @xmath125 : ( a ) linear decay of @xmath126 as a function of the number @xmath1 of generations ( see eq . [ transloc ] ) .
( b ) behavior of the slope @xmath127 as a function of the disorder strength @xmath0 : it vanishes linearly @xmath128 with @xmath129 ( see eq .
[ critipoole5 ] ) and @xmath6 ( see eq .
[ xiloc ] ) , title="fig:",height=226 ]    in the localized phase , one expects that the typical transmission @xmath122 defined by eq .
[ defttyp ] decays exponentially with the number @xmath1 of generations @xmath130 where @xmath131 represents the localization length that diverges at the delocalization transition @xmath132 .
we show on fig .
[ figtransloc ] our numerical results for the pool of size @xmath116 ( eq . [ pool5 ] ) : the exponential decay with @xmath1 of eq .
[ transloc ] is shown on fig .
[ figtransloc ] ( a ) for various disorder strength @xmath0 . the corresponding slope @xmath127 is shown as a function of @xmath0 on fig . [ figtransloc ] ( b )
: we find that this slope vanishes linearly in @xmath133 , in agreement with the exact result @xcite @xmath134 and in agreement with figure 7 of ref .
@xcite concerning a random - scattering model on the cayley tree .       of the infinite tree in the delocalized phase : ( a ) @xmath137 as a function of the disorder strength @xmath0 ( b ) same data in a log - log plot to measure the exponent of the essential singularity of eq .
[ essentialsingularitytransdeloc ] : @xmath138 as a function of @xmath139 : the slope is of order @xmath140 , title="fig:",height=226 ]   of the infinite tree in the delocalized phase : ( a ) @xmath137 as a function of the disorder strength @xmath0 ( b ) same data in a log - log plot to measure the exponent of the essential singularity of eq . [ essentialsingularitytransdeloc ] : @xmath138 as a function of @xmath139 : the slope is of order @xmath140 , title="fig:",height=226 ]    in the delocalized phase , the typical transmission remains finite in the limit where the number of generations @xmath1 diverges @xmath141 as shown on fig .
[ figtransdeloc ] , we measure the following essential singularity behavior of the typical transmission @xmath142    the presence of essential singularities in transport properties near the localization transition on the bethe lattice has been found in @xcite via the supersymmetric formalism ( see also @xcite where similar results have been obtained for the case where anderson tight - binding model is replaced by a non - linear @xmath23-model ) . in particular , eq .
71 of ref .
@xcite ) states that the leading critical behavior of the diffusion constant is given by : @xmath143 ( eq .
71 of ref .
we note that the exponent in this essential singularity is @xmath144 instead of the exponent of order @xmath145 that we measure .
the reason for this difference could be that the results of ref .
@xcite are based on the computation of the _ disorder - averaged two - point density - density correlation function _ ( see eq . 2 of ref .
@xcite ) , whereas our numerical results concern the _ typical value of landauer transmission _ , and not the disorder - averaged transmission which is expected to be governed by rare events ( see below the section [ distritrans ] ) .
also it is not clear to us what exactly represents the diffusion constant computed in @xcite , because it is known that on the bethe lattice the dynamics is not diffusive but ballistic in the delocalized phase @xcite , and that more generally the random walk on the bethe lattice is not diffusive because the tree - geometry induces an effective bias away from the origin ( see for instance @xcite and references therein ) . as a consequence , some implicit reinterpretation of the bethe lattice seems to underlie the statements of ref . @xcite
that makes difficult a precise comparison with our present numerical results .
if there exists some finite - size scaling in the critical region for the typical landauer transmission of the form @xmath146 the matching of our results in the localized phase ( see eq .
[ transloc ] ) and in the delocalized phase ( eq . [ essentialsingularitytransdeloc ] ) requires a finite - size correlation length exponent @xmath147 of order @xmath148 in another traveling / non - traveling phase transition studied in ref @xcite ( see the summary in appendix c ) , it has been obtained that the finite - size scaling exponent @xmath149 is determined by the relaxation rate towards the finite value in the non - traveling phase .
we have thus studied the relaxation length towards the finite value in the delocalized phase .
we find that our data for @xmath150 are compatible with the form @xmath151 where @xmath152 is a random stationary process as a function of @xmath1 .
we find that its autocorrelation function is exponential @xmath153 and we measure that the relaxation length @xmath154 diverges with an exponent @xmath155 of the order of the exponent @xmath147 of eq .
[ nufstrans ] .
we thus obtain that the critical properties are qualitatively similar to the critical properties described in ref @xcite ( see the summary in appendix c ) : the traveling phase is characterized by a velocity that vanishes linearly , but the finite size scaling is governed by the relaxation length towards the asymptotic finite value of the non - traveling phase .
exactly at criticality , we thus expects the following stretched exponential decay of the typical transmission @xmath156 where the exponent @xmath157 is related to the other exponents by ( see the scaling relations of eqs [ kapparelation ] and [ nutyprelation ] in appendix c ) @xmath158 from our previous estimate of the exponent @xmath159 , this would correspond to the numerical value @xmath160 we have not been able to measure this stretched exponential behavior exactly at criticality from our data , because a precise measure of the exponent @xmath157 would require to be exactly at the critical point .       of the probability distribution @xmath161 of the logarithm of landauer transmission : ( a ) in the delocalized phase ( here @xmath162 ) , the probability distribution @xmath161 does not move with @xmath1 , and vanishes with a discontinuity at the boundary @xmath163 .
( b ) in the localized phase ( here @xmath164 ) , the probability distribution @xmath161 moves with @xmath1 as a traveling wave of fixed shape .
, title="fig:",height=226 ]   of the probability distribution @xmath161 of the logarithm of landauer transmission : ( a ) in the delocalized phase ( here @xmath162 ) , the probability distribution @xmath161 does not move with @xmath1 , and vanishes with a discontinuity at the boundary @xmath163 .
( b ) in the localized phase ( here @xmath164 ) , the probability distribution @xmath161 moves with @xmath1 as a traveling wave of fixed shape .
, title="fig:",height=226 ]    up to now , we have only discussed the behavior of the typical transmission of eq .
[ defttyp ] as a function of @xmath1 and @xmath0 .
we now turn to the probability distribution of the logarithm of the transmission around its averaged value , i.e. we consider the distribution of the relative variable @xmath165 we find that as @xmath10 , this variable remains finite not only in the delocalized phase where @xmath166 is finite , but also in the localized phase where @xmath167 decays linearly in @xmath1 ( eq . [ transloc ] ) .
this means that in the localized phase , the probability distributions @xmath168 actually propagates as a traveling wave with a fixed shape around its moving center @xmath167 as shown on fig .
[ fighistotrans ] ( b ) .
this phenomenon has already been seen for random scattering models on the bethe lattice ( see fig .
8b of ref .
this is in contrast to the broadening with @xmath169 observed in low dimensions @xmath170 ( see eq . [ transd ] ) : the cayley tree thus corresponds to @xmath34 in eq .
[ transd ] .     of
@xmath171 as @xmath0 varies in the localized phase for various disorder strengths @xmath172 : the exponential decay as @xmath173 corresponds to slopes of order @xmath174 , @xmath175 , @xmath176 , @xmath177 .
( b ) width @xmath178 of the probability distribution of @xmath179 as a function of @xmath0 : it remains finite both in the localized phase and in the delocalized phase , but it presents a cusp singularity at criticality . , title="fig:",height=226 ]   of @xmath171 as @xmath0 varies in the localized phase for various disorder strengths @xmath172 : the exponential decay as @xmath173 corresponds to slopes of order @xmath174 , @xmath175 , @xmath176 , @xmath177 .
( b ) width @xmath178 of the probability distribution of @xmath179 as a function of @xmath0 : it remains finite both in the localized phase and in the delocalized phase , but it presents a cusp singularity at criticality .
, title="fig:",height=226 ]    the fact that the shape is fixed can be used numerically to measure more precisely the tails of this probability distribution by accumulating data over iterations : we show on fig .
[ figtranswidth ] ( a ) the histograms of @xmath27 obtained by this procedure for various disorder strength @xmath0 .
an essential property of this distribution is the right exponential tail @xmath180 in terms of the rescaled transmission @xmath181 ( see eq .
[ variableu ] ) , this corresponds to the power - law decay @xmath182 we find that the selected exponent @xmath183 slightly grows as the disorder strength @xmath0 decreases , from a value of order @xmath177 for the strong disorder @xmath184 towards a value of order @xmath174 near criticality ( see figure [ figtranswidth ] a ) .
the velocity of the traveling wave propagation of the whole distribution is actually determined by this power - law tail , as explained in detail in appendix a. in particular , the critical value @xmath185 has been predicted in @xcite ( see eq .
[ betacdemi ] in appendix a )
. an important consequence of the power - law tail of eq .
[ tailpowert ] with @xmath186 is that all integer moments @xmath187 of the landauer transmission will be governed by rare events .
finally , we show on fig .
[ figtranswidth ] ( b ) the width of the probability distribution of @xmath188 as a function of @xmath0 : it remains finite both in the localized phase and in the delocalized phase , but it presents a cusp singularity at criticality because in the delocalized phase , the distribution presents a singularity at finite distance from the typical value , as is clearly visible on fig . [ fighistotrans ] (
a ) : the transmission is bounded by @xmath189 , i.e. the histogram of @xmath188 presents a discontinuity at @xmath163 .
the landauer transmission studied in the previous section is of course the most appropriate observable to characterize the transport properties and to find the transition between the conducting / non - conducting phases .
however , one expects that these transmission properties that emerge when the disordered sample is linked to incoming and outgoing wires are related to the nature of eigenstates of the disordered sample in the absence of these external wires ( see fig .
[ figtree ] ) to determine whether a normalized eigenstate @xmath190 is localized or delocalized , the usual parameters are the inverse participation ratios ( i.p.r . ) @xmath191 of arbitrary power @xmath192 .
another important quantity to characterize the spatial extent of an eigenstate is its entropy @xmath193 the normalization condition yields the idendity @xmath194 , and the entropy corresponds to @xmath195 on an hypercubic lattice of size @xmath169 in dimension @xmath196 , containing @xmath197 sites , the delocalized phase is characterized by the following typical behaviors @xmath198 and @xmath199 these scalings are the same for pure homogeneous eigenfunctions with @xmath200 for all sites @xmath201 . at criticality , the scalings of the i.p.r .
involve a whole series of non - trivial exponents and the wave function is said to be multifractal ( see the review @xcite ) . finally in the localized phase , the i.p.r .
and the entropy are finite .    in the present section
, we discuss the statistical properties of these i.p.r . and
entropy for zero - energy eigenstates for the cayley tree of fig .
[ figtree ] as a function of the number @xmath1 of tree generations and of the disorder strength @xmath0 .     where each interior site has @xmath202 neighbors .
on the figure the tree ends at generation @xmath42 . , height=226 ]      we refer to @xcite where it is explained how eigenstates of finite trees can be constructed and how the density of states can be obtained .
note that on a tree , boundary sites dominate so that one needs a substraction procedure to obtain the appropriate bulk density of states . in the following , we are not interested into the density of states , which does not contain any information on the localized / delocalized nature of the spectrum .
we wish to study instead the spatial properties of eigenstates of zero - energy @xmath55 ( center of the band ) .
the schrdinger equation [ schodingerinside ] yields as before the recursion of eq .
[ riccatiinside ] for the riccati variables defined in eq .
[ riccati ] .
the difference with the scattering case is now in the boundary conditions @xmath203 that replace eq .
[ riccatifirst ] . as a consequence ,
the riccati variables are now real ( and not complex ) .
the energy @xmath55 will indeed be an eigenstate only if the schrdinger equation is also satisfied at the center of the finite tree at the center site @xmath204 since the on - site energy @xmath205 is a random variable drawn with some distribution , we may consider that we choose @xmath205 to satisfy eq .
[ riccaticenter ] to obtain a typical eigenstate of zero energy and to study its spatial properties .
let us first describe zero - energy eigenstate in the pure case to stress the peculiarities of the tree geometry .
in the pure case where all on - site energies vanish @xmath78 , the zero - energy eigenstate with radial symmetry reads @xmath206 where the root amplitude @xmath207 is determined by the normalization condition @xmath208 \label{normapur}\end{aligned}\ ] ] this state is the analog of the scattering state described in section [ puretrans ] .
in particular , one sees again an exponential decay of the wavefunction in this delocalized case which is very peculiar to the tree geometry : in the normalization condition of eq .
[ normapur ] , this exponential decay is exactly compensated by the exponential growth of the number of sites at generation @xmath209 , so that all generations @xmath210 carry exactly the same weight in the normalization condition , and the root weight @xmath211 vanishes as @xmath212 in the thermodynamic limit , where @xmath82 is the number of generations of the tree .    for this state , the usual inverse participation ratio of eq .
[ defiq ] for @xmath213 reads @xmath214 \label{ipppur}\end{aligned}\ ] ] more generally for @xmath215 , one obtains the following decay with the number of generations @xmath216 in terms of the total number of sites @xmath217 this corresponds to @xmath218 this behavior is thus very anomalous with respect to the corresponding decay as @xmath219 for the i.p.p .
of pure states in finite dimension @xmath196 ( see eq .
[ delociqdimd ] ) .
moreover , since the eigenfunction normalization yields the idendity @xmath194 , one sees that there exists some discontinuity as @xmath220 .
in particular , the entropy of eq .
[ defs ] which can usually be computed as the derivative of eq .
[ sderiiq ] should be computed directly here , and one obtains @xmath221 \label{entropiecaspurcalcul}\end{aligned}\ ] ] using eq [ normapur ] , one obtains the following behavior for large @xmath1 @xmath222 in terms of the total number of sites of eq .
[ calncaspur ] , the entropy is proportional to the logarithm of the number of sites @xmath223 that should be compared with the growth as @xmath224 for pure eigenstates in dimension @xmath196 ( see eq .
[ delocsdimd ] ) .    in conclusion , the tree geometry induces very anomalous scalings in @xmath225 for the logarithm of i.p.r . of pure states
[ ipppurcaln ] ) with respect to the case of finite dimension @xmath196 , whereas the entropy is well behaved in @xmath226 ( see eq .
[ entropiecaspurfin ] ) . in the disordered case ,
we expect to observe similar behaviors , whereas the i.p.r . and
the entropy will remain finite in the localized phase .      to compute the inverse participation ratio @xmath227 ( eq . [ defiq ] ) via recursion , one needs to introduce besides the riccati variable of eq .
[ riccati ] the auxiliary variables defined by @xmath228 where the sum is over all sites @xmath229 that are descendants of the site @xmath230 of the tree .
these variables satisfy the following recurrence inside the tree @xmath231
@xmath232 \label{reccq}\end{aligned}\ ] ] and the initial conditions at the boundaries @xmath233 at the central root , on needs to impose the normalization @xmath234 \label{normapsi2}\end{aligned}\ ] ] that determines the weight @xmath235 of the root in terms of the variables @xmath236 of the branches .    to compute the inverse participation ratio of parameter @xmath192 of eq .
[ defiq ] , one needs the variables @xmath237 together with the variables @xmath236 @xmath238 = \frac{\left [ 1 + \displaystyle \sum_{m=1}^{k+1 }   c^{(q ) } ( 1,m)\right]}{\left [ 1 + \displaystyle \sum_{m=1}^{k+1 }   c^{(1 ) } ( 1,m)\right]^q } \label{iqtree}\end{aligned}\ ] ]      to compute recursively the entropy of eq .
[ defs ] , one needs similarly to introduce the auxiliary variable @xmath239 where the sum over @xmath229 denotes the sum over the site @xmath230 and all its descendants : @xmath240 thus represents the entropy for the branch containing @xmath230 and its descendants .
the initial conditions at the boundaries are simply @xmath241 and the recursion inside the tree @xmath231 can be written as @xmath242 where the first term represents the weighted contribution of the branches entropies and where the second term represents the mixing entropy @xmath243 with the weights @xmath244 normalized to @xmath245 . at the central site of the tree ,
one uses the same formula but with @xmath246 branches instead of @xmath47 branches .        to study the statistical properties of eigenstates , we have used again the pool method explained in [ poolmethod ] .
the only difference is that for the transmission , we have followed the recursions for the complex riccati variables , whereas here we follow the recursions for the real riccati variable @xmath91 and for the auxiliary variables @xmath236 , @xmath247 and @xmath23 described above .       of a normalized zero - energy eigenstate : ( a ) evolution of the typical entropy @xmath248 with the number @xmath1 of generations :
it grows linearly in @xmath1 in the delocalized phase , whereas it remains finite in the localized phase .
( b ) asymptotic width @xmath249 of the relative variable @xmath250 as a function of the disorder strength @xmath0 : the width converges to @xmath251 ( as @xmath212 ) in the delocalized phase , whereas it is finite in the localized phase , title="fig:",height=226 ]   of a normalized zero - energy eigenstate : ( a ) evolution of the typical entropy @xmath248 with the number @xmath1 of generations : it grows linearly in @xmath1 in the delocalized phase , whereas it remains finite in the localized phase .
( b ) asymptotic width @xmath249 of the relative variable @xmath250 as a function of the disorder strength @xmath0 : the width converges to @xmath251 ( as @xmath212 ) in the delocalized phase , whereas it is finite in the localized phase , title="fig:",height=257 ]    we first consider how the eigenstate entropy defined in eq .
[ defs ] evolves with the number @xmath1 of generations . as shown on fig .
[ figentropy ] ( a ) , the typical value grows linearly in @xmath1 in the delocalized phase @xmath252 where the factor @xmath253 varies smoothly with @xmath0 and does not vanish continuously near the critical point . in this delocalized phase , we moreover find that the width of the relative variable @xmath254 decays to zero as @xmath10 @xmath255^{1/2}_{w < w_c }    \opsimeq_{n \to \infty } \frac{1}{n }   \label{entropywidthdeloc}\end{aligned}\ ] ] in the localized phase on the contrary , the typical value and the width of the relative variable @xmath254 remain finite as @xmath256 ( see fig .
[ figentropy ] ( a ) and ( b ) ) .     in the localized phase ( a ) @xmath257 as a function of the disorder strength @xmath0 near the critical point .
( b ) same data as a function of @xmath258 to measure the exponent of the power law of eq .
[ critientropy64loc ] : the slope is of order @xmath259.,title="fig:",height=226 ]   in the localized phase ( a ) @xmath257 as a function of the disorder strength @xmath0 near the critical point .
( b ) same data as a function of @xmath258 to measure the exponent of the power law of eq .
[ critientropy64loc ] : the slope is of order @xmath259.,title="fig:",height=226 ]    as shown on fig .
[ fig64loc ] , we find that the typical entropy of an infinite tree of the localized phase diverges with the following power law near criticality @xmath260 the entropy measures the size of the region where the weight @xmath261 is concentrated , whereas the localization length @xmath262 describes the far exponential decay of the transmission @xmath70 ( eq . [ transloc ] ) .
we thus conclude that in the localized phase , besides the localization length @xmath262 known since @xcite , there exists a larger diverging length scale @xmath263 that characterizes the size where the weight @xmath261 is concentrated .      as explained above , the root weight @xmath264 and of the i.p.r .
@xmath265 are determined by the auxiliary variables @xmath237 introduced in eq .
[ defcq ] .
we thus expect that the localized and delocalized phases correspond to the following behaviors for the auxiliary variables @xmath237    \(i ) in the localized phase , the auxiliary variables @xmath237 will remain finite random variables as the number of generations diverge @xmath10 .
then the root weight @xmath266 and the i.p.r . @xmath227 remain finite as @xmath10 .
\(ii ) in the delocalized phase , the auxiliary variables @xmath237 will instead grow exponentially with the number @xmath1 of generations with some lyapunov exponents @xmath267 defined by @xmath268 the corresponding typical behaviors of the root weight and of the i.p.r . then reads @xmath269 and @xmath270    we find numerically that the probability distributions of auxiliary variables @xmath271 move as traveling waves of velocity @xmath272 ( see eq .
[ lyapunovcq ] ) with fixed shape ( see figure [ fighistoc1c2 ] of appendix b where the fixed shape is shown for @xmath273 and @xmath213 ) .
we explain in appendix b how the lyapunov exponents @xmath272 can be determined via a tail analysis that yields the identity ( see eq .
[ selectionqidentity ] ) @xmath274 so that the logarithm of the i.p.r .
@xmath227 decays slower than linearly in @xmath1 , because the coefficient in eq .
[ lyapunoviq ] exactly vanishes .
physically , one could expect the decay to be of order @xmath275 as in eq .
[ ipppurcaln ] concerning pure states on the cayley tree , but we have not been able to measure the behavior of the i.p.r . , because within the pool method that we use , it turns out that the pool - dependent critical point is not the same for the variables @xmath236 and for @xmath247 ( see more details in section [ fspool ] ) , so that it does not seem easy with the pool method to extract reliable results concerning @xmath265 .    in the remaining of this section
, we thus focus on the statistical properties of the root weight @xmath235 that is determined by the auxiliary variable @xmath236 via eq .
[ normapsi2 ] .
the fact that the pool - dependent critical point for @xmath236 is the same as the critical point found for the landauer transmission can be understood from the tail analysis of the appendices a and b that involve exactly the same integral kernel .
: the figure shows that the slope @xmath276 vanishes linearly @xmath277 ( see eq .
[ xideloc ] with @xmath6 ) .
( b ) width @xmath278 of the relative variable @xmath279 in the limit @xmath256 as a function of the disorder strength @xmath0 : it remains finite both in the localized and delocalized phase , but it presents a cusp at the critical point @xmath118 . , title="fig:",height=188 ]   : the figure shows that the slope @xmath276 vanishes linearly @xmath277 ( see eq .
[ xideloc ] with @xmath6 ) .
( b ) width @xmath278 of the relative variable @xmath279 in the limit @xmath256 as a function of the disorder strength @xmath0 : it remains finite both in the localized and delocalized phase , but it presents a cusp at the critical point @xmath118 .
, title="fig:",height=188 ]    in the delocalized phase , the root weight @xmath266 presents the typical decay of eq .
[ lyapunovpsi0 ] .
we show on fig .
[ fig44deloc ] ( a ) that the corresponding lyapunov exponent @xmath276 vanishes linearly in the critical region @xmath280 i.e. we find the same exponent @xmath6 as for the divergence of the correlation length of the localized phase ( see eq . [ nueq1 ] ) .
this can be understood from the tail analysis presented in the appendices a and b.    we find that the width of the relative variable @xmath279 remains finite in the limit @xmath256 both in the localized phase and in the delocalized phase : its behavior as a function of the disorder strength @xmath0 presents a cusp at @xmath118 as shown on fig .
[ fig44deloc ] ( b ) .     of an infinite tree in the localized phase : ( a )
@xmath281 as a function of the disorder strength @xmath0 ( b ) @xmath282 as a function of @xmath258 to measure the exponent of the essential singularity of eq .
[ essentialsingularity44loc ] : we measure asymptotically a slope of order @xmath283.,title="fig:",height=226 ]   of an infinite tree in the localized phase : ( a ) @xmath281 as a function of the disorder strength @xmath0 ( b ) @xmath282 as a function of @xmath258 to measure the exponent of the essential singularity of eq .
[ essentialsingularity44loc ] : we measure asymptotically a slope of order @xmath283.,title="fig:",height=226 ]    in the localized phase , the root weight remains finite as @xmath10 . as show on fig .
[ fig44loc ] , we measure the following essential singularity @xmath284    as in section [ fsstrans ] , we now discuss the finite - size scaling in the critical region .
if there exists some finite - size scaling in the critical region for the typical root weight , the matching of our results in the delocalized phase ( see eq .
[ lyapunovpsi0 ] and [ xideloc ] ) and in the localized phase ( eq . [ essentialsingularity44loc ] ) requires a finite - size correlation length exponent @xmath285 of order @xmath286 we have performed the analysis described in section [ fsstrans ] and we measure that the relaxation length @xmath154 towards the finite value of eq . [ essentialsingularity44loc ] diverges with an exponent @xmath287 of the order of the exponent @xmath285 of eq .
[ nufszero ] .
we thus obtain that the critical properties are again qualitatively similar to the critical properties described in ref @xcite ( see the summary in appendix c ) : the traveling phase is characterized by a velocity that vanishes linearly , but the finite size scaling is governed by the relaxation length towards the asymptotic finite value of the non - traveling phase .
exactly at criticality , we thus expects the following stretched exponential decay of the typical transmission @xmath288 where the exponent @xmath289 is related to the other exponents by ( see the scaling relations of eqs [ kapparelation ] and [ nutyprelation ] in appendix c ) @xmath290 from our previous estimate of the exponent @xmath291 , this corresponds to the numerical value @xmath292 again , as explained after eq .
[ transrhonume ] , we have not been able to measure this stretched exponential behavior exactly at criticality from our data , because a precise measure of the exponent @xmath293 would require to be exactly at the critical point .
in summary , for the anderson model on the bethe lattice , we have studied numerically the statistics of the landauer transmission @xmath2 and the statistics of eigenstates at the center of the band @xmath55 , as a function of the disorder strength @xmath0 and the number @xmath1 of generations .
we have shown that both the localized phase and the delocalized phase are characterized by the traveling wave propagation of some probability distributions . in the text ,
we have presented detailed numerical results , and in appendices a and b , we have explained how the velocities of the traveling waves are determined by the tails , via the properties of the integral kernel introduced in @xcite .
the anderson transition then corresponds to a traveling / non - traveling critical point for these traveling waves , and the critical properties obtained are very similar to the traveling - wave phase transition studied in @xcite : ( i ) the finite value of the non - traveling phase presents an essential singularity ( ii ) the relaxation length towards this essential singularity determines the finite - size scaling in the critical region ( iii ) the finite - size correlation length exponent @xmath149 is different from the value @xmath6 that govern the vanishing of the velocity in the traveling - wave phase .
we thus hope that in the future , these properties for anderson localization on the cayley tree will be better understood by extending the methods that have been developed recently for the class of the fisher - kpp traveling waves ( see @xcite and references therein ) .
how results obtained on the bethe lattice are related to properties of anderson localization in finite dimension @xmath196 is of course a difficult question . during the history of localization ,
many values for the upper critical dimension have been proposed such as @xmath294 . even within the supersymmetric community , there seems to be different interpretations .
mirlin and fyodorov @xcite consider that the essential singularities that appear on the bethe lattice , are directly related to the exponential growth of the number of sites with the distance , and will become conventional power - law behaviors as soon as @xmath295 , so that the upper critical dimension is actually @xmath296 . on the contrary
, efetov @xcite argues that the results obtained for the bethe lattice are relevant to anderson transition in high dimension @xmath196 when reinterpreted within the so - called effective medium approximation. for the landauer transmission , the upper critical dimension @xmath297 can be defined as the dimension where the exponent @xmath298 concerning the width of the sample - to - sample distribution in the localized phase ( see eq.[transd ] ) vanishes @xmath299 .
this means that for @xmath300 , the probability distribution would travel as a traveling wave in the whole localized phase . to determine whether @xmath296 or @xmath301 , it would be thus interesting to understand whether the traveling wave propagation with a fixed shape is possible only on trees or whether it can also occur in sufficiently high dimension @xmath196 .    in low dimensions ,
the probability distribution of the logarithm of the transmission is known to broaden with some exponent @xmath302 ( see eq .
[ transd ] ) in the localized phase , whereas the probability distribution of the logarithm of i.p.p .
is known to shrink with the system size in the delocalized phase @xcite ( see also fig .
( 4a ) of @xcite where the same phenomenon occurs for the directed polymer in @xmath303 dimensions ) . however traveling wave propagations of the probability distribution @xmath304 of inverse participation ratios
@xmath227 have actually been found in finite dimensions but only _ exactly at criticality_. this phenomenon has been first obtained for the power - law random band matrix model @xcite and has been then observed for the anderson model in dimensions @xmath305 @xcite : the motions of the typical values @xmath306 with @xmath307 determine the multifractal spectrum , whereas the reduced variables @xmath308 keep fixed distributions presenting power - law tails with @xmath192-dependent exponents .
we have actually observed the same behavior for the directed polymer in @xmath303 dimensions exactly at the localization / delocalization transition ( see fig .
( 3a ) and ( 3b ) of ref .
@xcite ) . in all these cases
, it would be thus very nice to better understand the relations between the tails of the distributions and the motions of the typical values , and to identify the wave equation that underlies these traveling waves .
finally , our work raises once again the question whether anderson transitions are characterized by several length scales that diverge with various @xmath19 exponents . here for the cayley tree , besides the exponent @xmath6 that has been predicted for a long time @xcite , we have found other diverging lengths with the following physical meanings .
( i ) in the localized phase , besides the exponent @xmath6 associated to the divergence of localization length @xmath262 that describes the far exponential decay of the transmission @xmath70 ( eq . [ transloc ] ) , we have found that the critical behavior of the entropy is governed by another diverging length scale with @xmath309 ( eq . [ xisloc ] ) that characterizes the size where the weight @xmath261 is concentrated .
( ii ) in the delocalized phase , we have found that the landauer transmission reaches its asymptotic value after a length diverging as @xmath310 ( eq . [ nufstrans ] ) that governs the finite - size scaling of the transmission in the critical region . even if the precision of numerically measured critical exponents
can be always discussed , we feel nevertheless that our numerical results are not compatible with the existence of the single exponent @xmath6 .
note that for the directed polymer on the cayley tree , we have also found previously that the critical properties involve two exponents @xmath311 and @xmath312 @xcite . whether there exists various exponents @xmath19 for anderson model in finite dimension
@xmath196 has been debated for a long time .
the majority of papers on anderson localization seem in favor of a single @xmath19 ( see the reviews @xcite ) , but it seems to us that numerical papers actually study always the same type of observables . in particular , we have not been able to find numerical studies concerning the entropy of eigenstates . on the theoretical side , we are not aware of many theories in favor of various @xmath19 exponents , except ( i ) the supersymmetric studies of refs @xcite that involve two different diverging length scales , called respectively the localization length and the phase coherence length .
( ii ) the pseudo - delocalization transition of the random hopping model in dimension @xmath24 , that is characterized by the two exactly known correlation lengths exponents @xmath22 and @xmath311 ( see @xcite and references therein ) , as in the other models described by the same strong disorder fixed point ( see the review @xcite ) .
we feel that if localization models exhibit several @xmath19 exponents both on the cayley tree which represents some mean - field @xmath313 limit and in some models in dimension @xmath24 , the possibility of various @xmath19 for anderson transition in dimension @xmath314 should be reconsidered by studying in detail the statistics of various observables .
in this appendix , we translate the analysis of ref .
@xcite concerning the distribution of the self - energy in the localized phase for the traveling - wave propagation of the landauer transmission discussed in section [ trans ] .
the recursion of eq .
[ riccatiinside ] for the complex riccati variable @xmath315 reads more explicitly in terms of its real and imaginary parts @xmath316 with @xmath317-\infty,+\infty[$ ] and @xmath318 .    @xmath319    where the random on - site energies @xmath320 are drawn from the flat distribution of eq .
the landauer transmission of eq .
[ transresume ] then reads @xmath321    in the delocalized phase , this recursion of eq .
[ ricattiinsidexy ] is difficult to analyze because one has to find the joint distribution of the two finite variables @xmath322 that remain stable upon iteration . in the localized phase however ,
the problem is simpler @xcite as we now recall .      in the localized phase , if the imaginary part @xmath323 converge towards zero exponentially in @xmath87 as @xmath324 where @xmath325 and where @xmath326 remains a finite random variable upon iteration , one has then to study the simpler recurrence for large @xmath87 @xcite @xmath327 note that the form of eq .
[ ynscaling ] corresponds to a traveling wave of velocity @xmath328 for the variable @xmath329 @xmath330    the real part @xmath331 satisfies now a closed recurrence independent of the @xmath323 , and its stable distribution @xmath332 satisfies the closed equation @xmath333 \label{petoilex}\end{aligned}\ ] ] an important property of this distribution is that it presents the power - law tail @xmath334 because whenever one of the @xmath47 variables @xmath335 on the right handside of eq .
[ petoilex ] is close to @xmath251 , the variable @xmath336 of the left - handside is large with @xmath337 .
the stable joint distribution @xmath338 satisfies @xmath339 \delta \left [ y - \sum_{m=1}^k \frac { y_m } { e^{-v } x_m^2 } \right ] \label{ricattistab}\end{aligned}\ ] ]      the idea of @xcite is to look for the power - law tail in @xmath340 that is compatible with the recursion equation , i.e. one assumes @xmath341 with some exponent @xmath342 . note that , in the traveling wave language of eq .
[ lnynscaling ] , this is equivalent to look for solutions with an exponential tail @xmath343 for the variable @xmath344 .
so this corresponds to the usual exponential tail analysis of fronts @xcite , except for the following difference : in usual studies of propagation into unstable phases , it is the forward tail that has an exponential decay and that determines the velocity , whereas in our present case , it is the backward tail that determines the propagation .    in laplace
transform with respect to @xmath340 , the power - law decay of eq .
[ tail ] corresponds to the following singular expansion near the origin @xmath345 equation [ ricattistab ] becomes for the laplace transform @xmath346 @xmath347 \label{ricattistablaplace}\end{aligned}\ ] ] using the expansion of eq .
[ devlaplace ] , one obtains that the function @xmath348 has to satisfy the eigenvalue equation @xmath349 where the function @xmath350 represents the stationary distribution of the variable @xmath351 @xmath352 \label{fonctionq}\end{aligned}\ ] ] ( see eqs 6.5 and 6.6 in @xcite ) .
the tail analysis thus leads to the eigenvalue problem the integral kernel appearing in eq .
[ eqpropreax ] @xmath353 for each @xmath354 and @xmath0 , the integral kernel is positive , and thus one expects some continuous analog of the perron - froebenius theorem : the iteration of the integral kernel will converge towards a positive eigenvector @xmath355 that is associated to the maximal eigenvalue @xmath356 of the kernel .
for instance , for the special case @xmath357 , the solution is simply @xmath358 ( see eqs [ petoilex ] and [ fonctionq ] ) and @xmath359 except for this case @xmath357 , we are not aware of any explicit solution for the eigenvalue @xmath356 ( even for the simpler case where the distribution of the on - site energies @xmath360 is a cauchy law and where the distributions @xmath361 and @xmath350 are also cauchy laws @xcite ) . with eqs [ eqpropreax ] and [ kernel ] ,
one concludes that at fixed @xmath0 , each mode @xmath354 is associated to the velocity @xmath362 ( eq . [ ynscaling ] ) @xmath363      the selection of the tail exponent @xmath354 of eq .
[ tail ] and of the corresponding velocity @xmath364 of eq [ ynscaling ] usually depend on the form of the initial condition @xcite . in our present case , the initial condition is completely localized ( see eq .
[ riccatifirst ] ) @xmath365 in this case , one expects that the solution that will be dynamically selected @xcite corresponds to the tail exponent @xmath183 and to the velocity @xmath366 determined by the following extremization @xmath367_{\beta=\beta_{selec}(w ) } = \left [ \frac{1}{\beta^2 } \ln   \lambda_0(\beta , w )   - \frac{1}{\beta } \frac { \partial_{\beta}\lambda_0(\beta , w ) } { \lambda_0(\beta , w ) } \right]_{\beta=\beta_{selec}(w ) } \label{selection}\end{aligned}\ ] ] the critical point is then determined by the two conditions @xmath368      since the general discussion is rather obscured by the absence of an explicit expression for the eigenvalue @xmath356 of the kernel of eq .
[ kernel ] , it is useful to consider the following strong disorder approximation ( called upper limit condition in @xcite ) , where the recursion for the @xmath331 in eq .
[ ricattilinearized ] is simply replaced by @xcite @xmath369 the argument is that in the limit of very large @xmath0 , the distribution of @xmath336 has also a width of order @xmath0 , so that the neglected terms in eq .
[ ricattilinearized ] are of order @xmath370 .
of course the approximation of eq .
[ strongapprox ] is not very well controlled , because it concerns random variables , and it suppresses important correlations between the variables @xmath322 of eq .
[ ricattilinearized ] . nevertheless , it is useful to consider it before returning to the true recursions of eq .
[ ricattilinearized ] , because with eq .
[ strongapprox ] , the eigenvalue @xmath356 is replaced by simple expression @xcite ( see also section 7.2 of @xcite where the approximation of independence between @xmath322 is considered for the cauchy case ) @xmath371 the velocity then reads @xmath372 for the flat distribution @xmath360 of random site energies ( eq . [ flat ] ) , one thus obtains @xmath373 the function @xmath374 is defined on the interval @xmath375 : it flows towards @xmath376 in the limits @xmath377 and @xmath378 with the following behaviors @xmath379 and @xmath380 it has a single maximum at @xmath381 given by the extremum condition of eq .
[ selection ] , so that @xmath381 is actually independent of @xmath0 . for @xmath382 ,
the selected exponent is of order @xmath383 for any @xmath0 .
we note that is value is close to the value @xmath177 that we measure for the large disorder @xmath184 ( see figure [ figtranswidth ] a ) . in conclusion
, this strong disorder approximation allows to see explicitly how things work on a simple example , and seems to give a reasonable value of @xmath384 for very large @xmath0 .
we now return to the analysis of the full problem of eq .
[ ricattilinearized ] . here
we should say that we do not agree with the discussion of @xcite concerning the selection of the tail exponent @xmath354 :    \(i ) in @xcite , the authors conclude that within the localized phase , @xmath183 decreases from unity and reaches @xmath386 exactly at criticality ( see the text between eqs .
6.8 and 6.9 in @xcite ) .
\(ii ) in our numerical results of section [ trans ] we have found instead that the selected exponent @xmath183 is always smaller than @xmath144 : it slightly grows as the disorder strength @xmath0 decreases and has a value close to @xmath387 near criticality ( see figure [ figtranswidth ] a ) .     as a function of the tail exponent @xmath354 for various disorder strength @xmath0 : ( i ) in the localized phase @xmath125 ,
the velocity is extremum at a @xmath0-dependent value @xmath388 and it is positive @xmath389 ( ii ) at the critical point @xmath390 , the velocity is extremum at @xmath386 where it vanishes @xmath391 ( iii ) in the delocalized phase , there exists another traveling - wave moving in the other direction of velocity @xmath392 : this velocity is extremum at @xmath393 for any @xmath136 ( see appendix b for more details ) , height=226 ]    we propose the following argument to justify our finding @xmath385 .
we think that the tail analysis is actually well defined only on the interval @xmath375 for the following reasons .
the eigenvalue eq . [ kernel ] relates the behavior of @xmath348 at @xmath394 to the behavior of @xmath348 near the origin @xmath395 . for @xmath396 ,
the function @xmath397 will be finite in the region where @xmath398 is finite , i.e. the integration is dominated by the region @xmath399 and we obtain the power - law decay @xmath400 however the function @xmath348 has to be integrable at @xmath401 to obtain from eq . [ tail ] the probability of the variable @xmath340 alone at large @xmath340 @xmath402 according to eq .
[ axinfinity ] , the function @xmath348 is integrable only for @xmath403 i.e. we obtain that the tail analysis has a meaning only for @xmath404 . in the limit @xmath377 , we expect from eq .
[ lambda0beta0 ] the same behavior as in eq .
[ vapproxflatzero ] @xmath405 exactly at @xmath393 , the authors of @xcite have argued that @xmath406_{\beta=1/2 } = 0   \label{partialdemi}\end{aligned}\ ] ] for any @xmath0 , as a consequence of the symmetry @xmath407 coming from the consideration of the adjoint kernel ( see more details around eqs .
6.7 and 6.8 in @xcite ) . here in contrast to @xcite
, we think that the region @xmath408 is not physical because of eq .
[ axinfinity ] , but the condition of eq .
[ partialdemi ] is useful to understand why the critical point determined by eqs [ criticonditions ] corresponds to the tail exponent @xcite @xmath409 to summarize the selection mechanism , the shape of the velocity @xmath362 as a function of the tail exponent @xmath354 in shown on fig .
[ figvelocity ] for various disorder strength @xmath0 :    \(i ) in the localized phase @xmath125 , the velocity derivative is negative at @xmath393 : @xmath410_{\beta=1/2 } = \left [ \frac{1}{\beta^2 } \ln   \lambda_0(\beta , w > w_c )   \right]_{\beta=1/2 } = - \left [ \frac{1}{\beta } v(\beta , w > w_c ) \right]_{\beta=1/2 } < 0 \label{derineg}\end{aligned}\ ] ] the velocity is extremum at a @xmath0-dependent value @xmath388 and the corresponding selected velocity is positive @xmath389 .
\(ii ) at the critical point @xmath390 , the velocity is extremum at @xmath386 where it vanishes @xmath391 .
\(iii ) in the delocalized phase , the velocity derivative is positive at @xmath393 : @xmath411_{\beta=1/2 } = \left [ \frac{1}{\beta^2 } \ln   \lambda_0(\beta , w < w_c )   \right]_{\beta=1/2 } = - \left [ \frac{1}{\beta } v(\beta , w < w_c ) \right]_{\beta=1/2 } > 0 \label{deripos}\end{aligned}\ ] ] we explain in appendix b that in this delocalized phase , there exists another traveling - wave moving in the other direction with the velocity @xmath392 : this velocity @xmath412 is then extremum at @xmath393 for any @xmath136 ( see appendix b for more details ) .    as a final remark
, we believe that the solution exactly at @xmath118 is not valid anymore because of eq .
[ axinfinity ] , and that the appropriate treatment exactly at criticality should replace the finite velocity motion in @xmath413 assumed in eq .
[ lnynscaling ] by the form @xmath414 where the anomalous exponent @xmath415 has been discussed in the text ( see eqs [ transfss ] and [ transrhonume ] ) .      in conclusion ,
the analysis of @xcite is thus very close to the traveling wave analysis of the directed polymer model @xcite : the difference is that in the directed polymer , the variables @xmath331 are random variables independent of the @xmath323 and the selected @xmath384 can be simply obtained from the non - integer moments of @xmath336 ( see the section 7.2 in @xcite ) , whereas in the present localization model , the variables @xmath331 and @xmath323 are correlated and one has thus to solve the eigenvalue problem for the integral kernel of eq . [ kernel ] . in the following appendix
, we explain how the same ideas can be used in the delocalized phase .
as explained in section [ eigen ] , in the delocalized phase , we have considered the recurrence for the real riccati variable @xmath331 @xmath416 together with the recurrences of eqs [ reccq ] for the auxiliary variables @xmath237 introduced in eq .
[ defcq ] .      to make more visible
the similarities with the previous appendix , it is convenient to perform the change of variables @xmath417 to obtain the following recursions for these variables @xmath418    the recursion of eq . [ ricattix ] for @xmath331 alone will as before converge to some distribution @xmath332 satisfying eq .
[ petoilex ] . in the delocalized phase
, one expects that the variables @xmath237 will grow exponentially ( see eqs [ lyapunovcq ] ) .
we thus set @xmath419 where @xmath267 governs the exponential growth and where @xmath420 remains a finite random variable upon the iteration @xmath421 the similarity with the discussions of appendix a is now obvious ( see eqs [ ynscaling ] and [ ricattilinearized ] ) . assuming the power - law decay @xmath422 with @xmath423 , one finds that the function @xmath424 has to satisfy the eigenvalue problem @xmath425 i.e. it is the same integral kernel as in eq .
[ eqpropreax ] with the correspondences @xmath426 and @xmath427 , but now we are interested in the phase @xmath267 .      we have argued above ( see section [ beta0demi ] ) that the problem for the function @xmath424 is well defined only for @xmath428 , i.e. the exponent @xmath429 that governs the power - law of eq .
[ dtail ] is restricted to the interval @xmath430     of the variable @xmath431 of eq .
[ scalingdn ] for various disorder strength @xmath432 : the exponent @xmath433 of eq .
[ dtail ] is of order @xmath434 in the whole delocalized phase .
( b ) for @xmath435 , comparison of the probability distributions @xmath436 and @xmath437 of the variable @xmath438 for @xmath273 and @xmath213 : the exponents of eq .
[ dtail ] are respectively @xmath434 and @xmath439 in agreement with eq .
[ muselec].,title="fig:",height=226 ]   of the variable @xmath431 of eq .
[ scalingdn ] for various disorder strength @xmath432 : the exponent @xmath433 of eq .
[ dtail ] is of order @xmath434 in the whole delocalized phase .
( b ) for @xmath435 , comparison of the probability distributions @xmath436 and @xmath437 of the variable @xmath438 for @xmath273 and @xmath213 : the exponents of eq .
[ dtail ] are respectively @xmath434 and @xmath439 in agreement with eq .
[ muselec].,title="fig:",height=226 ]    denoting as before @xmath356 the maximal eigenvalue of the kernel of eq .
[ kernel ] , each mode @xmath440 is associated to the lyapunov exponent ( see eq .
[ ynscaling ] ) @xmath441 from eq .
[ lambda0beta0 ] , we have the following behavior for @xmath442 @xmath443 near the other boundary @xmath444 , the property of eq .
[ partialdemi ] , yields that the partial derivative with respect to @xmath440 is still negative at @xmath444 @xmath445_{\mu \to 1/(2q ) } = \left [ -\frac{1}{\mu^2 } \ln   \lambda_0 ( 1/2,w )   \right]_{\mu \to 1/(2q ) } < 0 \label{derimu}\end{aligned}\ ] ] for the special case @xmath273 , the curve @xmath446 is directly related via @xmath447 to the curve @xmath364 discussed in appendix a and shown on fig .
[ figvelocity ] ( see the delocalized case @xmath136 ) .
since the selected tail exponent @xmath448 has to extremize @xmath449 , we conclude that the selected tail exponents are given by the simple values @xmath450 in the whole delocalized phase .
this is in agreement with our numerical results .
we show on fig .
[ fighistoc1c2 ] ( a ) the probability distribution of @xmath451 of the variable @xmath452 of eq .
[ scalingdn ] for various disorder strength @xmath432 : the tail exponent @xmath433 of eq .
[ dtail ] remains of order @xmath434 . for @xmath213
, we find similarly that the selected tail exponent remains the same within the whole delocalized phase and that it is of order @xmath439 . on fig .
[ fighistoc1c2 ] ( b ) , we compare for @xmath435 the probability distributions @xmath451 and @xmath453 of the variables @xmath431 : the exponents of eq .
[ dtail ] are respectively @xmath434 and @xmath439 .    from eq .
[ muselec ] , we obtains that the lyapunov exponents satisfy the identities @xmath454 as explained in the text around eq .
[ lyapunoviqidentity ] , this identity is important to understand the decay of i.p.r . with the number @xmath1 of generations in the delocalized phase .
as explained in section [ poolmethod ] , the pool method consists in representing a probability distribution @xmath455 by a large number @xmath89 of random variables @xmath456 .
in particular this introduces a cut - off in the tail of @xmath455 around @xmath457 with @xmath458 . in the field of traveling waves ,
the presence of such a cut - off has been much studied ( see @xcite and references therein ) : the leading correction to the selected velocity is usually logarithmic in @xmath89 @xmath459 for our present problem , where the critical point corresponds to a vanishing velocity @xmath460 , we thus expect that the true critical point @xmath461 corresponding to @xmath462 will be shifted towards a pool - dependent pseudo - critical point @xmath463 corresponding to @xmath464 .
as explained in the text , all the traveling waves encountered in the present paper vanish linearly in @xmath133 .
we thus expect that the difference @xmath465 $ ] will also decay only logarithmically as @xmath466 . from the analysis presented in appendix
a and b , it is obvious that the pool dependent critical point for the traveling wave of the landauer transmission ( see appendix a ) and for the auxiliary variable @xmath236 or @xmath467 ( see appendix b ) are the same because they are defined by the same condition in terms of the eigenvalue @xmath468 of the integral kernel .
however , the other variables @xmath237 or @xmath469 that have the same true critical point @xmath461 , will not have the same pool - dependent critical point @xmath463 , because the constant @xmath470 in eq .
[ vpool ] will depend on @xmath192 .
this is indeed what we observe with our numerical computations with the pool @xmath116 : the pseudo - critical point @xmath463 for the variable @xmath247 is below the pseudo - critical point @xmath463 for the variable @xmath236 .
we observe that the difference between the two becomes smaller for the bigger pool @xmath115 . in conclusion , numerical studies based on the pool method are valid for observables that are related to a single traveling wave , but one can not study the critical properties of observables that depend on two distinct traveling waves that have different pseudo - critical points . this is why in the text , we have not been able to present reliable numerical results for the statistics of i.p.r .
@xmath265 that depends on both @xmath236 and @xmath247 .
in this appendix , we briefly recall the finite - size scaling properties of the traveling / non - traveling phase transition studied in ref .
@xcite , because these properties are useful to interpret our numerical results described in the text .    for a one - dimensional branching random walk in the presence of an absorbing wall moving at a constant velocity @xmath328 ,
the survival probability @xmath471 presents a phase transition at @xmath472 with the following critical behaviors ( see more details in @xcite ) :    \(i ) for @xmath473 , it converges exponentially in time towards a finite limit @xmath474 @xmath475 where the limit @xmath476 presents an essential singularity @xmath477 and where the relaxation time @xmath478 diverges as @xmath479        \(iv ) the critical region is described by some finite - size scaling form governed by the relaxation time @xmath478 that appears in eq .
[ qbelow ] ( and not by the relaxation time @xmath482 appearing in eq .
[ qabove ] ) @xmath485 where the scaling function @xmath486 has the following asymptotic behaviors . for @xmath487 , has the one has the power - law @xmath488 to recover the finite limit of eq .
[ qessential ] and one has the scaling relation @xmath489 for @xmath490 , one has the power - law @xmath491 to recover the exponential decay of eq .
[ qabove ] , and one has the scaling relation @xmath492 for the problem considered in @xcite all the exponents are exactly known @xmath493 , @xmath494 , @xmath495 and @xmath496 .    in our numerical results presented in the text concerning the traveling / non - traveling phase transitions that occur in anderson localization on the bethe lattice , we find very similar critical behaviors : the velocity of the traveling wave vanishes linearly with @xmath493 , and the finite - size scaling in the critical region is governed by the other exponent @xmath19 that appears in the relaxation towards the finite value of the non - traveling phase .
majumdar and p.l .
krapivsky , phys .
e 62 , 7735 ( 2000 ) ; s.n .
majumdar and p.l .
krapivsky , phys .
e 65 , 036127 ( 2002 ) s.n .
majumdar and p.l .
krapivsky , physica a 318 , 161 ( 2003 ) ; s.n .
majumdar , phys .
e 68 , 026103 ( 2003 ) .
j. prior , a.m. somoza and m. ortuno , phys . rev .
b 72 , 024206 ( 2005 ) ; a.m. somoza , j. prior and m. ortuno , phys . rev .
b 73 , 184201 ( 2006 ) ; a.m. somoza , m. ortuno and j. prior , phys .
99 , 116602 ( 2007 ) .
e. brunet and b. derrida , phys .
e 56 , 2597 ( 1997 ) ; e. brunet , b. derrida , j. stat .
103 , 269 ( 2001 ) ; e. brunet , b. derrida , phys .
e 70 , 016106 ( 2004 ) ; e. brunet , b. derrida , a. h. mueller , s. munier , europhys .
lett . , 76 , 1 ( 2006 ) ; e. brunet , b. derrida , a. h. mueller , and s. munier ; phys .
e 73 , 056126 ( 2006 ) ; e. brunet , b. derrida , a. h. mueller , and s. munier phys .
e 76 041104 ( 2007 ) .