Criticality in the Quantum Kicked Rotor with a Smooth Potential

# Criticality in the Quantum Kicked Rotor with a Smooth Potential

Rina Dutta and Pragya Shukla Department of Physics, Indian Institute of Technology, Kharagpur, India.
July 30, 2019
###### Abstract

We investigate the possibility of an Anderson type transition in the quantum kicked rotor with a smooth potential due to dynamical localization of the wavefunctions. Our results show the typical characteristics of a critical behavior i.e multifractal eigenfunctions and a scale-invariant level-statistics at a critical kicking strength which classically corresponds to a mixed regime. This indicates the existence of a localization to delocalization transition in the quantum kicked rotor.

Our study also reveals the possibility of other type of transitions in the quantum kicked rotor, with a kicking strength well within strongly chaotic regime. These transitions, driven by the breaking of exact symmetries e.g. time-reversal and parity, are similar to weak-localization transitions in disordered metals.

###### pacs:
PACS numbers: 05.45+b, 03.65 sq, 05.40+j

## I Introduction

.

The analogy of the statistical fluctuations of dynamical systems and disordered systems is well-known in the delocalized wave-regime (corresponding to the metallic limit in disordered systems and the classically chaotic limit in dynamical systems) and has been explained using random matrix theory as a tool casati1 (); izrp (); fyd (); thaha (). A similar analogy exists for fully localized regimes of the wavefunctions too (i.e. between the insulator limit of disordered systems and the integrable limit of dynamical systems) casati1 (); izrp (); fyd (). It is therefore natural to probe the presence/ absence of the analogy in partially localized/ critical regimes of these systems. Our analysis shows that, similar to the Anderson Hamiltonian ( as dimension), the quantum kicked rotor (QKR) undergoes a localization-delocalization transition in the classically mixed regime. We also find quantum phase transitions in its chaotic regime due to breaking of the symmetries e.g. time-reversal and parity in the quantum system. Similar to disordered systems, the symmetry-breaking transitions in the QKR occur due to weak-localization effects. Similar phase transitions due to symmetry-breakings have been seen in a few other complex systems too e.g. the ensembles of distinguishable spins pizo ().

The connection of the kicked rotor to the Anderson Hamiltonian has been known for several decades fgp (); izrp (); casati2 (); chir (). A recent work garcia () further explores the connection and shows that, for the non-analytic potentials in the QKR, the eigenstates show multifractality or power-law localization mj (); caste (); cu1 (); emm (), a behavior similar to the eigenstates of a dimensional Anderson Hamiltonian leer () at its critical point. Our study shows existence of the multifractal eigenstates in the QKR with smooth potentials too e.g. at specific parametric conditions. Furthermore, similar to a critical Anderson system, the multifractality in the QKR is accompanied by a critical level statistics (size-independent and different from the two ends of the transition), a necessary criteria for the critical behavior mj (). This indicates a much deeper connectivity of the kicked rotor to the Anderson Hamiltonian, not affected just by the nature of the potential or the dimension of system. As discussed here, the connection seems to be mainly governed by the ”degree of complexity” (measured by the complexity parameter discussed later) and may exist among a wider range of dynamical and disordered systems.

The present study is motivated by a recent analytical work psco (); psal () leading to a common mathematical formulation for the statistical fluctuations of a wide range of complex systems. The work, based on the ensemble-averaging, shows that the fluctuations are governed by a single parameter besides global-constraints on the system psco (); psal (). referred as the complexity parameter, turns out to be a function of the average accuracy of the matrix elements, measured in units of the mean-level spacing. The fluctuations in two different systems, subjected to similar global-constraints, are analogous if their complexity parameters are equal irrespective of other system-details.

The -formulation was recently used by us to find the Gaussian Brownian ensemble (GBE) me () analog and the power-law random banded matrix ensemble (PRBME) mir () analog of the Anderson system (for arbitrary ) psand (). However it can not directly be applied to find the QKR analog; this is, in principle, due to inapplicability of the ensemble averaging to dynamical systems. Fortunately it is possible to derive for dynamical systems by a semi-classical route, using the phase-space averages psnl (). The ”semi-classical” was used by us psnl (); psef () to map the statistics of the time-evolution operator of the QKR to the circular Brownian ensembles (CBE) me (). The lack of a suitable criteria for the critical statistical behavior prevented us earlier from a critical QKR-analysis. Our present work pursues the analysis by first analytically identifying the critical QKR-behavior using the ”semi-classical” ; the limit , with independent of the size , gives the critical points of transition. This is followed by a numerical analysis of an ensemble of the QKR which confirms critical behavior at the semi-classically predicted values. This in turn suggests a paradoxical validity of the ensemble averaging in dynamical systems at least of QKR type; (also indicated by the CBE-QKR mapping). A subsequent comparison of the ”semi-classical” to the ”ensemble-based” of the Anderson Hamiltonian (arbitrary ) gives us its QKR analog; the analogy is numerically confirmed too.

The statistical behavior in the bulk of the spectrum of a standard Gaussian ensemble is known to be analogous to a standard circular ensemble (for large matrix sizes) me (); apps (). Our work extends this analogy to their Brownian ensemble counterparts too, that is, between the GBE and the CBE (by the mapping GBE AE QKR CBE). This in turn indicates a connection among a wide-range of physical systems which are known to be well-modeled by the Anderson Hamiltonian, the kicked rotor and the Brownian ensembles of Gaussian and circular type.

The paper is organized as follows: The section II briefly reviews the basic features of the kicked rotor and the Anderson Hamiltonian required for our analysis; it also discusses the parametric conditions in the kicked rotor which can support critical points. The section III deals with the numerical confirmation of the critical level-statistics and the multifractality of the eigenfunctions at critical parametric conditions in the QKR and their comparison with a dimensional Anderson Hamiltonian. We conclude in section IV with a brief discussion of our main results and open questions.

## Ii Kicked Rotor and Anderson Hamiltonian

The kicked rotor and the Anderson Hamiltonian have been subjects of intense study in past and many of their details can be found in several references casati1 (); izrp (); casati2 (); been (); thaha (). However, for self consistency of the paper, we present here a few details required for later discussion.

### ii.1 Kicked Rotor and Complexity Parameter

The kicked rotor can be described as a pendulum subjected to periodic kicks (of time-period ) with Hamiltonian given as

 H=(p+γ)22+Kcos(θ+θ0)∞∑n=−∞δ(t−nT). (1)

Here is the stochasticity parameter, and are the time-reversal and parity symmetry breaking parameters in the quantum Hamiltonian while acting in a finite Hilbert space.

Integration of the equations of motion , between subsequent kicks e.g. and gives the classical map,

 pn+1 = pn+Ksin(θn+θ0)(mod2π), θn+1 = θn+pn+1(mod2π) (2)

The map is area-preserving and invariant under the discrete translation . It also preserves the time-reversal symmetry and the parity for all values of and . Thus the classical dynamics depends only on , changing from integrable () to near integrable () to large scale chaos ().

The quantum dynamics can be described by a discrete time-evolution operator izrp (); chir () where

 B=exp[−ikcos(θ+θ0)] (3) G=exp[−iT(p+γ)2/2ℏ]=exp[iτ(∂∂θ−iγℏ)2/2] (4)

Here , , and are the position and momentum operators respectively; has discrete eigenvalues, () due to periodicity of , (). The choice of a rational value for results in a periodicity also for momentum operator , () and therefore in discrete eigenvalues for , (). The quantum dynamics can then be confined to a two dimensional torus (with a Hilbert space of finite size with an integer). The classical analog of this model corresponds to the standard mapping on a torus of size in the momentum ; thus the classical limit is with and chir (); izrp ().

For the dynamical-localization analysis, it is useful to express the matrix in the momentum-basis izrp (); izrp (); psnl ():

 Umn=1Nexp[iτ(m−γ/ℏ)2/4+iτ(n−γ/ℏ)2/4]N1∑l=−N1exp[−ikcos(2πl/N+θ0)]exp[−2πil(m−n)/N] (5)

where with if is odd and if is even. It is clear that the properties of , eq.(1), are recovered in the infinite matrix size limit.

The quantum dynamics under exact symmetry conditions (, ) can significantly be affected by relative values of . It was first conjectured izrp () and later on verified psnl () that the statistical properties are governed by the ratio of the localization length to the total number of states or, equivalently, (as with as the diffusion constant). However, to best of our knowledge, the critical behavior at was not probed before. The other parameters playing a crucial role in the quantum dynamics are and , the measures of time-reversal and parity-symmetry breaking respectively (with and ). Note the change of or is a canonical transformation, thus leaving the classical Hamiltonian unaffected. The corresponding quantum dynamics, however, is affected as the quantum Hamiltonian acting in a finite Hilbert space may not remain invariant under a unitary transformation.

Following eq.(5), the multi-parametric nature of is expected to manifest itself in the statistical behavior of its eigenvalues (quasi-energies) and the eigenfunctions. However, as shown in psnl () using semi-classical techniques, the quasi-energy statistics of is sensitive to a single parameter and the exact symmetry-conditions. Under exact time-reversal symmetry () and partially violated parity () (taking , equivalently, , without loss of generality), we have psnl ()

 Λkr,t=θ20Nk24π2=N3θ20K264π4M2 (6)

Note, for , is essentially the same as the one conjectured in izrp () for scaling behavior of the spectral statistics. Similarly, for the strongly chaotic case () with only parity symmetry ( and no time-reversal (),

 Λkr,nt=γ2qN348π2 (7)

with . As eq.(6) indicates, in the strongly chaotic limit ; similarly, from eq.( 7), for ; the statistics in these cases can be well-modeled iz1 (); psnl () by the standard random matrix ensembles of unitary matrices, known as the standard circular ensembles e.g. circular orthogonal ensemble (COE), circular unitary ensemble (CUE) etcme (). The cases with a slow variation of , or (partial localization, partial time-reversal or parity-violation respectively) and finite size correspond to a smooth variation of or between and . The intermediate statistics for these cases psnl (); psef () can be well-described by the circular Brownian ensembles. The latter are the ensemble of unitary matrices, described as and characterized by a single parameter () and the exact system-symmetries me (); apps (). The perturbation belongs to a standard Gaussian ensemble of the Hermitian matrices e.g Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE) etc me (). The circular Brownian ensemble analog of a QKR is given by the condition

 Λkr=Λcbe. (8)

As mentioned above, the studies psnl (); psef () did not explore the infinite size limit of and its application for the critical behavior analysis; we discuss it in next section.

### ii.2 Critical behavior of Quantum Kicked Rotor

For a critical point analysis of the QKR, we search for the system conditions leading to the critical eigenvalue-statistics and the multifractal eigenfunctions in the infinite size limit (). The parametric conditions for the critical level statistics, characterized by a non-zero, finite in the limit , can be obtained from eqs.(6, 7). The analysis suggests the possibility of several continuous families of critical points, characterized by the complexity parameter and the exact symmetries; we mention here only the three main cases:

(i) : Eq.(6) in this case leads to a size-independent with . For small- (in the mixed regime), this corresponds to a localization delocalization transition under the time-reversal conditions (no parity symmetry) with a continuous family of critical points characterized by . The bulk-statistics here is analogous to a circular Brownian ensemble (see eq.(8)) with (due to a GOE type perturbation with of a Poisson matrix me (), see eq.(8)). The two ends of the transition in this case are the Poisson () and the COE ensemble ().

(ii) : Eq.(7) in this case gives with . For large (in strongly chaotic limit), a finite gives the critical parameter for the transition from a time-reversible to a time-irreversible phase (both phases delocalized); it can be referred as the weak-localization critical point. The end-points of the transition are the COE () and the CUE (), with the critical statistics given by an intermediate circular Brownian ensemble with (due to a GUE type perturbation , with , of a COE matrix me () ).

(iii) : Here eq.(6) gives with describing the critical point family for the transition from a parity-symmetric phase to a parity fully violated phase (both phases time-reversal). For in the mixed regime (), a variation of leads to the Poisson COE transition. For in the strongly chaotic regime (), the transition end-points are the 2-COE () and the COE (). The critical statistics for a specific is given by the intermediate Brownian ensemble with (due to a GOE type perturbation , with , of the Poisson ensemble if , or, the 2-COE ensemble if ).

As eqs.(6, 7) indicate, a size-independent can be obtained by other combinations of too. This suggests critical behavior in the symmetry-spaces other than those mentioned above.

The critical nature of the system for specific parametric conditions can further be confirmed by an analysis of the eigenfunction fluctuations. The studies of a wide range of systems (see mj (); caste (); cu1 (); emm () and the references therein) reveal the presence of strong fluctuations in the eigenfunctions near a critical point. The fluctuations can be characterized through the set of generalized fractal dimension or , related to the scaling of the moment of the wavefunction intensity with size caste (): . The multifractality of the eigenfunctions can also be analyzed through the spectrum of singularity strengths mj (), related to by a Legendre transformation: (see mj (); emm () for details). In section III, we numerically analyze both and to detect the multifractality of the QKR-eigenfunctions.

### ii.3 Anderson Hamiltonian and the Complexity Parameter

The Anderson model for a disordered system is described by a -dimensional disordered lattice, of size , with a Hamiltonian in the tight-binding approximation leer (). In the site representation, turns out to be a sparse matrix of size with the diagonal elements as the site energies and the off-diagonals given by the hopping conditions. For a Gaussian type on-site disorder (of variance and zero mean) and a nearest neighbor (n.n.) isotropic hopping with both random (Gaussian) and/ or non-random components, can be modeled by an ensemble (later referred as the Anderson ensemble or AE) with following density

 ρ(H,v,b)=Cexp⎡⎣−∑kH2kk/2ω−∑(k,l)=n.n.H2kl/2η⎤⎦∏(k,l)=n.nδ(Hkl−t)∏(k,l)≠n.n.δ(Hkl) (9)

with as the normalization constant. As discussed in psand (), the above ensemble can be rewritten as

 ρ(H,h,b)=Cexp[−∑k≤l(1/2hkl)(Hkl−bkl)2] (10)

where , with for pairs connected by the hopping, for all pairs representing the disconnected sites. The single parameter governing the spectral statistics (see eq.(19) of psand ()) can then be given as

 Λa,e(E,Y) = (|αw−αi|F2γ0)ζ2dL−d≈|αw−αi|γ0N(FItyp2)2, (11)

with

 αw=ln|1−γ0ω|+(z/2)ln[|1−2γ0η||t+δt0|] (12)

with as the number of the nearest-neighbors, and as an arbitrary constant. Further is the mean level density, as the localization length and as the typical inverse participation ratio: .

## Iii numerical Analysis

The objectives of our numerical analysis are two fold: (i) a search for the critical points of the quantum kicked rotor, and, (ii) a comparison of its fluctuation measures with those of a 3-dimensional Anderson ensemble. For this purpose, we analyze the following cases:

(i) QKR1: quantum dynamics time-reversal but parity broken: which gives . This case corresponds to the critical set (i) in section II. B and is analyzed for many sizes () to verify the critical behavior.

(ii) QKR2: quantum dynamics with both time-reversal and parity broken: . This case belongs to the set (ii) in section II.B and its critical nature is also confirmed by analyzing many values.

(iii) QKR3: : same as QKR1 but with which gives . We consider this case to verify the analogy with a Anderson ensemble.

(iv) QKR4: quantum dynamics time-reversal but parity broken: (case (iii) in section II.B). This is also analogous to the Anderson system mentioned above in the QKR3 case, notwithstanding the crucial changes in and for the two QKR cases.

To explore critical behavior, we analyze large ensembles of the matrices for both QKR1 and QKR2 for various matrix sizes ; the ensemble in each case is obtained by varying in a small neighborhood while keeping fixed. The chosen range give in the mixed regime for QKR1 () and in the chaotic regime for QKR2 (). Prior to the analysis, the quasi-energies (the eigenvalues of ) are unfolded by the local mean level density (, a constant due to repulsion and a unit-circle confinement of the quasi-energies). The figures 1,2 display the nearest-neighbor spacing distribution and the number variance , the measures of the short and long-range spectral correlations respectively, for the QKR1 and the QKR2. Note the curves in figure 1 are intermediate to the Poisson and the COE limits; the size independence implies their survival in the infinite size limit too. This indicates the QKR1 as the critical point of transition from a localized phase to the delocalized phase. Similarly the curves in figure 2, intermediate to the COE and the CUE limits, suggest the QKR2 as the critical point of transition from the time-reversed phase to the time-irreversible phase (both phases in the chaotic regime).

To reconfirm the critical nature of the QKR1 and the QKR2, we numerical analyze the moments of their local eigenfunctions intensity for various sizes. The results shown for in figures 3a,4a indicate the multifractal nature of the eigenfunctions; for small- ranges, shows a behavior (or ) with and for the QKR1 and the QKR2 respectively (see table 1 for the first few values of ). These results are reconfirmed by a numerical study of the spectrum displayed in figures 3b,4b. For this purpose, we use the procedure based on the evaluation of moments, described in emm () (using eq.(4) and eq.(10) of emm ()) which has the advantage of full control over the finite-size corrections. The numerical results for and for the QKR (see table 1) indicate a parabolic form of (also confirmed by the fits shown in figure 3b, 4b) and therefore a log-normal behavior of the local eigenfunction intensity for large -regions (see emm (); mj (); caste () for details). As shown in figures 3c, 4c, the tail-behavior of can be well-approximated by the function

 f(u′)=√a/2πeb+au′e−a(u′+c)2/2). (13)

with (see figures 3c, 4c for numerical values of a,b,c).

The bulk statistical behavior of the Hermitian matrices is known to be analogous to the unitary matrices me (); apps (). This, along with the single parametric formulation of the statistics of the Hermitian matrix ensemble psco (), suggests the analogy of the QKR-Anderson ensemble statistics if their parameters are equal (besides similar global constraints e.g. global symmetries) psal (). Our next step is the comparison of the fluctuation measures of a time-reversal Anderson case with the QKR3, a time-reversal system with partially localized wavefunction in the momentum space (dynamical localization). For this purpose, we analyze a cubic () Anderson lattice of linear size () with a Gaussian site disorder (of variance , and mean zero), same for each site, an isotropic Gaussian hopping (of variance and mean zero) between the nearest-neighbors with hard wall boundary conditions; these condition correspond to the critical point for a disorder driven metal-insulator transition psand ()). A substitution of the above values (with ) in eq.(12) gives . As shown by the numerical analysis in psand (); pswf (), and which on substitution in eq.(11) gives (with ).

For AE-QKR comparison, we analyze the ensembles of 2000 matrices with matrix size for the 3d AE case and for the QKR case. The energy dependence of (see psand (); psnl ()) forces us to confine our analysis to only of the levels near the band-center from each such matrix. The levels are unfolded by respective local mean-level density in each case (so as to compare the level-density fluctuations on a same density-scale) me (). Figure 5a shows the AE-QKR3-QKR4 comparison of ; the good agreement among the three curves verifies the -dependence of the spectral correlations. This is reconfirmed by figure 5b showing a comparison of the number variance. Note that for QKR4 is same as the theoretical analog given by the condition (see eqs.(6,11)). However for QKR3 show a small deviation from the theoretically prediction (). This may be due to eq.(6) being a poor approximation at the integrable-nonintegrable boundary . Note that the classical limits of QKR3 and QKR4 are different ( for QKR3, for QKR4), although their ”semi-classical” are same ( with for QKR3, for any for QKR4). As clear from eq.(2), marks the boundary between the integrable and mixed dynamics; corresponds to the mixed nature of the dynamics (see figure 5c).

As discussed in pswf (), the eigenfunction fluctuations of finite systems are influenced by two parameters, namely, system size as well as . To compare dependence of an eigenfunction measure, therefore, same system size should be taken for each system. Figure 6a shows the distribution of the local eigenfunction intensity, , for the and the QKR3. The close proximity of the two curves suggests as the parameter governing the local eigenfunctions intensity too; (we have verified the analogy also with QKR4). A comparison of the AE-QKR3 multifractality spectrum, shown in figure 6b, reconfirms their close similarity at least on the statistical grounds.

The numerical confirmation of the statistical analogy of QKR4-AE-QKR3 systems supports our claim regarding single parametric ()-dependence of the statistics besides global constraints. Note the latter are same for both QKR3 and QKR4 (i.e parity violated, time-reversal preserved and mixed dynamics) which results in their statistics intermediate to same universality classes, namely, Poisson and COE although the transition parameters are different in the two cases.

The QKR-AE analogy can be utilized to connect them to other complex systems too. In pswf (); psand (), we studied the AE-connection with the PRBME (described by , mir ()) and the GBE (, psand ()). For the AE case considered above, the PRBME and the GBE analogs for the spectral statistics turned out to be and ; these systems are therefore the spectral statistical analogs of QKR3, QKR4 as well as of a circular Brownian ensemble with (see case(i) of section II.B).

## Iv conclusion

To summarize, we studied the statistical analogy of two paradigmatic models of dynamical and disordered systems, namely, the quantum kicked rotor and the Anderson Hamiltonian, in the partially localized regimes. Our results indicate the existence of critical behavior in the classically mixed regime of the QKR with a smooth potential. This is qualitatively analogous to a disorder driven metal-insulator transition in the Anderson system; the quantitative analogy for their statistical behavior follows if their complexity parameters are equal. Our study also reveals the possibility of other transitions in the QKR e.g from a symmetry preserving phase to a symmetry fully violated phase. These transitions are analogous to similar symmetry breaking transitions in disordered metals e.g the Anderson Hamiltonian in the weak disorder limit in the presence of a slowly varying magnetic field.

As with the Anderson transition, the QKR transitions are governed by the complexity parameter too. However, contrary to the -derivation for the Anderson case by an ensemble route, for the QKR is derived by a semi-classical method. The semi-classical -formulation is also numerically verified for the QKR-ensembles. This indicates an equivalence of the ensemble-averaging and the phase-space averaging for the statistical analysis. This further lends credence to the single parametric formulation of the statistical behavior of complex systems, irrespective of the origin of their complexity. However it needs to be examined for other dynamical systems.

Research has indicated a multi-parametric dependence of the spectral-statistics at long energy scales of the dynamical systems berry (), originating in the level-density oscillations due to short periodic orbits. However, these studies are not at variance with our work. This is because the ”semi-classical” -derivation in psnl () is based on the assumed equivalence of the traces of the operators with their phase-space averages. The assumption may not be valid on short time-scales of the dynamics. One should also understand the exact role of the ensemble-averaging for the statistical analysis of dynamical systems.

## V acknowledgment

We would like to express our gratitude to Prof. M.V.Berry for many helpful suggestions.

## References

• (1) B.V. Chirikov, Phys. Rep. 52, 263 (1979); G.Casati and L.Molinary, Prog. Theor. Phys. Suppl. 98, 287, (1989).
• (2) F. Izrailev, Physics Reports, 196, 299, (1990).
• (3) Y.V.Fyodorov and A.D.Mirlin, Int. J. Mod. Phys.B, 8, 3795, (1994).
• (4) M. Thaha and R. Blumel, Phys. Rev. Lett. 72, 72, (1994).
• (5) I. Pizorn, T. Prosen, S. Mossmann and T.H.Seligman, New. J. Phys. 10, 023020, (2008).
• (6) S. Fishman, D. R. Grempel and R. E. Prange, Phys. Rev. Lett., 49, 509, (1982).
• (7) B.V.Chirikov, F.M.Izrailev and D.L.Shepelyansky, Physica D 33, 77, (1988).
• (8) G. Casati, B. V.Chirikov, I. Guarneri and F.Izrailev, Phys. Rev. E 48, R1613 (1993). G. Casati, L.Molinari and F.M.Izrailev, Phys. Rev. Lett. 64, 1851 (1990). K.Zyczkowski, M.KUs, M.Lewenstein and F.Izrailev, Phys. Rev. A 45, 811 (1992).
• (9) A. M. Garcia-Garcia and J. Wang, arXiv:0710.5852v1.
• (10) M. Janssen, Phys. Rep. 295, 1, (1998).
• (11) C.Castellani and L.Peliti, J. Phys. A 19, L429 (1986); F.Wegner, Nucl. Phys. B 280, 210 (1987).
• (12) E. Cuevas, Phys. Rev. B 68, 184206, (2003); Phys. Rev. B 68, 024206, (2003).
• (13) F. Evers, A. Mildenberger and A. D. Mirlin, Phys. Rev. B 64, 241303, (2001).
• (14) P. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985).
• (15) P.Shukla, Phys. Rev. E, 71, 026266, (2005); Phys. Rev.E 75, 051113 (2007); Phys. Rev. E 62, 2098, (2000).
• (16) P.Shukla, J. Phys. A, (2008).
• (17) M.L. Mehta,Random Matrices, (Academic Press, New York, 1991); C.E.Porter, Statistical Theory of Spectra: Fluctuations (Academic Press, New York, 1965); F.Haake, Quantum Signatures of Chaos (Springer, Berlin, 1991).
• (18) A.D.Mirlin, Y.V.Fyodorov, F-M Dittes, J Quezada and T H Seligman, Phys. Rev. E, 54, 3221, (1996).
• (19) P.Shukla, J. Phys.: Condens. Matter 17, 1653, (2005).
• (20) P. Shukla and A.Pandey, Nonliearity ; A. Pandey, R. Ramaswamy and P. Shukla, Pramana-J. Phys. India, 41, L75, (1993).
• (21) P. Shukla, Phys. Rev. E, 53, 1362, (1996).
• (22) A. Pandey and P. Shukla, J. Phys. A: Math. Gen. 24, 3907, (1991).
• (23) J. Tworzydlo, A. Tajic and C.W.J. Beenakker, arXiv: cond-mat/0405122; Phys. Rev. B 69, 165318 (2004).
• (24) F. Izrailev, Phys. Rev. Lett. 56, 541, (1986).
• (25) P. Shukla, Phys. Rev. E, 75, 051113, (2007).
• (26) M.V.Berry, Proc. R. Soc. London Ser. A 400, 229, (1985).

## Appendix A Figure Caption

Fig. 1 [Color online]. Spectral Measures of QKR1:

(a) the nearest-neighbor spacing distribution of the eigenvalues for various sizes , with inset showing behavior on the linear scale, (b) Number variance for various sizes.

The convergence of the curves for different sizes indicates scale-invariance of the statistics. The behavior is critical, being different from the two end-points, namely, Poisson and CUE statistics even in infinite size limit.

Fig. 2 [Color online]. Spectral Measures of QKR2: (a) Distribution of the nearest-neighbor eigenvalue spacings for various sizes, with inset showing behavior on the linear scale, (b) Number variance for various sizes.

Again the statistics being intermediate between COE and CUE, and, convergence of the curves for different sizes indicates its critical behavior.

Fig. 3 [Color online]. Multifractality of QKR1 :

(a) Fractal Dimension along with the fit with (good only for ). A fit for the large regime suggest following behavior: . (b) Multifractal spectrum for various sizes along with the parabolic fit with and . (c) Distribution with of the local intensity of an eigenfunction for QKR1. The solid line represent the function given by eq.(13) with , and (corresponding to an approximate log-normal behavior of ), a good approximation in tail-region as expected. The inset shows the behavior on linear scale.

Fig. 4 [Color on line]. Multifractality of QKR2: (a) Fractal Dimension along with the fit with (good only for ). (b) Mulifractal spectrum for various sizes along with a parabolic fit, of the same form as in figure 3.b with and . (c) Distribution of the local intensity of an eigenfunction for QKR2, with same as in figure 3c. The solid line represent the function with , and (corresponding to a log-normal behavior of ) which fits well in the tail region of . The inset shows the behavior on the linear scale.

Fig. 5 [Color online]. Comparison of spectral statistics of the Anderson ensemble with QKR3 and QKR4: (a) , with inset showing the linear behavior, (b) . Note QKR3 and QKR4 turns out to be close to the one suggested by the relation (giving for QKR3 and for QKR4). (c) Phase-space behavior of the classical kicked rotor at (see eq.(2))

Fig. 6 [Color online]. Comparison of the multifractality of the eigenfunctions of QKR3 with the Anderson ensemble: (a) local intensity distribution of an eigenfunction with same as in figure 3,4. Also shown is the function with , a good fit in the tail region of . The inset shows the behavior on linear scale. (b) multifractal spectrum : Note the QKR3 analog of the Anderson ensemble here is same as the one in figure 5. Again the fit has the same parabolic form as in figures 3,4, with and .

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