Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: this article will be about discussing the possibilities to characterize _ randomness _ and _ dynamical complexity _ in quantum mechanics and relating this issue to the questions of non - equilibrium statistical mechanics . we shall try to illustrate , mainly by presenting various numerical examples , a possible ` cyclist approach ' @xcite towards the _ quantum many - body problem _ which is inspired by experiences gained in studies of quantum and classical chaos of one or few particles ( see e.g. @xcite ) . solving the many body problem in quantum mechanics presents a major challenge from its very beginnings . and along this way , many ingenious important analytical and efficient numerical methods of solution have been developed , for example bethe ansatz @xcite , later generalized and interpreted as the quantum inverse scattering problem @xcite , real space renormalization group methods @xcite , quantum monte carlo @xcite , density matrix renormalization group ( dmrg ) @xcite , and more recently various quantum information - theoretic based time - dependent dmrg ( tdmrg ) @xcite . the ultimate aim of any of these methods is to efficiently find analytical or numerical approximation to the solution for some of the physical observables in the quantum many body problem , however many methods work only under some specific conditions which are not always well understood . for example , bethe ansatz and quantum inverse scattering work only for a small subset of problems which are _ completely integrable _ , and which may often have very non - generic non - equilibrium properties , such as e.g. ballistic transport at high temperatures @xcite . quantum monte carlo techniques represent a very successful set of numerical methods which can yield thermodynamic equilibrium averages for generic ( non - integrable ) systems , however they are practically useless for computing non - equilibrium , or time dependent quantities , such as transport coefficients . traditional dmrg method @xcite is provably very successful for computing accurate ground state expectation values of almost any physical observables of one - dimensional interacting systems . and tdmrg @xcite promised to extend the success of dmrg to time dependent physics . indeed , the first numerical experiments looked very promising , but after a closer look one may realize that they have all been applied to a rather special subset of interacting systems and to a rather special subset of initial states . for generic ( non - integrable ) interacting systems or for sufficiently complex initial states , tdmrg should fail to provide an efficient computation as shall be discussed below . in our view this is to be expected , and represents an intrinsic characteristic of quantum complexity of such system and should correspond to many body extensions of the phenomena of quantum chaos @xcite where spectra and eigenvector coefficients can be described by statistical ensembles of random matrices . in the past two , almost three decades there has been a lot of activity in the so called field of quantum chaos , or quantum chaology @xcite , where people were trying to understand the essential and significant features of quantum systems which behave chaotically in the classical limit . classical chaos can be defined in terms of positive ( algorithmic , kolmogorov - chaitin ) complexity of systems orbits . still , the question whether such definition of complexity can be in an intuitively sensible way translated to quantum mechanics failed to be answered in spite of many efforts . it seems that quantum systems of finite ( one or few body ) chaotic systems are generically more robust against imperfections @xcite than their classical counterparts , which is consistent with a simple illustration based on wave - stability of unitary quantum time evolution @xcite . nevertheless , it has been suggested that exponential sensitivity to initial conditions , the essential characteristics of classically chaotic systems , has many fingerprints in quantum mechanics . for example , in one of the pioneering works on quantum loschmidt echoes ( or fidelity decay ) , jalabert and pastawski @xcite have shown that in the semi - classical regime , decay of system s sensitivity to external perturbation , as defined by fidelity , is exponential with the rate which can be related to classical lyapunov exponents . however , this is only true in sufficiently semi - classical regime , where effective planck constant is smaller than effective strength of perturbation , and where fidelity can be essentially explained classically @xcite . on the other hand , in purely quantum regimes , quantum fidelity decays in completely different manner than the classical fidelity , and quite surprisingly displays slower decay for systems with stronger decay of temporal correlations @xcite .    throughout this paper we shall discuss various ergodic properties of a simple toy model of non - integrable interacting quantum many body system , namely kicked ising ( ki ) chain @xcite , and its time - independent version , which ( both ) undergo a transition from integrable to quantum chaotic regimes when the direction of the ( kicking ) magnetic field is changed . by simulating the dynamics of the model in terms of tdmrg we find that it can be performed efficiently only in the integrable regimes . being interested in statistical mechanics of the model we shall describe numerical experiments addressing the questions on the relationship between the onset of quantum chaos in ki model and its quantum mixing property , namely the nature of relaxation to equilibrium , and the properties of non - equilibrium steady state . on one hand we argue that the regime of quantum chaos essentially corresponds to the regime of quantum ergodicity and quantum mixing where diffusive transport laws are valid . on the other hand , we conjecture that the transition to non - ergodic regime may occur before the system parameters reach the integrable point ( even in thermodynamic limit ) , and that non - ergodic to ergodic transition can be characterized with order parameters which change discontinuously at a critical value of system parameter . we argue that the process of relaxation to equilibrium in quantum chaotic ( mixing ) case can be described in terms of quantum analogues of ruelle resonances @xcite . furthermore , we conjecture that the eigenvectors of this process possess a certain scale invariance which can be described by simple power laws . we also discuss a possibility of numerical calculation of quantum dynamical entropies @xcite in a non - trivial setting of ki model , and find , quite remarkably , that positivity of such dynamical entropies _ does not correspond _ to any other measures of quantum chaos , namely quantum dynamical entropies appear to be positive even in the integrable regimes . about two thirds of the material presented in this paper constitutes a review of a selection of recent results related to quantum chaos in many - body systems , with a flavor of quantum information . however , about one third of the material , constituting a major part of section [ sect : statmech ] , is new and original and has not yet been published before . the paper is organized as follows . in short section [ sect : definitions ] we review basic definitions of complete integrability and quantum chaos , in particular in the context of many - body systems . in section [ sect : model ] we introduce ki model which will serve us as a very convenient and efficient toy model to numerically demonstrate all the phenomena discussed in this paper . in section [ sect : losch ] we discuss the robustness of quantum systems to external perturbations and decoherence in open quantum systems , mainly as characterized by quantum loschmidt echoes and entanglement between the system and the environment . in section [ sect : mps ] we discuss the time efficiency of best ( known ) classical simulation , namely using tdmrg , of locally interacting 1d quantum systems and its possible dependence on integrability of the model . in section [ sect : transport ] we relate standard criteria of quantum chaos to normal ( diffusive ) versus anomalous ( ballistic ) transport and discuss a simulation of a quantum heat current in non - equilibrium steady state . in section [ sect : statmech ] we discuss a problem of quantum relaxation to equilibrium , i.e. the quantum mixing property , and quantitative characterization of quantum dynamical complexity . this section contains a large portion of original intriguing numerical results which support few interesting and perhaps surprising conjectures . in section [ sect : conc ] we conclude and discuss some open problems . the central issue of this paper is to verify and demonstrate to what extend the _ complete integrability _ affects non - equilibrium properties of quantum many - body systems and their dynamical complexities , and conversely , whether _ ( strong ) non - integrability _ generically coincides with the established criteria of _ quantum chaos_.    let us start by giving some established definitions ( see e.g. @xcite ) of the basic terms needed to understand the issue . a classical hamiltonian many - body system with @xmath0 degrees of freedom is _ completely integrable _ , if there exist @xmath0 functionally independent global smooth phase space functions ( integrals of motion ) which are mutually in involution , i.e. all poisson brackets among them vanish . in such a case global canonical transformation to canonical _ action - angle _ variables can be constructed , and dynamics can be explicitly solved - at least in principle . for locally interacting _ infinite _ systems ( @xmath1 ) analytic methods for an explicit construction of integrals of motion and canonical action - angle variables are known which usually go under a common name of _ inverse scattering method_. explicit solution by a mapping to an iso - spectral sch " rodinger problem in terms of inverse scattering technique is usually understood as a definition of _ complete integrability _ in such context .    definition of _ quantum complete integrability _ is less unique . nevertheless , algebraic , non - commutative versions of the inverse scattering technique exist and can be applied to some quantum interacting lattices in one - dimension , generalizing the famous bethe ansatz solution of the heisenberg spin @xmath2 chain , and this is used as the most general known definition of _ quantum complete integrability_. other versions of integrability for finite @xmath0 quantum systems have been proposed but they do not correspond to the integrability of the underlying classical limit , if the latter exists , so these ideas will not be considered in this paper .    it has not been generally accepted yet , though demonstrated in many occasions , that completely integrable quantum systems constitute only a small subset of physical models and posses many exceptional ( non - generic ) non - equilibrium dynamical properties , like for example anomalous transport at finite temperatures ( see e.g.@xcite . )    on the other extreme of ergodic hierarchy we have chaotic systems . in classical hamiltonian dynamics of few particles , chaos is best defined in terms of positive algorithmic ( kolmogorov ) complexity of systems trajectories , or equivalently , by exponential sensitivity to initial conditions . however , bounded quantum systems of finite number of particles can not be dynamically complex as their excitation spectrum is discrete , and hence the evolution is necessarily quasi - periodic ( or almost periodic ) . still , quite surprisingly , even for such systems certain dynamical properties are _ random _ and _ universal _ , if the underlying classical limit is sufficiently strongly chaotic @xcite . but genuine dynamical complexity may emerge in thermodynamic limit . however , there is still no completely satisfactory definition of dynamical chaos for infinite quantum systems . the commonly used working definition is the reference with the _ random matrix theory _ @xcite , namely the many - body quantum system is said to be quantum chaotic if its excitation spectrum or some other dynamical properties can be ( on certain energy , or time scale ) well described by ensembles of random hermitian matrices with appropriate symmetry properties . throughout this paper we shall use , either for illustration of theoretical results , or for numerical experimental studies , the following 1d locally interacting quantum lattice system , namely a chain of @xmath0 qubits , or spin @xmath2 particles , coupled with nearest neighbour ising interaction and periodically kicked with spatially homogeneous , but _ arbitrarily oriented _ magnetic field . in a suitable dimensionless units our model can be written in terms of a three - parameter periodic time - dependent hamiltonian @xmath3 @xmath4 is a unit - periodic dirac delta and @xmath5 are the usual pauli spin variables on a finite lattice @xmath6 , satisfying the commutation ( lie ) algebra @xmath7=\sum_\gamma 2\ii\varepsilon_{\alpha \beta\gamma } \delta_{jk } { \sigma^{\gamma}_{j}}$ ] . sometimes it will turn fruitful if we extend the set of pauli matrices by identity matrix and assign them a numerical superscript @xmath8 , namely @xmath9 . the finite chain will often be treated with _ periodic boundary conditions _ , @xmath10 , and sometimes the thermodynamic limit ( tl ) @xmath1 will be considered , in particular when we shall study the statistical mechanics of ( [ eq : hamiltonian ] ) in section [ sect : statmech ] . although kicked hamiltonian one - particle models have been very popular in the field of nonlinear dynamics and quantum chaos for decades , for example the chirikov s kicked rotator model @xcite , the use of kicked systems in quantum many body physics has been so far very limited . if one is not only interested in zero temperature ( ground state ) or low temperature physics , then as we shall try to demonstrate in this review , kicked many - body models like ( [ eq : hamiltonian ] ) provide simpler and clearer phenomenological picture of global dynamics than time - independent models . the main reason is that since energy is not a conserved quantity , the entire hilbert space of many - body configurations is accessible for non - integrable dynamics , and the notions of quantum ergodicity and mixing ( see e.g. ref.@xcite for definitions and further references ) can be defined more clearly than in the time - independent , autonomous setting . perhaps the first kicked interacting quantum lattice has been introduced in ref.@xcite . even if in traditional solid state physics such kicked dynamics would represent un - physically strong excitations , one has to realize that kicked quantum chains could be attractive options as benchmark models of quantum state manipulation and quantum computation in optical lattices . the _ kicked ising _ ( ki ) model ( [ eq : hamiltonian ] ) has been defined for the first time in ref.@xcite , generalizing the integrable ki model with transverse field introduced in ref.@xcite . clearly , as shown there @xcite for the case of _ transverse _ field , @xmath11 , the time - dependent model ( [ eq : hamiltonian ] ) can be considered _ completely integrable _ since it can be solved explicitly , for example by wigner - jordan - bogoliubov transformation to non - interacting spinless fermions , and a large class of its time correlation functions can be calculated explicitly . in addition , an infinite sequence of local conserved quantities ( integrals of motion ) can be constructed in such a case . there is another , trivial _ completely integrable _ regime of ki model , namely the case of _ longitudinal _ field , @xmath12 yet another , more non - trivial _ completely integrable _ regime of ki model is found when the magnitude of dimensionless field is a multiple of @xmath13 , namely @xmath14 , since then the magnetic kick can be considered as a multiple of @xmath15 rotation , and generated by a slightly obscure set of non - interacting spinless fermions .    however , in a general case of titled magnetic field when both components @xmath16 and @xmath17 are non - vanishing , and @xmath18 is non - integer , the model is non - integrable , and is conjectured to be not amenable to exact analytical methods . as discussed in the following sections , non - integrable ki model can display a variety of regimes according to the criteria of quantum chaos , quantum ergodicity and quantum mixing . in fact , recently the spectral statistics of ki model in strongly non - integrable regime has been studied in detail and random matrix behaviour has been clearly confirmed at short energy ranges @xcite .    due to kicked nature of interaction , the evolution propagator of ki model for one unit period of time ( the so - called floquet operator ) , starting just before the kick , can be constructed explicitly in terms of a time - ordered product @xmath19 the last line suggests a simple efficient quantum protocol to simulate ki model in terms of simple 1-qubit @xmath20 and 2-qubit @xmath21 quantum gates . if we write a compact kicked ising 2-qubit gate @xmath22 and introduce a left - to - right ordered operator product , namely @xmath23 , then we can write ki protocol as a simple string of @xmath24gates @xmath25    the ki model has a better known autonomous limit , namely time - independent ising chain in a tilted magnetic field @xmath26 which is again non - integrable unless the field is transverse @xmath27 , or longitudinal @xmath28 . one of the central questions about the dynamics of complex quantum systems is their robustness against small imperfections in the hamiltonian . while it is clear that due to unitarity the quantum evolution is stable against variation of initial states @xcite , it is not so clear how stable it is against variation of the hamiltonian , either being static or time - dependent - perhaps even noisy . let us write the hamiltonian as @xmath29 , where @xmath30 is the unperturbed hamiltonian , @xmath31 is a hermitian perturbation operator and @xmath32 is a small perturbation parameter . peres @xcite proposed the following measure of stability of quantum evolution : let us start from some fixed initial state @xmath33 , and write the time evolutions of this state generated with the unperturbed and perturbed hamiltonian , as @xmath34 and @xmath35 , respectively . then the stability is characterized by fidelity , _ i.e. _ the squared overlap between these two states @xmath36    there are two alternative interpretations of quantum fidelity : ( i ) one can interpret ( [ eq : fidelity ] ) as a _ quantum loschmidt echo _ , namely the probability that the initial state @xmath33 , which is propagated forward with perturbed evolution @xmath37 , and after that propagated backwards in time with unperturbed evolution @xmath38 , is again measured in the same ( initial ) state . alternatively ( ii ) is just the square modulus of expectation value of the unitary echo operator @xmath39 which is the quantum propagator in the _ interaction picture_.    there have been three main theoretical approaches to understanding of fidelity decay in quantum dynamical systems :      in a seminal work jalabert and pastawski derived quantum loschmidt echo for systems which posses well defined classical limit . they have shown that under certain conditions , namely that the perturbation strength is large enough - typically larger than appropriately scaled planck constant , and that initial states have certain classical interpretations - like coherent states , position states , etc , the quantum loschmidt echo decays exponentially @xmath40 with the rate @xmath41 which is _ perturbation independent _ and typically equals the lyapunov exponent of the underlying classical dynamics . more recently a completely classical interpretation of this so - called lyapunov decay has been given @xcite in terms of _ the classical loschmidt echo _ and a theory for exponents @xmath41 has been developed in terms of the full spectrum of lyapunov exponents for classical dynamical systems with few @xcite or many @xcite degrees of freedom .      in the opposite regime , where the scaled planck constant is bigger than the perturbation strength @xmath32 , one can use time - dependent quantum perturbation theory to second order to derive a simple linear response formula for fidelity decay @xcite @xmath42 in terms of 2-point time correlation function of the perturbation @xmath43 where @xmath44 is the perturbation operator in the interaction picture and @xmath45 is an expectation value in the initial state .    from this formula - which can be viewed as a kind of kubo - like linear response theory for dissipation of quantum information , an interesting conclusion can be drawn : fidelity decays faster for systems with slower decay of temporal correlations , or alternatively phrased , quantum system is more robust against external perturbations if it relaxes to equilibrium faster . one can actually go beyond the second order perturbation theory , and sum up the born series for fidelity to all orders in many specific situations @xcite . since we are here mainly interested in qubit ( spin 1/2 ) chains , we shall only review specific results for kicked ising chain @xcite . we shall discuss three different specific values of system parameters , in all three we take @xmath46 : ( a ) integrable regime of transverse field @xmath27 , ( b ) weakly non - integrable regime with @xmath47 , and ( c ) strongly non - integrable regime with @xmath48 . all the results , correlation functions and fidelity , are for random initial states , which can be interpreted as pure states of maximum quantum information , and averaging over random states is equivalent to a tracial state @xmath49 . we consider the evolution operator ( [ eq : floquet ] ) and perturb it by changing the magnetic field such that the perturbation is generated by the transverse component of the magnetization @xmath50 , namely @xmath51^t$ ] .    in the integrable case , or in general , in non - ergodic case , where the correlation function @xmath52 ( note that @xmath53 ) approaches a non - vanishing plateau value as @xmath54 , namely time averaged fluctuation of magnetization @xmath55 is non - vanishing @xmath56 , one can sum up born series to all orders , yielding a gaussian decay of fidelity @xmath57 the only assumptions are that @xmath58 is long enough for the time average in the definition of @xmath59 to converge , and short enough that fidelity is still above the finite size plateau value @xmath60 . note that @xmath59 can be considered as an analog of drude weight or charge stiffness in the linear response solid state transport theory ( see e.g. @xcite ) . in figs.([fig : corr],[fig : fid ] ) ( a , b ) one can observe the correlation plateaus of correlation functions for integrable and weakly non - integrable cases and respective good agreement with a gaussian decay of fidelity ( [ eq : gauss ] ) . note that in the integrable case the correlation function and the plateau @xmath59 have been calculated also analytically @xcite . note that for comparing different lattice sizes @xmath0 a size - scaled value of the perturbation strength @xmath61 , with @xmath62 , has been fixed rather than @xmath32 itself . one the other hand , for sufficiently strong integrability breaking , say in case ( c ) of ki model , the correlation function @xmath63 decays to zero in tl , which can be interpreted as quantum mixing behviour . it has been found that quantum mixing behaviour typically corresponds to random matrix ( quasi)energy level statistics , see e.g. ref.@xcite and the subsequent sections of the present paper . we shall call such behaviour the regime of quantum chaos . here again , born series for fidelity can be summed up to all orders , yielding an exponential decay @xcite @xmath64 where @xmath65 is a transport coefficient . the assumptions for validity of ( [ eq : ne ] ) are that @xmath58 is larger than a characteristic mixing time scale on which @xmath66 decays and short enough so that finite size effects in quantum correlation function @xmath66 does not yet affect the transport coefficient @xmath67 . this regime of fidelity decay is sometimes referred to as the fermi golden rule regime . again , as demonstrated in figs . ( [ fig : corr],[fig : fid])c the agreement of numerical data for ki model with the theory is excellent . note that decay time of fidelity scales as @xmath68 for ergodic and mixing dynamics , and as @xmath69 for non - ergodic dynamics , making the latter much more sensitive to small perturbations . complex quantum systems can often be well described by statistical ensembles of hamiltonians @xcite . assuming that @xmath30 and @xmath31 both belong to canonical ensembles of gaussian random matrices , one can evaluate ensemble averages of the fidelity and relate @xmath66 to the spectral form factor of the underlying random matrix ensemble . the perturbative ( linear response ) theory has been developed in ref . @xcite , see also ref .  @xcite . furthermore a non - perturbative ( super - symmetric ) averaging has been successfully applied to obtain exact expressions for average fidelity amplitude in the most interesting cases @xcite . rmt theory of fidelity has been successfully applied to chaotic quantum systems and even to several experimental situations @xcite .    following more applied philosophy , groups from toulouse and como have performed a series of numerical experiments @xcite analyzing the robustness of several reasonable models of quantum computer hardware under small imperfections , being either due to ( static ) unwanted inter - qubit coupling or due to stochastic ( noisy ) unwanted coupling to the environment , when performing quantum algorithms simulating various toy models of classical and quantum single - particle chaos @xcite , like for example quantum kicked rotator @xcite . the numerical results are in line with a general linear response theory , stating that static perturbations are in general more dangerous than noisy ones . similar thinking as in the previous subsection can be applied to perhaps even more fundamental problem of quantum physics , to the problem of decoherence @xcite . here we shall limit ourselves to an abstract unitary model of decoherence , where we treat a complete unitary evolution over two subsystems , a central system c , and an environment e , and then address a relevant information about the central system ( i.e. the part of the system which is of physical interest ) by partially tracing over the environment . such a discussion can indeed be followed with a close link to the problem of loschmidt echoes , by writing the total hamiltonian @xmath29 as an ideal ( unperturbed ) separable evolution @xmath70 perturbed by a small coupling @xmath31 between the system and the environment . we shall also assume that we start from initial pure state which is a product state @xmath71 we are interested in the properties of a generally mixed state of the central system obtained by partial tracing over the environment @xmath72\ ] ] then , under an ideal evolution , the state of the system remains a product state at all times , so the state of the central system remains pure and there is no decoherence . decoherence is usually characterized in terms of decaying off - diagonal matrix elements of @xmath73 in a suitable pointer state basis , for example in the eigenbasis of the central system hamiltonian @xmath74 . in fact for certain special forms of the perturbation operator @xmath31 , the magnitudes of off - diagonal matrix elements of @xmath73 can be literally written as fidelities , or quantum loschmidt echoes , in the evolution of the environment perturbed by system - environment coupling @xcite . in such framework , the lyapunov regime of fidelity decay corresponds exactly to the lyapunov growth of decoherence discussed by zurek and paz @xcite .    however , another interesting indicator of decoherence is the growth of _ bi - partite entanglement _ between the system and the environment . perhaps this notion is even more general since it does not depend on a particular choice of the pointer state basis . since the state of the universe ( central system + environment ) is always pure , the characterization of entanglement is easy , say in terms of von neuman entropy of @xmath73 , @xmath75 , or linear entropy @xmath76 which is a negative logarithm of purity @xmath77 . usually , the quantities @xmath78 and @xmath79 essentially give equivalent results - namely they can both be interpreted as the logarithm of an effective rank of the state @xmath73 , but the linear entropy , or purity , is more amenable to analytical calculations .    again considering chaotic models for the central system and the environment , miller and sarkar @xcite ( see also @xcite ) have been able to observe the lyapunov regime of entanglement growth , namely @xmath78 is for sufficiently strong perturbations @xmath32 found to grow linearly @xmath80 with the rate @xmath81 which is perturbation independent and given by ( the sum of positive ) classical lyapunov exponents of the unperturbed , separated ( sub)systems . on the other hand , in the purely quantum regime of small perturbation @xmath32 , again typically smaller than effective planck constant , one can use time - dependent perturbation theory and derive perturbation dependent entanglement entropy @xcite , namely the purity can be explicitly expressed as @xmath82 where @xmath83 is a particular integrated correlation function of the perturbation . in this regime we again find that quantum systems , and quantum environments , which display faster relaxation to equilibrium , are more robust against decoherence due to typical couplings . for studies of bi - partite entanglement , in particular in ki model system , see refs .  @xcite . there exists even closer connection to fidelity theory , namely one can prove a general inequality @xcite stating that purity is always bounded by the square of fidelity @xmath84 in other words : @xmath85 gives an upper bound for the growth of the linear entanglement entropy , @xmath86 .    in a slightly different context , _ bit - wise _ entanglement between a pair of qubits of a quantum register representing a time dependent quantum state , where the rest of the register is considered as an environment , has been demonstrated to be an indicator of quantum chaos @xcite and even a signature of classical chaos @xcite . in the theory of classical dynamical systems there is a fundamental difference between integrable and chaotic systems as outlined in section [ sect : definitions ] . chaotic systems , having positive algorithmic complexity , unlike the integrable ones , can not be simulated for arbitrary times with a finite amount of information about their initial states . computational complexity of individual chaotic trajectories is linear in time , however , if one wants to describe statistical states ( phase space distributions ) or observables of chaotic classical systems , up to time @xmath58 , exponential amount of computational resources is needed , typically @xmath87 where @xmath81 is the kolmogorov s dynamical entropy related to exponential sensitivity to initial conditions . but on the other hand , how difficult is it to simulate isolated and bounded quantum systems of many interacting particles using classical resources ? in analogy with the classical ( chaotic ) case , we might expect that the best classical simulation of typical quantum systems ( in tl ) is exponentially hard , i.e. the amount of computing resources is expected to grow exponentially with time .    even though there is no exponential sensitivity to initial conditions in quantum mechanics , there is a tensor - product structure of the many - body quantum state space which makes its dimension to scale exponentially with the number of particles , as opposed to linear scaling in the classical case , and due to presence of entanglement generic quantum time evolution can not be reduced to ( efficient ) classical computation in terms of non - entangled ( classical like ) states . here we propose the idea to use _ the computation complexity of best possible classical simulation of quantum dynamics , as a measure of quantum algorithmic complexity_. this section essentially reviews the article @xcite . we note that our proposal is different from existing proposals of quantum algorithmic complexity @xcite , namely we consider merely the classical complexity of ( best ) classical simulations of quantum states . in the sense of mora and briegel @xcite , quantum algorithmic complexity _ per unit time _ of initially simple time - evolving quantum states propagated by locally interacting hamiltonian @xmath88 is clearly vanishing , since an approximate quantum circuit which reproduces the state after time @xmath58 is a simple repetition of trotter - suzuki decomposition of @xmath89 . recently , a family of numerical methods for the simulation of interacting many - body systems has been developed @xcite which is usually referred to as time - dependent density - matrix - renormalization group ( tdmrg ) , and which has been shown to provide an efficient classical simulation of some interacting quantum systems . of course , it can not be proven that tdmrg provides the best classical simulation of quantum systems , but it seems that it is by now the best method available . simulations of locally interacting one - dimensional quantum lattices were actually shown rigorously to be efficient in the number @xmath0 of particles @xcite ( i.e. , computation time and memory resources scale as polynomial functions of @xmath0 at fixed @xmath58 , whereas here we are interested in the scaling of computation time and memory with physical time @xmath58 ( in tl , @xmath90 ) , referred to as _ time efficiency_.    in this section we address the question of time efficiency implementing tdmrg for a family of ising hamiltonian ( [ eq : hamaut ] ) which undergoes a transition from integrable ( transverse ising ) to non - integrable quantum chaotic regime as the magnetic field is varied . we mainly consider time evolution in operator spaces @xcite , say of density matrices of quantum states , or ( quasi ) local operator algebras . note that time evolution of pure states is often ill defined in tl @xcite . as a quantitative measure of time efficiency we define and compute the minimal dimension @xmath91 of matrix product operator ( mpo ) representation of quantum states / observables which describes time evolution up to time @xmath58 within fidelity @xmath92 . our main question concerns possible scaling of @xmath91 for different types of dynamics , and indeed we shall demonstrate a correspondence between , respectively , quantum chaos or integrability , and exponential or polynomial growth of @xmath91 .    the key idea of operator - valued tdmrg @xcite is to represent any operator in a matrix product form , @xmath93 in terms of @xmath94 matrices @xmath95 of fixed dimension @xmath96 . the number of parameters in the mpo representation is @xmath97 and for sufficiently large @xmath96 it can describe any operator . in fact , the minimal @xmath96 required equals to the maximal rank of the reduced super - density - matrix over all bi - partitions of the chain . the advantage of mpo representation lies in the fact that doing an elementary local one or two qubit unitary transformation @xmath98 can be done locally , affecting only a pair of neighboring matrices @xmath95 . we will illustrate evolution of ising chain ( [ eq : hamaut ] ) , with _ open _ boundary conditions , for two different magnetic field values : ( i ) an integrable ( regular ) case @xmath99 with transverse magnetic field , and ( ii ) non - integrable ( quantum chaotic ) case @xmath100 with@xmath101 tilted magnetic field . to confirm that @xmath102 , and @xmath103 , indeed represent generic quantum chaotic , and regular , system , respectively , we calculated level spacing distribution ( lsd ) of their spectra ( shown in fig .  [ fig : lsd ] ) . lsd is a standard indicator of quantum chaos  @xcite . it displays characteristic level repulsion for strongly non - integrable quantum systems , whereas for integrable systems there is no repulsion due to existence of conservation laws and quantum numbers . evolution by tdmrg proceeds as described in @xcite using an approximate trotter - suziki factorization , for some time step @xmath104 , of the evolution operator @xmath105 in terms of 2-qubit gates . and for each two qubit gates , the matrices @xmath95 can be updates using a singular value decomposition with some truncation error @xmath106 . interestingly , it has been found @xcite that the sum of all truncation errors up to time @xmath58 , denoted by @xmath107 ( provided @xmath104 is small enough so that the trotter - suzuki factorization error is negligible ) is simply proportional to _ fidelity error _ , namely @xmath108 where @xmath109 and numerical constant @xmath110 does not dependent on either @xmath104 , @xmath96 or @xmath0 . the central quantity we are going to study is @xmath111 which is the minimal dimension @xmath96 of matrices @xmath112 in order for the total truncation error @xmath107 to be less than some error tolerance @xmath113 , for evolution to time @xmath58 . in numerical experiments shown we take @xmath114 except for simulating thermal states of quenched hamiltonians where @xmath115 .      for a reference we start by investigating the evolution of pure states following the basic tdmrg algorithm @xcite . we studied time efficiency of simulation of _ pure states _ in schr " odinger picture , for which many examples of efficient applications exist , however all for initial states of rather particular structure , typically corresponding to low energy sectors of few quasi - particle excitations or to low dimensional invariant subspaces . treating other , typical states , e.g. eigenstates of unrelated hamiltonians , linear combinations of highly excited states , or states chosen randomly in the many - particle hilbert space , we found that , irrespectively of integrability of dynamics , tdmrg is _ not _ time - efficient , i.e. @xmath91 typically grows exponentially in time even in the integrable case of transverse field ( consistently with a linear growth of entanglement entropy  @xcite ) . numerical results are summarized in fig . [ fig : pure ] . we continue by discussing the time efficiency of operator - valued tdmrg method using mpos ( [ eq : mpo ] ) . let us first study the case where the initial operator is a local operator in the center of the lattice @xmath116 . in the integrable case time evolution @xmath117 can be computed exactly in terms of jordan - wigner transformation and toeplitz determinants @xcite , however for initial operators with infinite _ index _ [ sect.2 , 1st of refs.@xcite ] is half the number of fermi operators in jordan - wigner transformation of @xmath118 and is a conserved quantity for @xmath103 . ] like e.g for @xmath119 , @xmath1 , the evolution is rather complex and the effective number of terms ( pauli group elements ) needed to span @xmath117 grows exponentially in @xmath58 . in spite of that , our numerical simulations shown in fig . [ fig : local ] strongly suggest the linear growth @xmath120 for initial operators with infinite index . quite interestingly , for initial operators with _ finite _ index , @xmath121 saturates to a finite value , for example @xmath122 for @xmath123 , or @xmath124 for @xmath125 . in non - integrable cases the rank has been found to grow exponentially , @xmath126 with exponent @xmath127 which does _ not _ depend on @xmath113 , properties of @xmath128 or @xmath0 , for large @xmath0 . for @xmath129 we find @xmath130 .      in physics it is often useful to consider extensive observables , for instance translational sums of local operators , e.g. the hamiltonian @xmath88 or the total magnetization @xmath131 . as opposed to local operators , extensive initial operators , interpreted as w - like states in operator space , contain some long - range ` entanglement ' so one may expect that tdmrg should be somewhat less efficient than for local operators . indeed , in the integrable case we find for extensive operators with _ finite _ index that @xmath91 does no longer saturate but now grows linearly , @xmath120 , whereas for extensive operators with _ infinite _ index the growth may be even somewhat faster , most likely quadratic @xmath132 but clearly slower than exponential . in the non - integrable case , we again find exponential growth @xmath133 with the same exponent @xmath134 as for local initial observables . the results are summarized in fig . [ fig : transl ] .      in the last set of numerical experiments we consider time efficiency of the evolution of a thermal initial state @xmath135 under a sudden change of the hamiltonian at @xmath136 , namely @xmath137 . again , we treat two situations : in the first case we consider change after which the hamiltonian remains integrable , @xmath138 , while in the other case the change breaks integrability of the hamiltonian , @xmath139 . the initial state is prepared by means of imaginary time tdmrg from initial identity operator using the same mpo rank @xmath96 as it is later used for real time dynamics . we find , consistently with previous results , that at high temperature ( @xmath140 ) the rank @xmath91 grows very slowly , perhaps slower than linear , in the integrable case , and exponentially @xmath133 , in the non - integrable case . interestingly , at lower temperatures we find exponential growth in both cases , even in the integrable one . this is not unreasonable as the initial ( thermal ) state can be expanded in a power series in @xmath141 and the higher orders @xmath142 become less local with longer entanglement range as we increase the power @xmath143 . these results are summarized in fig . [ fig : rhot ] . in this section we would like to demonstrate the connection between quantum chaos , ( non)integrability and transport in non - equilibrium steady states of interacting quantum chains @xcite . within the linear response theory , the property of quantum mixing ( which is typically implied by quantum chaos , i.e. by the validity of random matrix level statistics @xcite , see also correlation functions in sect.[sect : losch ] ) , is typically synonymous for normal quantum transport since it implies _ finite _ kubo transport coefficients @xcite ( provided temporal correlation functions decay fast enough ) . however , here we would like to address the connection between quantum chaos and transport in far - from - equilibrium steady state , which may be beyond applicability of linear response . in particular , we are interested in the validity @xcite of the fourier law @xmath144 in quantum chains , relating the macroscopic heat flux @xmath145 to the temperature gradient @xmath146 .    to investigate this problem one has to deal with a finite open system connected to heat baths . here we consider an interacting quantum spin-@xmath2 chain ( [ eq : hamaut ] ) which exhibits the transition from integrability to quantum chaos as a parameter , e.g. the magnetic field , is varied . the standard treatment of this problem is based on the master equation , thus limiting numerical investigations to relatively small system sizes . by using this method , in an interesting paper @xcite , the decay of current correlation function in a model of non - integrable chain of quantum spins is computed . however , these results were not fully conclusive and the conclusions rely on linear response theory . also , in @xcite lindblad formalism was used to study the validity of fourier s law for different type of spin - spin interaction .    here we describe a different approach ( see ref.@xcite for details ) namely we follow the evolution of the system described by a _ pure _ state , which is _ stochastically _ coupled to an idealized model of heat baths . stochastic coupling is realized in terms of a local measurement at the boundary of the system and _ stochastic _ but _ unitary _ exchange of energy between the system and the bath . by this method we have been able to perform very effective numerical simulations which allow to observe a clear energy / temperature profile and to measure the heat current @xmath145 . again we consider ising spin chain in the magnetic field ( [ eq : hamaut ] ) of size @xmath0 , where the first and the last spin are coupled to thermal baths at temperatures @xmath147 and @xmath148 , respectively . in the non - integrable regime where the spectral statistics is described by rmt - the regime of quantum chaos - we found very accurate fourier law scaling @xmath149 , where @xmath0 is the size of the chain . in the integrable and near - integrable ( _ non - ergodic _ ) regimes instead , we found that the heat transport is ballistic @xmath150 . let us describe the numerical simulations . again we consider the autonomous model ( [ eq : hamaut ] ) at three particular cases : ( _ i _ ) _ quantum chaotic case _ @xmath151 at which lsd agrees with rmt and thus corresponds to the regime of quantum chaos , ( _ ii _ ) _ integrable case _ @xmath152 , at which lsd is close to the poisson distribution , and ( _ iii _ ) _ intermediate case _ @xmath153 at which the distribution lsd shows and intermediate character of weak level repulsion and exponential tail @xcite . let us now turn to study the energy transport in this model system . to this end we need to couple both ends of the chain of spins to thermal baths at temperatures @xmath154 . we have devised a simple way to simulate this coupling , namely the state of the spin in contact with the bath is statistically determined by a boltzmann distribution with parameter @xmath155 . our model for the baths is analogous to the stochastic thermal baths used in classical simulations @xcite and we thus call it a _ quantum stochastic thermal bath_. in the representation basis of @xmath156 the wave function at time @xmath58 can be written as @xmath157 where @xmath158 represents the _ up _ , _ down _ state of the @xmath159-th spin , respectively . the wave function at time @xmath58 is obtained from the unitary evolution operator @xmath160 . the interaction with the bath is not included in the unitary evolution . instead , we assume that the spin chain and the bath interact only at discrete times with period @xmath161 at which the states of the leftmost ( @xmath162 ) and the rightmost ( @xmath163 ) spins are stochastically reset . thus , the evolution of the wave function from time @xmath58 to time @xmath164 can be represented as @xmath165 where @xmath166 represents the _ stochastic _ action of the interaction with the left and right baths at temperatures @xmath167 and @xmath168 respectively . the action of @xmath169 takes place in several steps :    * ( _ i _ ) the wave function is first rotated by the angle to the eigenbasis of the components @xmath170 , @xmath171 of the edge spins along the field @xmath172 , that is @xmath173 . here , @xmath174 stands for the magnitude of the magnetic field . * ( _ ii _ ) a local measurement of the observables @xmath175 , @xmath176 is performed . then the state of the spins at the borders collapses to a state ( @xmath177,@xmath178 ) with probability @xmath179 , we put all coefficients @xmath180 with @xmath181 to zero . * ( _ iii _ ) the new state of the edge spins @xmath182 is stochastically chosen . after this action simulating the thermal interaction with the baths each of the edge spins is set to _ down _ , ( _ up _ ) state with probability @xmath183,(@xmath184 ) . if the new state of @xmath162 , or @xmath163 , is different than the one after step ( ii ) , then a simple unitary spin flip is performed to the wave - function . the probability @xmath185 depends on the canonical temperature of each of the thermal baths : @xmath186 * ( _ iv _ ) finally , the wave function is rotated back to the @xmath187 basis , @xmath188 . this completes the description of the interaction with the quantum stochastic bath . this interaction thus ( periodically ) resets the value of the local energy @xmath189 of the spins in contact with the baths . therefore , the value of @xmath161 controls the strength of the coupling to the bath . we have found that , in our units , @xmath190 provides an optimal choice . we have nevertheless performed simulations for other values of @xmath161 with qualitatively similar results . in particular , for weak couplings ( @xmath191 ) the heat conductivity does not depend on the coupling strength . note that our method does not correspond to a standard stochastic unraveling of a master equation for the density operator ( _ e.g. _ , in the lindblad form ) @xcite . however , we have tested that , using our prescription , averages over the ensemble of quantum trajectories or time averages of one given quantum trajectory are sufficient to reconstruct a density matrix operator @xmath192 that correctly describes the internal thermal state of the system in and out of equilibrium . for each run the initial wave - function @xmath193 of the system is chosen at random . the system is then evolved for some relaxation time @xmath194 after which it is assumed to fluctuate around a unique steady state . measurements are then performed as time averages of the expectation values of suitable observables . we further average these quantities over different random realizations of `` quantum trajectories '' .     out of equilibrium energy profile @xmath195 for the chaotic chain . the temperatures of the baths @xmath196 and @xmath197 , are both in the high temperature regime . results for chains of size @xmath198 ( crosses ) , @xmath199 ( open circles ) and @xmath200 ( solid circles ) are shown . the dashed line was obtained from a linear fit of the data for @xmath200 for the @xmath201 central spins . insets ( i ) and ( ii ) show the energy profile for the integrable and intermediate cases respectively , for @xmath202 . ]    in order to compute the energy profile we write the hamiltonian ( [ eq : hamaut ] ) in terms of local energy density operators @xmath203 : @xmath204 such that the total hamiltonian ( [ eq : hamaut ] ) can be written as @xmath205 apart from the boundary corrections . first we have performed equilibrium simulations in order to show that time averaged expectation values of the local energy density can be used to determine canonical local temperature . to this end we set the left and right baths to the same temperature @xmath155 . for low @xmath155 , the energy per site @xmath206 saturates to a constant which , together with the entire energy profile @xmath195 , is determined by the ground state . however , for larger @xmath207 , the energy profile is constant within numerical accuracy , and numerical simulations give @xmath208 , all results being almost independent of @xmath0 for @xmath209 . the numerical data for @xmath210 can be well approximated with a simple calculation of energy density for a two - spin chain ( @xmath211 ) in a canonical state at temperature @xmath155 , namely @xmath212 . therefore , if the temperatures of both baths are in high @xmath155 regime , then we can define the local temperature via the relation @xmath213 . we stress that equilibrium numerical data shown are _ insensitive _ to the nature of dynamics ( consistent with results of ref.@xcite ) , whether being chaotic , regular or intermediate . size dependence of the energy current in the chaotic chain with @xmath214 and @xmath197 . we show @xmath215 ( open circles ) and @xmath216 ( open squares ) . the dashed lines corresponds to @xmath217 scaling . in the inset , the size dependence of the energy current is shown for the integrable ( solid circles ) and the intermediate ( solid squares ) cases . ]    in fig . [ fig:3 ] we show the energy profile @xmath195 for an out of equilibrium simulation of the chaotic chain . in all non - equilibrium simulations , the temperatures of the baths were set to @xmath218 and @xmath219 . after an appropriate scaling the profiles for different sizes @xmath0 collapse to the same curve . more interesting , in the bulk of the chain the energy profile is in very good approximation linear . in contrast , we show that in the case of the integrable ( inset i ) and intermediate ( inset ii ) chains , no energy gradient is created which is a characteristic of ballistic transport . we now define the local current operators through the equation of continuity : @xmath220 =   - ( j_{n+1 } - j_n)$ ] , requiring that @xmath221 $ ] . using eqs . ( [ eq : h_local ] ) and ( [ eq : hamaut ] ) the local heat current operators are explicitly given by @xmath222 . in fig . [ fig:4 ] we plot @xmath215 as a function of the size @xmath0 of the system for sizes up to @xmath200 . the mean current @xmath145 is calculated as an average of @xmath223 over time and space @xmath159 . the energy difference was obtained from the energy profile as @xmath224 . two spins near each bath have been discarded in order to be in the bulk regime . since @xmath225 is an effective size of the truncated system , the observed @xmath226 dependence confirms that the transport is normal ( diffusive ) . moreover , also the quantity @xmath227 , where @xmath228 , shows the correct scaling with the size @xmath0 . on the other hand , in integrable and intermediate chains we have observed that the average heat current does not depend on the size @xmath229 , clearly indicating the ballistic transport . let us now come back to the kicked ising model ( [ eq : hamiltonian ] ) and try to consider some very elementary but fundamental questions considering its dynamics and non - equilibrium statistical mechanics . observing data of fig . [ fig : corr ] in section ( [ sect : losch ] ) one can conclude that the model perhaps displays an interesting order to chaos , or non - ergodicity to ergodicity transition when the integrability breaking parameter is increased . now we would like to inspect this transition more closely , and in particular understand the rate of relaxation to equilibrium in the ergodic and mixing case . we should stress right from the start that we are unable to prove any non - trivial statements about the model , but we can provide many suggestive numerical experiments which can be performed in a rather efficient way . we have learned in section [ sect : mps ] that it may be more fruitful to consider time evolution in the operator algebra spaces instead of in the spaces of pure states . let us go now a bit deeper into this subject .    for some related results on high - temperature relaxation in isolated conservative many - body quantum systems see e.g. refs.@xcite . time automorphism on unital quasi - local @xmath230 algebra ( see e.g. @xcite for introduction into the subject ) @xmath231 , @xmath232 of an infinite ki lattice , for one period of the kick , can be explicitly constructed by the following observations . formally , for any @xmath233 , @xmath234 , where @xmath235 is given by either ( [ eq : floquet],[eq : protocol],[eq : protocol2 ] ) . let @xmath236}$ ] , with @xmath237 , denote a finite , local @xmath238 dimensional algebra on a sub - lattice @xmath239 \subset \z$ ] , which is spanned by operators @xmath240 . it is straightforward to prove that dynamics is _ strictly local _ @xmath241 } \to { \cal a}_{[m-1,n+1]}. \label{eq : locality}\ ] ] in other words , the homomorphism @xmath242 ( [ eq : locality ] ) is a simple nontrivial example of a quantum cellular automaton as defined by schumacher and werner @xcite . @xmath231 can also be treated as a hilbert space with respect to the following inner product @xmath243 , where @xmath244 is a tracial state , @xmath245 for @xmath246}$ ] . this hilbert space can be in fact considered as a 1d lattice of 4-level quantum systems ( qudits with @xmath247 ) with the orthonormal basis @xmath248 labeled by an infinite sequence of 4-digits @xmath249 . restricting for a moment to a dimer lattice @xmath250}$ ] we can write the adjoint action of a 2-qubit gate @xmath251 ( [ eq : wgate ] ) in terms of a @xmath252 unitary matrix @xmath253 where , very explicitly @xmath254 . \quad\quad\ ] ] we should note that the map is unital @xmath255 , and that due to anti - unitary symmetry of ki model , the matrix @xmath256 is _ real_. we extend the map @xmath257 to entire algebra @xmath231 by @xmath258 , for any @xmath259}$ ] , @xmath260}$ ] . now , following the protocol ( [ eq : protocol2 ] ) we finish the construction of the time automorphism as a string of right - to - left ordered 2-qudit ( with @xmath247 ) gates @xmath261 explicitly , for any local observable @xmath262}$ ] , we have the following    * algorithm 1 * :    1 . set an initial vector : @xmath263 . 2 .   for @xmath264 : @xmath265 3 . the result is @xmath266    the algorithm produces exact result in @xmath267 multiplications and about the same number of additions . let @xmath268 } : { \cal a}_{\z } \to { \cal a}_{[m , n]}$ ] denote a _ linear orthogonal projector _ , satisfying @xmath268}(a\otimes e ) = ( \one|e)a$ ] if @xmath269}$],@xmath270}$ ] . let us define a truncated time evolution operator @xmath271 } = \ap_{[m , n ] } \at : { \cal a}_{[m , n ] } \to { \cal a}_{[m , n]}$ ] which we can actually implement on a computer with a finite memory register . the `` infinite - temperature '' time correlation function of ( traceless ) local quantum observables @xmath272}$ ] can be written as @xmath273 , where @xmath274 . an interesting question is , up to what time @xmath58 the @xmath275 can be computed numerically exactly with a finite computer register @xmath276 $ ] of @xmath277 qudits of size @xmath278 ? due to locality ( [ eq : locality ] ) of time homomorphism one can easily prove that @xmath279}^t a } , \quad { \rm for } \quad t \le 2l \label{eq : excorloc}\ ] ] hence the correlation functions are computable exactly up to time @xmath280 , and as we shall see later , truncated correlation function @xmath281 } a}$ ] often well approximates @xmath275 even at later times , or even its asymptotic decay . let us continue our discussion by considering time evolution for _ translationally invariant extensive _ ( tie ) observables . given some quasi - local observable @xmath282 we shall construct the corresponding tie observable by a formal mapping , @xmath283 in terms of lattice translation automorphisms @xmath284 , @xmath285 . the image of the entire quasi - local algebra under this mapping @xmath286 , is not a @xmath230 algebra , but it is a linear space which can be again turned into a hilbert space with the following inner product @xmath287}(x)}{\ap_{[-n , n]}(y)}\ ] ] where the domain of projector @xmath288}$ ] is extended to @xmath289 by continuity . orthonormal basis of @xmath289 is given by tie observables @xmath290 , for _ orders _ @xmath291 , and for uniqueness of notation , requiring @xmath292 . we shall interchangeably represent finite sequences of 4-digits with integers , @xmath293 . let @xmath294 be @xmath295 dimensional subspace spanned by tie observables @xmath296 with order @xmath297 , i.e. for @xmath298 having at most @xmath299 base-4 digits , so we have an inclusion sequence @xmath300 . since time and space automorphisms commute @xmath301 , one can immediately extend the time map onto the space of tie observables , @xmath302 by continuity . formally , we have @xmath303 . furthermore , locality ( [ eq : locality ] ) implies @xmath304 so we can write a simple adaptation of algorithm 1 for explicit construction of a time map of an arbitrary finite - order tie observable @xmath305 :    * algorithm 2 * :    1 . take the following pre - image of the tie observable @xmath306}$ ] , namely @xmath307 . 2 .   compute @xmath308 of @xmath309}$ ] according to * algorithm 1*. 3 .   transforming back to @xmath289 the result reads @xmath310    let us further define the _ natural truncations _ of tie space to order @xmath299 , @xmath311 as orthogonal projections @xmath312 , and truncated time evolution operators @xmath313 , which are naturally implemented on a computer by simply truncating overflowing coefficients @xmath314 . physically interesting question now concerns computation of time correlation functions between a pair of finite order ( say @xmath315 ) tie observables @xmath316 , namely @xmath317 , for example in fig . [ fig : corr ] we have shown the case of @xmath318 . as a consequence of locality ( [ eq : locality ] ) , and translational invariance , we find that the truncated evolution on @xmath294 , reproduces correlation functions exactly @xmath319 few remarks are in order : ( i ) truncated translationally invariant time evolution @xmath320 is perhaps more natural object to study than truncated local time evolution @xmath271}$ ] , for a simple reason that the truncation @xmath321 commutes with a shift @xmath322 , while @xmath268}$ ] does not . ( ii ) a space @xmath289 can be identified with translationally invariant linear functionals over @xmath231 , namely @xmath323 , @xmath324 , @xmath282 . we have @xmath325 and @xmath326 where hermitian adjoint maps @xmath327 and @xmath328 simply correspond to _ time reversed _ dynamics . ( iii ) convex subspace of positive translationally invariant functionals @xmath329 is an interesting invariant subspace of physical states , @xmath330 .      in classical mechanics of chaotic systems one typically observes that states ( phase - space densities ) develop small details in the course of time evolution at an exponential average rate . consequently , introducing a small stochastic noise of strength @xmath113 to perron - frobenius operator ( pfo , i.e. liouvillian propagator for discrete time dynamics ) , makes it non - unitary and shifts its spectrum inside the unit circle . typically , the effect of noise is equivalent to an ultraviolet cutoff - truncation of pfo - at the fourier scale @xmath331 , and often the leading eigenvalues of truncated pfo - the so - called ruelle resonances - remain frozen inside the unit circle in the limit @xmath332 ( or @xmath333 ) @xcite . for a general introduction to relaxation phenomena in classical hamiltonian dynamics see e.g. @xcite . let us now draw some some analogies with our quantum setting . we have seen that that the evolution @xmath334 somehow most closely resembles liouvillian evolution of classical hamiltonian dynamics . in hilbert space topology , operator @xmath334 is unitary and its spectrum lies on a unit circle , just like in the case of classical pfo . however , truncated @xmath335 matrices @xmath320 represent natural `` ultraviolet '' cutoff truncations for increasing orders @xmath299 . let us check numerically if some eigenvalues of these matrices remain frozen when @xmath336 .    indeed , as we demonstrate in fig . [ fig : s ] , we find several eigenvalues which converge as @xmath299 increases in the case of strongly non - integrable ( quantum chaotic ) case , with a gap between an eigenvalue of maximal modulus and the unit circle , whereas in the integrable case a set of @xmath299 eigenvalues touch the unit circle ( actually eigenvalue 1 is @xmath299-fold degenerate ) . numerical results suggest the following speculative conclusions . let @xmath337 be the converged ( frozen ) eigenvalues of @xmath320 , and @xmath338 , @xmath339 the corresponding _ right _ and _ left _ eigenvectors , respectively . then for arbitrary pair @xmath340 , the time correlation function can be expressed in terms of spectral decomposition ( see e.g. @xcite ) @xmath341 the above relation is the contribution of the point spectrum and is exact if the spectrum is pure - point . however , in classical cases one may quite typically have various singular components and branch cuts @xcite . note that the denominator @xmath342 is finite although both vectors should have infinite @xmath343 norm @xmath344 , @xmath345 , for any eigenvalue away from the unit circle , @xmath346 . there is a simple relation between the spectrum of pfo and the ergodic properties of dynamics : ( i ) if there is a _ spectral gap _ , i.e. there exists @xmath347 such that for all @xmath159 , @xmath348 , then dynamics is _ exponentially mixing _ , @xmath349 ( ii ) if some eigenvalues are on the unit circle , meaning that the corresponding eigenvector coefficients should be in @xmath343 , then the system is _ non - mixing _ since there are correlation functions which do not decay . ( iii ) if some eigenvalues are at @xmath350 then the system is _ non - ergodic _ since the correlation functions may have non - vanishing time - averages . if @xmath351 is a complete set of orthonormalized eigenvectors corresponding to eigenvalue @xmath350 , @xmath352 ( and note that since we are on the unit circle : @xmath353 ) then @xmath354 the latter ( iii ) happens in generic completely integrable quantum lattices , where @xmath351 correspond to an infinite sequence of conservation laws @xcite . furthermore , we have a strong numerical evidence that also in certain non - integrable quantum lattices @xcite , and also in ki model @xcite , one has a regime where few normalizable ( pseudo - local ) but not local ( like in integrable models ) conservation laws exist . this situation we call the regime of _ intermediate _ dynamics and is characterized by a non - vanishing _ stiffness _ @xmath355 signalling _ ballistic transport_.    in figure [ fig : c ] we compare the time autocorrelation function of the transverse magnetization @xmath356 computed in three different ways : ( 1 ) from exact time evolution @xmath357 on a finite lattice of length @xmath0 with periodic boundary conditions , ( 2 ) iteration of truncated ta matrix on infinite lattice @xmath358 and ( 3 ) asymptotics based on ( few ) leading eigenvalue resonance(s ) ( using formula ( [ eq : corr ] ) in terms of @xmath359 and @xmath360 . ) .    . gray level indicates the spectral gap @xmath361 of @xmath320 as a function of @xmath16 and @xmath17 at fixed @xmath362 . white regions correspond to @xmath363.,width=302 ]    we note that the leading eigenvalue and eigenvector of truncated ta @xmath320 can most efficiently be computed using our algorithm 2 as a key step of an iterative _ power - method_. in this way we were able to perform calculations of the leading ruelle resonances up to @xmath364 in contrast to full diagonalization of truncated matrices @xmath320 which were feasible only up to @xmath365 . let us observe the structure of the eigenvector coefficients @xmath366 corresponding to the leading eigenvalue . numerical results ( see fig .  [ fig : v ] , see also subsect.[sect : scaling ] later ) strongly suggest self - similar behaviour upon multiplying the code @xmath298 by @xmath367 which is a consequence of the fractal structure of the transfer matrix ( see illustration in ref.@xcite ) .     as a function of @xmath16 at @xmath368 at different truncation orders @xmath299 . in the inset we indicate a line of transition in 2d phase diagram . , width=151 ]     as a function of @xmath369 at @xmath370 at different truncation orders @xmath299 . in the inset we indicate a line of transition in 2d phase diagram . , width=151 ]    the stiffness @xmath371 and the spectral gap @xmath372 may be considered as _ order parameters _ , characterizing a particular kind of dynamical phase transition , namely the transition from non - ergodic dynamics - ordered phase , where @xmath355 for a typical @xmath373 which is not orthogonal to all conservation laws @xmath351 , i.e. normalizable eigenvectors of @xmath334 with eigenvalue @xmath350 . ] , or @xmath374 , to an ergodic and mixing dynamics - disordered phase , where @xmath375 , for all traceless @xmath324 , or @xmath376 . we know that ki model is non - ergodic in integrable regimes let us consider a fixed transverse field case @xmath11 , and start to switch on a small amount of longitudinal field @xmath16 . the interesting question is whether the transition happens for infinitesimal integrability breaking parameter @xmath16 in tl , or at a finite - _ critical _ field . in fig . [ fig : phasediagram ] , we fix @xmath362 and plot a two dimensional phase diagram of the spectral gap @xmath377 . it is clear that we have different behaviours in different regions of parameter space , for example we identify two : ( i ) type i transition : if the transverse field is roughly on the interval @xmath378 $ ] , then the spectral gap opens in the fastest possible manner which is allowed by a @xmath379 symmetry and the analyticity of the problem , namely @xmath380 . see fig . [ fig : typei ] . ( ii ) type ii transition : if the initial transverse field @xmath381 , or @xmath382 , then the gap opens up in a much more abrupt - perhaps a discontinuous way . we give an example by scanning the diagonal transition , i.e. putting @xmath369 and increasing @xmath16 from zero . numerical results , shown in fig.[fig : typeii ] can not be made fully conclusive , but they are not inconsistent with a conclusion that an abrupt transition to ergodic behaviour takes place at @xmath383 . there is another possibly interesting way of characterizing the transition , i.e. in terms of pseudo - local translationally invariant conservation laws @xcite . such conservation laws are the square normalizable elements @xmath384 , which are mapped onto themselves under the dynamics @xmath385 . we shall first make a non - trivial variational mpo ansatz for elements of @xmath289 , namely let us take an auxiliary vector space @xmath386 , where @xmath387 . then any operator @xmath388 , which is formally written in terms of mpo on the _ infinite _ spin chain @xmath389 where @xmath390 have the block matrix form ( for @xmath391 ) : @xmath392 represents a translationally invariant pseudo - local operator , i.e. an element of @xmath289 , provided that @xmath393 and @xmath394 , where @xmath395 is a spectral matrix norm and @xmath396 s stand for arbitrary matrices / vectors . of course , converse can not be generally true , not any element of @xmath289 can be written as mpo ( [ eq : timpo ] ) with finite @xmath96 , but still there are elements of the form ( [ eq : timpo ] , [ eq : mpodata ] ) which are not in @xmath294 for _ any _ finite @xmath299 . there exist a straightforward algorithm which performs time evolution on mpo data ( [ eq : mpodata ] ) , namely @xmath397 is also of the form ( [ eq : timpo],[eq : mpodata ] ) with dimension @xmath398 . we shall now make the following simple numerical experiment . let us fix @xmath96 , setting @xmath399 representing the simplest @xmath334-invariant subclass of ( [ eq : timpo],[eq : mpodata ] ) , and optimize ( maximize ) the fidelity - like quantity @xmath400 within this class of operators . let us write an operator which maximizes @xmath401 for a given @xmath96 as @xmath402 . note that @xmath403 gives a strong - topology measure of conservation of observable @xmath402 in one step of time evolution , so @xmath404 only for exact conservation laws . increasing @xmath96 may improve fidelity , if pseudo local conservation laws exist to which @xmath402 may converge , however in ergodic and mixing situation where no exact pseudo - local conservation laws exist , increasing @xmath96 should have no significant effect to fidelity @xmath405 . this is exactly what we observe in ki model following a line of type ii transition ( see fig .  [ fig : conslaw ] ) . numerical results of section [ sect : mps ] suggested that operator space entanglement measures can be used to characterize the complexity of time evolution , namely the minimal required rank @xmath406 of mpo ansatz is simply related to entanglement entropy of a time - evolving local observable , which is interpreted as a hilbert space vector . in order to make things as simple and precise as possible we go back to the time evolution automorphism @xmath242 over the quasi - local spin algebra @xmath231 . let us take some truncation order @xmath407 and consider the truncated map @xmath408}$ ] . starting with local operators on a single site @xmath409}$ ] this truncated map is exact up to time @xmath410 . writing time evolved observable at any instant of time as @xmath411 } a(0 ) =   \sum_{\ve{s } } a^{(t)}_{(s_{-n},\ldots s_n ) } { { | s_{-n},\ldots , s_n ) } } $ ] in terms of a ` wave - function ' @xmath412 , and partitioning a sub - lattice at @xmath413 , as @xmath414 =   [ -n , m-1 ] \cup [ m , n]$ ] , we can define @xmath415 reduced super - density matrix as @xmath416 since @xmath417 is interpreted as a ` _ pure state _ ' , namely a vector from @xmath418}\otimes{\cal a}_{[m , n]}$ ] , the operator space entanglement is most simply characterized either by von neuman entropy @xmath419 or linear entropy @xmath420 where @xmath421 ^ 2 $ ] is a purity of reduced super - density matrix .    in fig . [ fig : pursuperrho ] we plot the entanglement purity @xmath422  for close to symmetric bipartition @xmath423 where the entanglement is expected to be maximal  for three different characteristic cases of ki model , which will be in the following referred to as : _ quantum chaotic _ ( qc ) , @xmath424 , _ integrable _ ( in ) , @xmath425 , and _ non - ergodic _ ( ne ) non - integrable case , @xmath426 . we find that in qc case purity decreases exponentially @xmath427 , meaning @xmath428 , where the exponent @xmath127 is independent of the initial observable @xmath429 and asymptotically independent of @xmath299 . on the other hand , in in case @xmath422 does not decay at all so the resulting dynamical entropy @xmath430 , whereas in ne case @xmath422 decays slowly , likely slower than exponentially . the eigenvectors of @xmath334 corresponding non - unimodular eigenvalues seem to exhibit a certain scaling invariance ( fig . [ fig : v ] ) . here we would like to explore this property in a little bit more detail . for that purpose we again explore the map @xmath408}$ ] on the quasi - local algebra @xmath431}$ ] since the representation of dynamics is conceptually simpler ( algorithm 1 ) than dynamics @xmath432 on @xmath433 . in particular , it is worth to mention that if one traces out an additional qudit after each time step , then the dynamics is exact on a _ closed _ set of @xmath434 super - density matrices and has a simple explicit form in terms of a completely positive matrix map . ] @xmath435 t^\dagger_{n - m}\right\ } ,   \label{eq : rdynamics1 } \\ & & t_{r } = \prod_{0\le k \le r}^- \one_4^{\otimes k}\otimes w\otimes \one_4^{\otimes ( r - k ) } , \label{eq : rdynamics2}\end{aligned}\ ] ] namely no truncation is needed since at time @xmath436 we are describing an exact observable on @xmath437}$ ] . we write @xmath438 for tracing out the least significant qudit , and @xmath439 is an elementary @xmath440 projector . note that @xmath441 is just a matrix of @xmath271}$ ] in the canonical basis @xmath442 . it is rather trivial to exactly solve this dynamics for a small finite @xmath443 , however this does not yield physically very useful information about the ki dynamics . one would wish to study the correct tl by first taking @xmath444 , and only after that @xmath445 , however this task seems almost computationally intractable . still , we were able to make some modest numerical experiments exploring this question , suggesting that for sufficiently strong integrability breaking ( say qc case ) the asymptotic matrix @xmath446 has a remarkable scale invariance if we coarse - grain it by tracing over its @xmath440 blocks : @xmath447 where @xmath448 is some scaling factor . this seems to be true for both orders of the limits @xmath449 , @xmath444 , although better numerical results have been obtained for the ` incorrect ' limit , namely letting the number of iterations @xmath449 for a finite register size @xmath159 , and then checking the convergence of results with increasing @xmath159 .    before discussing numerical results we note another useful observation let us define and briefly study ki chain on a _ semi - infinite _ lattice @xmath450 $ ] , with the hamiltonian ( [ eq : hamiltonian ] ) for @xmath90 and open boundary condition on the left edge . now we consider a quasi - local algebra @xmath451 . the time automorphism is again strictly local @xmath452 } \to { \cal a}_{[0,n+1]}$ ] , and can be written as a semi - infinite product @xmath453 . in the definition of the truncated time map @xmath454 } = \ap_{[0,n ] } \at^+ : { \cal a}_{[0,n ] } \to { \cal a}_{[0,n]}$ ] truncation is needed only on the right edge . note that due to this property , simulation of local observables ( localized near the edge of the lattice ) is twice as efficient than on doubly - infinite lattice , meaning that with the same size of computer register one can exactly simulate for twice longer times . for example , computation of correlation functions @xmath455^t a}$ ] is exact @xmath456^t a } = \braket{b}{[\at^+_{[0,n]}]^t a},\quad { \rm for}\;\ ; t \le 2(n - q ) , \quad { \rm if}\;\ ; a , b\in { \cal a}_{[0,q]}.\ ] ] numerically inspecting the correlation functions of a simple local spin @xmath457 in fig . [ fig : cfsc ] we find clear asymptotic exponential decay for qc case , whereas for in and ne cases we find non - vanishing plateaus in the correlation function , i.e. non - vanishing stiffness @xmath458 which signals non - ergodicity and existence of local ( for integrable cases ) and pseudo - local ( for non - ergodic and non - integrable cases , e.g. ne ) conservation laws of @xmath459 . we note that the asymptotic correlation decay in non - integrable cases , say qc and ne , seems quite insensitive to increasing truncation order @xmath159 - indicating that the leading eigenvalues of @xmath454}$ ] remain frozen when increasing @xmath159 . note that the asymptotic exponents of correlation decay for a semi - infinite chain are not the same as for an infinite one , i.e. the point spectra of @xmath454}$ ] and @xmath460 are in general different , however we have some indications to believe that their phase diagrams should agree , namely ergodic regimes in the ki model on semi - infinite chain model are in one - to - one correspondence with ergodic regimes of the model on infinite chain .    however , the most remarkable feature of dynamics @xmath459 is the following . splitting the truncated semi - lattice as @xmath461=[0,m-1]\cup[m , n]$ ] and following the time evolution of an observable @xmath462}]^t a(0 ) =   \sum_{\ve{s } } a^{(t)}_{(s_{0},\ldots s_n ) } { { | s_0,\ldots , s_n ) } } $ ] in terms of a ` super - wave - function ' @xmath463 , one can again define the reduced super - density matrix as @xmath464 the dynamical equation for @xmath465 in terms of @xmath466 is _ exactly the same _ as for doubly - infinite lattice , namely eqs . ( [ eq : rdynamics1],[eq : rdynamics2 ] ) . hence also the conjecture ( [ eq : scaling ] ) on scaling invariance of @xmath446 should be the same for the two lattice topologies . however , numerical computations are much easier and thus the results are more suggestive for the semi - infinite case .     in log scale ( units are arbitrary ) versus the partial tracing index @xmath143 , for cases qc ( a ) , ne ( b ) and in ( c ) . diamonds , stars , and squares represent data for the semi - infinite ki lattice with truncation size @xmath467 , ( diamonds , stars ) and @xmath468 ( squares ) , all for @xmath469 . finite number of time steps @xmath470 ( diamonds ) , just before the absorbing boundary affects any of the data shown , is compared to steady state observable @xmath471 ( stars ) . in the non - integrable cases ( a , b ) , data are compared also with steady state @xmath471 simulation of the same local initial observable @xmath472 on the two - sided ( doubly - infinite ) ki lattice with truncation size @xmath473 ( @xmath364 ) , and @xmath474 ( triangles ) . , width=604 ]    let us now discuss some numerical results . we have always started from local initial operator @xmath472 . we took @xmath159 as large as allowed by existing computing resources , namely @xmath467 for the semi - infinite chain and @xmath473 for the infinite chain , and compared the data for asymptotic matrices @xmath475 ( in numerics @xmath58 has been chosen such that the results converged , typically @xmath476 ) with finite time data @xmath477 where time @xmath478 was set as large as allowed so that the data were still exact and no truncation was needed , typically @xmath479 . in all cases , numerical results were quite insensitive to small changes in truncation order @xmath159 . first we have computed the scaling of the principal matrix element , or the _ partial norms _ @xmath480 which , assuming ( [ eq : scaling ] ) , should asymptotically scale as @xmath481 ( see fig.[fig : scalsinorm ] ) . if asymptotic dynamics @xmath449 is determined by normalizable eigenvectors of @xmath459 , which necessarily correspond to uni - modular eigenvalues , then we should have @xmath482 . in qc case a clear scaling was observed for both topologies ( @xmath483 and @xmath484 ) with the same exponent @xmath448 , however the exponent was slightly different for @xmath475 and @xmath477 . in the cases of non - ergodic dynamics ( ne and in ) the results for two topologies were quite different . for semi - infinite topology we find very clearly that @xmath485 indicating that @xmath486 and even @xmath487 can be written as @xmath343 convergent sums of local operators .     for initial tracing index @xmath488 , and @xmath489 ( diamonds ) , @xmath490 ( stars ) , @xmath491 ( squares ) , @xmath492 ( triangles ) , for non - integrable cases of ki chain on semi - infinite lattice truncated at @xmath467 , namely for case qc ( a , b ) , and case ne ( d , e ) . finite time data at @xmath470 kicks are shown in ( a , d ) , while steady state observables @xmath471 are analyzed in ( b , e ) , all starting from initial observable @xmath493 . for comparison , asymptotic steady state @xmath471 data for two - sided , doubly infinite ki chain , truncated at @xmath494 , and with initial tracing index @xmath474 , are shown in ( c , f ) , namely for qc case ( c ) and ne case ( f ) . note that data in plot ( b ) are practically exactly overlapping.,width=604 ]    as a more quantitative test of conjecture ( [ eq : scaling ] ) we compare the upper - left ( ` most important ' ) @xmath252 block of the super density matrix scaled to a unit principal element @xmath495 . in fig.[fig : scalsi1d ] we plot the diagonal elements @xmath496 for different @xmath143 , while in fig.[fig : scalsi2d ] we plot 2d charts of the entire scaled density matrices @xmath497 . indeed we find for qc case that the matrix @xmath497 is practically insensitive to increasing truncation ( @xmath143 ) and to topology of the lattice ( semi - infinite versus infinite ) , both for finite time @xmath498 and ` steady - state ' @xmath471 . on the other hand , in non - ergodic cases , the scaling ( [ eq : scaling ] ) is typically broken . however , it seems to be observed in the steady state ( @xmath498 ) of ne case , which is ( in our setting ) probably a non - physical but still quite robust effect due to truncation ( a kind of absorbing boundary condition ) . - the grayness level is proportional to @xmath499 - for initial tracing index @xmath469 , and @xmath500 ( columns plots ) , for six different cases , qc , ne , and in of truncated semi - infinite ki lattice ( @xmath467 ) , at finite time @xmath470 , and asymptotic steady state @xmath471 ( row plots ) , always starting form initial observable @xmath493 . @xmath252 ` most important ' matrix elements are plotted with the matrix site @xmath501 at lower - left corner.,width=604 ]    summarizing this subsection , we conjectured that in the regime of quantum chaos reduced super - density matrices of time - evolved observables typically obey the scaling law ( [ eq : scaling ] ) with the exponent @xmath448 which only depends on global dynamics on quasi - local algebra of observables and not on a particular choice of initial observable . we expect that accurate numerical calculations of exponent @xmath448 would be possible within a certain quantum dynamic renormalization group scheme , however its precise formulation at present remains an open problem . in section [ sect : mps ] , further elaborated in subsection [ sect : opent ] , we have proposed an entanglement in operator space of quasi - local algebra as a possible new measure of quantum algorithmic complexity . here we would like to compare this briefly to some established proposals of quantum dynamical entropy , such as for example cnt entropy @xcite or af entropy @xcite , both being ideally suited for a quantum dynamical system formulated in terms of time automorphism @xmath242 , tracial invariant state @xmath502 , and quasi - local @xmath230 algebra @xmath231 . we shall here briefly review af entropy which is conceptually simpler . let us start by taking a set of say @xmath503 elements @xmath504 of @xmath231 which form a partition of unity , namely @xmath505 . following the dynamics up to integer time @xmath58 , @xmath506 , we construct a set of @xmath507 elements depending on a multi - index @xmath508 , namely @xmath509 from the homomorphism property and unitality of dynamics it follows that @xmath510 is also a partition of unity , @xmath511^ * a^{(t)}_{\ve{\alpha } } = \one $ ] . hence using an invariant state @xmath502 one can form a _ positive _ , _ hermitian _ , _ trace - one _ , @xmath512 matrix @xmath513^ * a^{(t)}_{\ve{\alpha}}\right).\ ] ] @xmath514 can clearly be interpreted as a density matrix pertaining to dynamically generated partition , and its von neumann entropy generated per unit time defines the af entropy @xmath515,\ ] ] where , strictly speaking , supremum over all possible generating partitions has to be taken . it is not surprising that practical evaluation of @xmath516 is impossible except for rather trivial cases , such as dynamics generated by shift automorphism @xmath517 @xcite . however , one can easily show that a related _ linear af entropy _ @xmath518 ^ 2 , \label{eq : aflin}\ ] ] is tractable much more easily , while the behaviour of @xmath519 and @xmath520 is likely to be similar in practice . the key observation is to write the purity @xmath521 ^ 2 =   \sum_{\ve{\alpha},\ve{\beta } } \rho^{(t)}_{\ve{\beta},\ve{\alpha } } \rho^{(t)}_{\ve{\alpha},\ve{\beta}}$ ] in terms of dynamics over a product algebra @xmath522 , with time automorphism @xmath523 and an invariant state @xmath524 . to this product structure we add a linear map @xmath525 depending on a generating partition @xmath526 @xmath527 if @xmath528 is a unit element in @xmath231 then @xmath529 is a unit element in @xmath530 , and purity can be written using a transfer - matrix - like approach as @xmath531^t ( \one \times \one)\right\ } .\ ] ] the idea can be worked out in detail for a _ complete _ generating partition of a local sub - algebra @xmath532}$ ] of size @xmath533 . there it turns out that , due to locality of dynamics ( [ eq : locality ] ) , the resulting purity is _ independent _ of @xmath534 , i.e. the supremum ( [ eq : aflin ] ) is already achieved by a rather modest partition of @xmath535 elements , and the map @xmath536 can be _ factored _ into a direct product of two maps acting separately on _ two independent copies _ of a _ semi - infinite lattice_. leaving out some technical details of derivation , the final result reads as follows . let us consider a time dependent @xmath537 matrix @xmath538 , representing a state on @xmath539 , with initial value @xmath540 , and dynamics given by the following completely positive matrix map @xmath541}\left\{t_{t+1 }   \left [ e_{00}\otimes \one_4\otimes q(t ) \otimes e_{00 } \right ]   t^\dagger_{t+1}\right\ } , \label{eq : qdynamics}\ ] ] where @xmath542 } r)_{j , k } = \sum_{s=0}^{15 } r_{s + 16j , s+16k } $ ] traces out two least important qudits , and unitary time evolution matrix @xmath543 is given in ( [ eq : rdynamics2 ] ) . then the purity , and linear dynamical entropy ( lde ) , are simply given as @xmath544 ^ 2 $ ] , @xmath545 , respectively . the asymptotic increase of @xmath79 per unit time yields the linear af entropy ( [ eq : aflin ] ) . again , for practical calculations it is convenient to consider truncation of matrices @xmath538 after each iteration ( [ eq : qdynamics ] ) , say at dimension @xmath278 . in fig . [ fig : smb ] we plot lde for different cases of ki dynamics , and we observe that lde is always clearly growing linearly @xmath546 , so the af entropy is always strictly positive , even in non - ergodic ( ne ) and integrable ( in ) cases , and that the results are robust and stable against changing the truncation order @xmath299 . perhaps this finding appears surprising , but one has to bear in mind that af and cnt entropies can be positive even for simpler dynamics , such as quasi - free flows on @xmath230 algebras .    for _ ergodic _ dynamical systems on @xmath230 algebras one can use shannon - mcmillan - breiman ( smb ) theorem ( for classical smb theorem see e.g. @xcite , and @xcite for a possible quantum extension ) , which states that for _ typical _ sequences @xmath547 , _ multi - time correlation function _ ( mtcf ) should decay exponentially @xmath548 where the exponent @xmath81 , which should be essentially independent of @xmath549 , is equal to a dynamical entropy of the map @xmath242 with respect to an invariant state @xmath502 . for an interesting application of smb theorem in the context of quantum dynamical chaos see ref . @xcite . in our numerical experiments we considered truncated dynamics @xmath408}$ ] , writing truncation order as @xmath550 , and computed two kinds of mtcf : ( i ) for a uniform sequence @xmath551 , where @xmath552 ( or @xmath553 in which case the truncated lattice was placed as @xmath554 $ ] , hence @xmath555 ) , we estimated mtcf ( [ eq : mctf ] ) as @xmath556 } a_1 \at_{[-n , n ] }   \cdots \at_{[-n , n ] } a_1}$ ] . ( ii ) for a random sequence @xmath547 , we took a complete generating partition @xmath557 on @xmath558}$ ] with @xmath559 elements , for @xmath560 or @xmath561 , and computed an average square modulus of mtcf @xmath562 } a_{\alpha_1 } \at_{[-n , n ] }   \cdots \at_{[-n , n ] } a_{\alpha_{t-1}}}|^2 \rangle}}$ ] over 20 randomly sampled sequences @xmath547 . in both cases ( i , ii ) we have found a rather good agreement between @xmath563 and lde @xmath79 for the case of quantum chaotic dynamics ( fig .  [ fig : smb ] ) . for average random mtcf we found rather good agreement even for non - ergodic cases ( ne and in ) , however mtcf for the uniform sequence exhibited big oscillating fluctuations there , indicating simply that smb theorem does hold for non - ergodic dynamics of ki model . summarizing , it seems that af entropy , even though it very cleanly generalizes the concept of kolmogorov - sinai classical dynamical entropy to quantum dynamical systems , can not be used as an indicator of quantum chaos or an indicator of computational complexity of quantum dynamical systems . nevertheless , it has to be stressed that quantum dynamical entropies can be related to the notion of quantum algorithmic complexity , mentioned in section [ sect : mps ] , by a quantum version of brudno theorem established in ref . in the present article we have reviewed several possible approaches to describe dynamic instabilities , relaxation phenomena , and computation complexity in the simple model of one - dimensional non - integrable locally interacting quantum many body dynamics . we have argued that the essential non - equilibrium statistical properties of quantum dynamical systems - in the absence of idealized external baths - may be crucially related to the integrability of the system , or in the complementary case , to the existence of the regime of quantum chaos . the general flavor which remains after such studies is that the non - integrable quantum many body problem at high temperature will preclude any exact and complete solution by its very nature . still it is hoped that a more complete ergodic theory of such systems could be developed , allowing for example for exact calculation of relaxation rates , scaling exponents of resonance eigenvectors , etc . it has been shown in many cases that integrable quantum systems have anomalous non - equilibrium statistical mechanics at high temperature , for example they exhibit ballistic transport . this can easily be understood as being a consequence of existence of ( an infinite sequence of ) exact local conservation laws which prevents quantum ergodicity , similarly as existence of canonical action variables in classical integrable systems prevents classical ergodicity . an interesting open question is the following : how strong integrability breaking perturbation is needed , in generic cases , to break all the exact conservation laws and yield normal ( diffusive ) transport ? in other words : a quantum kam theory is needed ! numerical experiments shown in this paper suggest an interesting possibility , namely that in some cases a finite , critical perturbation strength is required . i would like to thank g. casati , t. gorin , c. mejia - monasterio , t. h. seligman , and m. nidari , for collaborations on related projects and acknowledge stimulating discussions with f. benatti , j. eisert , s. fishman and p. zoller . this work has been supported by the grants p1 - 0044 , and j1 - 7347 of slovenian research agency . 9                  g.  vidal , phys . rev .  lett . * 91 * , 147902 ( 2003 ) ; _ ibid . _ * 93 * , 040502 ( 2004 ) ; s.  r.  white and a.  e.  feiguin , phys .  rev .  lett . * 93 * , 076401 ( 2004 ) ; a.  j.  daley , c.  kollath , u.  schollw " ock and g.  vidal , j.  stat .  mech . * 4 * , p04005 ( 2005 ) ; g.  vidal , cond - mat/0605597 . d.  gobert , c.  kollath , u.  schollw " ock and g.  sch " utz , phys . e * 71 * , 036102 ( 2005 ) ; u. schollw " ock and s.  r.  white , in g.  g. batrouni , and d.  poilblanc ( eds . ) : effective models for low - dimensional strongly correlated systems , p.155 , aip , melville , new york ( 2006 ) . f.  verstraete , j.j . garc ' ia - ripoll and j.i . cirac , phys.rev.lett . * 93 * , 207204 ( 2004 ) ; m.  zwolak and g.  vidal , phys.rev.lett . * 93 * , 207205 ( 2004 ) ; a.e . feiguin and s.r . white , phys.rev.b * 72 * , 220401(r ) ( 2005 ) . r.  sch " affer , h .- j . st " ockmann , t.  gorin , t.  h.  seligman , phys . lett . * 95 * , 184102 ( 2005 ) ; r.  sch " affer , t.  gorin , h .- j . st " ockmann , t.  h.  seligman , new j. phys . * 7 * , 152 ( 2005 ) ; t.  gorin , t.  h.  seligman and r.  l.  weaver , phys . rev . e * 73 * , 015202(r ) ( 2006 ) . o.  bratelli and d.  w.  robinson , _ operator algebras and quantum statistical mechanics i , ii _ , ( springer - verlag , 1979 , 1981 ) ; b.  simon , _ the statistical mechanics of lattice gases _ , ( princeton university press , 1993 ) .
Abstract:

Answer: the problem of characterizing complexity of quantum dynamics - in particular of locally interacting chains of quantum particles - will be reviewed and discussed from several different perspectives : ( i ) stability of motion against external perturbations and decoherence , ( ii ) efficiency of quantum simulation in terms of classical computation and entanglement production in operator spaces , ( iii ) quantum transport , relaxation to equilibrium and quantum mixing , and ( iv ) computation of quantum dynamical entropies . discussions of all these criteria will be confronted with the established criteria of integrability or quantum chaos , and sometimes quite surprising conclusions are found . some conjectures and interesting open problems in ergodic theory of the quantum many problem are suggested . _ invited review article for special issue of journal of physics a on quantum information _