Critical parameters for the one-dimensional systems with long-range correlated disorder

Critical parameters for the one-dimensional systems with long-range correlated disorder

Yi Zhao, Suqing Duan, and Wei Zhang 111Author to whom any correspondence should be addressed, Email: Institute of Applied Physics and Computational Mathematics, P.O. Box 8009(28), Beijing 100088, P. R. China

We study the metal-insulator transition in a tight-binding one-dimensional (1D) model with long-range correlated disorder. In the case of diagonal disorder with site energy within and having a power-law spectral density , we investigate the competition between the disorder and correlation. Using the transfer-matrix method and finite-size scaling analysis, we find out that there is a finite range of extended eigenstates for , and the mobility edges are at . Furthermore, we find the critical exponent of localization length () to be . Thus our results indicate that the disorder strength determines the mobility edges and the degree of correlation determines the critical exponents.

72.15.Rn, 71.30.+h, 73.20.Jc

I Introduction

Anderson’s localization theory points out that in a system with uncorrelated diagonal disorder, all one-electron states are spatially localized with an exponentially decaying envelope with a characteristic localization length , when the disorder strength is larger than a critical value.1 () From the scaling theory,2 () it is also well-known that an infinitesimal disorder can cause localization of all states in one and two dimensional systems in the thermodynamic limit. In recent years, it was found that spatial correlation of disorder potential played an important role in the nature of long-range charge transport of low-dimensional systems. In the presence of short-range correlated disorder, there exists a set of discrete resonant energy levels of extended states.3 (); 4 (); 5 (); 6 (); 7 (); 8 (); 9 (); 10 (); 11 (); 12 () For example, a single extended eigenstate was found in the random dimer model,4 (); 8 (); 12 () which was verified by measuring microwave transmission on the semiconductor superlattices with intentional correlated disorder.13 () These models with short-range correlated disorder do not possess a true disorder induced metal-insulator transition in the thermodynamic limit. More recently, a one-dimensional (1D) disordered model with long-range correlation considered by de Moura and Lyra has arisen a great interest, as it can result in a continue band of extended states under appropriate conditions.14 () The transition from localization to delocalization occurred at the critical points was examined experimentally in a single-mode waveguide with inserted correlated scatters.15 () More recent experiments in ultra-cold atoms/Bose-Einstein condensate push the research in this direction further. bec ()

In the disordered systems, the critical parameters, such as the critical disorder strength , mobility edge , critical exponents , (defined by the dependence of localization length on or near the critical points or ) are of particular importance and interest, since they determine the nature of the localization-delocalization transitions and the universal properties of disordered systems. Even for three-dimensional (3D) Anderson model which has been studied for more than forty years, the accurate determination of the critical parameters is still under recent investigations,16 () where the critical parameters for different lattices are calculated and it was found that is close to . While for 3D systems with scale-free disorder, it was found that and the extended Harris criterion is obeyed. 16b ()

In spite of many studies3 (); 4 (); 5 (); 6 (); 7 (); 8 (); 9 (); 10 (); 11 (); 12 (); 13 (); 14 (); 15 (); 17 (); 18 (); 19 (); 20 () on the 1D systems with correlated disorder, many of them are in a perturbative level and there are still several important questions to be answered, such as the positions of the mobility edges and the value of the critical exponent . The determination of these critical parameters and the comparison with the results of 3D models mentioned above will definitely deepen our understanding on disordered systems, especially the dependence on dimensionality and correlation. Moreover, the correlated disordered system is a kind of system which is intermediate between ordered system and pure disordered system. The investigation of the universal properties of this type system and comparison with quasi-periodic system is of fundamental importance. In addition, the accurate values of these critical parameters (mobility edges) have important applications in charge transport in low dimensional (correlated) disordered systems, such DNA molecules,21 (), Bose-Einstein condensate in 1-D optical lattices with speckle potential.bec (). In the present paper, we calculate these critical parameters for the 1D systems with the correlated diagonal disorder. And we also resolve some inconsistencies in the literatures.14 () One of the advantages of our study is that we are able to investigate localization-delocalization transition and obtain the accurate critical parameters in 1D system with large size and the boundary effects is negligible, unlike the 3D system where all the sizes in three dimensions can not be very large due to the computational limitation.

We first investigate the localization properties of the 1D tight-binding model with long-range correlated diagonal disorder (the hopping constant is set to be unity). The diagonal on-site energies are distributed in with being the strength of disorder. The site energies have a power-law spectral density with . The function is the Fourier transform of the two-point correlation function and is related to the wavelength of the undulations on the random parameter landscape by .14 () We find that the critical disorder width at fixed energy for . Moreover, we can determine the positions of the mobility edges , which separate the extended and localized energy eigenstates. The effective energy band width of the extended states has the linear relationship for . Thus, a phase diagram in the space is obtained. We also discuss the -dependence of the critical exponent of the localization length of the eigenstates. Interestingly, we find that the disorder strength determines the mobility edges and the degree of correlation determines the critical exponents. Our results also indicate that the nature of a correlated disordered system is somehow between that of the pure random (without any correlation) and pseudo-random (quasi-periodic) systems. Our approach is non-perturbative and our results are hard to be obtained by perturbative calculations.

The organization of the paper is as follows. In the next section, we introduce the 1D disordered systems with long-range correlation and the basic approach in our calculation. In section III, we present our results of the mobility edge and critical exponents. The paper is summarized in section IV.

Ii Model and the Basic approach

We consider a 1D tight-binding Hamiltonian with long-range correlated disorder


where is the Wannier state localized at site with on-site energy , and is the nearest-neighbor hopping amplitude. We set for simplicity. The on-site energies can be constructed by the Fourier filtering method23 (); 24 (); 25 () as follows: (i) generate a random sequence with a Gaussian distribution; (ii) get its Fourier components using the fast-Fourier transformation method; (iii) establish a new sequence by the relation ; (iv) construct the sequence , which is the inverse Fourier transform of ; (v) adjust the scale of the sequence reaching to . It can be checked that the disordered on-site energy is long-range correlated with the power spectrum .26 () The exponent characterizing the correlation reflects the roughness of the energy sequence. The larger the value of is, the smoother the energy landscape is.

The Schrödinger equation for the wavefunction amplitudes is


Here, is the eigenenergy. Our interest focuses on the critical behavior in the thermodynamical limit. A useful way to find the critical points is based on the finite size analysis of the normalized localization length :26 () If increases as the system size , which implies that would exceed when is large enough, the electron stays on a extended state; Otherwise, the state is localized, because is a decreasing function of and smaller than unity in the thermodynamic limit. So the dependent of can give the information of the localization-delocalization transition. Thereby, we can use this quantity to determine the critical disorder strength at fixed energy and the critical energies at fixed disorder strength.

In the regime of weak disorder, it was found17 () that the localization length for a state of energy (not near band center and band edge) is determined by the correlation function of the disorder potential


where , , , . For our disordered system with long range correlation with the spectral density , the localization length has the normalized form , for ; , for . Thus we see that in the weak disorder regime the states are localized for and is a critical value for the appearance of extended states.

To determine the full phase diagram and study the critical parameters, we use the transfer matrix method. Using the two-component vector , where the superscript denotes the transpose, Eq. (2) can be written in a recursive form




We obtain the transfer-matrix equation , with . The localization length at energy is the inverse of the Lyapunov coefficient , which is the largest eigenvalue of the limiting matrix .27 () The large numbers appeared in the calculations have been took care of by dividing a large number and multiplying it at the end of calculation. The reorthogonalization method reorth () has also been used and the same critical parameters are obtained. All the values of in this paper are based on a geometrical mean of disorder configurations.

Iii Results and discussions

iii.1 Critical disorder strength and mobility edges

Figure 1: The normalized localization length as a function of the system size for fixed . is taken as and in and , respectively.

We perform the calculations of the normalized localization length for different values of and energies. Figure 1 shows plots of the normalized localization length versus for , and , respectively. The system size ranges from to . We notice that is monotonically decreasing with for ; while the picture is different for : the values of rise with increasing for , whereas the values of decline with increasing for . The size independence of at indicates a critical point of a continuous phase transition.

In Fig. 2, we display the dependence of on the disorder strength for the state of with typical values of . One may observe that there emerges an intersection point for different system size when , which corresponds to the critical point. This critical point doesn’t exist for . From the inset of Fig. 2, we could explicitly note the different dependence of on the system size in the vicinity of . It tells us the state with is localized for and is extended for and .

Figure 2: The normalized localization length as a function of the disorder width at fixed for systems with different size. and . Insert: The detail behavior of near the critical point .
Figure 3: The normalized localization length as a function of at fixed for . Detail behavior of near the critical point is shown in the inset.

Table I(a) lists the critical values of for various energies. It can be seen that the critical disorder strength is independent of whenever . At the same time, we find a simple relationship between and as:


For the special state at band center with energy , the critical value , equals to the bandwidth. This conclusion for state was also reached by H. Shima et al.26 ()

Using the same method, we also obtain the critical energies at fixed , as shown in Fig. 3. The results are given in Table I(b) for . The critical energies/the mobility edges can be determined by the following equations:


The eigenstates between and are delocalized. Therefore, the effective bandwidth of the delocalized states is (for )


Figure 4 shows the phase diagrams in the and plane. The critical value is also consistent with the perturbative calculation discussed above. We emphasize that mobility edges or the effective bandwidth is only dependent on the disorder strength as long as . Also, in the limit of , the effective bandwidth becomes which is the well-known result in the 1D ordered crystal lattice. It is interesting to compare with the pseudo random systems with quasi-periodic on-site energy , Q an irrational number, . In this model, the critical value and the mobility is also at , independent of .28 () The linear relations between the mobility edge (critical disorder strength ) and disorder strength (energy ) are simple and neat, further more they have a direct consequence on the critical exponents as will be discussed in the next subsection.

I(a) Critical values of disorder strength for different and energies








2.50 1.000(1) 2.001(1) 3.001(2) 4.001(2) 3.001(2) 2.001(2) 1.000(2)
3.00 1.0000(4) 2.0000(5) 3.0001(5) 4.0001(6) 3.0001(6) 2.0001(6) 1.0001(4)
3.50 1.0000(1) 2.0000(1) 3.0000(1) 4.0000(1) 3.0000(1) 2.0000(1) 1.0000(1)
4.00 1.0000(1) 2.0000(1) 3.0000(1) 4.0000(1) 3.0000(1) 2.0000(1) 1.0000(1)

I(b) Critical values of energies for different and disorder strength

W=3.0 W=2.0 W=1.0

2.50 -0.5004(6) 0.5003(5) -1.0001(4) 1.0002(5) -1.5001(5) 1.5001(5)
3.00 -0.5000(2) 0.5001(2) -1.0000(3) 1.0000(2) -1.5000(1) 1.5001(2)
3.50 -0.5000(1) 0.5000(1) -1.0000(1) 1.0000(1) -1.5000(1) 1.5000(1)
4.00 -0.5000(1) 0.5000(1) -1.0000(1) 1.0000(1) -1.5000(1) 1.5000(1)
Table 1: The values of critical disorder strength and critical energies (mobility edges) obtained by the coordinates of the intersection points in the cases . The errors for the last digits are listed in the parentheses.
Figure 4: Phases diagrams in the and plane. The delocalized and localized states are separated by and . Here . We take in plane as an example.

iii.2 Critical exponents based on finite-size Scaling Analysis

Figure 5: Scaling function of for . Data from various chains, whose size ranges from to , overlap each other. As an example, we give the optimized values of and in the case of .

II(a) Critical values for disorder strength and exponent for different and energy

E=0.0 E=0.5 E=1.0 E=1.5

2.25 4.0015(5) 2.15(5) 3.0013(3) 2.13(4) 2.0008(7) 2.12(6) 1.0002(6) 2.12(6)
2.50 4.0007(5) 1.86(4) 3.0007(3) 1.83(3) 2.0004(3) 1.84(3) 1.0003(3) 1.81(3)
2.75 4.0003(3) 1.65(2) 3.0003(2) 1.65(2) 2.0001(2) 1.66(2) 1.0000(2) 1.66(3)
3.00 4.0002(1) 1.51(1) 3.0001(1) 1.50(1) 2.0001(1) 1.50(1) 1.0000(1) 1.50(2)
3.25 4.0000(1) 1.40(1) 3.0000(1) 1.40(1) 2.0000(1) 1.40(1) 1.0000(1) 1.41(2)
3.50 4.0000(1) 1.31(1) 3.0000(1) 1.30(1) 2.0000(1) 1.31(1) 1.0000(1) 1.31(1)
3.75 4.0000(1) 1.24(1) 3.0000(1) 1.24(1) 2.0000(1) 1.24(1) 1.0000(1) 1.25(2)
4.00 4.0000(1) 1.20(2) 3.0000(1) 1.19(1) 2.0000(1) 1.19(1) 1.0000(1) 1.20(2)
4.25 4.0000(1) 1.13(2) 3.0000(1) 1.13(1) 2.0000(1) 1.13(1) 1.0000(1) 1.14(1)
4.50 4.0000(1) 1.10(2) 3.0000(1) 1.10(1) 2.0000(1) 1.10(2) 1.0000(1) 1.11(2)

II(b) Critical values for energy and exponent for different and disorder strength

W=3.0 W=2.0 W=1.0

2.25 0.5007(4) 2.15(4) 1.0004(3) 2.14(3) 1.5003(4) 2.12(4)
2.50 0.5003(2) 1.84(2) 1.0002(2) 1.85(2) 1.5001(2) 1.84(3)
2.75 0.5001(1) 1.66(1) 1.0001(1) 1.65(1) 1.5000(1) 1.66(1)
3.00 0.5001(1) 1.51(1) 1.0001(1) 1.51(1) 1.5000(1) 1.52(1)
3.25 0.5000(1) 1.40(1) 1.0000(1) 1.41(1) 1.4999(1) 1.43(2)
3.50 0.5000(1) 1.30(1) 1.0000(1) 1.31(1) 1.5000(1) 1.32(1)
3.75 0.5000(1) 1.24(1) 1.0000(1) 1.24(1) 1.5000(1) 1.25(1)
4.00 0.5000(1) 1.19(1) 1.0000(1) 1.19(1) 1.5000(1) 1.19(2)
4.25 0.5000(1) 1.14(1) 1.0000(1) 1.14(1) 1.5000(1) 1.15(1)
4.50 0.5000(1) 1.10(1) 1.0000(1) 1.10(2) 1.5000(1) 1.11(2)
Table 2: The values of critical exponent and critical points and got from the fitting with Eq. (10) for . The errors are obtained from bootstrap method with 95 confidence.

In the light of the finite-size scaling analysis, the normalized localization length for finite system may be related to the localization length for infinite system by the following scaling law:27 ()


where the parameter varies as or in the vicinity of the critical point. Unlike three dimensional systems with usually small values of scaling variable (i.e. the size of cross section), in our one-dimensional systems, the irrelevant scaling exponents ander () are of little effect due to scaling variable with very large values, i.e., the size of the system.

Firstly, we expand the scaling function in polynomial form:


Here defined as or , is a nonlinear combination of and with the parameters and to be determined. Secondly, based on Eq. (10), we can obtain the values of and by fitting which we got from previous calculation. The expansion in the scaling function must be carried up to , that makes sure the degree of confidence reaches to . As shown in Fig. 5, we have an excellent scaling curve, accompanied by the optimized values of and . The critical values for obtained from finite-size scaling analysis are identical to those obtained in section III(A). The values of critical exponent and their errors estimated from the fitting procedure are given in Table II. It is clear that the critical exponent is independent on or . With increasing the correlation (), decreases to . One may note that . Actually it is a direct consequence of the linear relation between the disorder strength and mobility edges. From and Eqs. (6) and (7), one finds that . Thus one has . It is interesting to compare our results with those in 3D systems. For 3D Anderson model (without correlation) is close , 16 () while in the 3D system with scale free disorder, . 16b ()

Then we find the dependence of on follows the rule:


with is approximately and as seen in Fig. 6. We have found that for infinite large correlation of . It is interesting to note that the case with infinite large long-range correlation () looks like the quasiperiodic system, since the position of the mobility edges and the critical exponents (1) in both systems are the same. 28 () In another limit situation , , which is not easy to be obtained from direct calculation due to the critical slowing down phenomena. This result is apparently different from at mentioned in Ref. 14 (). Our results are based on finite size scaling. The data show very good scaling behavior as shown in Fig. 5. The critical exponents and mobility edges are obtained simultaneously by using fitting formula Eq. (10). The mobility edges obtained in this way are identical to those obtained in section III(A). Therefore our results are reliable. It is interesting that the mobility edges are determined only by the disorder strength and the critical exponents are determined only by the degree of correlation .

Figure 6: The -dependence of the averaged critical exponent (for different values of and ) and the fitting curve from Eq. (11) with and . Inset: as a function of , where the straight slope is set to .

Iv Conclusion

In summary, we have investigated the universal localization properties of 1D tight-binding model where the disorder is long-range correlated with a power-law spectral density , . We have obtained the complete phase diagram and neat results for the critical parameters. There is a finite range of extended eigenstates, and the effective bandwidth decreases linearly with increasing for . The positions of the mobility edges separating localized and extended states depend only on the disorder strength whenever . Using the finite-size scaling analysis, we find the critical exponent . In particular when and when . Compared to the 3D uncorrelated Anderson model, the positions of the mobility edges are apparently different. And the critical exponent is no longer equal to 16 (); ander (). So, the transition in the 1D disordered system with correlation is in a new universal class. The existence of mobility edge and dependence of the critical exponent indicate that the nature of a correlated disordered system is somehow between that of the pure random (without any correlation) and pseudo-random (quasi-periodic) systems. Our work sheds some light on the metal-insulator transitions.

V Acknowledgements

This work is supported in part by the National Natural Science of China under No. 10574017, 10744004, 10874020, National Fundamental Research of China grant under No. 2006CB921400 and a grant of the China Academy of Engineering and Physics.


  • (1) P. W. Anderson, Phys. Rev. 109, 1492 (1958).
  • (2) E. Abrahams, P. W. Anderson, D. G. Licciardello, and T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979).
  • (3) J. C. Flores, J. Phys. Condens. Matter 1, 8471 (1989).
  • (4) D. H. Dunlap, H. L. Wu, and P. Phillips, Phys. Rev. Lett. 65, 88 (1990).
  • (5) H.-L. Wu and P. Phillips, Phys. Rev. Lett. 66, 1366 (1991).
  • (6) S. N. Evangelou and D. E. Katsanos, Phys. Lett. A 164, 456 (1992).
  • (7) A. Bovier, J. Phys. A 25, 1021(1992).
  • (8) S. N. Evangelou and A. Z. Wang, Phys. Rev. B 47, 13126 (1993).
  • (9) J. C. Flores and M. Hilke, J. Phys. A 26, L1255 (1993).
  • (10) M. Hilke, J. Phys. A 27, 4773 (1994).
  • (11) F. C. Lavarda, M. C. dos Santos, D. S. Galvão, and B. Laks, Phys. Rev. Lett. 73, 1267 (1994).
  • (12) J. Heinrichs, Phys. Rev. B 51, 5699 (1995).
  • (13) V. Bellani, E. Diez, R. Hey, L. Toni, L. Tarricone, G. B. Parravicini, F. Domínguez-Adame, and R. Gómez-Alcalá, Phys. Rev. Lett. 82, 2159 (1999).
  • (14) F. A. B. F. de Moura and M. L. Lyra, Phys. Rev. Lett. 81, 3735 (1998).
  • (15) U. Kuhl, F. M. Izrailev, A. A. Krokhin, and H.-J. Stöckmann, Appl. Phys. Lett. 77, 633 (2000).
  • (16) G. Roati, C. D’Errico, L. Fallani, M. Fattori, C. Fort, M. Zaccanti, G. Modugno, M. Modugno and M. Inguscio Nature 453, 895 (2008); J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer and A. Aspect, Nature 453, 891 (2008).
  • (17) A. Eilmes, A. M. Fischer, and R. A. Römer, Phys. Rev. B 77, 245117 (2008).
  • (18) M. L. Ndawana, R. A. Römer, M. Schreiber, Europhys. Lett. 68, 678-684 (2004).
  • (19) F. M. Izrailev and A. A. Krokhin, Phys. Rev. Lett. 82, 4062 (1999).
  • (20) S. Russ, J. W. Kantelhardt, A. Bunde, and S. Havlin, Phys. Rev. B 64, 134209 (2001).
  • (21) Wei Zhang and Sergio E. Ulloa, Phys. Rev. B 69, 153203 (2004).
  • (22) Wei Zhang and Sergio E. Ulloa, Phys. Rev. B 74, 115304 (2006).
  • (23) P. Carpena, P. Bernaola-Galan, P. C. Ivanov, and H. E. Stanley, Nature (London) 418, 955(2002); 421, 764 (2003); D. Holste, I. Grosse, and H. Herzel, Phys. Rev. E 64, 041917 (2001); Wei Zhang and Sergio E. Ulloa, Microelectronics Journal 35, 23 (2004).
  • (24) D. Saupe, in The Science of Fractal Images, edited by H.-O. Peitgen and D. Saupe (Springer, New York, 1988); J. Feder, Fractals (Plenum Press, New York, 1988).
  • (25) C.-K. Peng, S. Havlin, M. Schwartz, and H. E. Stanley, Phys. Rev. A 44, R2239 (1991).
  • (26) S. Prakash, S. Havlin, M. Schwartz, and H. E. Stanley, Phys. Rev. A 46, R1724 (1992).
  • (27) H. Shima, T. Nomura, and T. Nakayama, Phys. Rev. B 70, 075116 (2004).
  • (28) B. Kramer and A. McKinnon, Rep. Prog. Phys. 56, 1469 (1993).
  • (29) A. MacKinnon and B. Kramer, Z.Phys. B 53, 1 (1983).
  • (30) S. Das Sarma, S. He, and X. C. Xie, Phys. Rev. Lett. 61, 2144 (1988); Phys. Rev. B 41, 5544 (1990).
  • (31) K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82, 382 (1999).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description