Mobility edges in one-dimensional finite-sized models with large quasi-periodic disorders

Qiyun Tang(湯起蕓) and Yan He(賀言)

College of Physics,Sichuan University,Chengdu 610064,China

Keywords: quasi-periodic disorders,mobility edges,Aubry-Andre model

1.Introduction

The phenomena of localization[1]have been extensively studied in condensed matter physics for many years.For threedimensional(3D)systems,a threshold of disorder strength exists, above which the Anderson localization will take place.For generic band structures,a closely related notion is the mobility edge[2]which is the critical energy separating the localized and extended energy eigenstates.Since the direct observation of mobility edges in a 3D system is quite challenging,attention has been turned to low-dimensional systems.However, the number of dimensions is a very important factor for Anderson localization.It is argued by the scaling theory[3,4]that for one-dimensional(1D)systems even infinitesimal small disorders will make all eigenstates completely localized.Thus,one cannot find any mobility edges in the 1D band structure.

About forty years ago, the low-dimensional quasiperiodic systems with correlated disorders were proposed and captured a lot of attention.The most famous example is the so-called Aubry-Andre (AA) model.[5,6]In essence, the AA model is just a 1D tight-binding model with incommensurate or quasi-periodic on-site potentials.By Fourier transformation,one can show that the hopping terms and potential terms of the AA model can transfer into each other.By equating the hopping terms with the potential terms,one can see that a selfdual symmetry exists in the AA model.It is easy to see that this self-dual point will separate between the localized and extended states.Unfortunately,one still cannot find any mobility edges in the standard AA model.

In order to make mobility edges possible in the AA-type models, Das Sarma and co-workers[7-9]introduce the slowly varying quasi-periodic potentials into the AA model.This new class of 1D model was shown to support the mobility edges and has been extensively studied ever since.The success of slowly varying quasi-periodic model stimulated a lot of more works on the mobility edges in 1D quasi-periodic systems.[10-13]It is also found that the mobility edge can appear in many different types of modified AA model such as in Refs.[14-19].During the same time, the mobility edges of many different models have been studied.[20-23]A variety of novel physical phenomena induced by the mobility edge and methods to determine the mobility edges are also explored extensively.[24-30]In addition, the research of mobility edge has been further extended to the domain of non-Hermitian systems.[31-37]

In this paper, we present another version of modified AA model, which can also accommodate mobility edges.In this modified AA model,we simply replace the irrational frequencyαwith another much smaller irrational frequency such asα/MwhereM ～102.Due to this much smaller irrational frequency, our model has a quite long quasi-period, which is roughly a few percent of the system size.Here we emphasize that the system size must be kept finite.We will call this model the“AA model with large quasi-periodic disorders”.We will present a detailed study of mobility edges of this large quasiperiodic model with disorders in on-site potentials or in hopping coefficients or in both.We will see that, in each case,their mobility edges are very similar to the AA model with slow-varying quasi-periodic disorders.

We can determine the location of these mobility edges by the so-called energy-matching method that we have proposed in our previous works.[38]In essence,we can approximate the above quasi-periodic disordered model by an ensemble of different periodic models.Then the region of extended states can be approximately obtained by the overlaps of energy bands of all these periodic models.We will demonstrate this method works quite efficiently in determining the mobility edges for different versions of AA models with large quasi-periodic disorders.In the end of this paper, we will present qualitative arguments in support of the existence of mobility edges in the models with large quasi-periodic disorders.

The rest of this paper is organized as follows.In Section 2,we introduce the two types of quasi-periodic disordered models,which will be studied closely.Then we present a brief review of the energy matching method in Section 3.In Section 4, we will apply this method to determine the mobility edges of various disordered models and also use detailed numerical calculations to verify them.In Section 5,we provide a qualitative explanation of the existence of mobility edges in large quasi-periodic disordered models.At last, we briefly conclude in Section 6.

2.The large quasi-period disordered model and slowly varying disordered model

In this paper, we will mainly focus on the modified AA models with very small irrational frequencies or large quasiperiodic disorders.The Hamiltonian of this type of models can be summarized as

Herewandμare the amplitude of the disorder andilabels the lattice site.αandβare certain order 1 irrational numbers,which determine the quasi-periodic behavior.In particular, if one setα=0, then the model only contains disorders in potential terms.Similarly, one can also put the disorders in the hopping term or in both hopping and potential terms.Hereφis some fixed phase angle,which can tune the relative disorder between the potential terms and hopping terms.Most importantly,we also introduce a large integerMin the denominators,which is used to control the order of magnitude of irrational frequencies.It is easy to see that the length of quasi-period ofwiandμiis given by

which is quite long when compared to the system size.In this sense,we call it the long quasi-period model.

Throughout the whole paper, we will often compare the models of large quasi-periods disorders with the models of slowly varying quasi-periodic disorders.As we discussed in the introduction, the model with slowly varying quasi-period is well known to support the mobility edges.Its Hamiltonian can be expressed as

One usually assumes the exponentvto satisfy 0

Due to the characteristic thatw′iandμ′iapproach constants in the largeilimit, one can give a heuristic argument to determine the mobility edge of the slowly varying quasiperiodic model in certain cases.In this method, the asymptotic constancy ofw′iandμ′iis considered to be an important condition for the generation of mobility edges.The large quasi-periodic model does not have this property that disorders gradually approaching constants, but it still produces mobility edges similar to the slowly varying quasi-periodic models.This shows that the asymptotic constancy ofw′iandμ′iis not a necessary condition for the existence of mobility edges.

There is a simple physical picture that can roughly explain why the mobility edges can exist in the slowly varying disordered models.For smalli,the model of Eq.(4)will favor the localized states due to the disorders.On the other hand,for largei, the model of Eq.(4) becomes uniform and favors the extended state.It is the competition of these two effects that gives rise to the mobility edges in the slowly varying disordered models.In Section 5, we will provide some simple arguments to explain the appearance of mobility edges in the large quasi-periodic disordered models.

3.Introducing the energy-matching method

The“energy-matching method”is an effective method to determine the mobility edge.The basic idea of this method is to approximate the quasi-periodic disordered models by an ensemble of periodic models.For convenience,we label these periodic models bya=1,...,N.Then their Hamiltonian can be written as

Here periodic boundary conditions such ascN+1=c1are imposed andwaandμado not depend on the lattice site.For each ?Ha, we can diagonalize the Hamiltonian in the momentum space to find the extended eigenstates with the following energy bands:

Although we cannot diagonalize the quasi-periodic disordered model analytically, one can expect intuitively that the eigenstate of a given energyEis extended, if this energyEis located inside in the energy band of all the periodic model ?Ha.Therefore,the energy region of the extended state of the model of Eq.(1)will be the overlap of all energy bands of Eq.(7)

In a previous paper,[38]we have verified that the energymatching method can give the correct location of mobility edges in this case.

In the next section,we will apply this method to calculate the mobility edges of the models with large quasi-period and the models with slowly varying disorders in several different examples.We will verify that the numerical results of mobility edges are consistent with the result of energy-matching method.In the meantime,the above two types of models will give rise to the same mobility edges.

4.Numerical verification of mobility edges

In this section, we choose the inverse participation ratio(IPR)and Lyapunov exponent,[39-41]which are common physical quantities of disordered systems, as the indicators to distinguish between extended states and localized states.These indicators can help us to numerically verify the accuracy of the mobility edges, which are obtained by the energy-matching method in the previous section.

The inverse participation ratio (IPR) of then-th normalized eigen-wave functionψnis defined asIt is well-known that the Lyapunov exponent is the inverse of the localization length.Therefore, one expects the non-zero Lyapunov exponentγ/=0 for localized states.On the other hand,the Lyapunov exponent will be almost zero for extended states.

4.1.The AA model with disordered potential

As the first example, we consider the case withα=0 andβ=(√5-1)/2, which means the disorders are only in the potential.This is also the original version of the AA model.In this case,wais just a constantw.According to energy-matching method, the lower bound of the region of extended states is determined by the maximum value of the bottom of each energy band of ?Ha, which can be denoted asEbottom=μa-2(t+w).Similarly, the upper bound of extended states is determined by the minimum value of the top of each energy band of ?Ha,which isEtop=μa+2(t+w).It is easy to see that the resulting energy region of extended states is

Asμis tuned,the above upper and lower bounds generate two lines of mobility edges separating extended states from localized states.

To verify the above mobility edges,we plot the eigenenergy of the AA model with large quasi-period disordered potential of Eq.(1)as a function ofμin Fig.1(a).The parameters used in this calculation are listed in the figure caption.The color of each point represents the IPR value of the eigenstates.One can see that the boundary of the dark region matches the mobility edges that we have decided above by the energymatching method.Similarly, we plot the eigenenergy of the slowly varying AA model Eq.(4)with disordered potential in Fig.1(b).The parameters ofα,βin the slowly varying AA model are the same as those in the AA model with the large quasi-period.We also assumew′=wand the slowly varying exponent asv=0.5.It can be seen from the figure that the mobility edges of the two types of models are the same, and they are both consistent with the calculated results of Eq.(12).

In order to make more quantitative observations, in Figs.1(c)-1(f),we also plot the IPR values and Lyapunov exponent as a function of eigenenrgy for a few selected points asμ=0.9,1.6,2.5.While figures 1(c)and 1(d)show the results of the AA model with a large quasi-period, and figures 1(e)and 1(f) show the results of slowly-varying AA model.One can see that the IPR drops from the order of magnitude of 10-2to 10-4when the eigenenergy across certain critical values,and the Lyapunov exponent will also change from non-zero to zero.This clearly signals the transition from localized states to extended states.These figures also support the conclusion that the mobility edges of the two types of models are both consistent with formula(12).

Fig.1.(a)Eigenenergy of the model of Eq.(1)as a function of μ with α =0,β =(√5-1)/2,w=0.5 and M=200.(b)Eigenenergy of the model of Eq.(4)as a function of μ′ with w′ =w and v=0.5, other parameters are the same as(a).The colors in panels(a)and(b)represent the IPR of each eigenstate.The red and green curves represent the mobility edges.(c) and (e) The IPR and Lyapunov exponent of Eq.(1) as a function of eigenenergy for (w,μ)=(0.5,0),(0.5,1.6),(0.5,3.2), other parameters are the same as (a).(d) and (f) The IPR and Lyapunov exponent of Eq.(4)with w′=w and v=0.5,other parameters are the same as(c).

4.2.The AA model with the same disorders in hopping and potential

Now we turn to a more complicated model withβ=αandφ=0.In this case,μiandwihave the same disorder.For this model, we can also use the energy-matching method to determine the mobility edge.The lower bound of the energy region of Eq.(8) is determined by the extreme valuesE=μ-2(t+w)andE=-μ-2(t-w).Depending on the value ofμ, we will take the larger one as the lower bound.Whenμ>2w, clearly, the lower bound of the energy region of the extended state isEmin=μ-2(t+w).One the other hand,ifμ<2w,we haveEmin=-μ-2(t-w)as the lower bound.

Similarly,we can also get extreme valuesE=μ+2(t+w)andE=-μ+2(t-w)for the upper bound of the energy region of the extended state.Obviously,regardless of the value ofμandw, the energyE=-μ+2(t-w) is smaller.So the upper bound of the energy region of the extended states isEmax=-μ+2(t-w).In summary,whenβ=α,we find the condition for extended states is given by

To verify the above mobility edges,we show the eigenenergy of the AA model with large quasi-periodic disorders in both the hopping and potential in Fig.2(a).The parameters used are listed in the figure caption.Again,the color of each point represents the IPR value of the eigenstates.Similarly,we also plot the eigenenergy of the slowly varying AA model with the same disorders in Fig.2(b).The values of the parameters of the slowly varying AA model are the same as those in the AA model with the large quasi-period,but with the varying exponentv=0.5.It can be seen from these figures that the mobility edges of the two types of models are both consistent with the calculated results of Eq.(13).In Figs.2(c)and 2(d), the IPR values as a function of eigenenergy are shown for a few points withμ=0.9,1.6,2.5.Whereas panel(c)shows the results of AA model with large quasi-period and panel(d)shows the results of slowly-varying AA model.It provides a more accurate verification of the mobility edges of Eq.(13).

Fig.2.(a)Eigenenergy of the AA model Eq.(1)as a function of μ with α =β =(√5-1)/2, w=0.5 and M=200.(b)Eigenenergy of the slowly varying AA model Eq.(4)as a function ofμ′ with w′=w and v=0.5,other parameters are the same as(a).The colors in panels(a)and(b)represent the IPR of each eigenstate.The red and green curves represent the mobility edges.(c)The IPR of the AA model as a function of eigenenergy for(w,μ)=(0.5,0.9),(0.5,1.6),(0.5,2.5),other parameters are the same as(a).(d)The IPR of the slowly varying AA model with w′=w and v=0.5,where other parameters are the same as(c).

4.3.The AA model with different disorders in hopping and potential

As the last example, we consider the case with different disorders in hopping and potential.This can be easily achieved by settingφ/=0.In this case,wjandμjare not the same disorder, but they are also not completely independent.To determine the mobility edges, we still consider an ensemble of periodic models labeled byawithwa=wcos(2παa/M)andμa=μcos(2πβa/M-π/2)=μsin(2πβa/M).Here we set the phase difference to beφ=-π/2.Since in the periodic modelHa,waandμaare only constants, we can simplify the Eq.(8)to the following

hereR±ais given by

To verify the above mobility edges,we plot the eigenenergy of the above model with color indicating its IPR value in Fig.3(a).The parameters are given in the figure caption.As a comparison, we also display the eigenenergy of the slowly varying AA model of Eq.(4)with the same type of disorders in Fig.3(b).It can be seen from these figures that the mobility edges of the two types of models are consistent with the calculated results of Eq.(16).In Figs.3(c)and 3(d),we show the IPR values as a function of eigenenergy for a few selective points ofμ=0,1,2.The panels(c)and(d)correspond to the results of the AA model with large quasi-period and slowlyvarying AA model,respectively.It provides more quantitative evidence that supports the mobility edges of Eq.(16).

Fig.3.(a)Eigenenergy of the AA model Eq.(1)as a function ofμ with α=β =(√5-1)/2,w=0.5,φ =-π/2 and M=200.(b)Eigenenergy of slowly varying AA model Eq.(4)as a function ofμ′ with w′=w and v=0.5,other parameters are the same as(a).The colors in panels(a)and(b)represent the IPR of each eigenstate.The red and green curves represent the mobility edges.(c)The IPR of the AA model as a function of eigenenergy for (w,μ)=(0.5,0),(0.5,1),(0.5,2), and other parameters are the same as (a).(d) The IPR of the slowly varying AA model with w′=w and v=0.5,and other parameters are the same as(c).

5.Qualitative explanations of the emergence of mobility edge with large quasi-period

In the original AA model,there is a duality point located atμ=2t.It has been proved rigorously[42]that the eigenstates are all localized forμ> 2tand are all extended forμ<2t.Because of this, no mobility edge exists for a fixedμin the original AA model.Our modification of AA model is to make the irrational frequency very small asα →α/M.The appearance of mobility edges in the large quasi-period models seems to contradict the above-mentioned theorem.We should emphasize that the observed mobility edges only show up for finite-sized systems.It is the combined effects of the large quasi-period and finite system size that make the mobility edge appear.

Here we want to show the effects ofMon the model of Eq.(1).In Fig.4,we plot the eigenenergy of the AA model of Eq.(1)with colors indicating its IPR value for several differentMvalues.The parameters are given in the figure caption.In panel(a),we haveM=1 which is exactly the original AA model.One can clearly see that the duality point located atμ=2 separates the extended and localized states.Therefore,for a fixedμ, there are no mobility edges.From panels (b)to (c), we gradually increaseMfrom 60 to 120.Then for a fixedμ,there appear more and more localized states.But the location of the mobility edge is still not very clear.When we reachM=180,the mobility edge becomes quite sharp and its location is the same as the mobility edge of the slowly varying disordered the AA model.These figures clearly show the relationship between the values ofMand the appearance of mobility edges.

Fig.4.Eigenenergy of the AA model Eq.(1) as a function of μ with α =0, β =(√5-1)/2, w=0.The colors represent the IPR of each eigenstate.From(a)to(d),M=1,60,120,180,respectively.

To see the effects ofMmore clearly, we propose another quantity that can reflect the degree of localization of each eigenstate.In the energy-matching method,we approximate the large quasi-period model by a lot of periodic models.For a given eigenenergyE, we can find the number of periodic models whose energy bands containE.The eigenstates become more localized if this number becomes smaller.In Fig.5, we plot the number of periodic models as a function of eigenerngyEforμ=1.5 andμ=2.5.For the left and right panels,M=50 and 200, respectively.One can see that forμ=1.5,the number of periodic models rapidly decreases from the band center to the band edges.For largerM, this number decreases even faster.This indicates that the eigenstates near the band edges are essentially localized forμ<2.

Fig.5.Thenumber ofperiodic modelsasa func√tionof eigenenergyE forμ=1.5(blue dots)and 2.5(red dots).From left to right,M=50,200,respectively.Other parameters are α =0,β =(5-1)/2,M=50.

Fig.6.The IPR value as a function of eigenenergy E forμ=1.5(red dots)and 2.5(green dots).From left to right,N=1000,10000,respectively.Other parameters are α =0,β =(5-1)/2,M=50.

The mobility edges of the AA model for a fixedμare also observed in Ref.[43].In this paper,it is argued that the original AA model can be thought as a 1D projection of the 2D Hofstadter model with incommensurate magnetic flux density.From this point of view, the localization of the eigenstate of the AA model near the band edge can be understood by semiclassical calculations of the 2D Hofstadter model with small magnetic flux density.We refer readers to Ref.[43]for more detailed explanations.These results are consistent with what we observed.In our paper, we also consider the generalized AA models and find similar results.

To see the effects of system sizeN,we plot the IPR value as a function of eigenenergyEforN=103and 104in Fig.6.Here the IPR is calculated forμ=1.5 and 2.5.The firstμis smaller than the criticalμc=2 and the other is greater than it.One can see that if the system size is smaller, the IPR of states with these twoμbecome closer.One can also see that for fixedμ=1.5, the IPR also decreases a lot as one moves from the band edges to the band center.This shows there is some criticalEforμ<2, which separates the extended and localized states.This feature is more observable for a finitesized system.Because in this case,the IPR difference betweenμ<2 andμ>2 is not that large.

6.Conclusion

In this paper, we have shown that the AA-type model with the large quasi-period can also support mobility edges.By inserting the large quasi-periodic disorders in the hopping coefficients, the potentials or both, we demonstrate that the mobility edges in these examples are almost the same as the corresponding slowly varying disordered models.The precise locations of the mobility edges are deduced by the energy matching method.We argued that the appearance of mobility edges for a fixedμin the large quasi-periodic disordered models is due to the combined effects of long quasi-period and finite system size.The longer quasi-period makes the IPR difference between the band edges and band center bigger.On the other hand, the smaller system size makes the IPR difference across the dual point smaller.Combining the above two effects,the mobility edges for a fixedμbecome more apparent for a model with a large quasi-period and finite system size.

Acknowledgments

Project supported by the National Natural Science Foundation of China(Grant No.11874272)and Science Specialty Program of Sichuan University(Grant No.2020SCUNL210).