the metal - insulator transition induced by disorder in the problem of anderson localization remains poorly understood from a theoretical point of view . on the one hand , it is true that the prediction for the existence of a metal - insulator transition induced by disorder has been verified by numerical simulations . on the other hand ,
the critical exponents at the metal - insulating transition measured numerically still can not be extracted from the field theories which are believed to embody the metal - insulator transition.@xcite the paradigm of this unpleasant situation is the plateau transition in the lowest landau level of the integer quantum hall effect.@xcite the development of reliable methods to compute analytically critical properties at a disorder - induced metal - insulator transition remains to this date an open problem .
the simplest known example of disorder - induced criticality in the context of anderson localization was solved by dyson in 1953.@xcite according to dyson , the global density of states ( dos ) in the thermodynamic limit and close to the band center is anomalous when a single quantum particle hops with a random amplitude between the nearest - neighbor sites of a chain excluding any other form of disorder , in short , the one - dimensional random hopping problem .
subsequent works showed that dyson s spectral anomaly is related to a diverging localization length,@xcite an anomalous decay of the envelope of wave functions,@xcite and an anomalous probability distribution of the resistance upon approaching the band center.@xcite    the coalescence of the concept of universality imported from the theory of critical phenomena with the symmetry classification imported from random matrix theory has led to a classification of diffusive regimes and disorder - induced metal - insulator transitions in terms of 10 universality classes that are uniquely characterized by the intrinsic symmetries preserved by the disorder ( assumed weak ) for any given dimensionality of space.@xcite according to this classification , dyson s singular dos is a trademark of the chiral universality classes .
the chiral universality classes are realized in problems of anderson localization for which the hamiltonian anticommutes with some unitary operator for all realizations of the disorder and in which case the dos is known to be singular at the band center in zero,@xcite quasi - one,@xcite and two dimensions.@xcite even though the chiral universality classes are critical at the band center , they have been far less studied than the orthogonal , unitary , and symplectic universality classes , which we will refer to as the standard universality classes .
for example , the probability distribution of the so - called mesoscopic local density of states ( ldos ) has been known since 1989 from the work of altshuler and prigodin for a single chain with weak on - site disorder,@xcite whereas it has not yet been computed for the one - dimensional random - hopping problem .
the first aim of this paper is to fill this gap .
if a diffusive regime exists , perturbative techniques can be used to describe the crossovers between universality classes.@xcite however , little is known quantitatively on the crossover between two different universality classes , here defined in terms of nonlinear - sigma models , say , from the diffusive to the localized regimes .
this is even true in one dimension , which is by far the most studied laboratory for anderson localization as was recently illustrated by a flurry of works on the validity of one - parameter scaling in one dimension.@xcite the absence of a diffusive regime renders the concept of universality classes ambiguous in one dimension in that symmetry alone does not specify a universality class in one dimension .
for example , the ratio of the localization length to the mean free path in a wire of finite width interpolates smoothly between its two limiting values for fully preserved or completely broken time - reversal symmetry , respectivley , as a function of a weak magnetic field.@xcite similarly , fine tuning of microscopic parameters is required to achieve delocalization ( i.e. , diverging localization length ) of quasiparticles at the fermi energy in a dirty superconducting wire of finite thickness with both broken spin - rotation and time - reversal symmetries whereas delocalization becomes generic in the thick wire scaling limit.@xcite one can nevertheless define a _ symmetry class _ in one dimension by demanding that single - parameter scaling holds .
lack of universality in one dimension can then be understood as the fact that there is no scaling limit for which single - parameter scaling becomes a generic property of an ensemble of random microscopic hamiltonians of a given symmetry .
symmetry classes in one dimension can not be construed as an enumeration of stable or unstable fixed points of some putative effective field theory .
computing the crossover between symmetry classes in one dimension nevertheless remains a well - defined problem .
time - reversal symmetry can not be broken for a spinless particle constrained to move in a one - dimensional and simply connected world , in which case the only possible crossover between symmetry classes takes place between the chiral and standard classes .
the second aim of this paper is to compute the crossover between these two symmetry classes in one dimension for the probability distributions of the conductance and of the ldos from the ballistic to the localized regimes assuming a weak disorder at the microscopic level . in that regard ,
our results are complementary to the analysis made by schomerus and titov of the validity of one - parameter scaling in one dimension in which they computed exactly the first four cumulants of the logarithm of the conductance to leading order in the ratio of the length of the chain to the localization length assuming a weak disorder at the microscopic level.@xcite we thus propose an approximate ansatz for the crossover whereby the approximation consists of assuming that the phase and amplitude of the reflection coefficient separate . by comparing the cumulants of the lyapunov exponent of the transfer matrix computed with our approximate ansatz to the ones computed with the large - deviation ansatz of schomerus and titov deep in the localization regime , we deduce that our ansatz captures the first cumulant , but fails with the higher cumulants .
the paper is organized as follows .
we formulate the one - dimensional problem of anderson localization as a continuous schrdinger equation with a relativistic kinetic energy in sec .
[ model ] .
a functional renormalization group equation is derived for the probability distribution of the reflection coefficient in sec .
[ functional renormalization group equations ] .
exact solutions to these coupled functional renormalization group equations are given for the chiral and the standard classes , as well as an approximate solution that describes the crossovers between these two classes in sec .
[ sec : solutions to the functional renormalization group equations ] .
the accuracy of the approximate solution for the crossover regime is estimated in sec .
[ sec : validity of the separation ansatz for the crossover regime ] .
the probability distribution of the ldos is computed exactly in the chiral and standard classes whereas it is computed approximately in the crossover regime between these two classes in sec .
[ sec : probability distribution of the ldos ] .
we conclude in sec .
[ sec : conclusion ] .
( 220,150)(-120,-130 ) ( -80 , 20)(1,0)140 ( -80,-10)(1,0)140 ( 60,-10)(0,1)30    ( -40 , 5 ) ( 0,10)-180 ( -40 , 15)(-1 , 0 ) 20 ( -40,-5)(-1 , 0 ) 20 ( -120 , 17)@xmath0 ( -95 , 5)@xmath1 ( -110 , -5)@xmath2    ( -70,-10)(0,1)30    ( -70,20)(10,0)10(1,-1)30    ( -70,10)(1,-1)20 ( -70,0)(1,-1)10    ( 30,20)(1,-1)30 ( 40,20)(1,-1)20 ( 50,20)(1,-1)10    ( -80,-30)(1,0)140 ( -80,-60)(1,0)140 ( 60,-60)(0,1)30    ( -40,-45 ) ( 0,10)-180 ( -40,-35)(-1 , 0 ) 20 ( -40,-55)(-1 , 0 ) 20 ( -120 , -33)@xmath3 ( -95,-45)@xmath1    ( -70,-60)(0,1)30    ( -70,-30)(10,0)10(1,-1)30    ( -70,-40)(1,-1)20 ( -70,-50)(1,-1)10    ( 30,-30)(1,-1)30 ( 40,-30)(1,-1)20 ( 50,-30)(1,-1)10    ( -80 , -80)(1,0)150 ( -80,-110)(1,0)150    ( -40 , -95 ) ( 0,10)-180 ( -40 , -85)(-1 , 0 ) 20 ( -40,-105)(-1 , 0 ) 20 ( -120 , -83)@xmath4 ( -95 , -95)@xmath1    ( -70,-110)(0,1)30 ( 60,-110)(0,1)30    ( -70,-80)(10,0)10(1,-1)30    ( -70,-90)(1,-1)20 ( -70,-100)(1,-1)10    ( 30,-80)(1,-1)30 ( 40,-80)(1,-1)20 ( 50,-80)(1,-1)10    ( -70,-120)(1,0)130 ( 60,-120)(-1,0)130 ( -5 , -130)@xmath5
in this paper we want to investigate the statistical properties of the ldos @xmath6 and dimensionless conductance @xmath7 for a single quantum particle , which is ( i ) restricted to a simply connected and strictly one - dimensional ( 1d ) geometry and ( ii ) subjected to a static and weak random environment ( disorder ) that interpolates smoothly between the standard and chiral symmetry classes . condition ( i ) implies that time - reversal symmetry is preserved for every realization of the disorder . in 1d and assuming a metallic ground state in the absence of the disorder , the ground state consists of two distinct fermi points .
plane waves at the fermi points are called left- and right - movers , respectively . in the standard symmetry class , the channels for forward and backward scatterings induced by an on - site disorder potential in the basis of left- and right - movers are equally likely up to nonrandom oscillatory factors . in the chiral symmetry class ,
the disorder potential is off - diagonal in the basis of left- and right - movers .
presuming weak disorder in condition ( ii ) allows us to choose a kinetic energy that is a first - order differential operator , for it only makes sense to linearize the spectrum about the two fermi points when the disorder potential is weak .
we thus choose to model the dynamics of the single quantum particle and the static random environment by the hamiltonian    [ eq : def our anderson localization problem ] @xmath8 where we have set the fermi velocity @xmath9 and @xmath10 to one .
the boundary conditions that will be imposed in the sequel are depicted in fig .
[ fig : geometries ] .
they all obey condition ( i ) .
we have denoted by @xmath11 the unit @xmath12 matrix and by @xmath13 the @xmath14 pauli matrices .
the static random environment is represented by three independent potentials @xmath15 , @xmath16 that we choose , for mathematical convenience , to be supported on the interval @xmath17 , and to be white - noise and gaussian distributed with vanishing mean , ( disorder averaging is denoted by @xmath18 ) @xmath19    the three variances @xmath20 carry the dimensions of inverse length and need not be equal .
a finite value @xmath21 can also be easily accommodated in our formalism.@xcite in this case , the clean system is insulating as @xmath22 opens up a dimerization gap at the fermi energy and the random environment mimics to a first approximation the static fluctuations of the phonons responsible for the dimerization sufficiently close to the peierls transition.@xcite the problem of anderson localization in the continuum defined by eq .
( [ eq : def our anderson localization problem ] ) is a coarse - grained version of a single particle hopping on a simply - connected chain with a uniform nearest - neighbor hopping amplitude subjected to weak and independent random fluctuations of an on - site potential and of the nearest - neighbor hopping amplitude .    the generic symmetry obeyed by eq .  ( [ eq : def our anderson localization problem a ] ) is time - reversal invariance , @xmath23 if the condition @xmath24 is satisfied , or , for that matter , any equivalent condition obtained by a redefinition of the pauli matrices through a position independent su(2 ) rotation around @xmath25 , eq .
( [ eq : def our anderson localization problem a ] ) has an additional chiral symmetry in that it anticommutes with @xmath26 , @xmath27 consequently , all nonvanishing eigenvalues of @xmath28 come in pairs with opposite signs .
the origin of the chiral symmetry of the continuum model for any given realization of the disorder is a sublattice symmetry of an appropriate lattice regularization . a microscopic random hamiltonian realizing the chiral symmetry is , for example , a tight - binding hamiltonian on a hypercubic lattice with strictly vanishing on - site energies and random hopping matrix elements restricted to nearest - neighbor sites , in short , the random hopping problem .
the ldos , dos , conductance , shot - noise power , etc .
, can all be extracted from the elements of the @xmath14 scattering matrix @xmath29 at the energy @xmath30 that connects incoming ( `` i '' ) and outgoing ( `` o '' ) states @xmath31 the amplitude for incoming and outgoing waves to the left- and right - hand side of the disordered region is denoted by @xmath32 here .
the matrix elements @xmath33 are the reflection coefficients whereas the matrix elements @xmath34 are the transmission coefficients .
alternatively , one can work with the @xmath14 transfer matrix @xmath35 at the energy @xmath30 which follows from the scattering matrix ( [ eq : def s matrix ] ) through the basis transformation implied by @xmath36 conservation of probability dictates that the scattering matrix is unitary or , equivalently , that the transfer matrix is pseudo - unitary    @xmath37    conservation of probability also implies the existence of the polar decompositions    [ eq : polar decomposition s+m ] @xmath38    where @xmath39 , @xmath40 , @xmath41 , @xmath42 are independent complex numbers with @xmath43 and @xmath44 .
the polar decomposition is not unique .
for example , the sign of @xmath45 can be absorbed into a redefinition of @xmath39 , @xmath40 , @xmath41 , @xmath42 .
the variable @xmath45 , when restricted to the halfline @xmath46 , has the geometrical interpretation of a radial coordinate on a riemannian manifold.@xcite the time - reversal symmetry of the hamiltonian implies the transformation laws    @xmath47    which enforce the constraints @xmath48 and @xmath49 or , equivalently , @xmath50 . the chiral symmetry of the hamiltonian implies the transformation laws    @xmath51    which enforce the constraints @xmath52 and @xmath53 or , equivalently , @xmath54 , @xmath55 , @xmath56 .
then , at the band center @xmath57 with the condition ( [ eq : def chiral limit for g s ] ) , the scattering matrix becomes hermitian .    under the assumptions that ( i ) the disorder is weak and ( ii ) @xmath39 , @xmath40 , @xmath41 , @xmath42 are all
independently and uniformly distributed on the unit circle in the complex plane , it was essentially shown in ref .   that the probability distribution @xmath58 of @xmath59 obeys the fokker - planck equation @xmath60 in the geometry of fig .
[ fig : geometries]c .
the mean free path @xmath61 is some function of @xmath62 , @xmath63 , @xmath64 , and @xmath65 . in a more general context of multi channel quantum wires assumption ( ii )
is known as the isotropy assumption,@xcite whereas the fokker - planck equation ( [ eq : fp std limit ] ) is a special case of the dorokhov - mello - pichard - kumar ( dmpk ) equation.@xcite a consequence of assumption ( ii ) is one - parameter scaling as encoded by the fokker - planck equation ( [ eq : fp std limit ] ) , i.e. , the probability distribution of @xmath66 and , consequently , of @xmath67 depends on the single dimensionless parameter @xmath68 .
it is evident that assumption ( ii ) breaks down at the band center @xmath57 and with the chiral condition ( [ eq : def chiral limit for g s ] ) because @xmath69 and @xmath70 are then both real valued . correspondingly , the fokker - planck equation obeyed by the probability distribution @xmath58 of @xmath71 is different from eq .
( [ eq : fp std limit ] ) and given by the diffusion equation @xmath72 in the geometry of fig .
[ fig : geometries]c.@xcite the mean free path @xmath61 is some function of @xmath64 . again , the fokker - planck equation ( [ eq : fp chiral limit ] ) encodes one - parameter scaling as the probability distribution of @xmath66 and , consequently , of @xmath67 depends on the single dimensionless parameter @xmath68 .
there are no known multiple - parameter scaling equations for the ldos , dos , conductance , etc . , which describe the crossover between the two limiting cases described by eqs .
( [ eq : fp std limit ] ) and ( [ eq : fp chiral limit ] ) for values of @xmath5 ranging from the ballistic ( @xmath5 smaller than the mean free path ) to the localized regime ( @xmath5 larger than the mean free path ) to the best of our knowledge ( see refs .   for discussions of two - parameter scaling ) .
we fill this gap in the remainder of the paper .
( 240,160)(0,-140 ) ( 0,10)(a ) ( + 10,-10)(1,0)220 ( + 10,-50)(1,0)220 ( 50,-50)(0,1)40 ( 190,-50)(0,1)40 ( 120,-60)(-1,0)70 ( 120,-60 ) ( 1,0)70 ( 117,-70)@xmath5 ( 40,-50)(0,1)40 ( 20 , 0 ) ( 1,0)20 ( 70 , 0)(-1,0)20 ( 40 , 5)@xmath73 ( 50,-20)(1,-1)30 ( 50,-30)(1,-1)20 ( 50,-40)(1,-1)10 ( 50,-10)(10,0)11(1,-1)40 ( 190,-20)(-1,1)10 ( 190,-30)(-1,1)20 ( 190,-40)(-1,1)30 ( 40,-50)(0,10)4(1,1)10 ( 50,-50)(0,10)4(-1,1)10 ( 0,-80)(b ) ( 10,-100)(1,0)220 ( 10,-140)(1,0)220 ( 40,-140)(40,0)2(0,1)40 ( 160,-140)(40,0)2(0,1)40 ( 80,-140)(-10,0)4(-1,1)10 ( 40,-100)(10 , 0)4(1,-1)10 ( 40,-100)(0,-10)4(1,-1)10 ( 80,-140)(0 , 10)4(-1,1)10 ( 200,-140)(-10,0)4(-1,1)10 ( 160,-100)(10 , 0)4(1,-1)10 ( 160,-100)(0,-10)4(1,-1)10 ( 200,-140)(0 , 10)4(-1,1)10 ( 10,-113)@xmath74 ( 10,-130)@xmath75 ( 208,-113)@xmath76 ( 208,-130)@xmath77 ( 85,-113)@xmath78 ( 85,-130)@xmath79 ( 130,-113)@xmath80 ( 130,-130)@xmath81 ( 55,-122)@xmath82 ( 175,-122)@xmath83
we shall show in secs .   [
sec : solutions to the functional renormalization group equations ] and [ sec : probability distribution of the ldos ] that the statistical properties of the dimensionless conductance @xmath7 and of the ldos @xmath6 follow from the knowledge of the statistical properties of the reflection coefficient @xmath84 which has been decomposed into its square modulus @xmath85 and phase @xmath86 . for a disordered wire of length @xmath87 as depicted in fig .
[ fig : composition law and rg step ] , the reflection coefficient @xmath88 is related to the entries of the scattering matrix @xmath89 for the slice of length @xmath73 and to the entry @xmath90 of the scattering matrix @xmath91 for the much longer segment of length @xmath5 by the composition law @xmath92 when the width @xmath73 of the slice is much larger than the lattice spacing @xmath93 but much smaller than the mean free path @xmath61 , here defined by @xcite @xmath94 we infer the continuous langevin process    [ eq : langevin eq for r and phi ] @xmath95    by using the relations @xmath96 which are valid up to first order in the disorder potentials @xmath97 , and by expanding the right - hand side of eq .
( [ eq : composition law reflection amplitude ] ) to the same order .
we note in passing that the mean free path ( [ eq : def ell ] ) is simply given by @xmath98 in the born approximation . for later convenience
, we have continued the energy @xmath30 to the upper part of the complex plane , @xmath99 @xmath100 .
the continuous langevin process ( [ eq : langevin eq for r and phi ] ) can also be formulated as the fokker - planck equation obeyed by the joint probability distribution function @xmath101 . using the standard methods of chapter 3 in ref .
say , one finds@xcite    [ eq : master equation in terms of r and phi ] @xmath102 where @xmath103 and @xmath104    we have introduced the short - hand notations @xmath105 observe that the dependence on the polar angle @xmath106 arises solely from @xmath107 and @xmath108 in the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) .
thus , provided initial and boundary conditions are also periodic on the interval @xmath109 , the joint probability distribution @xmath101 is periodic on the interval @xmath109 even though the phase @xmath106 is originally defined modulo @xmath110 .
the periodicity @xmath111 is a consequence of the langevin process ( [ eq : langevin eq for r and phi ] ) being invariant under @xmath112 and @xmath113 . a dimerization gap @xmath114 produces the changes    @xmath115    the polar decomposition ( [ eq : polar decomposition s+m ] ) suggests the change of variable @xmath116 in terms of which the langevin process ( [ eq : langevin eq for r and phi ] ) becomes    [ eq : langevin eq for x and phi ] @xmath117    in turn , the joint probability distribution function @xmath118 obeys the fokker - planck equation    [ eq : master equation in terms of x and phi ] @xmath119 where @xmath120 and @xmath121    the langevin process ( [ eq : langevin eq for x and phi ] ) or the fokker - planck equation ( [ eq : master equation in terms of x and phi ] ) are invariant under @xmath122 and @xmath113 .    equations ( [ eq : master equation in terms of r and phi ] ) or ( [ eq : master equation in terms of x and phi ] ) are , from a conceptual point of view , the main result of this paper as they determine under what conditions one - parameter scaling holds in a 1d weakly disordered quantum wire .
for example , eq .  ( [ eq : master equation in terms of x and phi ] ) encodes the functional renormalization group flow of the joint probability distribution of the radial
coordinate @xmath123 and phase @xmath106 of the reflection coefficient @xmath1 .
if it is the property that the functional renormalization group flow of the probability distribution of @xmath66 depends solely on the dimensionless ratio @xmath68 aside from @xmath66 , which is understood as one - parameter scaling for @xmath66,@xcite then two examples of one - parameter scaling can be constructed from eq .
( [ eq : master equation in terms of x and phi ] ) as follows .
first , the fokker - planck ( dmpk ) equation ( [ eq : fp std limit ] ) follows from insertion of the ansatz @xmath124 into eq .
( [ eq : master equation in terms of x and phi ] ) with @xmath125 and from integration over @xmath126 of the resulting equation .
second , the fokker - planck ( dmpk ) equation ( [ eq : fp chiral limit ] ) follows from insertion of the ansatz @xmath127 \mathcal{x}(x;l ) \hphantom{aa}\end{aligned}\ ] ] with @xmath128 into eq .
( [ eq : master equation in terms of x and phi ] ) and from integration over @xmath106 of the resulting equation .
one - parameter scaling is not the rule for generic initial and boundary conditions of eq .
( [ eq : master equation in terms of x and phi ] ) and for generic values of @xmath65 and @xmath20 due to the absence of a diffusive regime in 1d.@xcite this is most clearly seen by the fact that the localization length @xmath129 , as defined by the typical dependence of the conductance @xmath67 on @xmath130 , is finite , @xmath131 , for the standard symmetry class , but diverges in the chiral limit @xmath128 .
hence , the localization length must depend on a second microscopic parameter aside from @xmath61 for generic values of @xmath65 and @xmath20 .
equations ( [ eq : master equation in terms of r and phi ] ) and ( [ eq : master equation in terms of x and phi ] ) describe the full crossover from the chiral to the standard class with no restrictions on the values of @xmath66 and @xmath5 .    for any given absorption @xmath132 , the parameter space @xmath133 with @xmath134 and @xmath135 of eqs .
( [ eq : master equation in terms of r and phi ] ) and ( [ eq : master equation in terms of x and phi ] ) is three - dimensional and is depicted in fig .
[ fig : parms ] .
remarkably , the parameter @xmath136 always enters eqs .
( [ eq : master equation in terms of r and phi ] ) and ( [ eq : master equation in terms of x and phi ] ) as a prefactor to @xmath137 or @xmath138 and conversely @xmath137 or @xmath138 are always multiplied by @xmath136 . as we shall see in sec .
[ sec : solutions to the functional renormalization group equations ] this property turns out to be crucial in our study of the crossover between the chiral and standard symmetry classes .
equation ( [ eq : master equation in terms of r and phi ] ) was derived in the standard class by abrikosov , melnikov , and kumar , and by rammal and doucot among others starting from the continuous nonrelativistic schrdinger equation.@xcite the large @xmath66 limit @xmath139 of eq .
( [ eq : master equation in terms of x and phi ] ) was derived by schomerus and titov when @xmath125 and @xmath140 starting from a langevin process for the eigenfunctions of the discrete schrdinger equation on a chain .
@xcite    ( 180,170)(0,10 ) ( + 60,+70 ) ( 0,1)100 ( 0,+55 ) ( 4,1)180 ( 30,+90 ) ( 3,-2)120    ( 20,60 ) ( 0,1)40 ( 100,80 ) ( 0,1)40 ( 80,20 ) ( 0,1)70 ( 120,30 ) ( 0,1)70 ( 160,40 ) ( 0,1)70    ( 60,15 ) ( 4,1)120    ( 20,100 ) ( 4,1)80    ( 80,90 ) ( 4,1)80    ( 95,10 ) ( -3,2)90 ( 85,90 ) ( 3,-2)90 ( 30,160 ) ( 3,-2)90    ( 80,80)(1,1)13 ( 80,75)(1,1)20 ( 80,70)(1,1)26 ( 80,65)(1,1)32 ( 80,60)(1,1)40 ( 80,55)(1,1)46 ( 80,50)(1,1)53 ( 80,45)(1,1)59.5 ( 80,40)(1,1)66.5 ( 80,35)(1,1)73 ( 80,30)(1,1)80 ( 80,25)(1,1)80 ( 80,20)(1,1)80    ( 160,95)(-1,-1)73.5 ( 160,90)(-1,-1)66.5 ( 160,85)(-1,-1)59.5 ( 160,80)(-1,-1)53 ( 160,75)(-1,-1)46 ( 160,70)(-1,-1)40 ( 160,65)(-1,-1)33 ( 160,60)(-1,-1)26.5 ( 160,55)(-1,-1)19.5 ( 160,50)(-1,-1)13    ( 30,95)(1,-1)15 ( 30,100)(1,-1)30 ( 30,105)(1,-1)45 ( 30,110)(1,-1)60 ( 30,115)(1,-1)75 ( 30,120)(1,-1)90 ( 30,125)(1,-1)90 ( 30,130)(1,-1)90 ( 30,135)(1,-1)90 ( 30,140)(1,-1)90 ( 30,145)(1,-1)90 ( 30,150)(1,-1)90 ( 30,155)(1,-1)90 ( 30,160)(1,-1)90    ( 120,95)(-1,1)15 ( 120,90)(-1,1)30 ( 120,85)(-1,1)45 ( 120,80)(-1,1)60 ( 120,75)(-1,1)75    ( 170,105)@xmath136 ( 48,165)@xmath141 ( 150,15)@xmath65 ( 65,11)@xmath142 ( 13,50)@xmath143 ( 58,60)@xmath144 ( 93,70)@xmath145 ( 62,113)@xmath146 ( 21,65)ch1 ( 102,85)ch2 ( 21,105)kw1 ( 102,125)kw2 ( 20,60 ) ( 100,80 ) ( 20,100 ) ( 100,120 )
exact and approximate solutions to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) are derived in this section .
we will consider three regimes of parameter in fig .
[ fig : parms ] .
the first regime corresponds to the chiral symmetry class .
the second regime corresponds to the standard symmetry class .
the third regime describes a specific path connecting the two symmetry classes .    the chiral symmetry class is given by the condition    [ eq : def chiral regime parameters ] @xmath147 which defines the band center of the energy eigenvalue spectrum of hamiltonian ( [ eq : def our anderson localization problem a ] ) when @xmath148    i.e. , when hamiltonian ( [ eq : def our anderson localization problem a ] ) anticommutes with @xmath26 as in eq .  ( [ eq : def chiral symmetry for h ] ) .
the standard symmetry class occurs whenever    [ eq : strong+weak def of std limit ] @xmath149 or @xmath150 or @xmath151    since any one of these three conditions decouples @xmath152 from @xmath106 in eq .
( [ eq : master equation in terms of r and phi ] ) and is consistent with a stationary probability distribution of @xmath106 which is uniformly distributed . in this paper
we shall only consider the crossover between the chiral and standard symmetry classes along the line @xmath153    we must naturally choose the initial and boundary conditions to be imposed on the solutions to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) .
our choice is motivated by the computation of the probability distribution of the ldos that we will carry out in sec.[sec : probability distribution of the ldos ] .
two geometries depicted in fig .
[ fig : bc for ldos ] will be considered in sec.[sec : probability distribution of the ldos ] .
the relevant geometries for the probability distribution of the reflection coefficient @xmath1 that enter the computation of the ldos are depicted in fig .
[ fig : geometries ] . in figs .
[ fig : geometries]a and [ fig : geometries]b the disordered quantum wire is closed on the right - hand side and connected to a reservoir on the left - hand side . in fig .
[ fig : geometries]a conservation of probability holds . in fig .
[ fig : geometries]b conservation of probability does not hold due to the presence of an absorption encoded by a finite imaginary part @xmath154 of @xmath30 . in fig .
[ fig : geometries]c the disordered quantum wire is open at both ends and flux conservation ( no absorption ) is assumed .
we choose the ideal initial condition @xmath155 and @xmath156 when @xmath157 for the semiopen and fully open geometries of fig .
[ fig : geometries ] , respectively . in the chiral symmetry class ,
we demand that the phase of @xmath1 is either 0 or @xmath111 with probability 1/2 initially . in the standard symmetry class ,
we demand that the phase of @xmath1 is uniformly distributed in the interval @xmath158 initially . in view of the periodicity of the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) we always impose initial conditions with a periodicity of @xmath111 .
we shall now construct exact solutions to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) in the chiral and standard symmetry classes , respectively , for the three geometries of fig .
[ fig : geometries ] .
we shall also solve the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) in the crossover regime ( [ eq : crossover regime weak if a ] ) for the geometry of fig .
[ fig : geometries]a .
we could not construct exact solutions to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) in the crossover between the chiral and standard symmetry classes for the geometries of figs .
[ fig : geometries]b and [ fig : geometries]c .
we shall construct , however , an approximate solution to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) in a crossover regime that interpolates between the exact solutions for the geometries of figs .
[ fig : geometries]b and [ fig : geometries]c .
the accuracy of this approximate solution will be discussed in sec .
[ sec : validity of the separation ansatz for the crossover regime ] .      in geometry [ fig : geometries]a , an incoming plane wave from an ideal lead is perfectly reflected by the disordered quantum wire .
hence , the reflection coefficient @xmath1 of the disordered quantum wire must be a pure phase , @xmath159 , for all @xmath160 .
we thus insert the separation - of - variable ansatz    @xmath161    into the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) that , after integration over @xmath152 , reduces to @xmath162 \frac{\partial\phi}{\partial\phi } \nonumber \\ & & + \left [   2\zeta^{\ } _ { 0 } + 1 + \zeta \cos(2\phi ) \right ] \frac{\partial^{2}\phi}{\partial\phi^{2}}. \label{eq : fp if closed at one end and no absorption}\end{aligned}\ ] ]    by assumption , we have switched off the absorption @xmath125 .
our choice of initial condition shall depend on the proximity to the parameter regimes ( [ eq : def chiral regime parameters ] ) and ( [ eq : strong+weak def of std limit ] ) but shall always be periodic on the interval @xmath109 .      in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) the fokker - planck equation ( [ eq : fp if closed at one end and no absorption ] ) reduces to    @xmath163\frac{\partial^{2}\phi}{\partial\phi^{2 } } \nonumber \\ \label{eq : fp if closed at one end and no absorption ch}\end{aligned}\ ] ]    and has the normalized , stationary , and periodic solution @xmath164    whenever the phase @xmath106 is @xmath144 or @xmath111 with probability @xmath165 initially .
the same conclusion could have been anticipated by inspection of the langevin process ( [ eq : langevin eq for r and phi ] ) or from the fact that the chiral symmetry class ( [ eq : def chiral regime parameters ] ) imposes the condition that the scattering matrix is hermitian , see eq .
( [ eq : chiral symmetry s ] ) .      in the standard symmetry class ( [ eq : strong def of std limit ] ) , @xmath107 and @xmath108 drop out from the fokker - planck equation ( [ eq : fp if closed at one end and no absorption ] ) :    @xmath166    its normalized , stationary , and periodic solution is @xmath167    whenever the phase @xmath106 of @xmath1 is uniformly distributed initially .
more generally , the solution to eq .
( [ eq : fp if closed at one end and no absorption std ] ) converges to the uniform probability distribution ( [ eq : solution std limit semi - open geometry ] ) as soon as @xmath168 for any given initial probability distribution .
this is also true for the standard symmetry class ( [ eq : yet another def of std limit ] ) and ( [ eq : weak def of std limit ] ) as @xmath169 and @xmath170 are ineffective in these limits .    when @xmath171 and @xmath172 the stationary probability distribution of @xmath1 need not be uniform in @xmath158 . for example , consider the case @xmath173 , @xmath174 , @xmath175 , for which the fokker - planck equation ( [ eq : fp if closed at one end and no absorption ] ) reduces to    [ eq : band center anomaly ] @xmath176 \frac{\partial^{2}\phi}{\partial\phi^{2 } } \label{eq : band center anomaly a}\end{aligned}\ ] ] and has the stationary and normalized solution@xcite @xmath177 with a periodicity of @xmath111 . here
@xmath178 denotes the value of the gamma function for @xmath179 .
equation ( [ eq : band center anomaly b ] ) is a special case of the stationary solution @xmath180    to eq .
( [ eq : fp if closed at one end and no absorption ] ) with @xmath181 as shown by titov.@xcite in the chiral limit @xmath182 and @xmath183 , this solution becomes unnormalizable and we have to use eq .
( [ eq : sum of delta functions ] ) as a stationary solution instead .      for simplicity
, we consider the crossover regime ( [ eq : crossover regime weak if a]).@xcite the chiral symmetry class ( [ eq : def chiral regime parameters ] ) is reached when @xmath181 . the standard symmetry class ( [ eq : weak def of std limit ] ) is reached when @xmath184 . in the crossover regime , the fokker - planck equation ( [ eq : fp if closed at one end and no absorption ] ) reduces to @xmath185 a stationary solution to eq .
( [ eq : chiral 1 fp with no absorption and r=1 ] ) is constructed from the solution @xmath186 valid for @xmath187 and on the interval @xmath109 , by periodic extension to the interval @xmath158 .
the constant @xmath188 is chosen so that the solution ( [ eq : sol cross if geo c ] ) is properly normalized .
the solution interpolates between eqs .
( [ eq : sum of delta functions ] ) and ( [ eq : solution std limit semi - open geometry ] ) as @xmath189 and @xmath190 , respectively , as it should be . equation ( [ eq : sol cross if geo c ] ) is an important intermediary result as we will use it to construct an approximation to the probability distributions for the conductance and ldos in the crossover regime that interpolate between the exact chiral and standard limiting probability distributions .      in geometry [ fig : geometries]b , an incoming plane wave from an ideal lead
is not perfectly reflected by the disordered quantum wire due to a finite absorption @xmath191 .
hence , the reflection coefficient @xmath1 of the disordered quantum wire is not a pure phase anymore , @xmath192 .      in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) the phase @xmath106 and squared magnitude @xmath152 of the reflection coefficient @xmath1
must necessarily separate since @xmath193 is necessarily real valued , @xmath194 \mathcal{r}(r;t ) .
\hphantom{aa } \label{eq : separation variable chiral geometry b}\end{aligned}\ ] ] after insertion of eq .
( [ eq : separation variable chiral geometry b ] ) into the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) one finds that @xmath195 obeys    @xmath196\mathcal{r } \nonumber \\ & & + \big[2\omega r + 3 ( 1-r)(1 - 3r ) \big ] \frac{\partial \mathcal{r}}{\partial r } \nonumber \\ & & + 2r(1-r)^{2}\frac{\partial^{2 } \mathcal{r}}{\partial r^{2 } } \label{eq : fp chiral geometry b}\end{aligned}\ ] ]    with the stationary solution @xmath197    the constant @xmath188 is chosen so that the solution ( [ eq : solution r for closed wire with absorption ( chiral ) ] ) is properly normalized ( this is always possible for @xmath198 ) .
the stationary solution ( [ eq : solution r for closed wire with absorption ( chiral ) ] ) has a square - root singularity at @xmath199 , an essential singularity at @xmath200 , and is normalizable for @xmath201 $ ] .
we shall see below that the probability for @xmath152 to be in the vicinity of 0 is enhanced in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) compared to the standard symmetry class ( [ eq : strong+weak def of std limit ] ) .      in the standard symmetry class ( [ eq : strong+weak def of std limit ] )
the phase @xmath106 and squared magnitude @xmath152 of the reflection coefficient @xmath1 can reasonably be taken to separate in view of the initial condition of a uniformly distributed phase @xmath106 , @xmath202 after insertion of the separation - of - variable ansatz ( [ eq : separation variable std geometry b ] ) into the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) one finds that @xmath195 obeys    @xmath203\mathcal{r } \nonumber \\ & & + \big [   2\omega r + ( 1-r)(1 - 5r )   \big ]   \frac{\partial \mathcal{r}}{\partial r } \nonumber\\ & & + r(1-r)^{2}\frac{\partial^{2 } \mathcal{r}}{\partial r^{2 } } \label{eq : fp standard limit closed both sides with abs}\end{aligned}\ ] ]    with the stationary and normalized solution @xmath204    the stationary solution ( [ eq : solution r for closed wire with absorption ( std ) ] ) has an essential singularity when @xmath200 . the stationary solution ( [ eq : solution r for closed wire with absorption ( std ) ] )
was found by pradhan and kumar in ref .  .
when the initial condition on the probability distribution of @xmath106 is not that of a uniform distribution , it is presumed that the deviations of the solution to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) from the separation - of - variable ansatz ( [ eq : separation variable std geometry b ] ) with the stationary solution ( [ eq : solution r for closed wire with absorption ( std ) ] ) are short transients .
we do not know of an analytical verification of this assumption .
there is no compelling reason to believe that the squared amplitude @xmath152 and the phase @xmath106 of the reflection coefficient @xmath1 separate in the crossover regime .
however , short of an explicit exact solution to the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) in the crossover regime , we shall nevertheless pursue a strategy relying on an approximation built on a separation - of - variable ansatz that allows us to interpolate between the chiral and standard symmetry classes .
one possible alternative to the separation - of - variable ansatz that we tried is a perturbative expansion about the standard symmetry class ( [ eq : strong+weak def of std limit ] ) using a fourier expansion of the joint probability distribution for @xmath152 and @xmath106 .
however , this approach has the drawback that it can not be expected to interpolate all the way to the chiral symmetry class ( [ eq : def chiral regime parameters ] ) in view of the nature of the singularities characterizing the limiting probability distributions .
for this reason we will not present the perturbative approach in this paper .
our starting point is the separation - of - variable ansatz    [ eq : separation - of - variable ansatz in geometry b ] @xmath205 we further assume that the probability distribution function for the phase @xmath106 shows a quick relaxation to its stationary solution , i.e. , @xmath206 is assumed to be @xmath207-independent , @xmath208    insertion of eq .
( [ eq : separation - of - variable ansatz in geometry b ] ) into the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) followed by an integration over @xmath106 yields    [ eq : separation ansatz fp for r geometry b ] @xmath209 \mathcal{r } \nonumber\\ & & + \big [   2\omega   r + ( 1-r)(1 - 5r ) \nonumber \\ & & \hphantom{aaaaaaa } -2\ , \zeta \alpha ( 1-r)(1 - 2r ) \big ] \frac{\partial\mathcal{r}}{\partial r } \nonumber\\ & & + ( 1-\zeta \alpha ) r(1-r)^{2 } \frac{\partial^{2}\mathcal{r}}{\partial r^{2 } } \label{eq : resulting ode for r if separation ansatz}\end{aligned}\ ] ] where the real constant @xmath210 is the ratio of two fourier expansion coefficients , @xmath211 a stationary solution to eq .
( [ eq : resulting ode for r if separation ansatz ] ) is @xmath212 } { \left .
^{-\zeta \alpha/(1-\zeta \alpha ) } \left(1-r\right)^{2/(1-\zeta \alpha)}}. \label{eq : solution for r if separation ansatz}\end{aligned}\ ] ]    this solution reduces to the known results in the standard [ eq .
( [ eq : solution r for closed wire with absorption ( std ) ] ) ] and chiral [ eq .
( [ eq : solution r for closed wire with absorption ( chiral ) ] ) ] symmetry classes . in the standard symmetry class ( [ eq : strong+weak def of std limit ] )
the probability distribution of the phase is uniform @xmath213 and hence @xmath214 , whereas in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) it is a sum of the two delta functions @xmath215 , @xmath216 , @xmath217 .
there remains considerable arbitrariness in the choice of @xmath218 because many different @xmath218 share the same @xmath210 . in practice
we will choose the stationary solution ( [ eq : sol cross if geo c ] ) since we then recover the true solution in the semiopen geometry of fig .
[ fig : geometries]a as @xmath219 .
we postpone to sec .
[ sec : validity of the separation ansatz for the crossover regime ] the discussion of the accuracy of the approximation implied by the separation - of - variable ansatz ( [ eq : separation ansatz fp for r geometry b ] ) with the choice ( [ eq : sol cross if geo c ] ) for the stationary probability distribution of @xmath106 .      in geometry
[ fig : geometries]c , an incoming plane wave from an ideal lead is partially reflected and partially transmitted by the disordered quantum wire . following abrikosov in ref .
we trade the squared magnitude @xmath152 of the reflection coefficient @xmath1 for the resistance ( inverse of conductance )    @xmath220 \label{eq : def varrho}\end{aligned}\ ] ]    under which @xmath221    obeys the fokker - planck equation    [ eq : master equation in terms of varrho and phi ] @xmath222 where @xmath223 and @xmath224      in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) insertion of the separation - of - variable ansatz @xmath225 \mathcal{q}(\varrho;t)\end{aligned}\ ] ] into the fokker - planck equation ( [ eq : master equation in terms of varrho and phi ] ) yields    [ eq : fp+initial geometry c chiral ] @xmath226 with the initial condition @xmath227    the normalized solution to eq .
( [ eq : fp+initial geometry c chiral ] ) has the integral representation    [ eq : sol if chiral geo c ] @xmath228 with the auxiliary function and variable @xmath229 , \qquad t^{\prime}:=2t.\end{aligned}\ ] ] the integrations over @xmath230 and @xmath231 can be performed in closed form , @xmath232 .
\nonumber\\ & & \end{aligned}\ ] ]    by shifting the minimum of the resistance @xmath233 from 1 to 0 , @xmath234 , we recover the same asymptotics for the probability distribution of @xmath235 computed by stone and joanopoulos in ref .   with a log - normal distribution of the nearest - neighbor hopping amplitudes of the discrete schrdinger equation
when @xmath236 or @xmath237 .
the change of variable @xmath238 yields the probability distribution @xmath239 \nonumber\\ & & \end{aligned}\ ] ] for the conductance @xmath7 in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) .
the same result follows from solving the diffusion ( dmpk ) equation ( [ eq : fp chiral limit ] ) and performing the change of variables @xmath240.@xcite      in the standard symmetry class ( [ eq : strong+weak def of std limit ] ) insertion of the separation - of - variable ansatz @xmath241 into the fokker - planck equation ( [ eq : master equation in terms of varrho and phi ] ) yields    [ eq : abrikosov fp ] @xmath242 with the initial condition @xmath243    as shown by abrikosov in ref .  , the normalized solution to eq .
( [ eq : abrikosov fp ] ) has the integral representation @xmath244 ^ 2/4 t } % \frac{1}{4
t } % \left [ % \ln % \frac{\varrho - u}{u(1-u ) } % \right]^{2 } %   } } { [ u(1-u)(\varrho - u)]^{1/2 } } d_0(u;\varrho ) % \ln % \left [ % \frac{\varrho - u}{u(1-u ) } % \right ] \nonumber \\ & = & \sqrt{\frac{4}{\pi t^{3}}}\ , \int_{y^{\ } _ { \varrho}}^{+\infty } { d}y\ ,    \frac{\,y\ , e^{-(t/4)-(y^{2}/t ) } } { \big ( \cosh^{2}y-\varrho \big)^{1/2 } } \label{eq : sol if std geo c}\end{aligned}\ ] ] with @xmath245 $ ] and @xmath246.\ ] ] the change of variable @xmath238 yields the probability distribution @xmath247 with @xmath248 $ ] for the conductance @xmath7 in the standard symmetry class ( [ eq : strong+weak def of std limit ] ) .      as was already the case in sec .
[ subsubsec : crossover regime b ] we could not solve the fokker - planck equation ( [ eq : master equation in terms of varrho and phi ] ) exactly in the crossover regime .
we thus try the separation - of - variable ansatz @xmath249 we are again assuming that the probability distribution @xmath218 is @xmath207 independent and we made the specific choice for @xmath218 given by eq .
( [ eq : sol cross if geo c ] ) . insertion of eq .
( [ eq : separation variable crossover ansatz geometry c ] ) into the fokker - planck equation ( [ eq : master equation in terms of varrho and phi ] ) followed by an integration over the phase @xmath106 yields    [ eq : separation fp+initial cond crossover c ] @xmath250 with the initial condition @xmath251    the normalized solution to eq .
( [ eq : separation fp+initial cond crossover c ] ) has the integral representation    [ eq : solution cross if separation var geo c ] @xmath252 \nonumber \\ & & \times \int_{0}^{1}{d}u\ , \frac{(\frac{1}{u}-1)^{\zeta\alpha/(1-\zeta\alpha ) } }       { [ u(1-u)(\varrho - u)]^{1/2 } } % % u^{-\frac{1}{2}-\frac{\zeta \alpha}{1-\zeta \alpha } } % % ( 1-u)^{-\frac{1}{2}+\frac{\zeta \alpha}{1-\zeta \alpha } } % % ( \varrho - u)^{-\frac{1}{2 } } \nonumber \\ & & \times \int_{0}^{1}{d}v\ , \frac{e^{-d^{2}(u , v;\varrho)/4t^{\prime } } }       { v^{1/2}(1-v)^{1/(1-\zeta \alpha ) } } d(u , v;\varrho ) \nonumber\\ & & \end{aligned}\ ] ] with the auxiliary function and variable @xmath229 , \quad t^{\prime}:=(1-\zeta \alpha)t .
\hphantom{aa}\end{aligned}\ ] ]    as it should be , we recover eq .
( [ eq : sol if chiral geo c ] ) when @xmath253 and eq .
( [ eq : sol if std geo c ] ) when @xmath254 .
the approximate probability distribution @xmath255 for the conductance @xmath7 in the crossover regime is obtained from eq .
( [ eq : solution cross if separation var geo c ] ) through @xmath256 the dependence on @xmath207 for the mean and variance of the conductance in the crossover regime computed from eq .
( [ eq : prob dist g cross geo c ] ) are depicted in fig .
[ fig : conductance ] . in the chiral symmetry class @xmath217 , @xmath257 , the mean and variance decay like @xmath258 . for @xmath259 and @xmath254 ,
the mean and variance decay exponentially with @xmath207 .
( a )     the mean ( a ) and variance ( b ) of the conductance computed from the probability distribution ( [ eq : prob dist g cross geo c ] ) as a function of @xmath160 for @xmath260 and @xmath217 . , width=283 ]    ( b )     the mean ( a ) and variance ( b ) of the conductance computed from the probability distribution ( [ eq : prob dist g cross geo c ] ) as a function of @xmath160 for @xmath260 and @xmath217 . ,
we have constructed an approximate solution to the fokker - planck equations ( [ eq : master equation in terms of r and phi ] ) and ( [ eq : master equation in terms of varrho and phi ] ) in the crossover regime ( [ eq : crossover regime weak if a ] ) by assuming the separation of the squared magnitude @xmath152 from the phase @xmath106 of the reflection coefficient @xmath1 . in this section
we shall investigate the accuracy of this approximation to the crossover regime ( [ eq : crossover regime weak if a ] ) for very large values of the length @xmath5 of the disordered quantum wire , i.e. , @xmath5 much larger than the localization length @xmath129 , which will be defined below .
to this end , we shall compare the cumulants of the radial coordinate @xmath66 computed from the separation - of - variable ansatz with the same cumulants computed from a large - deviation ansatz introduced by schomerus and titov in refs .   and .
we begin from eq .
( [ eq : large x limit fp x phi all derivatives to right ] ) without absorption and in the crossover regime .
we seek the asymptotic behavior when @xmath261 of the solutions to eq .
( [ eq : large x limit fp x phi all derivatives to right ] ) in the crossover regime .
we try the separation - of - variable ansatz @xmath262 where @xmath218 is the stationary and normalized solution ( [ eq : sol cross if geo c ] ) .
insertion of eq.([eq : separation of varaibale ansatz geo c large x crossover ] ) into eq .
( [ eq : large x limit fp x phi all derivatives to right ] ) with subsequent integration over @xmath106 yields    @xmath263    equation ( [ eq : fp x if large t and crossover ] ) is nothing but the diffusion equation with a constant drift proportional to the amount of chiral symmetry breaking @xmath264 , i.e. , it is solved by @xmath265^{2 }       }       { \left ( 1-\zeta \alpha\right)t } \right ]       } { \sqrt{\pi \left(1-\zeta \alpha\right ) t } } \label{eq : sol x if large t and crossover}\end{aligned}\ ] ]    normalized to the real line @xmath266 . in the approximation ( [ eq : sol x if large t and crossover ] ) , the only nonvanishing cumulants are the first and second cumulants ,    [ eq : cumulants x in cross geo c ] @xmath267    provided @xmath268 ( i.e. , away from the chiral symmetry class ) , the random variable @xmath66 becomes self - averaging in the thermodynamic limit .
it then makes sense to identify the localization length @xmath129 through @xmath269 which , together with eq .
( [ eq : first cumulant x cross geo c ] ) , gives @xmath270 in the standard symmetry class ( [ eq : strong+weak def of std limit ] ) @xmath254 and the localization length is twice the mean free path @xmath271 as is well known.@xcite in the chiral symmetry class ( [ eq : def chiral regime parameters ] ) @xmath253 and the localization length defined by eq .
( [ eq : def xi cross geo c ] ) is not a finite number anymore but a random variable as the right - hand side of eq .
( [ eq : first cumulant x cross geo c ] ) vanishes whereas the right - hand side of eq .
( [ eq : second cumulant x cross geo c ] ) remains finite .
for the crossover regime ( [ eq : crossover regime weak if a ] ) , the constant @xmath210 implicitly depends on @xmath65 .
the asymptotic behavior of @xmath210 for small @xmath272 can be evaluated as , with the probability distribution ( [ eq : sol cross if geo c ] ) , @xmath273 since a logarithmic divergence occurs both in the @xmath66 and @xmath274 integrals in the numerator as we take @xmath275 , while it occurs only in the @xmath274 integral in the denominator around @xmath276 .
it is then natural to identify @xmath277 with the typical localization length@xcite @xmath278 when approaching the chiral symmetry class ( [ eq : def chiral regime parameters ] ) through the crossover regime ( [ eq : crossover regime weak if a ] ) .
we have also computed the dependence of @xmath210 on @xmath65 by evaluating eq .
( [ eq : alpha vs varepsilon ] ) numerically . in fig .
[ fig : alpha ] @xmath279 is plotted as a function of @xmath65 in the crossover regime ( [ eq : crossover regime weak if a ] ) .
we observe that @xmath279 asymptotically behaves as @xmath280 in the vicinity ( @xmath281 ) of the chiral symmetry class ( [ eq : def chiral regime parameters ] ) .
numerical evaluation of @xmath279 as a function of @xmath65 in the crossover regime ( [ eq : crossover regime weak if a ] ) .
the integral @xmath210 defined in eq .  (
[ eq : def alpha ] ) is evaluated numerically for the stationary distribution of the phase ( [ eq : sol cross if geo c ] ) .
the line representing @xmath282 is a guide to the eyes .
, width=283 ]      schomerus and titov in refs . and devised a systematic method to compute all the cumulants of the conductance and the ldos to leading order in @xmath283 . in practice , they carried out explicitly the computation of the first four cumulants .    to draw a connection to their work
we note that whereas the parametrization @xmath284 of the magnitude of the reflection coefficient is the natural one from a geometrical point of view according to eqs .
( [ eq : polar decomposition s+m ] ) , ( [ eq : fp std limit ] ) , and ( [ eq : fp chiral limit ] ) , it is    @xmath285    that behaves statistically as the logarithm of the conductance @xmath7 ( or the logarithm of the ldos @xmath6 according to the borland conjecture , see ref .  ) in the regime defined by the condition @xmath286    we denote by @xmath287 the joint probability distribution of @xmath230 and @xmath106 , and introduce @xmath288 .
we infer from eq .
( [ eq : large x limit fp x phi all derivatives to right ] ) that @xmath287 obeys the fokker - planck equation    [ eq : large u limit fp u phi all derivatives to right ] @xmath289 where @xmath290    in the large @xmath230 limit and in the absence of absorption .
the fact that @xmath291 can be written as a total derivative plays an essential role when solving , for any given purely imaginary @xmath292 and any given positive integer @xmath293 , the eigenvalue equation    [ eq : large deviation ansatz ] @xmath294 f^{\ } _ k \label{eq : large deviation ansatz a}\end{aligned}\ ] ] for the eigenvalue @xmath295 and the normalized eigenfunction @xmath296 of @xmath106 with periodicity @xmath111 that results from insertion of the large - deviation ansatz @xmath297    into eq .
( [ eq : large u limit fp u phi all derivatives to right ] ) .    for large @xmath298 , the integration over @xmath292 is dominated in eq .
( [ eq : large deviation ansatz b ] ) by the region of size @xmath299 close to the origin @xmath300 .
assume that for a nonvanishing @xmath292 of order @xmath301 , all eigenvalues @xmath295 are real valued , of descending order , and that the largest eigenvalue @xmath302 is separated from @xmath303 by a finite gap .
consequently , truncation of the summation over @xmath293 to the single term @xmath304 on the right - hand side of eq .
( [ eq : large deviation ansatz b ] ) produces an exponentially small error @xmath305 .
\label{eq : large deviation ansatz c}\end{aligned}\ ] ]    next , assume the expansions @xmath306 where the expansion coefficient @xmath307 encodes the @xmath308th cumulant of @xmath230 , @xmath309 insertion of the expansion ( [ eq : expansion mu s and f s ] ) into the eq .
( [ eq : large deviation ansatz a ] ) for @xmath304 can be solved iteratively for @xmath310 . for @xmath311 , eq .
( [ eq : large deviation ansatz a ] ) reduces to @xmath312 whose normalized solution is nothing but the stationary and normalized solution to eq .
( [ eq : fp if closed at one end and no absorption ] ) . for @xmath313 , eq .
( [ eq : large deviation ansatz a ] ) reduces to @xmath314 the expansion coefficient @xmath315 is obtained from integrating eq .
( [ eq : ode order n=1 ] ) over @xmath106 as a result of @xmath291 being a total derivative ,    [ eq : solving for 1 cumulant u ] @xmath316 where @xmath317 [ eq : computation alpha00 ]    after comparison of eqs .
( [ eq : relating cumulants u and mun0 ] ) and ( [ eq : solving for 1 cumulant u ] ) on the one hand with eqs .
( [ eq : first cumulant x cross geo c ] ) and ( [ eq : def alpha ] ) on the other hand , we conclude that the first cumulant of the logarithm of the conductance is computed exactly to leading order in @xmath318 within the separation - of - variable ansatz for the crossover regime ( [ eq : crossover regime weak if a ] ) . as a corollary , the crossover ( [ eq : typical localization length in cross geo c ] ) for the typical localization length
is given exactly by the separation - of - variable ansatz .
carrying on the iteration to compute the expansion coefficients @xmath307 , however , we find that the separation - of - variable ansatz breaks down for @xmath319 as it does not agree anymore with the large - deviation ansatz from schomerus and titov in refs . and .
( 210,160)(-210,-120 ) ( -190 , 20)(1,0)190 ( -190,-10)(1,0)190 ( 0,-10)(0,1 ) 30 ( -190,-10)(-2,-1)10 ( -190 , 20)(-2 , 1)10    ( -135 , 5)(0 , 10 ) 180 ( -135 , 15 ) ( 1 , 0 ) 20 ( -135,-5 ) ( 1 , 0 ) 20 ( -128 , 5)@xmath320    ( -65 , 5 ) ( 0,10)-180 ( -65 , 15)(-1 , 0 ) 20 ( -65,-5)(-1 , 0 ) 20 ( -80 , 5)@xmath321    ( -100,-10)(0,1)30 ( -103,-20)@xmath274 ( -190,-30 ) ( 1,0)190 ( 0,-30)(-1,0)190 ( -105,-40)@xmath5    ( -210,35)@xmath0    ( -200 , -60)(1,0)200 ( -200,-90)(1,0)200 ( -200,-90)(0,1 ) 30 ( 0,-90)(0,1 ) 30    ( -135,-75)(0 , 10 ) 180 ( -135,-65 ) ( 1 , 0 ) 20 ( -135,-85 ) ( 1 , 0 ) 20 ( -128,-75)@xmath320    ( -65,-75 ) ( 0,10)-180 ( -65,-65)(-1 , 0 ) 20 ( -65,-85)(-1 , 0 ) 20 ( -80,-75)@xmath321    ( -100,-90)(0,1)30 ( -103,-100)@xmath274 ( -200,-110 ) ( 1,0)200 ( 0,-110)(-1,0)200 ( -105,-120)@xmath5    ( -210,-50)@xmath3
in this section , the probability distributions @xmath322 and @xmath323 of the ldos @xmath6 in the geometries of figs .  [
fig : bc for ldos]a and [ fig : bc for ldos]b , respectively , are calculated based on the probability distribution function of the reflection coefficient @xmath192 obtained in sec .
[ sec : solutions to the functional renormalization group equations ] .      for any given realization of the disorder potentials
@xmath97 , the ldos @xmath324 is defined by @xcite @xmath325 where the matrix elements @xmath326 in real space of the green function are obtained from solving @xmath327 g(y , z)= \delta(y - z)\sigma^{\ } _ { 0 } , \nonumber\\ & & e=\mathcal{e}^{\prime}+{i}\mathcal{e}^{\prime\prime } , \qquad \mathcal{e}^{\prime\prime}\geq0 , \nonumber\\ & & ( \eta:\hbox { positive infinitesimal } ) \label{eq : def green fct}\end{aligned}\ ] ] with the boundary conditions associated to the geometries in fig .
[ fig : bc for ldos ] . as long as @xmath274 is far away from the ends of the disordered quantum wire , say in the middle , the statistical properties of @xmath324 should be independent of @xmath274
. however , the statistical properties of @xmath324 depend strongly on the boundary conditions at the ends of the disordered quantum wire.@xcite for completeness , we will understand under the chiral limit of the ldos ( [ eq : def meso ldos ] ) the property that the green function ( [ eq : def green fct ] ) anticommutes with @xmath26 , i.e. , @xmath328 whereas we will understand under the standard limit of the ldos ( [ eq : def meso ldos ] ) the conditions on the green function ( [ eq : def green fct ] ) that    [ eq : def std limit from green fct ] @xmath329 or @xmath330 or @xmath331    schomerus _ et al .
.   expressed the probability distribution of the ldos in terms of the probability distributions of the reflection coefficient @xmath320 ( @xmath321 ) for outgoing plane waves reflected from the region to the left ( right ) of @xmath274 for the continuous nonrelativistic schrdinger equation in the standard symmetry class ( [ eq : def std limit from green fct ] ) .
here we extend their analysis to the continuous relativistic hamiltonian ( [ eq : def our anderson localization problem a ] ) . a crucial observation in ref .
is that , for weak disorder , the green function can be expanded in terms of free scattering states within a small interval of order of the mean free path @xmath61 since one can assume that there are no impurities in this interval .
the solution to eq .
( [ eq : def green fct ] ) for @xmath274 and @xmath292 belonging to such an impurity - free interval can be written as a linear combination of @xmath332 , @xmath333 , @xmath334 , @xmath335 , except for the discontinuity at @xmath336 ,    @xmath337    where @xmath338 , @xmath339 , @xmath340 is the heaviside function that vanishes when @xmath341 and equals one when @xmath342 , and we have momentarily reinstated the fermi velocity @xmath9 and @xmath10 in the relativistic dispersion relation @xmath343 .
the relative amplitudes among @xmath344 as well as the factor @xmath345 are determined by taking into account multiple scattering from the left and right boundaries of an interval free of impurities .
combining eqs .
( [ eq : def green fct ] ) and ( [ eq : plane wave expansion green ] ) reproduces in the relativistic limit the mesoscopic relation @xmath346 between the ldos at @xmath274 and the reflection coefficients @xmath320 and @xmath321 from the disordered region to the left and right of @xmath274 , respectively , derived by schomerus _
et al . _  in ref .
for the nonrelativistic continuous schrdinger equation . in the nonrelativistic case , the diagonal matrix element of the green function in real space gives the microscopic local density of states .
the microscopic local density of states exhibits @xmath347 oscillations.@xcite the ldos ( [ eq : ldos ] ) is obtained after smearing out these oscillations .
the effect of the relativistic approximation ( [ eq : def our anderson localization problem a ] ) is to remove the @xmath347 oscillations from the outset .
henceforth let us measure the ldos @xmath324 in units of @xmath348 , the dos of a clean wire .    the probability distribution function @xmath349 of the ldos @xmath324 defined by eq .
( [ eq : ldos ] ) is now simply given by    @xmath350    where @xmath351 ( @xmath352 ) is the length of the segment of the disordered quantum wire to the left ( right ) of the point @xmath274 where the ldos is measured .
the total length of the disordered quantum wire is evidently given by ( see fig .  [
fig : bc for ldos ] ) @xmath353 the dependence of the probability distribution for the ldos on the boundary conditions enters through the dependence of the probability distributions @xmath354 and @xmath355 on the boundary conditions corresponding to the geometries of fig .
[ fig : geometries].@xcite this dependence on boundary conditions is implied by the dependence on @xmath356 or @xmath357 of the probability distribution of the two reflection coefficients @xmath320 and @xmath321 , respectively .
when a disordered quantum wire of finite length is closed , one must broaden the single - particle energy levels beyond the mean level spacing to make sense of the energy dependence of the mesoscopic ldos .
this is achieved here by introducing a finite absorption @xmath358 , in which case the dependence on @xmath359 or @xmath360 becomes immaterial .
the absorption @xmath358 is switched off when a disordered quantum wire of finite length @xmath5 is connected to a perfect lead .
we consider first the geometry depicted in fig .
[ fig : bc for ldos]a and then the geometry depicted in fig .
[ fig : bc for ldos]b .    ( a )     ( a ) three traces for the probability distribution of the ldos are plotted with @xmath361 ( from the bottom up at @xmath362 ) so as to interpolate between the chiral and standard symmetry classes ( [ eq : def chiral limit from green fct ] ) and ( [ eq : def std limit from green fct ] ) , respectively .
each trace is obtained from numerical integration of eq .
( [ eq : distribution for the ldos ] ) for the semi - open geometry of fig .
[ fig : bc for ldos]a with @xmath363 .
( b ) three traces for the probability distribution of the ldos are plotted with @xmath361 ( from the bottom up at @xmath364 ) so as to interpolate between the chiral and standard symmetry classes ( [ eq : def chiral limit from green fct ] ) and ( [ eq : def std limit from green fct ] ) , respectively .
each trace is obtained from numerical integration of eq .
( [ eq : distribution for the ldos ] ) for the closed geometry of fig .
[ fig : bc for ldos]b with @xmath365 .
we have verified that numerical integration of eq .
( [ eq : distribution for the ldos ] ) agrees with a representation of eq .
( [ eq : distribution for the ldos ] ) in terms of elementary functions when possible .
, width=283 ]    ( b )     ( a ) three traces for the probability distribution of the ldos are plotted with @xmath361 ( from the bottom up at @xmath362 ) so as to interpolate between the chiral and standard symmetry classes ( [ eq : def chiral limit from green fct ] ) and ( [ eq : def std limit from green fct ] ) , respectively .
each trace is obtained from numerical integration of eq .
( [ eq : distribution for the ldos ] ) for the semi - open geometry of fig .
[ fig : bc for ldos]a with @xmath363 .
( b ) three traces for the probability distribution of the ldos are plotted with @xmath361 ( from the bottom up at @xmath364 ) so as to interpolate between the chiral and standard symmetry classes ( [ eq : def chiral limit from green fct ] ) and ( [ eq : def std limit from green fct ] ) , respectively .
each trace is obtained from numerical integration of eq .
( [ eq : distribution for the ldos ] ) for the closed geometry of fig .
[ fig : bc for ldos]b with @xmath365 .
we have verified that numerical integration of eq .
( [ eq : distribution for the ldos ] ) agrees with a representation of eq .
( [ eq : distribution for the ldos ] ) in terms of elementary functions when possible .
, width=283 ]      to calculate the probability distribution of the ldos in a disordered quantum wire connected to a reservoir on the left and closed on the right , we need ( i ) the joint probability distribution @xmath366 for the square @xmath367 of the magnitude of the reflection coefficient @xmath320 and its phase @xmath368 on the segment of length @xmath351 to the left of the point at which the ldos is measured and ( ii ) the probability distribution @xmath369 of the phase @xmath370 of the reflection coefficient @xmath321 on the segment of length @xmath352 to the right of the point at which the ldos is measured .
the calculation of @xmath369 is outlined in sec.[subsec : disordered wire closed on the right - hand side without absorption ] .
the calculation of @xmath366 is outlined in sec.[sec : disordered wire opened at both ends without absorption ] .      the probability distribution of the ldos in the chiral symmetry class ( [ eq : def chiral limit from green fct ] ) is ( see the appendix [ app sec : proof ] )    [ eq : log normal distribution ldos if ch in geo semi - open ] @xmath371 this is the probability distribution of a log - normal distributed random variable with @xmath372    the average ldos is @xmath373 which grows exponentially with @xmath374 . for comparison , the probability distribution of the ldos in the standard symmetry class ( [ eq : def std limit from green fct ] ) , which was obtained in refs .   and
, is also log - normal ,    [ eq : log normal distribution ldos if std in geo semi - open ] @xmath375 , \nonumber\\ & & \\ & & \left\langle \ln\nu \right\rangle^{\ } _ { l^{\ } _ { \mathrm{l}}}= -\frac{l^{\ } _ { \mathrm{l}}}{\ell } , \\ & & \left\langle \left ( \ln\nu + \frac{l^{\ } _ { \mathrm{l}}}{\ell } \right)^2 \right\rangle^{\ } _ { l^{\ } _ { \mathrm{l } } } = \frac{2l^{\ } _ { \mathrm{l}}}{\ell}.\end{aligned}\ ] ]    incidentally , the average ldos is not affected by the disorder in the standard symmetry class , @xmath376 .
the probability distributions for the ldos in the chiral and standard symmetry classes are depicted in fig .
[ fig : distribution of the ldos]a . in the thermodynamic limit @xmath377 , the logarithm of the ldos is self - averaging in the standard symmetry class
whereas this is not the case in the chiral symmetry class .
this difference could have been anticipated from the identification @xmath378 , valid when @xmath379 , with the radial coordinate @xmath66 that obeys the fokker - planck equations ( [ eq : fp std limit ] ) and ( [ eq : fp chiral limit ] ) , respectively .
numerical evaluation of the integrals on the right hand side of eq .
( [ eq : distribution for the ldos ] ) is depicted in fig .
[ fig : distribution of the ldos]a for the crossover regime ( [ eq : crossover regime weak if a ] ) within our approximation which is encoded by the separation - of - variable ansatz ( [ eq : separation variable crossover ansatz geometry c ] ) , ( [ eq : separation fp+initial cond crossover c ] ) , ( [ eq : solution cross if separation var geo c ] ) , and ( [ eq : sol cross if geo c ] ) .      to calculate the probability distribution of the ldos in a disordered quantum closed at both ends with a finite absorption @xmath380 , we need the joint probability distributions @xmath381 and @xmath382 . we can safely ignore the dependence on @xmath383 in the regime @xmath384 as is done in sec .
[ subsec : disordered quantum wire closed on the right - hand side with absorption ] .
the probability distribution of the ldos in the chiral symmetry class ( [ eq : def chiral limit from green fct ] ) is@xcite ( see appendix [ app sec : proof ] )    [ eq : distribution ldos if ch in geo closed ] @xmath385^{2 } \frac{1}{\nu } \exp\left [ -\frac{\omega(1+\nu^{2})}{2\nu } \right ] k^{\ } _ { 0 } \left ( \frac{\omega|1-\nu^{2}|}{2\nu } \right ) \label{eq : distribution ldos if ch in geo closed a}\end{aligned}\ ] ] where @xmath386 , @xmath310 are the modified bessel functions and has the asymptotics @xmath387    for comparison , the probability distribution of the ldos in the standard symmetry class ( [ eq : def std limit from green fct ] ) , which was obtained in refs .   and ,    @xmath388       }       { \sqrt{u-(1+\nu^2)/(2\nu ) }       } \left [ u k^{\ } _ 0\left(\omega\sqrt{u^2 - 1}\right ) + \sqrt{u^2 - 1 } k^{\ } _ 1\left(\omega\sqrt{u^2 - 1}\right ) \right],\end{aligned}\ ] ]    has the asymptotics @xmath389    the probability distributions for the ldos in the chiral and standard symmetry classes are depicted in fig .
[ fig : distribution of the ldos]b .
spectral weight from the tails @xmath390 and @xmath391 of the probability distribution in the standard symmetry class is redistributed around @xmath392 in the chiral symmetry class .
consequently , in the limit @xmath219 , the mean value of the mesoscopic ldos in the chiral symmetry class becomes @xmath393 ^ 2 } \int^\infty_1{d}\nu\left(1+\frac{1}{\nu^2}\right ) \nonumber \\ & & \times \exp\left[-\frac{\omega}{2}\left(\nu+\frac{1}{\nu}\right)\right ] k_0\left(\frac{\omega}{2}\left(\nu-\frac{1}{\nu}\right)\right ) \nonumber \\ & { \buildrel \omega\to0   \over = } & \frac{1}{(\omega/2)[\ln(\omega/2)]^2}.\end{aligned}\ ] ] this result is identical to the smeared dyson singularity of the dos in the chiral symmetry class , @xmath394 ^ 2}.\end{aligned}\ ] ]      numerical evaluation of the integrals on the right hand side of eq .
( [ eq : distribution for the ldos ] ) is depicted in fig .
[ fig : distribution of the ldos]b for the crossover regime ( [ eq : crossover regime weak if a ] ) within our approximation which is encoded by the separation - of - variable ansatz ( [ eq : separation - of - variable ansatz in geometry b ] ) , ( [ eq : separation ansatz fp for r geometry b ] ) and ( [ eq : sol cross if geo c ] ) .
the probability distribution of the mesoscopic local density of states ( ldos ) @xmath6 for a strictly one - dimensional problem of anderson localization with chiral symmetry ( the chiral symmetry class ) was computed in closed form for two simply connected geometries assuming a weak disorder . as is the case when the chiral symmetry is maximally broken ( the standard symmetry class ) , the probability distribution of @xmath6 strongly depends on boundary conditions . in a semiopen geometry we found that the probability distribution of @xmath6 is log - normal with a vanishing mean as opposed to log - normal with a finite mean in the standard symmetry class . in a closed geometry with absorption we found that the probability distribution of @xmath6 has a double - peak structure whereby the second peak turns out to be a logarithmic singularity at @xmath392 .
we verified that the smeared dyson singularity is reproduced by the mean value of the ldos when the absorption is sufficiently small .
furthermore , we found one exact and proposed two approximate solutions to the functional renormalization group equation obeyed by the joint probability distribution for the squared modulus and the phase of the reflection amplitude of a finite and possibly dissipative wire with semiopen and open boundary conditions , respectively , from which we could extract the approximate crossover of the ldos @xmath6 and conductance @xmath7 between the chiral and standard symmetry classes . we could show that our approximation is exact for the first cumulant but fails to describe higher cumulants of the logarithms of @xmath6 and @xmath7 in the crossover regime .
we are indebted to p.  w.  brouwer for an initial collaboration on this project and for his generous feedback at later stages .
we would also like to thank h.  schomerus and m.  titov for useful comments .
this work was supported in part by japan society for the promotion of science ( s.r . ) .
in the appendix we prove eqs .
( [ eq : prob dis ldos chiral semi - open ] ) and ( [ eq : distribution ldos if ch in geo closed a ] ) .
we make use of the fact that the probability distributions for the squared modulus and phases of the reflection coefficients @xmath395 factorize in the chiral symmetry class . we have chosen here to represent @xmath396 by @xmath397 . we thus need to compute                                              t.  nagao and k.  slevin , j.  math .
phys .  * 34 * , 2075 ( 1993 ) ; j.  j.  m.  verbaarschot and i.  zahed , phys .
lett .  * 70 * , 3852 ( 1993 ) ; s.  hikami and a.  zee , nucl .
b * 408 * , 415 ( 1993 ) ; a.  v.  andreev , b.  d.  simons , and n.  taniguchi , nucl .
b * 432 * , 487 ( 1994 ) .
l.  i.  deych , d.  zaslavsky , and a.  a.  lisyansky , phys .
lett .  * 81 * , 5390 ( 1998 ) ; l.  i.  deych , a.  a.  lisyansky , and b.  l.  altshuler , _ ibid .
_ * 84 * , 2678 ( 2000 ) ; l.  i.  deych , m.  v.  erementchouk , a.  a.  lisyansky , and b.  l.  altshuler , _ ibid .
_ * 91 * , 096601 ( 2003 ) .
an alternative derivation of the fokker - planck equation ( [ eq : master equation in terms of r and phi ] ) is the following .
first , a discrete langevin process with the discrete time step `` @xmath73 '' is written down by expanding @xmath409 in eq .
( [ eq : composition law reflection amplitude ] ) up to second order in @xmath97 .
second , a continuous fokker - planck equation for the probability distribution @xmath410 is derived from the discrete langevin process using the methods of appendix a3.5 of ref .  , say . here , @xmath411 denotes averaging over @xmath97 restricted to the thin slice of length @xmath73 .