the critical behavior of time - independent systems , especially in terms of the spectral statistics and the eigenstate statistics , has attracted great attention . on the spectrum side ,
hofstadter s butterfly spectrum of the harper model has been a paradigm for critical spectral statistics , representing a multifractal spectrum @xcite of a system exactly on the metal - insulator transition point . on the eigenstate side , mainly through studies in time - independent models , such as the power - law random banded matrix ( prbm ) model @xcite and the standard anderson tight - binding model ( tbm ) @xcite , it has been well - established that for a system on a metal - insulator transition point or the anderson transition point @xcite , its eigenstates show clear fractal features .
this background of understanding the critical behavior of time - independent systems motivated our interest in the critical behavior of periodically driven systems .
below we first introduce recent related studies of critical floquet spectra , and then briefly describe the motivation and the results of this work .
it is well known that the floquet ( quasi - energy ) spectrum of a delta - kicked version of the harper model also displays hofstadter s butterfly patterns @xcite .
interestingly , though the kicked harper model ( khm ) can be classically chaotic , its spectrum , due to its fractal nature , does not follow the bohigas - giannoni - schmit conjecture @xcite at all .
this makes the khm not only a fruitful model for gaining new insights into the issue of quantum - classical correspondence in classically chaotic systems , but also an intriguing model to study critical spectral statistics .
indeed , for quite a long time , studies of fractal floquet spectra were largely restricted to the khm and its variants @xcite . in a proposal to experimentally realize hofstadter s butterfly floquet spectrum in cold - atom laboratories , wang and gong @xcite recently demonstrated that hofstadter s butterfly floquet spectrum can be synthesized by use of a double - kicked cold - atom rotor system @xcite under a quantum resonance condition .
@xcite then showed that the butterfly floquet spectrum of the cold - atom system studied in refs .
@xcite is equivalent to that of the standard khm if and only if one system parameter takes irrational values .
in addition to motivating a cold - atom realization of critical floquet spectra of periodically driven systems , refs . @xcite seem to have offered a general strategy for synthesizing critical floquet spectra in driven systems .    using an approach extended from refs .
@xcite , recently wang and gong @xcite showed that the floquet spectra of a class of driven su(2 ) systems also display butterfly patterns and multifractal properties that are characteristics of highly critical spectra .
this establishes a completely different class of critical driven systems without a connection with the khm context .
interestingly , the driven su(2 ) model in ref .
@xcite can be understood as a simple extension of the well - known kicked top model ( ktm ) @xcite in the quantum chaos literature .
because the ktm has just been experimentally realized in a cold @xmath0cs system @xcite , it can be expected that a critical driven su(2 ) system may also be experimentally realized using the collective spin of a @xmath0cs atomic ensemble .
an alternative experimental realization may be based on a driven two - mode bose - einstein condensate @xcite , which represents a strongly self - interacting driven system .
given the above - mentioned class of driven quantum systems with critical floquet spectra , it becomes necessary to study the behavior of the associated floquet eigenstates . theoretically speaking ,
because driven su(2 ) systems always have a finite number of floquet eigenstates , the eigenstate analysis becomes much easier than in the khm , with the latter necessarily involving an infinite number of eigenstates for a fractal floquet spectrum .
a careful investigation of the floquet eigenstates over the entire spectrum will help to better understand the critical behavior in time - dependent systems in general .
experimentally speaking , information about the eigenstate statistics may be more directly accessible to measurements than a fractal spectrum .    to analyze the critical behavior of the floquet eigenstates in driven su(2 ) systems
, we adopt the same approach as in previous studies of time - independent systems .
that is , we shall numerically examine the fluctuations of the eigenstates @xcite .
the eigenstate fluctuations can be characterized by a set of inverse participation ratios ( ipr ) : @xmath1 where @xmath2 is the index of the eigenstates , @xmath3 represents one eigenstate under investigation , and @xmath4 are the basis states . for convenience
we focus on the ipr @xmath5 ( i.e. , @xmath6 ) . by analogy to critical eigenstate behavior in time - independent systems , we expect that @xmath5 scales anomalously with the hilbert space dimension @xmath7 as @xmath8 where @xmath9 is a fractal dimension of a particular eigenstate @xmath10 . but
is there also a unique fractal dimension @xmath11 for the average behavior of all the floquet eigenstates , for example , via the slope of the averaged @xmath12 , denoted @xmath13 , versus @xmath14 ? to that end , we shall examine if , as the system gets closer to the thermodynamic limit ( @xmath15 ) , the distribution of @xmath16 shows signs of a scale - invariant form @xcite . in other words , whether the distribution function of @xmath12 , denoted @xmath17 $ ] , only shifts as @xmath7 varies .
certainly , the system under our study has only a finite size @xmath7 . in time - independent anderson - transition studies using the prbm or the tbm ,
it was conjectured that the variance of @xmath12 , denoted @xmath18 , scales with @xmath7 as @xmath19 with @xmath20 , @xmath21 , and @xmath22 being three adjustable parameters @xcite . for a @xmath23-dimensional system on the anderson transition point
, it was shown that @xmath22 is related to @xmath11 by @xmath24 where @xmath25 equals @xmath26 or @xmath27 depending upon whether or not the system has time - reversal symmetry @xcite .
as one main task of this work , we shall examine if these results for time - independent systems still hold for critical floquet eigenstates .
furthermore , we hope to see how the criticality of the eigenstates of unitary operators differs from the criticality of the eigenstates of self - adjoint operators .
results along this direction will also be relevant to recent investigations on the  unitary anderson model " @xcite , the thue - morse sequence generating multifractal eigenstates of the quantum baker s map @xcite , the one - parameter model of quantum maps showing multifractal eigenstates @xcite , as well as recent experimental and theoretical studies of anderson transition in kicked - rotor systems @xcite .
we now briefly summarize the main findings of this work .
for the driven su(2 ) systems studied here , we consider two different parameter regimes : in one regime the floquet spectra display clear butterfly patterns , and in the other regime , the butterfly patterns of the floquet spectra have dissolved due to increased strength of the driving fields . for both regimes , we find that @xmath28 $ ] is not as smooth as observed in the tbm or prbm , indicating some non - universal features in dynamical systems .
the @xmath28 $ ] for cases with dissolved butterfly patterns is however smoother .
for either regime , it is found that the ensemble average @xmath29 does scale linearly with @xmath14 , with the slope of the @xmath30 vs @xmath14 curve clearly defining the fractal dimension @xmath11 for all the eigenstates .
we also find it possible to fit the variance of @xmath31 by eq . , but with the exponent @xmath22 given by @xmath32 instead ( with @xmath33 ) ; i.e. , a factor of _ two _ is missing from the denominator as compared with eq .
( 4 ) for time - independent critical systems . to further understand this difference
, we propose a random matrix model , which we call  power - law random banded unitary matrix " ( prbum ) model . by tuning the parameters of the prbum
, the @xmath11 value associated with the prbum can be varied .
more interestingly , we observe that the variance @xmath18 of the prbum also follows eq .
( [ var_goe ] ) , with the exponent @xmath22 again given by eq .
( [ dq_coe ] ) .
this suggests that our findings about the floquet eigenstate statistics based on driven su(2 ) systems do reflect some general aspects of critical floquet eigenstates .
this paper is organized as follows . in sec .
[ sec3 ] , we present detailed results of the eigenstates statistics in our driven su(2 ) models , with or without time - reversal symmetry . in sec .
[ sec4 ] , we introduce the prbum to represent a class of critical floquet operators , discuss the statistics of the eigenstates of prbum , and then compare the associated results with those found in actual dynamical systems . in sec .
[ sec5 ] , we study the eigenstate statistics of the standard kicked top model @xcite , which represents a classically chaotic , but non - critical , driven system .
section [ sec6 ] concludes this work .
the focus in ref . @xcite is on the fractal spectral statistics . here , using the same model we study the statistics of the floquet eigenstates .
the first floquet operator under study is given by @xmath34 where @xmath35 are angular momentum operators satisfying the su(2 ) algebra and @xmath36 is the conserved total angular momentum quantum number that defines a @xmath37-dimensional hilbert space .
readers can refer to ref .
@xcite for detailed descriptions and motivations of this model .
this model is also called as a  double - kicked top model " ( dktm ) in ref . @xcite .
eigenstates of the @xmath38 operator are denoted as @xmath39 , with @xmath40 .
states @xmath41 will be chosen as our representation for eigenstate analysis . to analyze the floquet eigenstates ,
it is necessary to express the floquet operator in symmetric basis states , a procedure that block - diagonalizes the floquet matrix . on the one hand
, this will simplify our analysis ; on the other hand , this is necessary for the sake of comparison between an actual dynamical system and the prbum model proposed below .    the dktm floquet operator @xmath42 in eq .
( [ floquet_coe ] ) has a parity unitary symmetry @xmath43 where @xmath44 .
this symmetry can be used to block diagonalize the @xmath42 matrix into two disconnected sub - matrices associated with either odd - parity or even - parity subspaces . without loss of generality we only present below results for the @xmath36-dimensional odd - parity subspace . besides the parity symmetry
, @xmath42 also has a time - reversal anti - unitary symmetry @xmath45 , with @xmath46 where @xmath47 is the complex conjugation operator . to explore the implication of this time - reversal symmetry for the eigenstate statistics
, we shall also consider a variant of @xmath42 , i.e. , @xmath48 evidently , @xmath49 differs from @xmath42 only in the last factor , i.e. , @xmath50 in @xmath42 is replaced by @xmath51 . because of this difference , we call @xmath42 in eq .
( [ floquet_coe ] ) the @xmath52 model and call @xmath49 the @xmath53 model .
it is easy to check that the latter does not have the parity symmetry or the time - reversal symmetry . for the @xmath53 model , which can not be reduced to any block - diagonal form
, we examine the eigenstates of the full floquet matrix .    for both cases
we define a dimensionless system parameter , @xmath54 this choice of @xmath55 being @xmath56 times the golden mean is to ensure that the resulting floquet eigenstate statistics is indeed representative of driven systems with fractal floquet spectra . as detailed below , we consider _ two _ different regimes for the product @xmath57 . in the first regime defined by @xmath58 ,
the floquet spectra show clear butterfly patterns ; in the second regime defined by @xmath59 , the butterfly spectra have dissolved , with fractal dimensions of the spectra increased @xcite .
this is a time - reversal symmetric system .
because dyson s circular ensemble of random unitary matrices @xcite with time - reversal symmetry is called  circular - orthogonal - ensemble " ( coe ) , we regard the @xmath52 model as an example of critical coe statistics .
figure [ jx - jxdistalphaj1.0](a ) shows the distributions of the _ logarithm _ of the ipr @xmath5 , denoted @xmath28 $ ] , for different @xmath36 .
it is seen that the distribution function @xmath28 $ ] is not as smooth as that observed in early anderson - transition studies @xcite .
nevertheless , it is clear that as @xmath36 increases , the left tail of @xmath28 $ ] systematically shifts to the left direction associated with more negative @xmath12 .
the profile of @xmath28 $ ] , though somewhat changes as @xmath36 increases , does maintain its main features as @xmath36 increases .
due to these features that are similar to early findings for the critical eigenstates in time - independent systems , it can be expected that the average of @xmath12 will show a scaling behavior with @xmath60 . as shown in fig .
[ jx - jxdistalphaj1.0](b ) , this is indeed the case .
therein , @xmath30 , obtained by averaging @xmath12 over all eigenstates ( in the odd - parity subspace ) , displays an excellent linear behavior with @xmath61 . from the slope of the fitting line in fig .
1(b ) , we are able to obtain the fractal dimension @xmath62 .    the distribution profile @xmath28 $ ] in fig .
1(a ) is seen to display rich features , with significant fluctuations and multiple notable peaks .
qualitatively , this reflects that our system is an actual dynamical system and hence the underlying rich dynamics will manifest itself through some non - universal statistical features .
related to this observation we also note that in our calculations , all the floquet eigenstates are treated equally and all of them are used for averaging .
this is in contrast to the common procedure in analyzing time - independent critical systems , where only those energy eigenstates in a certain small energy window around zero eigenvalue are included to examine the distribution of @xmath12 @xcite .
the justification for including all floquet states in our analysis is as follows : the quasi - energy spectra lie on a unit circle and hence all states with different eigenphases on the unit circle should be treated on equal footing . to double
check this understanding , we have also taken windows of different widths centered around _
zero _ value of the eigenphase and then calculate the distribution of @xmath31 .
no improvement in the smoothness of @xmath63 is found .
rather , we obtained similar distribution of @xmath12 with clear fluctuations .
it is also tempting to connect the non - universal features of @xmath63 with the phase space structures of the underlying classical limit .
however , such a perspective , which calls for a good understanding of quantum - classical correspondence in critical systems , is unlikely to succeed because the classical limit of our dynamical model is completely chaotic @xcite .    in fig .
[ jx - jxdistalphaj1.0](c ) , we plot @xmath64 $ ] vs @xmath60 ( filled circles ) , where @xmath65 is the variance of @xmath16 and @xmath20 is a fitting parameter , whose value is found by fitting our data points with the empirical formula given in eq .
( [ var_goe ] ) . as seen in fig .
[ jx - jxdistalphaj1.0](c ) , the fitting is reasonably good , yielding that @xmath66 $ ] scales as @xmath67 , with @xmath68 [ @xmath11 is already determined by the fitting in fig .
[ jx - jxdistalphaj1.0](b ) ] , @xmath69 , and @xmath70 .
despite obvious fluctuations around the fitting curve , the result in fig .
[ jx - jxdistalphaj1.0](c ) suggests that the tool borrowed from traditional anderson - transition studies for time - independent systems can be still useful here . furthermore ( probably more interestingly ) , the fitting in fig .
[ jx - jxdistalphaj1.0](c ) also unexpectedly reveals a big difference from what can be expected from eq .
( [ dq_goe ] ) with @xmath33 and @xmath71 : here @xmath68 instead of @xmath72 .
therefore , an intriguing difference between time - independent critical systems and periodically driven critical systems is observed here .     for the @xmath52 model , with @xmath73 and @xmath74 in the representation of odd - parity basis states defined in the text .
the size of the floquet matrix ensemble is important for numerical simulation . in order to construct the necessary ensemble
, we consider a range of @xmath75 , i.e. , @xmath76 , yielding respectively @xmath77 and @xmath78 matrices for the different values @xmath36 .
( b ) @xmath79 , the mean value of @xmath12 averaged over all floquet eigenstates , as a function of @xmath60 .
the slope of the fitting line gives @xmath80 .
( c ) logarithm of @xmath66 $ ] as a function of @xmath60 , where @xmath36 is the dimension of the odd - parity hilbert subspace
. filled circles are our numerical results for the @xmath52 model , and the solid line is the fitting of the numerical results using the empirical formula given in eq .
( [ var_goe ] ) with @xmath81 , @xmath82 and @xmath83 .
the plotted variables here and in all other figures are dimensionless.,width=283,height=264 ]      as mentioned above , for this parameter regime the butterfly patterns in the floquet spectra have dissolved almost completely .
we present the associated eigenstate statistics in fig .
[ jx - jxdistalphaj10.0 ] . in fig .
[ jx - jxdistalphaj10.0](a ) , we show the distribution profile of @xmath12 for different @xmath36 .
in contrast to the previous case shown in fig .
1(a ) , @xmath84 $ ] is now much smoother ( essentially only one peak is left ) . from the same panel , we also see a systematic left - shift of the distribution function as @xmath36 increases .
this systematic left - shift leads to an evident linear behavior of the average value of @xmath12 as a function of @xmath61 , as shown in fig .
[ jx - jxdistalphaj10.0](b ) .
the slope of the fitting line in fig .
[ jx - jxdistalphaj10.0](b ) gives the fractal dimension @xmath85 . comparing this result with that in fig .
1(b ) , one sees that though the fractal dimension of the floquet spectra increases due to increasing @xmath57 @xcite , the fractal dimension of the associated eigenstates may decrease .    in fig .
[ jx - jxdistalphaj10.0](c ) , we examine the variance of @xmath12 as a function of @xmath60 ( again , for the odd - parity subspace ) .
same as in fig .
1(c ) , we fit our results with the empirical formula given in eq .
( [ var_goe ] ) .
the fitting in fig .
[ jx - jxdistalphaj10.0](c ) is better than that in fig .
1(c ) , consistent with the fact that the distribution of @xmath12 is quite smooth here .
the fitting in fig .
[ jx - jxdistalphaj10.0](c ) gives @xmath86 , and @xmath83 , where the value of @xmath11 is found in fig .
[ jx - jxdistalphaj10.0](b ) .
the finding that @xmath22 is not equal to @xmath72 but @xmath11 again strengthens our early observation from fig
. 1 .     for the @xmath52 model with @xmath59
is presented for the odd - parity subspace .
other parameters are the same as in fig .
[ jx - jxdistalphaj1.0 ] .
( b ) @xmath87 is plotted as a function of @xmath60 .
the slope of the fitting line gives @xmath88 .
( c ) logarithm of @xmath66 $ ] as a function of @xmath60 , where @xmath36 is the dimension of the odd - parity hilbert subspace .
filled circles are our numerical results for the @xmath52 model , and the solid line is the fitting of the numerical results using the empirical formula given in eq .
( [ var_goe ] ) , with @xmath89 , @xmath90 and @xmath83.,width=283,height=264 ]       for the @xmath53 model , with @xmath91 and @xmath92 .
the numbers in the brackets are the size of the floquet matrix ensemble .
the dimension of the hilbert space is given by @xmath93 .
( b ) @xmath79 , the ensemble mean value of @xmath12 averaged over all floquet eigenstates , as a function of @xmath14 .
the slope of the curve of @xmath79 vs @xmath14 gives @xmath94 .
( c ) logarithm of @xmath95 $ ] vs @xmath96 .
filled circles are numerical results for the @xmath53 model , and the solid line is the fitting of the numerical results using the empirical formula given in eq .
( [ var_goe ] ) , with @xmath97 , @xmath98 , and @xmath99.,width=302,height=283 ]    to verify if our findings above are general , we now turn to the @xmath53 model [ eq . ( [ floquet_cue ] ) ] . due to the lack of time - reversal symmetry here , this case can be regarded as an example of critical  circular - unitary - ensemble " ( cue ) statistics .
all the eigenstates of the floquet operator @xmath49 will be considered .      for this regime where the butterfly patterns of the floquet spectra can be clearly seen , fig .
[ jx - jydistalphaj1.0](a ) displays the distribution of @xmath12 for different hilbert space dimension @xmath93 .
analogous to the previous case with time - reversal symmetry , @xmath28 $ ] displays interesting fluctuations . as @xmath7 increases ,
@xmath28 $ ] undergoes changes in its profile , shifts its left tail , but also maintains many features . in fig .
[ jx - jydistalphaj1.0](b ) we obtain again a nice linear scaling behavior of @xmath79 with @xmath14 . from the slope of the linear scaling , we obtain the fractal dimension @xmath100 .
this @xmath11 value is different from that for the @xmath52 model with the same values of @xmath57(note that the spectral statistics for the @xmath53 model also differs from that for the @xmath52 model @xcite ) .
same as in fig .
[ jx - jxdistalphaj1.0](c ) , in fig .
[ jx - jydistalphaj1.0](c ) we study the variance of @xmath12 [ now denoted @xmath18 ] as a function of @xmath96 , using the fitting formula given in eq .
( [ var_goe ] ) .
the fitting , though with clear fluctuations , yields that @xmath95 $ ] scales as @xmath102 , with @xmath103 , @xmath104 , and @xmath105 [ @xmath11 value obtained from fig .
[ jx - jydistalphaj1.0](b ) ] .
remarkably , though eq .
( 4 ) with @xmath33 and @xmath106 ( because of the lack of time - reversal symmetry ) predicts @xmath107 , here we have @xmath105 instead .
the important common feature shared by the @xmath53 model and the @xmath52 model is thus the missing of a factor of 2 in the numerically obtained @xmath22 value as compared with the empirical formula for time - independent critical systems .
this interesting finding also supports the use of eq .
( [ var_goe ] ) as a tool for understanding floquet eigenstate statistics . our numerical observations here will be further strengthened by a random matrix model .       for the @xmath53 model with @xmath59
is presented for different hilbert space dimension @xmath93 .
other parameters are the same as in fig .
[ jx - jydistalphaj1.0 ] .
( b ) @xmath79 is plotted as a function of @xmath96 .
the slope of this linear curve gives @xmath108 .
( c ) logarithm of @xmath95 $ ] vs @xmath96 .
filled circles are our numerical results for the @xmath53 model , and the solid line is the fitting of the numerical results using the empirical formula given in eq .
( [ var_goe ] ) , with @xmath109 , @xmath110 and @xmath99.,width=283,height=264 ]    just like the @xmath52 model , in this regime the butterfly patterns of the floquet spectra have dissolved .
the statistical properties of the floquet eigenstates are shown in fig .
[ jx - jydistalphaj10.0 ] . in fig .
[ jx - jydistalphaj10.0](a ) , the distributions of @xmath12 is seen to be much smoother than those seen in fig .
this is somewhat expected from our early findings in the @xmath52 model .
figure [ jx - jydistalphaj10.0](b ) shows a linear scaling of @xmath13 vs @xmath96 , with its slope giving @xmath108 . in fig .
[ jx - jydistalphaj10.0](c ) , we study the variance of @xmath12 as a function of @xmath96 , as compared with the empirical formula given in eq .
( [ var_goe ] ) : the fitting with the empirical formula is excellent , yielding @xmath111 , and @xmath99 , where the value of @xmath11 is determined in fig .
[ jx - jydistalphaj10.0](b ) . once again
, here we find @xmath99 instead of @xmath112 [ as suggested by eq .
( 4 ) with @xmath106 ] .
in studies of time - independent critical systems , the prbm model at criticality @xcite has proved to be fruitful .
the prbm is an ensemble of random hermitian matrices whose matrix elements @xmath114 are independently distributed gaussian random numbers with mean @xmath115 and the variance satisfying @xmath116^{-1}. \label{varh_ij}\ ] ] the case @xmath117 represents the critical point and @xmath118 is a parameter characterizing the ensemble .
a straightforward interpretation of this model is that it describes a one - dimensional sample with random long - range hopping , with the hopping amplitude decaying as @xmath119 .
motivated by our results above for critical floquet states , we aim to propose a class of random unitary matrices , whose floquet eigenstate statistics can show some general aspects of critical statistics and can be used to shed some light on actual dynamical systems .
our natural starting point for generating such random unitary matrices are the hermitian prbm .
to generate a random unitary matrix from a hermitian matrix in the prbm ensemble , we employ the algorithm by mezzadri , whose original motivation is to generate cue random matrices @xcite from general gaussian random matrices . for the sake of completeness ,
we have presented a description of mezzadri s algorithm in appendix [ app1 ] . for our purpose ,
that is , to generate a critical random unitary matrix , we first set the starting point of mezzadri s algorithm as a prbm ensemble at the critical point ( @xmath120 ) .
we then generate an ensemble of random unitary matrices ( denoted @xmath121 ) of the cue class .
significantly , because of the use of prbm as the input for mezzadri s algorithm , we find that the variance of the matrix elements @xmath122 thus obtained also satisfies a power - law , i.e. , @xmath123^{-1}. \label{var_unitary}\ ] ] here the parameter @xmath124 is a common prefactor of the matrix elements , which can be determined by the unitary condition .
the parameters @xmath125 and @xmath126 in eq .
( [ var_unitary ] ) depend on the parameters @xmath127 and @xmath128 of the prbm used .
as three computational examples , panels ( d)-(f ) of fig .
[ var_prbum ] present the dependence of @xmath129 $ ] upon @xmath130 , for three ensembles of random unitary matrices we generated , with sizes @xmath131 , and @xmath132 . if the scaling of @xmath133 with @xmath134 is indeed a power law as described by eq .
( [ var_unitary ] ) , then one should see a linear dependence of @xmath129 $ ] in @xmath130 .
this is indeed the case in figs .
[ var_prbum](d)-(f ) .
note that the deviations in figs .
[ var_prbum](d)-(f ) from the fitting straight lines at very large values of @xmath130 are due to two trivial reasons .
first , for very large @xmath134 , the value of @xmath133 is vanishingly small and hence @xmath135 $ ] becomes extremely large , thus yielding large fluctuations .
second and more importantly , for a fixed matrix size , if @xmath134 is very large , then the available number of matrix elements become insufficient for good statistics . indeed , as the matrix size increases from @xmath136 to @xmath137 , it is seen from figs .
[ var_prbum](d)-(f ) that the validity window of the linear fitting gradually extends to larger values of @xmath130 .     on @xmath130
is plotted , where @xmath138 represents a matrix element at the @xmath139th row and @xmath140th column . here
the prbm ensemble as the input is set at the critical point @xmath120 with the parameter @xmath141 .
panels ( a)-(c ) are for prbum - coe , with the dimension @xmath142 and @xmath132 , respectively . the fitting function of eq .
( [ var_unitary ] ) ( solid lines ) gives @xmath143 , @xmath144 , and @xmath145 .
panels ( d)-(f ) are for prbum - cue , with the dimension @xmath146 and @xmath132 , respectively . the fitting function of eq .
( [ var_unitary ] ) ( solid lines ) yields @xmath147 , @xmath148 , and @xmath149.,width=283,height=238 ]    the random unitary matrices generated in the above manner , with their matrix elements satisfying the power - law scaling of eq .
( [ var_unitary ] ) , are defined as  power - law random banded unitary matrix " of the cue type ( prbum - cue ) .
as detailed in appendix a , one can then generate prbum of the coe type ( prbum - coe ) via @xmath150 . as shown in panels ( a)-(c ) of fig .
[ var_prbum ] , the variance of the matrix elements of prbum - coe also obeys eq . ( [ var_unitary ] ) , with different values of @xmath125 and @xmath126 .    to check whether the prbum - coe and prbum - cue ensembles show critical statistics
, we analyzed their eigenstates , especially in terms of the distribution and the scaling of @xmath12 .
it is found that as we tune the parameter @xmath128 of the prbm used in the algorithm , the resulting fractal dimensions @xmath11 can be also tuned continuously .
for example , the @xmath11 value of prbum can be made close to that of our driven su(2 ) models .
in particular , at @xmath141 , we obtain @xmath151 for prbum - coe and @xmath152 for prbum - cue , yielding @xmath153 and @xmath154 , respectively .
these two @xmath11 values are quite close to the @xmath11 values of the @xmath52 and @xmath53 models found in fig . 1 and fig .
3 . below we describe these findings in detail .
obtained for prbum - coe , with the matrix dimension @xmath155 and @xmath156 .
the numbers in the brackets give the size of the ensemble .
( b ) same as in fig .
[ jx - jxdistalphaj1.0](b ) and fig .
2(b ) , yielding @xmath157 .
( c ) same as in fig .
[ jx - jxdistalphaj1.0](c ) and fig .
2(c ) , but with @xmath36 replaced by @xmath7 .
the fitting curve gives @xmath158 , @xmath159 , and @xmath68.,width=283,height=264 ]    this random unitary matrix ensemble is intended to model a critical floquet operator with time - reversal symmetry .
the results for prbum - coe generated from prbm with @xmath160 are shown in fig .
[ coe_dist ] . in fig .
[ coe_dist](a ) , we show the distributions of @xmath12 for different values of the matrix dimension @xmath7 ( which is the counterpart of @xmath36 in the @xmath52 model ) , with all the eigenstates of the prbum - coe ensemble considered .
in contrast to the @xmath52 dynamical model with a small @xmath57 [ see fig .
1(a ) ] , @xmath84 $ ] here displays very smooth behavior .
figure [ coe_dist](b ) depicts a nice linear relation between @xmath161 and @xmath14 .
the slope of the straight line in fig .
[ coe_dist](b ) gives the fractal dimension @xmath157 , a value close to that in the @xmath52 model with @xmath162 .
as in fig .
[ jx - jxdistalphaj1.0](c ) , fig .
[ coe_dist](c ) shows the fitting of the variance of @xmath12 with @xmath7 , using eq .
( [ var_goe ] ) .
interestingly , the values of the fitting parameters are found to be @xmath163 , @xmath164 , both are similar to those determined in fig . 1(c ) .
more interestingly , this fitting shows that @xmath95 $ ] scales as @xmath102 , with @xmath68 .
this supports our finding in fig .
[ jx - jxdistalphaj1.0](c ) and fig .
we have also studied other cases of prbum - coe using other prbm as the input of mezzadri s algorithm .
for example , we find that if the parameter @xmath128 is set at @xmath165 , then the @xmath11 of the prbum - coe ensemble is around 0.24 , which is close to the @xmath11 value previously found in the @xmath52 model with @xmath166 .
these results clearly support our use of prbum - coe to illuminate the critical eigenstate statistics in the @xmath52 model .
this ensemble aims to model a critical floquet operator without time - reversal symmetry .
all eigenstates of an ensemble of prbum - cue matrices are used for our statistical analysis . for @xmath160 , fig .
[ cue_dist](a ) displays @xmath28 $ ] versus @xmath12 , showing again a smooth dependence .
figure [ cue_dist](b ) shows the corresponding @xmath79 versus @xmath167 , which yields the fractal dimension @xmath168 . in fig .
[ cue_dist](c ) , we fit the dependence of @xmath169 $ ] in @xmath96 , yielding @xmath95 \sim n^{-\gamma}$ ] , with @xmath105 [ instead of @xmath170 predicted by eq .
this also confirms our early observations in the @xmath53 model .
the values of the fitting parameters are found to be @xmath171 and @xmath172 , which are close to what we found in fig .
[ jx - jydistalphaj1.0](c ) .
we have also checked that if we perform analogous calculations for @xmath173 , then the @xmath11 value for the prbum - cue ensemble will be close to that found in fig .
4(b ) . given these results ,
we are led to the conclusion that prbum as proposed above do share some general aspects with periodically driven systems having critical eigenstate statistics .
obtained for prbum - cue , with the matrix dimension @xmath174 and @xmath175 .
the numbers in the brackets give the size of the ensemble .
( b ) same as in fig .
[ jx - jydistalphaj1.0](b ) and fig .
4(b ) , yielding @xmath168 .
( c ) same as in fig .
[ jx - jydistalphaj1.0](c ) and fig .
4(c ) , the fitting gives @xmath176 , @xmath172 , and @xmath105.,width=283,height=264 ]
for the standard classically chaotic kicked top model , for @xmath177 and @xmath92 .
the numbers in the brackets are the size of the floquet matrix ensemble . in constructing the ensembles we have considered a range of @xmath75 , i.e. , @xmath178 .
( b ) distributions of @xmath12 for the standard dyson s coe matrices , with the same matrix dimension as in the kicked top model and the same ensemble size .
( c ) analogous to fig .
[ jx - jxdistalphaj1.0](b ) and fig .
7(b ) , the scaling behavior of @xmath79 vs @xmath61 is shown .
open circles are numerical results for the kicked - top model , crosses are numerical results associated with coe random matrices , and the solid curve represents the theoretical prediction from the random matrix theory .
the scaling shows that @xmath179 in the standard kicked top model , which is dramatically different from our observations made from the double - kicked top model.,width=302,height=264 ]    finally , as a numerical  control " experiment , we study the floquet eigenstate statistics of the standard kicked top model .
this will help appreciate the difference between a normal driven system and a critical driven system , both of which can have a chaotic classical limit
. consider then the following floquet operator for the standard kicked top model @xcite , @xmath180 which is just the last two factors of eq .
( [ floquet_coe ] ) , with the same parity symmetry and time - reversal symmetry as the @xmath52 model .
in addition , we set the parameter @xmath181 at the same value as given in eq .
( [ parameter1 ] ) .
we construct a statistical ensemble by considering a range of @xmath75 , i.e. @xmath178 ( with chaotic classical limits ) .
we carry out the floquet eigenstate statistics in the odd - parity subspace , whose dimension is @xmath36 .
because the classical limit is found to be chaotic , we compare the statistics with that associated with dyson s coe matrices in random matrix theory ( rmt ) .
figure [ haake_top](a ) and ( b ) compare @xmath28 $ ] associated with @xmath182 with that obtained from coe matrices , for different @xmath36 . the difference between the actual dynamical system and the coe can hardly be seen .
figure [ haake_top](c ) depicts @xmath183 as a function of @xmath61 , with the results of the kicked top ( open circles ) almost on top of those of coe matrices ( crosses ) . the solid line in fig .
[ haake_top](c ) represents the theoretical curve for @xmath184 obtained from rmt , i.e. , @xmath185 .
the agreement between numerical coe results , analytical rmt result , and the kicked top system as a classically chaotic dynamical system is almost perfect . from the curve shown in fig .
[ haake_top](c ) , it is clear that @xmath11 here is unity and as such the system does not show critical behavior .
this non - critical behavior indicates that the floquet states of the kicked top model are essentially random states , a feature fundamentally different from our double - kicked top system that has a butterfly spectrum and critical statistics in the floquet eigenstates .
it is also interesting to note that in fig .
[ haake_top](a ) and ( b ) , as @xmath36 increases , @xmath28 $ ] becomes narrower and develops higher peaks .
this is an indication that , unlike the critical cases studied above , @xmath63 for the standard kicked top model approaches a dirac - delta type singular function with zero width ( i.e. @xmath186 ) as @xmath36 increases .
in this numerical study we have examined the statistics of the floquet eigenstates of a recently proposed double - kicked top model with multifractal floquet spectra . following the methodologies used in studies of anderson transition in time - independent systems
, we have shown that the floquet eigenstates associated with multifractal floquet spectra also display critical behavior .
in particular , we focus on the distribution of @xmath12 and examine how the quantity @xmath13 averaged over all states scales with the dimension of the hilbert space @xmath7 .
it is shown that @xmath187 scales linearly with @xmath14 , with the slope of this linear scaling giving the fractal dimension @xmath11 of the floquet eigenstates .
the values of @xmath11 are found to be far from unity ( as a comparison , we showed that similar analysis for a standard kicked top with a chaotic classical limit yields @xmath179 ) , constituting strong evidence that the floquet eigenstates are fractal and hence lying between localized and delocalized states .
though we have worked on @xmath5 only , we note that similar analysis can be done for @xmath188 defined in eq .
( [ pq ] ) .
one may then define a generalized fractal dimension @xmath189 and further establish the multifractal nature of the floquet eigenstates .
the variance of @xmath12 , denoted @xmath18 for a hilbert space of dimension @xmath7 , is also examined . in anderson - transition studies with prbm , @xmath18 is known to scale as @xmath102 with @xmath190 for one - dimensional systems , where @xmath191 for a system with ( without ) time - reversal symmetry . by contrast , in our critical driven system , @xmath18 is seen to scale similarly , but with @xmath192 .
this reflects an interesting difference between time - dependent systems and time - independent systems .
indeed , eigenstates of prbm are to model those of critical hermitian operators , whereas floquet eigenstates of a critical driven system should be understood in terms of critical unitary operators . to justify this understanding ,
we have introduced a random unitary matrix ensemble called prbum , with the variance of the matrix elements of the unitary matrices following a power - law distribution .
we show that the eigenstates of prbum share many critical statistical features with the double - kicked top model .
most important , the variance of @xmath12 of prbum does scale as @xmath193 , which is the same as in the double - kicked top model as a critical driven system .
we hence anticipate that this scaling property of the variance of @xmath31 may be general in critical driven systems .
these results complement the spectral results in ref .
@xcite and should motivate further mathematical and theoretical studies in critical driven systems .
j.w . acknowledges support from national natural science foundation of china ( grant no.10975115 ) , and j.g . is supported by the nus start - up fund ( grant no .
r-144 - 050 - 193 - 101/133 ) and the nus  yia " ( grant no .
r-144 - 000 - 195 - 101 ) , both from the national university of singapore .
this is a simple and numerically stable algorithm to generate the cue matrices from an ensemble of complex random matrices @xmath194 , whose elements are gaussian distributed random numbers with mean zero and variance _ unity_. in particular , applying the gram - schmmidt ortho - normalization method to the columns of an arbitrary complex matrix @xmath195 , one can factorize @xmath195 as : @xmath196 where @xmath197 is a unitary matrix and @xmath198 is an invertible upper - triangular matrix . one can easily prove that the above factorization is not unique . because of this non - uniqueness , the random unitary matrices @xmath199 are not distributed with haar measure @xcite , i.e. , the @xmath199 matrices are not uniformly distributed over the space of random unitary matrices .
fortunately , this factorization can still be made unique by imposing a constraint on the @xmath198 matrices . by some group
theoretical arguments , it was shown @xcite that if one finds a factorization such that the elements of main diagonal of @xmath198 become real and strictly positive , then @xmath199 matrices would be distributed with haar measure and hence form cue . following these results ,
the major steps of mezzadri s algorithm are the following .
first , we start with an @xmath200 complex gaussian random matrix @xmath195 .
second , we factorize @xmath195 by any standard @xmath201decomposition routine such that @xmath202 .
third , we create a diagonal matrix @xmath203 where @xmath204 are the diagonal elements of @xmath198 . as a final step ,
we define @xmath205 and @xmath206 . by construction ,
the diagonal elements of @xmath207 are always real and strictly positive , and as such @xmath208 would be distributed with haar measure and can be used to form the desired cue .
the symmetric coe matrices can be constructed from the cue matrices in a very simple manner . in particular ,
let @xmath121 be a member of the cue generated above , then it can be shown that @xmath150 will be a member of coe .
for the generation of prbum advocated in this work , we propose to replace @xmath195 in the first step by a member in the prbm ensemble that models anderson transition .
though there is no mathematical theory for our procedure , the uniformly distributed eigenphases ( not shown here ) of our prbum ensemble thus generated suggest its uniform distribution .
see also some related papers : r. lima and d. shepelyansky , phys .
* 67 * , 1377 ( 1991 ) ; r. artuso _ et .
al _ , _ ibid _ * 69 * , 3302 ( 1992 ) ; r. ketzmerick , k. kruse , and t. geisel , _ ibid _ * 80 * , 137 ( 1998 ) ; i. i. satija , phys .
e * 66 * , 015202 ( 2002 ) ; j. b. gong and p. brumer , phys . rev . lett .
97 , 240602 ( 2006 ) .
q. zhang , p. hnggi , and j. b. gong , , 053607 ( 2008 ) ; new journal of physics , * 10 * , 073008 ( 2008 ) .
q. xie and w. hai , eur .
j. d * 33 * , 265 ( 2005 ) ; m.p .
strzys , e.m .
graefe , and h.j .
korsch , new j. phys .
* 10 * , 013204 ( 2008 ) ; j. b. gong , l. molina - morales , and p. hnggi , , 133002 ( 2009 ) .