chaos in ordinary language means disorder , randomness , absence of law , and unpredictability .
it played the major role in philosophy of ancient greeks .
already empedocles ( cca.490430 b.c . )
viewed the real world , , as a combination of the world of perfect order , ( today we may tentatively translate it as ) , and the world of complete disorder , . since chaos seemed to be too elusive for any kind of rigorous description , the attention of science has for long been focused only to sphairos .
this has changed over the past decades when numerous explorations have shown that chaos is not as lawless as originally thought . at present , the term _ chaos theory _ belongs to the common speak @xcite .    in both classical and quantum physics @xcite
is built on very sophisticated mathematics .
it sticks on the original idea of randomness , but at the same time it creates a new type of universality in the representation of reality .
the universality of chaos is based on ergodic character of motions and implies high efficiency of the statistically oriented description . since majority of systems in nature
are chaotic , at least to a considerable extent , this type of description is often the ultimate one .
the routes to more deterministic descriptions may exist but they are usually impassable .
chaos thus belongs to very fundamental subjects in physics .    among the physical systems that materialize the signatures of chaos in nature , an important place is held by atomic nuclei .
already the fact that even more than 80 years after the discovery of their composition we are not able to present a satisfactory microscopic theory of nuclear structure is a good reason to think that nuclei are indeed chaotic objects par excellence .
it is not an accident that nuclei were at dawn of the field named quantum chaos .
however , atomic nuclei not only show well recognized signs of chaos , they also provide a variety of models that can be used to better understand the rules of chaos in a more general context . such a path is followed in this text .
our aim is not to draw an exhaustive map of the overlap between physics of chaos and physics of nuclei
such reviews exist @xcite and we may only recommend to read them .
instead , we want to make a passage through selected topics of classical and quantum chaos using models taken from the description of nuclear collective dynamics .
before diving into the ocean of chaos , we have to briefly introduce the models we will be playing with .
these are the old ( and the good ) geometric collective model ( below abbreviated as gcm ) and the newer and more sophisticated interacting boson model ( ibm ) of nuclei .
both these models attempt to describe nuclear collective motions  vibrations and rotations  assuming only the quadrupole type of nuclear deformations .
the number of degrees of freedom associated with both gcm and ibm is five .
two degrees of freedom correspond to the quadrupole deformation of the nucleus in the body - fixed frame ( one can think of an ellipsoid with a fixed volume ) ; their dynamics therefore represents vibrational motions .
the other three degrees of freedom correspond to the orientation of the whole nucleus in the laboratory frame  they describe rotations . in the following we will mostly consider only the non - rotational case , i.e. , fix the angular momentum ( which is an integral of motions ) to zero . in this case , we will be left just with two vibrational degrees of freedom , opening a direct link to other @xmath0 systems famously studied in the context of chaos .
we will see that nuclear vibrations represent a rather rich specimen of chaos in two dimensions .
the gcm , introduced already in 1952 by a.bohr ( for the review see @xcite ) , treats the nucleus as structureless drop described by quadrupole shape variables @xmath1 and rotation euler angles @xmath2 .
these 5 generalized coordinates are associated with a particular parametrization of the rank-2 tensor @xmath3 ( real , symmetric , and traceless in the cartesian form ) , describing the quadrupole deformation of the nucleus in the laboratory frame .
the deformation tensor @xmath3 can be introduced via the expansion of the nuclear radius @xmath4 to spherical harmonics @xmath5 , or alternatively via multipole moments of the nuclear mass ( charge ) distribution @xmath6 . in the first case
, the @xmath7 spherical component of the deformation tensor ( @xmath8,@xmath9,@xmath10 ) is associated with the coefficient ( up to complex conjugation ) at the @xmath11 term of the radius expansion . in the second case ,
the deformation tensor is obtained via the following integral ( up to a scaling factor ) of the mass ( charge ) density : @xmath12 .
any variables describing the shape of the nucleus must be invariant under rotations of the nucleus as a whole .
there are just two independent building blocks for all rotational invariants made of the quadrupole deformation tensor  namely the quadratic and cubic couplings of the quadrupole tensor to zero total angular momentum .
the shape variables @xmath13 and @xmath14 can be defined as a certain parameterization of these two elementary couplings:^{(\lambda)}$ ] of rank-@xmath15 and -@xmath16 spherical tensors @xmath17 and @xmath18 to the resulting spherical tensor of rank @xmath19 ( the definition is analogous to the elementary angular momentum coupling ) .
scalar product of spherical tensors is defined for @xmath20 as @xmath21^{(0)}$ ] . ]
@xmath22^{(0)}=\tfrac{1}{\sqrt{5}}\ \beta^2 \ , \qquad [ [ \alpha^{(2)}\times\alpha^{(2)}]^{(2)}\times\alpha^{(2)}]^{(0)}=-\sqrt{7}\ \beta^3\cos 3\gamma \label{bega } \,.\ ] ] alternatively , @xmath13 and @xmath14 can be derived from the eigenvalues of @xmath3 , which read as @xmath23 $ ] , @xmath24 , @xmath25\bigr)$ ] .
these values are the only components of the deformation tensor in the intrinsic frame ( if ) of the nucleus ( the frame connected with the principal axes of the quadrupole deformation ) , where the deformation tensor is diagonal .
the transformation between the if and laboratory frame ( lf ) is expressed in terms of the euler angles @xmath2 , describing the orientation of the deformed nucleus in lf .    ) with the same choice of parameters . adapted from ref.@xcite . ]
note that @xmath1 can be considered as radial coordinates in the cartesian plane @xmath26 . while the radius @xmath13 quantifies the degree of deformation of the nucleus ( @xmath27 for spherical shape ) , the angle @xmath14 parameterizes the type of the deformed shape .
the plane is divided into six equivalent angular sectors , which follow from an inherent discrete symmetry of the problem , given by the dynamical equivalence of three possible orientations of the deformed shape in the if ( along @xmath28 ) and two deformation types ( prolate , oblate ) .
the @xmath29 , @xmath30 , and @xmath31 axially symmetric prolate or oblates shapes are located at @xmath14 equal to @xmath32 , @xmath33 , and @xmath34 , respectively , while the intermediate @xmath14 values correspond to the related triaxial shapes .
note that these conclusions are consistent with the above - given form of the deformation tensor in the intrinsic frame . to become more familiar with the formalism ( which , as one has to agree , looks rather tough as the first sight )
, the reader is encouraged to consult specialized literature , e.g. , ref.@xcite .
the plane of deformation coordinates with the basic sample of shapes is depicted in fig.[begafi ] .
the momenta associated with the coordinates @xmath13 and @xmath14 are denoted as @xmath35 and @xmath36 , respectively . in the non - rotational case , as noted above , the euler angles @xmath2 are fixed , rendering zero values of the associated rotational momenta .
the resulting vibrational hamiltonian therefore contains just the shape variables @xmath1 and the corresponding momenta @xmath37 .
our analysis is performed mostly with the hamiltonian which contains purely quadratic momentum term of the kinetic energy , and quadratic , cubic and quartic coordinate terms of the potential energy : @xmath38}_{t_{\rm vib}}+\underbrace{a\beta^2+b\beta^3\cos 3\gamma+c\beta^4}_{v_{\rm vib}}\ , .
\label{hgcm}\ ] ] @xmath39 and @xmath40 are adjustable parameters satisfying physical constraints @xmath41 ( kinetic energy is positive ) and @xmath42 ( the potential is confining ) .
the geometric hamiltonian ( [ hgcm ] ) is closely related to the famous hnon - heiles hamiltonian , which was originally introduced in the context of stellar dynamics in galaxies and became a paradigmatic example of chaos in the classical domain @xcite .
the hnon - heiles hamiltonian is a special case of eq.([hgcm ] ) with @xmath43 . in this contribution ( as in ref.@xcite )
we claim that the inclusion of the confining @xmath44 term generates even much richer structures than those obtained in the @xmath43 case .
the coexistence of regular and chaotic types of classical motions generated by the full geometric hamiltonian at high energies is illustrated in fig.[begafi ] .
the two trajectories shown in this figure were calculated for the same model parameters and energy , but represent very dissimilar species of collective motions .
the regular orbit corresponds to a moderate , highly organized vibration , localized in a relatively small region around an axially symmetric equilibrium shape ( above the potential well to which the nucleus with given parameters drops in the low energy limit ) .
in contrast , the chaotic orbit represents some large - amplitude vibrations , which rather disorderly wander through small , medium and large oblate deformations with all three intrinsic orientations
.    a generalization of the simple hamiltonian ( [ hgcm ] ) can be achieved by considering higher - order terms in kinetic or potential energy .
we consider two extensions of the kinetic energy , particularly @xcite : @xmath45 \,,\qquad t_{\rm vib}^{(ii)}=\frac{1}{2m(1+\kappa\beta^2)}\left[p_\beta^2+\left(\frac{p_\gamma}{\beta}\right)^2\right ] \
, , \label{halt}\ ] ] both based on some specific deformation - dependent effective mass ( @xmath46 is an additional model parameter ) .
such or similar extensions of the effective mass are relevant for the description of nuclear rotational bands since they modify the behavior of nuclear moments of inertia . here we consider them to probe the sensitivity of the vibrational measures of chaos to the general form of the hamiltonian kinetic term .
another extension of the above hamiltonian is the inclusion of rotational terms , containing the intrinsic components of angular momentum and the corresponding moments of inertia .
this however leads beyond the 2d description and makes the analysis much more difficult .
we follow this line only in rather a restricted manner  considering only rotations around a fixed axis @xmath31 @xcite .
this can be done by introducing a third collective coordinate , angle @xmath47 , and the associate momentum @xmath48 , which coincides with the angular momentum of the nucleus around @xmath31 in both laboratory and intrinsic frames .
the quadrupole deformation tensor @xmath49 in the laboratory frame is obtained from the diagonal form @xmath50 in the intrinsic frame ( see above ) via rotation by angle @xmath51 around the @xmath31-axis .
this identifies @xmath47 as a variable connected with the if @xmath52 lf transition .
the rotational energy reads as : @xmath53 which has to be added to the vibrational hamiltonian in eq.([hgcm ] ) to get the total hamiltonian @xmath54 .
details can be found in ref.@xcite .
an important feature of the geometric hamiltonian is related to its scaling properties  the fact that dynamics of the system must be invariant under scale transformations affecting the units of relevant quantities , i.e. , energy , coordinate , and momentum .
let us return now to the bare vibrational hamiltonian ( [ hgcm ] ) .
since the scale transformations of this form can be absorbed in parameters @xmath55 , we deduce that some choices of these parameters must be dynamically equivalent .
a more detailed treatment @xcite shows that the quantity @xmath56 is invariant under all scale transformation and represents the only essential parameter of the bare vibrational model in the classical case .
the value of the mass parameter @xmath40 is relevant only in the quantum case , weighting the absolute density of quantum states in the discrete energy spectrum ( for @xmath57 we get the classical limit with infinite state density ) .
this effective parameter reduction implies a crucial simplification of the numerical analysis of the geometric model .
the quantization of the nuclear geometric model is achieved by replacing individual terms in the gcm hamiltonian by the corresponding operators in the hilbert space of wavefunctions @xmath58 .
this task is mathematically rather involved and was completed with a substantial contribution of the patron of this school , marcos moshinsky  see , e.g. , refs.@xcite and many others .
the basic form of the vibrational kinetic energy reads as follows : @xmath59 note that this operator acts in the space where the scalar product @xmath60 is given by an integration over the collective coordinates with a specific measure @xcite .
the expression of the potential @xmath61 is the same as in the classical case , see eq.([hgcm ] ) .
since we are mostly interested in the non - rotational motions , the rotational energy will not be included into the quantized hamiltonian .
examining eq.([tqgcm ] ) one immediately notices the differences from the kinetic energy of a particle in a 2d space , which in polar coordinates @xmath62 reads as @xmath63 .
the absence of euler angles in the description of vibrational motions suggests a possibility to quantize the vibrational energy in the standard 2d manner , i.e. , using the formula : @xmath64 note that this implies a change of the definition of scalar product and the necessity to impose the condition @xmath65 to keep link with the 5d quantization , where this condition is naturally satisfied .
in addition , we choose one of two possible angle - reflection symmetries of the wavefunctions , namely @xmath66 .
we therefore construct three different quantizations of the vibrational collective hamiltonian , all of them yielding the same classical limit : the 2d - even , 2d - odd , and the 5d quantizations ( in the 5d case , the reflection symmetry is fixed ) . from the viewpoint of nuclear theory , only the 5d option is correct , but from the viewpoint of chaos theory ( which is in our focus here ) we have the freedom to probe all three possibilities .
we will use this freedom to check the influence of the quantization scheme on the signatures of quantum chaos . for details , see ref.@xcite .    let us move to the alternative description of nuclear collective motions with the aid of the interacting boson model .
the ibm , introduced in 1970 s by f.iachello and a.arima ( for a review see @xcite ) , does a similar job as the gcm , but in terms of an ensemble of bosons with angular momentum 0 ( @xmath67-bosons ) and 2 ( @xmath68-bosons ) .
these bosons are commonly interpreted as an approximation for nucleon pairs in an even - even nucleus ( their total number @xmath69 is usually fixed at a half of the number of valence nucleons ) , but at the same time they give rise to quadrupole degrees of freedom needed to describe basic collective motions of nuclei .
the @xmath67- and @xmath68-bosons live in a hilbert space associated solely with their spin degree of freedom ( such space has a finite dimension ) and interact with each other through two - body interactions .
for instance , a popular form of the hamiltonian is composed of the @xmath68-boson number operator @xmath70 and a bosonic quadrupole operator @xmath71 as follows : @xmath72)^{(2)}}_{\hat{q}_\chi}\cdot\underbrace{(d^\dag s+s^\dag\tilde{d}+\chi[d^\dag\times\tilde{d}])^{(2)}}_{\hat{q}_\chi } \ , , \label{hibm}\ ] ] where we used @xmath68-boson annihilation operators @xmath73 having proper transformation properties of a rank-2 tensor .
the parameters @xmath74 and @xmath75 , weighting both terms in eq.([hibm ] ) , are more or less arbitrary , while @xmath76 from the quadrupole operator varies within interval @xmath77 for prolate shapes and @xmath78 $ ] for oblate shapes .
note that the @xmath70 term of the above hamiltonian is of the 1-body nature and the term @xmath79 generates specific 2-body interactions .
the link between the bosonic hamiltonian  la ( [ hibm ] ) and the geometric hamiltonian of the type ( [ hgcm ] ) goes via the use of coherent states in the ibm .
for instance , the glauber type of coherent states with an average total number of boson equal to @xmath69 is given by : ) onto a subspace with a sharp boson number , yielding a bosonic condensate state . ]
@xmath80\ket{0 } \label{glau } \,,\ ] ] where @xmath81 stands for a normalization factor to ensure @xmath82 .
the complex variables @xmath83 ( with @xmath84 ) which parameterize the coherent state , have transformation properties of a rank-2 spherical tensor and can be associated with a linear combination of the deformation tensor @xmath3 and the tensor of associated momenta @xmath85 . performing the same gymnastics as in the gcm case ( for details
see ref.@xcite ) , one obtains the five generalizad collective coordinates @xmath86 and the corresponding momenta .
the general bosonic hamiltonian @xmath87 is then translated into a function of these canonical variables by evaluating its expectation value in the respective coherent state @xmath88 . in the @xmath89 limit , this function coincides with the classical hamiltonian @xmath90 associated with @xmath87 .     and @xmath91 , using two truncation dimensions of the hamltonian matrix expressed in the 2d oscillator basis , namely @xmath9239905 and 59080 . panel ( a ) shows a quantum measure of chaos @xmath93 ( where @xmath94 is so - called brody parameter of the nearest - neighbor spacing distribution ; for explanation see below ) calculated for spectra determined by a numerical diagonalization in both truncated bases .
panels ( b ) and ( c ) show the absolute and relative differences between individual level energies , @xmath95 and @xmath96 , respectively .
it is seen that at a certain energy the two calculations yield diverging results .
only the fraction of the spectrum below the vertical dashed line ( corresponding to less than a half of the smaller dimension ) is reproduced correctly . adapted from ref.@xcite . ]    for instance , hamiltonian ( [ hibm ] ) yields the following classical limit @xcite : @xmath97\sqrt{1-\frac{t+\beta^2}{2 } } \nonumber\\ & & \!\!\!\!\!\!\!\!\!\!\!\!\!\!\ !
-\frac{2\chi^2 b}{7}\left[\left(\frac{t+\beta^2}{2}\right)^2-p_\gamma^2\right ] \ , , \qquad\qquad\qquad\qquad\qquad\qquad t = p_\beta^2+\left(\frac{p_\gamma}{\beta}\right)^2 \label{hibmclas } \,.\end{aligned}\ ] ] it needs to be stressed , with regard to the boundedness of the ibm energy spectrum , that the coordinates and momenta appearing in the coherent - state parameterization are restricted only to finite domains , namely : @xmath98 $ ] , @xmath99 $ ] , and @xmath100 $ ] .
thus the square root in eq.([hibmclas ] ) is always well defined .
we may conclude that the ibm classical hamiltonians are closely related to the gcm hamiltonians , discussed above , though by far not identical with them .
the most important difference between both models lies in the kinetic energy terms , which have a much more complex form in the ibm case . on the other hand ,
the potential energy terms , which can be extracted from the hamiltonian by setting all momenta to zero , are rather similar and yield analogous dynamical features .    in conclusion of this section
we say a few words on numerical computations .
classical trajectories for both gcm and ibm are determined with the aid of the fourth - order runge - kutta integration of the classical equations of motions .
some details can be found in ref.@xcite .
quantum spectra are obtained through a numerical diagonalization of the respective quantum hamiltonian .
calculations within the ibm ( which has just a finite dimensional hilbert space ) are naturally done in the complete basis composed of states with different numbers @xmath67- and @xmath68-bosons , coupled to definite angular momenta .
the calculations in gcm ( infinite dimensional space ) has to be performed in a truncated basis .
the results of these must be checked for convergence .
let us streess that due to the restriction to the subset of the quantum spectrum with angular momentum @xmath101 , we perform all gcm calculations in two dimensions .
in general , only a certain fraction of the gcm states obtained by the numerical diagonalization within a truncated 2d basis can be employed in the analyses . as illustrated in fig.[egcm ] , these converged states yield unambiguous and reliable results .
in physics , we like most to deal with integrable systems .
these are the systems which  roughly speaking  have a sufficiently large degree of symmetry to make the dynamical equations solvable .
more precisely @xcite , the full integrability of a classical hamiltonian system with @xmath102 degrees of freedom requires the following three conditions to be simultaneously satisfied : ( a ) there exists @xmath102 integrals of motions @xmath103 ( @xmath104 ) , represented as functions in the @xmath105-dimensional phase space which have zero poisson bracket with the hamiltonian , @xmath106 , ( b ) all these integrals are mutually , i.e. , have vanishing poisson brackets with each other , @xmath107 @xmath108 ( in quantum mechanics , we would say : they are compatible ) , and ( c ) the @xmath105-dimensional gradients @xmath109 are all linearly independent in every point of the phase space .
these conditions guarantee that the hamiltonian ( as well as any other conserved quantity ) can be written as a function of the given integrals , @xmath110 , and that there exists a canonical transformation @xmath111 to some new variables ( so - called action - angle variables ) in terms of which the solution of the hamilton equations becomes rather trivial ( the transformation itself being usually highly nontrivial ) .    for example , the dynamics of the geometric collective model and the interacting boson model is integrable if the corresponding hamiltonian does not depend on variable @xmath14 .
then the momentum @xmath36 becomes an integral of motion , in addition to the hamiltonian @xmath112 itself .
it is easy to verify that all three conditions for integrability in dimension @xmath0 are fulfilled in this situation : ( a ) and ( b ) are obvious , and ( c ) follows from the fact that the gradient of @xmath112 is a nonvanishing normal to the energy surface , while the gradient of @xmath36 is a nonvanishing tangent  therefore , they always define two perpendicular directions .    as we see from eqs.([hgcm ] ) and ( [ hibmclas ] ) , the gcm and ibm hamiltonians are independent of @xmath14 ( hence integrable ) for @xmath113 and @xmath114 , respectively . the @xmath36-dependent integral of motion
is usually written in a less trivial form @xcite : @xmath115 ( where @xmath116 are components of the angular momentum @xmath117 in the if ) , which defines a quantity in the nuclear context called .
note that the rigorous proof of integrability in the full dimension @xmath118 ( including rotational degrees of freedom ) requires a construction of the fifth integral of motions , which is additional to @xmath112 , @xmath119 , @xmath120 , and @xmath121 .
the ibm hamiltonian ( [ hibm ] ) has an extra realization of integrability , situated at the points @xmath122 of the parameter space ( @xmath75 remains an arbitrary scaling parameter ) .
it can be shown that this case of integrability of the ibm hamiltonian is connected with its underlying su(3 ) dynamical symmetry @xcite .
( red and blue dots ) and @xmath123 ( black dot ) . the two cases , differing just by a small displacement of the lighter body , show rather different sets of trajectories
. exponentially diverging distance between both alternative positions of the lighter body is shown in the inset . adapted from ref.@xcite . ]     for 3-body system ) .
panel ( a ) shows evolution of initially close trajectories for times @xmath124 units , in which the deviation is not large .
panel ( b ) displays the evolution for times @xmath125 units , the two final positions represented by the blue and red dots ( with the corresponding shapes drawn aside ) .
note that long - time deviation between the two vibrational orbits ( given in the inset ) reaches a saturation regime , which follows from the boundedness of the system for any finite energy . ]    already from these examples it must be obvious that integrability is a rather rare spice , therefore that by far not all the systems are integrable .
the full realization of this simple fact and its essential consequences came surprisingly late  only in the last decade of the 19th century .
it happened independently and almost simultaneously in the works of h.poincar ' e and a.lyapunov , who can undoubtedly be named the founders of chaos .
poincar ' e s story , in particular , can be used as an exemplary illustration of the random ways in which science makes its progress @xcite .
poincar ' e , who in the preceding years spent a lot of effort to solve the three body problem , submitted an account of his results to a scientific competition , organized in 1885 to celebrate the birthday anniversary of king oscar ii of sweden .
his work won the competition , but an assistant of one of the famous mathematicians in the jury noticed some unclear points in the derivations .
the thing was even much more complex than poincar ' e thought . in 1890
, he submitted a new version of his treatise , in which the mathematical notion of chaos  defined in the sense of instability of motions
was born . in one of the later writings poincar
' e explains : _ _ an example of this behavior is given in fig.[body3 ] .
the new treatise was published at poincar ' e own expenses , exceeding the reward of the competition , which may serve as a warning to those who think that science must always be profitable .
lyapunov came to similar conclusions ( concerning the instability of motions ) in 1892 in his doctoral thesis @xcite .
the description of chaos remained for long a rather exotic domain of abstract mathematics .
in integrable systems , the transformation to action - angle variables @xmath126 ensures that motions in the phase - space , described by hamilton equations @xmath127 and @xmath128 , are localized on manifolds , which are topologically equivalent to tori .
these motions are therefore highly organized .
poincar ' e himself developed a method how the tori can be visualized for systems with @xmath0 .
the method is based on cutting the phase space by a planar section and marking passages of individual trajectories through the section  today we call these pictures poincar maps .
systems results from the requirement that each point of the planar section is crossed by a single trajectory .
indeed , the point in the section is characterized by 2 variables ( usually a coordinate @xmath129 and the associated momentum @xmath130 ) , another variable is fixed by the section constraint ( often of the form @xmath131 ) , and the fourth variable ( e.g. , the momentum @xmath132 ) can  for hamiltonian systems  be computed from energy conservation .
therefore , we have determined all 4 variables in the 4-dimensional phase space . for @xmath133 systems with no integrals of motions additional to @xmath112
the poincar ' e map would have to be more than 2-dimensional , which is rather impractical . ] but what happens when the integrability is destroyed ? certainly , some of the tori must die , but what is the exact scenario and an underlying mechanism ?
plane in the phase space ; details of the hamiltonian can be found in ref.@xcite .
the concentric curves correspond to deformed phase - space tori .
all panels show examples of the poincar ' e - birkhoff fixed - point theorem .
courtesy of m.macek . ]    in 1912 , poincar ' e conjectured a theorem , in which he partly attacked this problem . the theorem was proven in 1913 and finalized in 1925 by g. birkhoff and became known as the poincar ' e - birkhoff fixed - point theorem .
consider an @xmath0 torus characterized by a pair of frequencies @xmath134 with ratio @xmath135 .
if @xmath136 is an irrational number , the trajectory never closes a loop
the motion is only quasi - periodic .
in contrast , if @xmath136 is rational , all trajectories on the given torus are periodic  after @xmath137 crossings of the poincar ' e surface they return to the initial point .
the theorem of poncar ' e and birkhoff is focused on such a rational tori , or more precisely , on the remnants of these tori after the destruction of integrability by a general perturbation .
they showed that each periodic torus leaves an even number of fixed points in the poincar ' e map  a half of these fixed points is stable ( elliptic ) , the other half is unstable ( hyperbolic ) .
this was the first indication that the transition from the full regularity of an integrable system to a complete chaos of a strongly perturbed system is a highly non - trivial ( and usually rather fascinating ! ) process .
figure  [ ibmsec ] shows some examples of poncar ' e maps generated by the ibm hamiltonian with @xmath101 ( hence @xmath0 ) .
alternating islands of elliptic and hyperbolic fixed points are clearly visible .    several decades later , an essential and famous result related to the order - to - chaos transition in classical hamiltonian systems was formulated and proven by a.kolmogorov ( 1954 ) , v.arnold ( 1963 ) , and j.moser ( 1962 ) .
it is nowadays known as the kam theorem @xcite .
what does it tell us ?
it attacks the problem of what exactly happens if an integrable system is subject to an increasingly strong non - integrable perturbation . which of the tori of the original system die , for a given strength of perturbation , and which survive in a deformed form ?
it turned out that the rational tori , such as those considered above , disappear immediately , leaving only the above described remnants . on the other hand ,
some of the irrational tori are rather they may survive ( in a deformed form ) even rather strong perturbations .
the clue for the torus to survive is that it must be .
this superficial statement can be quantified by means of the number theory , through the convergence properties of rational approximations of the number @xmath136 ( the ratio of torus frequencies for an @xmath0 system ) .
the tori characterized by irrational ratios @xmath136 ( slowly converging rational sequences ) are more resistant than the tori with ratios ( fast converging rational sequences ) .
some people say that this law demonstrates the higher immunity of ugliness in the world , but sometimes even the most terrible irrational orbits look really pretty .    ) with @xmath138 .
value of @xmath75 is given in each column .
the maps are constructed for energy @xmath139 of the potential local maximum at @xmath27 ( here , all sections are contracted to a single point with @xmath140 ) .
the upper and lower rows show the @xmath141 and @xmath142 sections , respectively .
each map visualizes 1000 intersections of 100 trajectories with the plane of the section .
adopted from ref.@xcite . ]
several examples of complex poncar ' e maps generated by the gcm vibrational hamiltonian ( [ hgcm ] ) are shown in fig.[gcmsec ] .
particularly the maps in the middle column ( @xmath143 ) show numerous fine structures , which must be caused by rather sophisticated regular orbits .
they appear in the system in spite of a large strength @xmath75 of the hamiltonian term violating the full integrability at @xmath113 . however , a surprising observation is that most of these orbits can not be identified with the survived tori from the @xmath113 system since the maps for lower values of @xmath75 ( c.f . ,
@xmath144 ) contain no trace of such orbits .
one may perhaps think of a of the given hamiltonian to an unknown integrable hamiltonian beyond the parameter space of the gcm .
note that also the @xmath145 maps contain some suspicious fine structures , whose origin can not be explained via the only known integrable regime of the model .
nuclear vibrations  even in their most simplistic description  are highly nontrivial type of motion !    to classify the overall degree of chaos for transitional systems such as those exemplified in fig.[gcmsec ] , one needs to divide the whole phase space into the part filled with regular orbits , and the part filled with chaotic orbits .
since each point of the phase space is crossed by a single trajectory , this division is unique , and one can introduce the degree of classical chaos through the following quantity , below shortened as the : @xmath146 \label{freg } \,.\ ] ] here , @xmath147 stands for the total of the fixed - energy hypersurface @xmath148 in the @xmath105-dimensional phase space , while @xmath149 is the occupied by regular orbits .
the fraction @xmath150 is therefore a quantity between 0 ( for fully chaotic systems ) and 1 ( for completely regular systems ) .
corresponding to the given orbits is shown in panel ( e ) , where we notice a clear difference between regular and chaotic cases .
this is verified more systematically in panel ( d ) , which shows @xmath151 at three different times for 2000 orbits initiating along the line @xmath152 of section ( b ) .
we observe that chaotic domains yield @xmath153 decreasing to negligible values , while in regular domains @xmath153 fluctuates around much higher values . adapted from ref.@xcite . ]
it is clear that any numerical estimate of the regular fraction ( [ freg ] ) for a system between order and chaos must be imperfect .
its precision depends primarily on the sampling density of the phase space by test trajectories .
each of these trajectories is either regular ( stable ) or chaotic ( unstable ) , and the regular fraction can be directly approximated by the fraction @xmath154 .
the convergence of this estimate to the exact value of @xmath150 is guaranteed ( assuming a uniform coverage of the phase space by initial points ) in the limit @xmath155 .
however , a truly hard computational problem is how to unmistakably determine the stability or instability of each individual orbit .
this requires to compute simultaneously the evolution of several orbits @xmath156 $ ] ( both @xmath29 and @xmath157 stand for @xmath102-dimensional vectors ) which initiate in a close vicinity of the given reference orbit @xmath158 $ ] .
if the maximal lyapunov exponent for the deviation @xmath159 of these orbits is positive , the reference orbit is chaotic , if the lyapunov exponent is zero , the orbit is regular . is defined through the long - time phase - space distance @xmath160 between a given orbit and an arbitrary accompanying orbit : @xmath161 .
this yields an exponential growth of the maximal distance for long times : @xmath162 .
fof instance , the lyapunov exponent for the trajectory in fig.[ginst ] , determined however only from a finite time interval of the linear growth of @xmath163 ( where @xmath47 stands for the distance in the coordinate plane ) , is @xmath164 .
] applicability of this asymptotic - time criterion is however hindered in bounded systems , where the deviation of different orbits has an upper limit for any finite energy @xmath165 , yielding therefore all lyapunov exponents formally zero ( cf . the time dependence of the deviation of the vibrational orbits in fig.[ginst ] ) .
several ways to bypass this computational problem has been proposed in the literature .
one of these ways  the one based on so - called alignment indices @xcite  was used in our gcm and ibm calculations . the procedure resorts to evaluating only the relative phase - space deviations of trajectories instead of absolute ones : @xmath166 , where @xmath167 . the relative deviation is clearly a unity vector in the phase space , @xmath168 , and only the mutual orientation of two or more such vectors is meaningful .
let us take a pair of randomly selected orbits , both of them accompanying the reference orbit @xmath169 , and calculate relative deviations @xmath170 and @xmath171 .
if the reference orbit is unstable , the relative deviations tend to align in either parallel or antiparallel directions  along the line of the maximal lyapunov exponent .
this means that the smaller of the lengths @xmath172 ( this quantity is called and denoted below as @xmath153 ) converges to zero for long enough times .
in contrast , for a stable reference orbit @xmath153 keeps oscillating within the definition interval @xmath173 $ ] .
an illustration of the efficiency of this method in the gcm is presented in fig.[sali ] .     for the gcm hamiltonian ( [ hgcm ] ) at @xmath139 as a function of control parameter @xmath74 and/or @xmath75 for @xmath174 .
three computation methods are compared by distinct curves .
the methods are distingished by initiating the trajectories on the @xmath141 or @xmath142 sections ( two curves denoted as @xmath150 ) , or in the whole phase space ( curve @xmath175 ) .
the horizontal axis is split into two parts : the first with @xmath176 and @xmath177 $ ] ( the scale invariant parameter @xmath178 $ ] ) , the second with @xmath179 and @xmath180 $ ] ( hence @xmath181 $ ] ) . adopted from ref.@xcite . ]    our procedure for calculating the regular fraction @xmath182 for a given gcm or ibm vibrational hamiltonian ( @xmath0 ) at fixed energy @xmath165 proceeds in the following steps : ( i ) we choose either the @xmath141 or @xmath142 section of the phase space and calculate the domain in this section that can be reached by trajectories with the given energy .
( ii ) we cover the accessible domain with tiny rectangular cells of the same area ; total number of these cells is @xmath183 .
( iii ) we perform a repeated computation in which an arbitrary point inside any cell of the accessible domain becomes an initial point of a calculated orbit ( note that all 4 canonical variables are fixed by selecting a point in the section and energy ) .
the stability of this orbit is determined by the alignment index technique .
if the orbit is stable , the initial cell as well as all the cells visited by the orbit during the evaluated time interval are counted and the resulting number is incremented to variable @xmath184 ( initially set to zero ) .
it can happen that a single cell is visited by both regular and chaotic orbits ; in that case the value incremented to @xmath184 for this cell is the properly weighted average between 0 and 1 .
( iv ) the calculation is repeated until all cells of the accessible domain are either initial for a trajectory or are visited by at least one of these trajectories .
( v ) we estimate the regular fraction as @xmath185 .
this method has an advantage of a systematic coverage of the phase space  or more precisely , its 2d section  with trajectories . on the other hand ,
it is clear that the resulting value of the regular fraction may depend on the phase - space section selected .
the efficiency of our procedure and its inherent ambiguities in the gcm are illustrated by fig.[freg0 ] , where we compare the values of @xmath182 ( for @xmath139 and variable control parameter @xmath75 and/or @xmath74 ) calculated for both @xmath141 and @xmath142 sections , as well as by directly counting the stable and unstable orbits generated without reference to any phase - space section .
it is seen that all three methods give approximately the same results . the differences between the three curves in fig.[freg0 ]
can be declared to be an internal uncertainty of the method at the selected level of resolution ( number @xmath186 ) .
figure [ freg0 ] , apart from its methodological aspect , has also a very interesting physical content .
we see that for moderate values of parameter @xmath75 the regular fraction quickly drops from value @xmath187 at @xmath113 to @xmath188 at @xmath189 . in view of the kam scenario , this is a perfectly expectable behavior since @xmath75 weights the nonintegrable perturbation in the hamiltonian ( [ hgcm ] ) .
however , at @xmath190 the regularity starts growing again !
in particular , in the subsequent region up to @xmath191 we observe a series of fine - structured peaks in the @xmath150 dependence . or @xmath192 is limited to energies @xmath193 and results from a specific shape of the potential . ] the increased regularity in this parameter domain
is connected with rather complex poincar maps , such as those shown in the middle column of fig.[gcmsec ] .
quite surprisingly , the peaks in fig.[freg0 ] follow a certain if expressed in scale invariant parameter @xmath194 ( the four local peaks appear at integer multiples of the value @xmath194 corresponding to the first peak ) .
the origin of this behavior is still unknown .
the complete dependence of @xmath150 on the gcm hamiltonian parameters and energy can be found in ref.@xcite .     for gcm hamiltonians with extended kinetic term @xmath195 ( panel a ) and @xmath196 ( panel b ) , see eq.([halt ] ) .
dependences for several values of parameter @xmath46 are given for both types of kinetic term .
the other parameters are the same as in fig.[freg0 ] ( which corresponds to the @xmath197 case in both panels ) . adapted from ref.@xcite . ]    in fig.[freg0ext ] , we show the regularity plots for generalized gcm hamiltonians with the kinetic energy terms according to eq.([halt ] ) .
one can point out the opposite prevailing influence of type i and ii extensions of the kinetic energy on the regular fraction : while the type i extension with growing @xmath46 tends to reduce @xmath150 , the type ii extension seems to rather increase it .
note that the influence of the rotational term ( [ trot ] ) on the regularity plot was analyzed in ref.@xcite .
let us note that similar structures as in figs.[freg0 ] and [ freg0ext ] are observed also in the interacting boson model , where people often discuss about a so - called @xcite .
it needs to be stressed that the general ibm hamiltonian is more sophisticated than the gcm one and therefore can not be characterized by a single scale invariant parameter such as @xmath194 .
in particular , the common parameterization ( [ hibm ] ) yields two essential parameters , namely @xmath76 and the proportion between @xmath74 and @xmath75 , so that the peak series from fig.[freg0 ] becomes a kind of in the regularity landscape .
the ridge follows an arc shape in the parameter plane and , as in the gcm , it is most pronounced for energies just above the @xmath27 local maximum of the potential .
the physics behind both phenomena is probably rather similar .
there is a fundamental problem to define chaos in quantum mechanics since linearity of the theory excludes any chance to introduce an exponential sensitivity to initial conditions .
consider two normalized initial states @xmath198 and @xmath199 differing by @xmath200 with @xmath201 and @xmath202 , @xmath203 .
the evolution of both states reads : @xmath204 where all scalar products remain conserved due to unitarity of the evolution operator . as a consequence , the squared distance of both states , @xmath205 , is independent of time .
quantum hilbert space does not permit any kind of butterfly - wing effect .
a different type of instability , the one which remains relevant even in the quantum domain , was pointed out by a.peres in 1984 @xcite . instead of the dependence on initial conditions he proposed to consider the dependence of quantum dynamics on details of the hamiltonian ( tiny changes of interaction parameters or small external perturbations ) .
although the resulting concept of so - called plays an important role in many branches of quantum physics , including quantum information applications , it does not provide
unfortunately!a selfsufficient definition of chaos in general quantum systems @xcite .
according to present understanding , chaos is firmly anchored in classical physics .
quantum physics can only explore a multitude of its distinct consequences ( signatures ) on the quantum level of description .
quantum chaos in this sense is therefore not a , but rather a research field located somewhere in the highlands along the border between quantum and classical physics .
m.berry , who undertook some pioneering expeditions to these wild territories , speaks about instead of quantum chaos @xcite .    for bounded quantum systems
 those with discrete energy levels  the most important signatures of chaos were identified in statistical properties of energy spectra .
the relevant mathematics was created already in the 1950 s , long before the quest for quantum chaos started @xcite .
it was when e.wigner looked for an achievable type of physical description related to long sequences of neutron resonances  plentiful and rather complex excited states of compound atomic nuclei .
wigner decided to give up the exact description  prediction of the properties of all individual resonances  and proposed to strictly separate smooth ( predictable ) and fluctuating ( statistical ) parts of the corresponding quantities .
while the smoothened dependences ( averages , which reflect only sort of properties of the system ) can be described in terms of some reasonable simplifying models , the fluctuations ( that depend on all details of the dynamics ) represents just random , unpredictable component of the problem . in the 1980 s it became clear that quantum signatures of chaos lie right in these fluctuating properties @xcite .    as an example , consider a discrete energy spectrum described by the state density @xmath206 .
to predict the position of every individual level @xmath207 might be too hard ( or impossible  especially if the system is as complex as a nucleus .
however , one can perhaps estimate a smoothened density @xmath208 , where @xmath209 stands for a suitable smoothening function ( e.g. , the gaussian of width @xmath210 ) .
indeed , the smoothed density is closely related to the size @xmath211 of the energy shell in the multidimensional phase space in units associated with the planck constant , which , in principle , can be evaluated with the aid of suitable techniques ( like the fermi gas approximation in the case of a nucleus ) .
in contrast , the remaining component of the state density , its fluctuating part @xmath212 , depends on tiny details of the hamiltonian and can not be predicted by simple models .
nevertheless , statistical features of @xmath212 reflect the degree of chaos in the classical limit of the system .    to extract the chaos - related information , one usually performs a transformation of the spectrum called @xcite .
it can be written in the following way : @xmath213 the unfolded sequence @xmath214 has a constant average density ( equal to one ) and therefore materializes statistical properties of the spectrum in the most explicit form .
it turns out that mutual correlations between the members of this sequence show some universal behavior , which depends on whether the classical counterpart of the system is regular or chaotic .
surprisingly , in regular ( integrable ) systems the unfolded sequences of energies with the same quantum numbers ( such as angular momentum or parity ) exhibit only weak or virtually no correlations , making the visual appearance of spectra rather disordered .
in contrast , chaotic systems yield strong correlations and produce seemingly more regular spectra .
( classical regular fraction ) , @xmath93 ( the brody parameter adjunct ) , and @xmath215 ( a shifted exponent of the power spectrum ) for the gcm .
panels ( a ) and ( b ) show results for @xmath216 and for @xmath217 , respectively .
the mass parameter was chosen such that the displayed energy intervals contain 30 ( panel  a ) and 40 ( panel  b ) thousand of quantum levels , calculated with the 2d quantization option ( see the text ) .
panels ( c ) and ( d ) show separately @xmath93 and @xmath215 for all three different quantization schemes ( the hamiltonian parameters are the same as in the respective left panel ) .
we conclude that  within the error bars ( given by the fitting methods of the respective quantities in finite samples of levels in a vicinity each energy)the correlation measures show no dependence on the quantization . ]
the correlations in unfolded spectra of chaotic systems are of both short- and long - range types .
the short range correlations are usually expressed through the statistical distribution of nearest neighbor spacing @xmath218 constructed for sequences of levels with the same quantum numbers .
this distribution reveals whether the system exhibits the so - called , a phenomenon indicating the absence of additional integrals of motions .
we use the following formulae @xcite : @xmath219 where the brody form with @xmath220 interpolates between the limiting cases @xmath221 ( poisson ) and @xmath222 ( wigner ) .
note that the symbol @xmath223 stands for the normalization factor and @xmath224 for a coefficient ensuring that @xmath225 .
the poissonian distribution , which describes an uncorrelated spectrum , is typical for integrable systems .
it implies the absence of level repulsion ( since @xmath226 for @xmath227 ) , which is related to the existence of additional integrals of motions that eliminate the corresponding offdiagonal matrix elements of the hamiltonian .
the wigner formula includes the level repulsion given by @xmath228 for small spacings and captures properties of chaotic systems ( with generally nonvanishing offdiagonal hamiltonian matrix elements ) . )
is only an approximation of the distribution predicted by the random matrix theory .
it applies in system symmetric under the time reversal . for other systems ,
the formula is modified to include even stronger ( @xmath229 ) level repulsion @xcite . ] system in between order and chaos are characterized by a transitional brody form with the parameter @xmath94 somehow related to the overall degree of chaos . to be more specific
, we expect a correlation ( generally nonlinear ) between values @xmath230 .
associated with the gcm unfolded energy spectrum at a single energy @xmath165 .
model parameters : @xmath231 , @xmath232 .
panel ( a ) displays the number @xmath137 of levels in the energy interval @xmath233 corresponding to the shortest periodic orbit ( period @xmath234 ) as a function of energy @xmath165 .
the value @xmath235 sets the minimal frequency @xmath236 for the evaluation of the power spectrum .
panel ( b ) shows in the log - log scale the power spectrum @xmath237 for a sample of 2048 levels in a vicinity of the energy indicated in panel ( a ) by the vertical line .
the slope of the cloud of points in the above @xmath238 ( vertical line in panel  b ) yields the exponent in the @xmath239 dependence .
the power spectrum is determined by two different methods : ( i ) by considering all 2048 levels ( black dots ) and ( ii ) by grouping the levels into 32 bunches of 64 levels and calculating an average power spectrum over the bunches ( green dots ) .
the corresponding fits of the slope , indicated by the full ( i ) and dashed ( ii ) tilted lines , yield values @xmath240 and @xmath241 . note that in fig.[qchagcm ] we used the method ( ii ) . ]
long range correlations in spectra of chaotic systems are also very strong
. they can be expressed in different ways , for instance via the variance of the number of unfolded levels in an interval of length @xmath242 , or equivalently via a so - called @xmath243 statistics @xcite .
an alternative approach @xcite is based on the fourier analysis of a sequence of statistical quantities @xmath244 .
the sequence apparently satisfies @xmath245 , expressing fluctuations of the cumulative number @xmath246 of unfolded levels around the straight line @xmath247 .
the fourier transformed sequence reads as @xmath248 $ ] , with @xmath40 being the number of levels included into the analysis , and yields a power spectrum @xmath249 .
the power spectrum contains information on the frequencies @xmath102 present in the @xmath250 ( number @xmath137 being identified with a discrete time ) and represents a statistical quantity associated with @xmath250 .
it turns out that the average of the power spectrum can be well approximated by the formula @xmath239 , where @xmath251 for strongly correlated spectra of chaotic systems , and @xmath252 for the uncorrelated poissonian spectra of regular systems @xcite .
hence we establish a link @xmath253 .
statistical properties of spectra were extensively studied for both the geometric collective model and the interacting boson model @xcite .
figure  [ qchagcm ] shows some of the gcm results , namely a comparison of classical regular fraction with the corresponding measures of short - range and long - range spectral correlations for two choices of the hamiltonian parameters .
a clear correspondence between @xmath150 and the dependences of @xmath254 and @xmath255 is observed , which confirms the above given statement concerning the link between chaos , on the classical level , and statistical properties of spectra , on the quantum level .
this , as some authors call it , has been already tested in numerous systems , but mostly in the regime of a weak energy dependence of the regular fraction @xcite . in particular , the most influential analyses have been performed with 2d billiard systems , which show a constant value of @xmath150 , totally independent of @xmath165 . here
we verify the bohigas conjecture in a more complex form ( though still based on a 2d model ) , demonstrating that even rapid changes of the classical regularity are closely followed by both short- and long - range quantum correlation measures .
moreover , the comparison is performed for different quantization schemes and shows that the correlation properties are quantization - independent .
these properties must therefore reflect only the features related to the classical limit of the system , in agreement with the bohigas conjecture .
we want to emphasize that the rapid variability of the gcm regular fraction with energy not only creates a desirable framework for testing the quantum chaos assumptions , but also represents a considerable computational challenge ! the statistical analysis of the spectrum can only be performed piece by piece  on samples of levels situated in a close vicinity of energy @xmath165 , where @xmath165 varies along the whole spectrum . on the one hand ,
the samples have to be sufficiently large to allow for a reliable statistical evaluation , on the other hand , the energy interval @xmath95 spanned by each sample needs to be sufficiently narrow with respect to the scale that determines local variations of @xmath182 . as an example , we illustrate in fig.[onef ] the fourier analysis of a single sample of levels around one particular energy @xmath165 , outlining the procedures capable to determine the local power - spectrum exponent @xmath256 .
in the rest of this section we describe an interesting visual approach to quantum chaos , which will complement the cute poincar ' e maps presented above in connection with classical chaos .
the method was proposed by a.peres @xcite ( also in 1984 , which seems to be _ annus mirabilis _ for quantum chaos ) and its practical application is limited to systems with two degrees of freedom .
it is based on the fact that the time average @xmath257 of an arbitrary quantity @xmath258 represents an automatic integral of motion , regardless of any specific dynamical properties of the system .
classically , the average is defined through infinite - time integration over the whole classical orbit , and therefore assigns a unique value @xmath259 to any point of the phase space ( the same value is assigned to all points visited by a single trajectory ) .
the integrals of motions @xmath257 do not generally satisfy the conditions for integrability ( any system has an infinite number of such integrals ) , but , as shown below , they provide a valuable analyzing tool for systems between order and chaos .    ) with @xmath260 [ the o(6 ) dynamical symmetry limit ] for @xmath261 bosons .
symbol @xmath194 stands for the o(5 ) seniority , while the other symbols denote casimir invariants of the su(3 ) and u(5 ) algebras , which are involved in the overall u(6 ) dynamical algebra of the ibm .
the lower row shows the lattices @xmath262 of mean values of the respective quantities @xmath258 , while the upper row shows the corresponding lattices of the variances @xmath263 .
courtesy of m. macek . ]     in the 5d quantization with @xmath264 ( upper rows ) and @xmath265 ( lower rows ) .
parameters @xmath74 and @xmath75 of the potential are indicated above the panels , @xmath174 , and @xmath266 . ]    in the quantum case , the time average is defined by the following formula : @xmath267\,\hat{p}\,\exp[-\tfrac{i}{\hbar}\hat{h}t ] \,,\qquad [ \hat{\bar{p}},\hat{h}]=0 \label{per } \,,\ ] ] where @xmath268 is an operator corresponding to the original observable .
it is easy to verify that the commutator of the time - average operator @xmath269 with the hamiltonian @xmath87 vanishes in the @xmath270 limit , irrespective of the concrete choice of @xmath268 .
quantum expectation values of both observables @xmath268 and @xmath269 are the same for the energy eigenstates of the system , @xmath271 , and it is therefore rather straightforward to create a 2d picture in which the discrete energies @xmath272 are plotted against the corresponding expectation values @xmath273 . for the resulting lattice of points we use the term .
it can be constructed for any observable @xmath268 .
peres lattices for systems with @xmath0 can be seen as quantum counterparts of the classical poincar ' e maps .
it turns out that geometric regularity ( irregularity ) of the lattice signals regular ( chaotic ) dynamics of the system on the classical level .
in particular , if the system is classically integrable , any integral of motion must be a function ( suppose a function ) of the canonical actions , hence @xmath274 and @xmath275 . in transition to the quantum regime , @xmath276 and @xmath277
become approximate multiples of the planck action and form a rectangular grid in the @xmath278 plane .
as the peres lattice in the @xmath279 plane is related to this @xmath278 grid by a sooth mapping , it will certainly be somehow distorted , but one expects that it will still look orderly . for non - integrable systems , on the other hand , there is no reason for the peres lattice to be regular . for given energy @xmath165
tend to be contracted to a single value due to ergodicity of the system . ] in this case , any partial arrangement of the lattice is a signature of fractional regularity of the dynamics . it should be noted that various peres lattices can in general be combined with each other or smoothly transformed without affecting regularity or chaoticity of the patterns present in the lattice .
for instance , the lattice of @xmath272 versus @xmath263 , which displays the variance of the quantity @xmath258 in the hamiltonian eigenstates , is just a combination of lattices for @xmath280 and @xmath258 .
some examples of peres lattices for the above ibm and gcm hamiltonians are given in figs.[peribm ] and [ pergcm ] .
the first of them , fig.[peribm ] , shows several incompatible lattices calculated for an integrable ibm hamiltonian with the o(6 ) dynamical symmetry .
the o(5 ) seniority , cf .
eq.([senior ] ) , whose average and variance lattices are drawn in the first column , represents a conserved quantity for any ibm hamiltonian ( [ hibm ] ) with @xmath114 , as can be directly verified through the null variance of the seniority in all hamiltonian eigenstates .
the other quantities are not integrals of motions , they nevertheless generate perfectly regular lattices in both average values and dispersions .
this is in agreement with integrability of the hamiltonian employed .
figure  [ pergcm ] depicts several peres lattices for the gcm , using eigenstates with @xmath101 in the 5d quantization .
the lattices were constructed using only the average values of the respective peres operator , the corresponding variances are not given this time .
two choices of the peres operator are presented : the first one is the @xmath101 seniority @xmath281 , which expresses an operator @xmath282 in the plane @xmath283 of the deformation coordinates ( in the 5d quantization ) .
the corresponding averages are denoted as @xmath284 and their lattices are shown in the first row of each panel .
the second choice of the peres operator coincides with the cubic term of the gcm potential , @xmath265 , which represents the nonintegrable perturbation of the @xmath113 integrable hamiltonian .
the averages , denoted @xmath285 , are seen in the second row of each panel .
peres lattices in fig.[pergcm ] are given for several points in the parameter plane @xmath286 , following a path which crosses all inequivalent gcm potentials differing by the scale invariant parameter @xmath194 .
at first , we fix @xmath176 and change @xmath287 , visualizing the onset of chaotic vibrations in deformed nuclei and their partial regularization at larger values of @xmath75 .
next we keep @xmath179 and change @xmath288 . in this way , we make a transition between two qualitatively different types of the potential  the one corresponding to deformed equilibrium shapes of the nucleus and the one associated with spherical nuclei ( the @xmath289 potentials have three degenerate minima at @xmath290 , while those with @xmath74 above a certain critical value @xmath291 have just a single global minimum at @xmath292 ) . and @xmath290 minima of the potential become degenerate and swap is at @xmath293 . ] finally , we fix @xmath294 and decrease @xmath295 , coming back to the integrable regime , but on the spherical side of the parameter plane .
the most chaotic lattices are the ones in panels ( c ) and ( e ) .
the first corresponds to @xmath296 , which coincides with the first minimum of @xmath150 at @xmath139 in fig.[freg0 ] .
both peres lattices in this panel can be compared with the energy dependences of @xmath182 and with the other measures of chaos shown in panels ( a ) and ( c ) of fig.[qchagcm ] .
we may notice that the regions with increased regularity in fig.[qchagcm ] correspond to the occurrence of regular parts of the lattice in fig.[pergcm](c ) .
the chaotic lattice in panel ( e ) corresponds to @xmath297 , which is a point situated near the between spherical and deformed shapes .
many other examples of peres lattices and their detailed comparison with the dependence @xmath182 can be found in refs.@xcite and @xcite .    from the examples shown here and in the literature it can be inferred that peres lattices significantly complement the standard methods used to investigate chaotic properties of quantum systems .
in fact , they even enable one to go beyond the standard methods .
this is so because quantum states can be _ individually assigned _ to distinctly regular or distinctly chaotic parts of the lattice , which is in contrast to the statistical approach that can only attribute some averaged regular / chaotic attributes to sufficiently large ensembles of states .
constructing a suitable peres lattice for a given quantum system with known spectral correlations is therefore similar to analyzing the classical dynamics with a given @xmath150 by means of poincar ' e sections .
two remarks on the applicability of the method are in order : first , to unambiguously deduce whether a given state is in a regular or chaotic domain of the spectrum , it is usually useful to construct peres lattices for several operators .
a comparison of several lattices helps particularly in situations when an overlap of two or more regular parts in a single lattice produces a seemingly chaotic mesh .
second , the method in the above - presented form is applicable only in systems with @xmath0 .
consider , for example , an integrable system with @xmath298 . to ensure a one - to - one mapping from the space of canonical actions @xmath299
, one would have to consider a @xmath300-dimensional lattice ( @xmath102-dimensional in general ) .
this is , however , rather difficult with only 2d physical literature being available . these unalterable limitations of us , poor 3d creatures , will probably always hinder our understanding of higher than 2d systems .
we are arriving to the end of our short tour .
its purpose was to show some known concepts of chaos theory by means of simple models taken from the nuclear context .
was the marriage of chaos with nuclei successful ?
we believe it was .
we think that the models like those used here provide an excellent environment for general investigations of chaos in both classical and quantum incarnations  environment , which is sufficiently complex to embrace diverse nontrivial phenomena , but at the same time simple enough to allow an intelligible description . in analogy with the billiard models , that have been extensively used for such investigations in the past decades , the nuclear models can be used as models with two degrees of freedom .
in contrast to the billiard models , however , they exhibit very rich dependence of chaos measures on both energy and hamiltonian parameters , which opens some new directions of potential research .
* geometrical chaos : attempts have been made to develop techniques that would be able to _ predict _ an approximate value of the regular fraction @xmath150 for an arbitrary choice of the system s parameters and energy . an interesting approach is based on the geometrical method , which reduces the dynamics of the system to free motions of particles in curved space @xcite . a related technique resorts to evaluating the convexity / concavity of the energy contours of the model potential @xcite .
it can be shown that in the gcm ( and/or ibm ) the geometrical method is only partially successful @xcite .
pre - estimation of the degree of chaos in a given system therefore remains one of the major tasks of chaos theory .
* chaos versus symmetry : people keep trying to establish a link of order and chaos with some generalized concepts of symmetry . if such a link exists and is universal , it would interpret partial regularity , as observed in various systems of nature , through some imperfect , fractional dynamical symmetries possessed by these systems .
this led to the ideas of so - called partial dynamical symmetry @xcite and quasi dynamical symmetry @xcite .
although both gcm and ibm are particularly well suited for this kind of research @xcite , its success in connection with chaos is so far only partial .
* implications of chaos : one can step from origins to consequences of chaos , investigating various chaos - assisted phenomena in specific systems . for example , the coexistence of regular and chaotic modes of dynamics turned out to be relevant for transitions between spherical and deformed shapes of nuclei @xcite .
regular dynamics was also shown to play a remarkable role in the emergence of rotational bands in highly excited nuclei @xcite , which suggests that the physics underlying the general phenomenon of adiabatic separation of intrinsic and collective motions in many - body systems might be closely related to chaos theory .
m. gutzwiller , _ chaos in classical and quantum mechanics _ , springer , new york , 1990 .
weidenmller and g.e .
mitchell , _ rev .
phys . _ * 81 * , 539 ( 2009 ) .
mitchell , a. richter , and h.a .
weidenmller , _ rev .
phys . _ * 82 * , 2845 ( 2010 ) .
rowe and j.l .
wood , _ fundamentals of nuclear models .
foundational models _ , world scientific , singapore , 2010 .
p. cejnar , p. str ' ansk ' y , and m. macek , _ nuclear physics news _ * 21 * , 22 ( 2011 ) .
m. h ' enon and c. heiles , _ astron .
j. _ * 69 * , 73 ( 1964 ) .
p. cejnar and p. str ' ansk ' y , _ phys .
lett . _ * 93 * , 102502 ( 2004 ) .
p. str ' ansk ' y , _ ph.d .
thesis _ , charles university , prague , 2009 .
p. str ' ansk ' y , m. kurian , and p. cejnar , _ phys .
_ c * 74 * , 014306 ( 2006 ) .
e. chacn , m. moshinsky , and r.t .
* 17 * , 668 ( 1976 ) .
e. chac ' on and m. moshinsky , _ j. math
. phys . _ * 18 * , 870 ( 1977 ) .
str ' ansk ' y , p. hruka , and p. cejnar , _ phys .
_ e * 79 * , 04620 ( 2009 ) .
str ' ansk ' y , p. hruka , and p. cejnar , _ phys .
_ e * 79 * , 066201 ( 2009 ) .
f. iachello and a. arima , _ the interacting boson model _ , cambridge university press , cambridge , 1987 .
r. l. hatch and s. levit , _ phys .
_ c * 25 * , 614 ( 1982 ) .
m. macek , p. str ' ansk ' y , p. cejnar , s. heinze , j. jolie , and j. dobe , _ phys .
_ c * 75 * , 064318 ( 2007 ) .
arnold , _ mathematical methods of classical mechanics _ , springer , new york , 1989 . j. barrow - green , _ poincar ' e and the three body problem _ ( history of mathematics , v. 11 ) , am .
a.m. lyapunov , _ the general problem of the stability of motion _ ,
taylor & francis , london , 1992 .
m. macek and a. leviatan , _ nuclear structure and dynamics 12 _ , ed .
r. bijker , _ aip conf .
proc . _ * 1488 * , 441 ( 2012 ) .
skokos , _ j. phys . _ a * 34 * , 10029 ( 2001 )
skokos , ch .
antonopoulos , t.c .
bountis , and m.v .
vrahatis , _ j. phys . _ a * 37 * , 6269 ( 2004 ) . y. alhassid and n. whelan , _ phys
_ * 67 * , 816 ( 1991 ) .
n. whelan and y. alhassid , _ nucl .
phys . _ a * 556 * , 42 ( 1993 ) .
a. peres , _ phys
_ a * 30 * , 1610 ( 1984 ) .
t. gorin , t. prosen , t.h .
seligman , and m. nidari , _ phys .
_ * 435 * , 33 ( 2006 ) . m. berry , _ physica scripta _ * 40 * , 335 ( 1989 ) . c.e .
porter ( ed . ) , _ statistical theories of spectra : fluctuations _ , academic press , new york , 1965 .
o. bohigas , m.j .
giannoni , and c. schmit , _ phys .
_ * 52 * , 1 ( 1984 ) .
t. guhr , a. mller - groeling , and h.a .
weidenmller , _ phys .
rep . _ * 299 * , 189 ( 1998 ) .
a. rela~ no , j.m.g .
g ' omez , r. a. molina , j. retamosa , and e. faleiro , _ phys .
_ * 89 * , 244102 ( 2002 ) .
a. peres , _ phys
_ * 53 * , 1711 ( 1984 ) .
l. horwitz , y.b .
zion , m. lewkowitz , m. schiffer , and j. levitan , _ phys .
lett . _ * 98 * , 234301 ( 2007 ) .
str ' ansk ' y _ et al .
_ , to be published .
a. leviatan , _ prog .
phys . _ * 66 * , 93 ( 2011 ) .
rowe , p. rochford , and j. repka , _ j. math .
_ * 29 * , 572 ( 1988 ) .
m. macek , j. dobe , and p. cejnar , _ phys .
_ c * 80 * , 014319 ( 2009 ) .
m. macek , j. dobe , and p. cejnar , _ phys .
_ c * 82 * , 014308 ( 2010 ) .
berezovoj , yu.l .
bolotin , and v.a .
cherkaskiy , _ phys .
lett . _ a * 323 * , 218 ( 2004 ) .
a. leviatan and m. macek , _ phys .
_ b * 714 * , 110 ( 2012 ) .
m. macek , j. dobe , p. str ' ansk ' y , and p. cejnar , _ phys .
* 105 * , 072503 ( 2010 ) .