I

II

III

! " # $% $# # %

& $ ! " # $% $#

& ' ( !

! % # '

# % ! )* %

+ # !!

, %+ ' #

! # #

'# % - ! " # $

% $# % & ' !

.!! ! %

!% , . /, "

# $% $# 0 $ (

!% , !# % !

'.# !% % # /#

% 1 !!, #

'# 2 & !

+ !% %

% ! % #

" '

% % # 2 ' ULB -% ! 0 '

3# &

% ! ' $% -.#

(ULB) 0! 4% ! % ! + ! & 3# & $#

-# '% %/ ULB

* 2005 '% - % % !# ! '

$% -.# Pierre Gaspard % # ( !% % !% % %+ #

2# % # % % ,

+ ! % '

! (56 ,

2# ! #

!# ANIC 4 ' $% -.#

IV

'# % ! "'* !!% 7 3 ! %

! 0 3 ,!#

$, 3% '

% % ! , % %

! (8 7'*"9& ' $! ! '% '!!

** :'$'' **; '

V

VI

1. Chaos in a Near Integrable Hamiltonian Lattice ,Rothos V. M, Antonopoulos Ch. & Drossos L., 2002, Int. J. Bif. Chaos, 12, 8, 1743 1754

1. Geometrical Properties of Local Dynamics in Hamiltonian Systems: the GeneralizedAlignment Index (GALI) Method ,Skokos Ch., Bountis T. & Antonopoulos Ch., submitted at February 2006

2. Stability of Simple Periodic Orbits and Chaos in a Fermi Pasta Ulam Lattice ,Antonopoulos Ch. & Bountis T., 2006, Physical Review E, 73, 056206

3. Chaotic Dynamics of N degree of Freedom Hamiltonian Systems ,Antonopoulos C., Bountis T. & Skokos C., 2006, Int. J. Bif. Chaos, 16, 6, 1777 1793

4. Detecting Order and Chaos in Hamiltonian Systems by the SALI Method ,Skokos Ch., Antonopoulos Ch., Bountis T. C. & Vrahatis M. N., 2004, J. Phys. A,37, 6269 6284

5. How does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? ,Skokos Ch., Antonopoulos Ch., Bountis T. C. & Vrahatis M. N., 2003, Prog. Theor.Phys. Supp., 150, 439 443

1. SALI: An Efficient Indicator of Chaos with Application to 2 and 3 Degrees of FreedomHamiltonian Systems ,Antonopoulos Ch., Manos A. & Skokos Ch., 2004. In Proceedings of the 1st Interna-tional Conference From Scientific Computing to Computational Engineering , 1st IC SCCE, ed. Tsahalis D. T., Patras Univ. Press, Vol. III, 928 933

2. Smaller Alignment Index (SALI): Determining the Ordered or Chaotic Nature of Orbitsin Conservative Dynamical Systems ,Skokos Ch., Antonopoulos Ch., Bountis T. C. & Vrahatis M. N., 2003. In Proceedingsof the Conference Libration Point Orbits and Applications , eds. Gomez G., Lo M. W.& Masdemont J. J., World Scientific, 653 664

VII

1. Estimation of the Dynamic Aperture of Symplectic Mappings Using Evolutionary Al-gorithms ,Ch. G. Antonopoulos, I. G. Petalas, T. C. Bountis & M. N. Vrahatis, 2006. InProceedings of Volos 18th Summer School and Conference , ed. T. C. Bountis.

2. A Fast and Reliable Method for Distinguishing Regular from Chaotic Motion in Hamil-tonian Systems ,Antonopoulos Ch., Skokos Ch., Bountis T. C. & Vrahatis M. N., 2003. In RecentAdvantages in Mechanics and Related Fields , Volume in Honour of Prof. Goudas C.L., Univ. of Patras, 325 333

3. Smaller Alignment Index (SALI): Detecting Order and Chaos in Conservative Dynam-ical Systems ,Skokos Ch., Antonopoulos Ch., Bountis T. C. & Vrahatis M. N., 2002. In Proceed-ings of the 4th GRACM Congress on Computational Mechanics , ed. Tsahalis D. T.,Vol. IV, 1496 1502

1. Detecting Irregular Orbits in Gravitational NBody Simulations,Gemmeke, J. F. & Portegies Zwart, S. F. & Kruip, C. J. H., 2006, Communications inNonlinear Science and Numerical Simulation (2006), doi:10.1016/j.cnsns.2006.09.012

2. The Resonant Structure of Jupiters Trojan AsteroidsI. LongTerm Stability and Dif-fusion,Robutel, P. & Gabern, F., 2006, Monthly Notices of the Royal Astronomical Society,Vol. 372, Issue 4, pp. 1463 1482

3. SelfConsistent Models of Cuspy Triaxial Galaxies with Dark Matter Haloes,CapuzzoDolcetta R. & Leccese L. & Merritt D. & Vicari A., 2006, arXiv:astroph/0611205,pdf version: http://de.arxiv.org/abs/astroph/0611205, Preprint submitted to ApJ, As-trophysics

4. Lyapunov Indices With Two Nearby Trajectories in a Curved Spacetime,Wu X. & Huang T. Y. & Zhang H., 2006, Physical Review D Particles, Fields,Gravitation and Cosmology, Vol. 74, Issue 8, Article number 083001

5. The Resonant Structure of Jupiters Trojan Asteroids I: Long Term Stability and Dif-fusion,Robutel P. and Gabern F., 2006, MNRAS, In press

VIII

6. Quasiperiodic and Chaotic Discrete Breathers in a Parametrically Driven System With-out Linear Dispersion,Maniadis P. & Bountis T., 2006, Physical Review E, Vol. 73, pp. 046211

7. Detecting Irregular Orbits in Gravitational N Body Simulations,Gemmeke J. F., Portegies Zwart S. F. & Kruip C. J. H., 2006, arXiv:astroph/0607343 v114 July 2006, pdf version: http://arxiv.org/abs/astroph/0607343, Preprint submittedto Communications in Nonlinear Science and Numerical Simulation

8. Sensitivity tools vs. Poincare sections,Barrio R., 2005, Chaos Solitons & Fractals, Vol. 25, No. 3, 711 726

9. On the Integrability and Chaos of an N = 2 Maxwell Chern Simons HiggsMechanical Model,de Assis L. P. G., Helayel Neto, Haas F. & Nogueira A. L. M. A., 2005, TechnicalReport no: CBPF NF 013/05, Centro Brasileiro de Pesquisas Fisicas CCP/CBPF,Rio de Janeiro, Brasil, May 2005, arXiv:hepth/0505159 v1 18 May 2005, pdf version:http://arxiv.org/abs/hepth/0505159

10. Several Diagnostic Indexes for Orbital Chaos,Wu X. & Huang T. Y., 2005, Progress in Astronomy, Vol. 23, No. 4, pp. 318 330,December 2005, ISSN 1000 8349

1. Barred Galaxies: Studying the Chaotic and Ordered Nature of Orbits,Manos T. and Athanassoula E., 2006, In Proceedings of the 7th Astronomy Conferenceof the Hellenic Astronomical Society (in press), astroph/0601331

2. Chaos and the Dynamical Evolution of Barred Galaxies,Manos T. and Athanassoula E., 2005, In Proceedings of the 5th International ConferenceThe fabulous destiny of galaxies: bridging past and present, Marseille, France (inpress), astroph/0510823

3. Detecting Chaotic and Ordered Motion in Barred Galaxies,Manos T. and Athanassoula E., 2005, In SF2A2005: Semaine de l AstrophysiqueFrancaise, eds. Casoli F., Contini T., Hameury J. M. and Pagani L., EDPSciences,Conference Series, 631 632

1. Caos e integrabilidade em teorias com supersimetria,de Assis L. P. G., 2005, Centro Brasileiro de Pesquisas Fisicas CBPF, Rio de Janeiro(in Portuguese)

IX

= ". 2 '# 0.# =>>? " , % "

@ $ #

1. Detecting Chaotic Orbits in N Body Simulations,Gemmeke J. F., 2005, Faculteit der Natuurwetenschappen, Wiskunde en Informatica,Universiteit van Amsterdam, The Netherlands

= 2 3 2# & $# -# '

% =>>A 3% '% " ',

" # $ $#

X

III

!"#$ VI

% &'&( %

) "*+, (" - ./#0# ,#1"#$ #-( ,"-( 2

=B 3 & 2# @ C

== 2# 3 $% B>

==B 9 B>

=== 7!+ 2# 3 B=

==D ! Poisson BA

=D 9 2# @ 3# & B?

=DB Poincare ', "# 2# & BE

=A 9 % KAM BF

=? % " 9# 2# 3# & BG

=E ' Lyapunov 3# & BC

=EB $!! BC

=E= $!! =D

=F ( 3% '! :SALI; =A

=G ( 0 3% '! :GALI; =E

=GB 9 0 3% '! :GALI; =C

=G= DB

=GD ' $# GALI 2 " DB

=GA 9 6! GALI ! 9! " DE

=C % $ 3# & AA

XI

XII 3

=CB 9 $% 9%+ AE

=C= 9 $% AF

=CD ' 0# 2# & ?>

=B> % Krein ??

=B>B 3! 0 2# & ??

=BB & '% )9 ?C

=BBB '!! & 1 * ?C

=BB= & '! 6 EB

=BBD '% )9 ED

=B= ,% E?

4 ,#,-( - #-( ,"-( "*+,', !"5', ) - 4 /"+,

*/ 67

DB ( &, SALI ! " EF

D= ( &, SALI ! 2 " F>

DD $ $# ". 2 FF

DA ( &, SALI 7 2 $ GD

D? &! SALI ) 3 . 2 9! G?

D?B " > 1 B : 1 )'; " ) G?

DE ',! 0 3% '! :GALI; GC

DEB 2 3 - '% GC

DF & CG

8 ,#,-( - #-( ,"-( "*+,', !"5', #**+, /"+, 9

*/ %:%

AB $ 7 " ' B>B

ABB # $# 7 FPU B>A

AB= 0% 0% " " B>?

ABD ' " " SPO1 FPU B>F

ABA 7 OPM FPU BB>

AB? ', ! $ 7 BEC BBD

ABE ' "# IPM OPM BEC BBE

ABF 7 SPO2 FPU BBF

A= $7 @ Lyapunov 2# $# B==

A=B % $7 9 2 3 B==

A== 6! ! 4% ' B=A

A=D @ Lyapunov )9 B=G

A=A & @ Lyapunov 2 $ BD>

A=? ( " SPO2

',# FPU BD=

3 XIII

AD ',! 0 3% '! :GALI; $!FPU BDA

AA & BDE

; 5/ - 5# "1*--. 1-#, - /"# %8%

?B '!! BAB

?= 3,.% ! BAD

?D 3 & BA?

?A $ !% BAE

?? $ BAG

??B 2 1 BAG??= 4 1 B?>??D -% & B?B

?E & B?A

6 (" %;7

EB $

HI 3!% 0 & " SPO1B?FE= $

-HI 9 $% Hill $ " SPO1 B?CED $

0HI & Fourier $ &

Q(u) SPO1 BE>EA $

3HI '. 0 (Wedge Product) BEDE? $

'HI ( & . GALI2 SALI BE?

?*#&

" 2005 # 50 .!! Fermi, Pasta Ulam (FPU) JDBK Los Alamos ($ " % , ! ! N = 32 % % + & % & %

! % !

% % ! %

%

# !% % .#

% ,% % &

&! Fermi, Pasta Ulam # % ! . # ! % %

! # % %

! ) .. ' %

% % )

.. : BB;

2# N = 32 # % % ! .# 2 2#

H =12

Nj=1

x2j +N

j=0

(12(xj+1 xj)2 + 13(xj+1 xj)

3 +14(xj+1 xj)4

)= E :BB;

% ! )9 = 0 = 0 FPU 1 :

; # = 0 = 0 FPU 1 :

;

B

)

:BB; 32 ! , .# %

xj(t) = xj+12xj+xj1+((xj+1xj)2(xjxj1)2

)+((xj+1xj)3(xjxj1)3

), j = 1, . . . , N.

:B=;

% !%

! % !

! % % & Maxwell Boltzmann Gibbs 18 # ( ! ! # ! % ! .% # !#

! %

& BBI & ! 32 % % Fermi,Pasta Ulam # % BC?? Los Alamos ($ " % % # , :normal mode q = 1; % !

$ ! . Fermi, Pasta Ulam % .

! % ! # ! ! E #

% % , , % , !% #

4 1 5 L N :normal modes; % ! BB = = 0 JA EG AG DCK

9 , % , :%; + , :recurrence time; 9 , % , %+

! % ! @ !

! ! E : % ! ; :BB; ,

! % ,

:B=; " # %I 0% ! !

! ! # %

. # ,M

" , ,# FPU % /

!

# # !% ! !

! % % ! % % #

4

! .! , ,% Izrailev Chirikov JACK 1966 % %+

,# FPU %+ ! / # # # " 1979 ! ! # / ! ! % 2#

# # % J=>K

! %

,# ,% % ,

.% % % : BB ; % ! ! . %

! J=C =GK " , Flach et al. JD=K .% % % % ! ! FPU % !% ! #

:$7; % ,

% % % ' 50 %, / Fermi, Pasta Ulam , ! , + %+

JBBK ( %

# , ! 2#

N E E/N = )' % # #

! , Lyapunov 9 ,% 2# %! # % #

,% ,

%

&

" = , % % !! % !# %

!

, ' +

2# I " 3% '! :SALI; 0 3% '! :GALI; $!, % %

!# # I 2#

Poincare Lyapunov % $ #" , & % %

'! 6 % !

%! # % & 2#

& , D +

, 2# %! #

% " # ! ,!

8

SALI GALI . ! , .!! % !# % Ao ?o ,

$+ ! , # #

! !% % !#

!,% " !, ,! ,

0 3% '! GALI # 2# 2 3 # % )9 # % 2# %! # %

& A , % %

# 2# # #

% % FPU Bose Einstein Condensation (BEC) '# % 2# %! # % # # !

! # , &

% , ,

% ! &

% 2#

' , A # $7 SPO2 FPU 1 : :BB;) % % % % % #

% % #, !

!% J=CK JD=K ! / . %+ #, SPO2 : B=;N % .# #

,# FPU k

# % ! :k 4; # !

SPO2 % /I k = 2(N + 1)/3

9 ! # ,% ! % %

SPO2 : B=; % ! N2 % N1 ! % $7 % + SPO1 % % % ! ' %

% $7 ! # SPO2 SPO1 , SPO2 !% !

! %

OPM(=Out of Phase Motion) 2# FPU BEC !

E N . % ! E/N &! !%+ ,

# Lyapunov 2# FPU BEC # OPM ! ! % % -% Lyapunov !!%+ ,

Li L1ei/N :BD;

;

& B=I ( ! ! % ! E2u/N # %

SPO2 2# :BB; # ! % N2 DeLuca et al. J=CK Flach et al. JD=K

! 2.76 3.33 % i = 1, 2, . . . ,K(N) K(N) 3N/4 " # % %I$# $7 !

% % , !

# ,M

& BDI &! , Lyapunov !# # $7 SPO1 SPO2 ! E = 2.62 $7 % % ( . !#

# % 102

0 # ! #

$7 # , % % ! % . !

6

% % ! % !% % ! Lyapunov 1 N1 !% % , Lyapunov $7 %

% %

BD &! , , % BD % .

! % , ! FPU #, $7 ) ..

2# :BD; ! Lyapunov Li, i = 1, . . . , N !%+ #

hKS =

i:Li>0

Li :BA;

% + % Kolmogorov Sinai :hKS; ,+ & , % %+ !# # !

! # # Lyapunov -% ! FPU BEC %Kolmogorov Sinai .

# %N # % :extensive quantity; &

& A ,% ! %

# , % % ,

% % & !

% N 0 ! % FPU N = 5 SPO1 !% % # !% , ! % , Poincare :x1, x1;! ! x3 = 0 $7 SPO1 % # ! ! ! !#

# , # Poincare

O# % %! % ! %

% + .! # !, " ! % % % %

,%+ % ! Poincare

" , ? 2#

# % $+ %

+ 3% '! SALI 3 ! % # % % /# !# & # % % #I

$ % ,! % # , # .

% % # % :

; M

7

0 / # !

# # . ! %

0 3% '! GALI $

@

!

"# $

% & ' ( " )SALI* + (" )GALI*! , #

% # Poincare#

Lyapunov#

' , ! - '# #

" .#

/ ' - '

!

! ,# # !#

# # % % +! . !

.%

0 ! n # # % %

# ! t % !# #

# % (x1(t), x2(t), . . . , xn(t)) n 1

'% # Rn " # :!# n

# #; !% D # Rn D Rn % +

C

%:

9 % !# t !,

,# .# :&3'; # . ,

xi(t) dxi(t)dt

= Fi(x1(t), . . . , xn(t)), i = 1 . . . n :=B;

% % !#! 9

.# %+ % , !

t % x1, . . . , xn D

9 Fi,i = 1 . . . n % ! !

xi %+ D '% : .C1; # ! (x1(0), x2(0), . . . , xn(0)) =(1, 2, . . . , n) D &3' :=B; % +

xi(t) = Si(t, i), i = 1 . . . n, t 0

!, .. # !# " %

! # !

# , D

V = (x1)(x2) . . . (xn) :==;

.. 0 # !%

:=B;

d

dtV = (V )

ni=1

xixi

= (V )( F ) :=D;

% # , D

F = 0 :=A;

% % F % # F !

" % % % t ,%+

, .# % % ,%+

)' 2# % % n = 2N # !, % !

H !

q =

!" # # %%

(q1, . . . , qN ) )$* p = (p1, . . . , pN )

H : RN RN R, H = H(q, p).

(

H + Hamilton 2#

.# % !, qi pi % ! Hamilton

qi =H

pi, pi = H

qi:=?;

! i = 1 . . . N

9 .# :=?; !, % ,

= H()

= [1, . . . , 2N ]T ,

a = qa (a N) a = paN (N < a 2N) =

[ 1

, . . . ,

2N

]T.

9 % 2N 2N % [0N ININ 0N

],

0N % N N % IN % % N N % ,

.# Hamilton :=?; % 2#

H = H(q, p)%

dH

dt=

2Ni=1

{H

pipi +

H

qiqi

}= 0 :=E;

% : %;

I1 = H(q, p) % N 1 # %

Ii = Ii(q, p) :=F;

dIidt = 0, i = 2, 3, . . . , N % % .

# .

% . ! Poisson % % [Ii, Ij ] = 0, i, j = 1, 2, . . . , N

! & % , ! #

Liouville Arnold % ! % N 1 % :

%)

,#Pfoliate; # , % # 2# ! ! %

#

! "

0 # 2#

%+ # ! 7!

+

L : RN RN RL(q, q) = T (q) U(q), :=G;

q = (q1, . . . , qN ) q = (q1, . . . , qN ) % %

# % T (q) = 12

Ni=1miq

2i % mi % + i

% U(q) % .# Euler Lagrange ! ! qi % ,# .# JAGK[GOLDSTEIN, ter Haar]

d

dt

(L

qi

)=L

qi, i = 1 . . . N. :=C;

7!+

L %+ (q, q) (q, p) H Legendre

L : RN RN RN RN

L(q, q) = (q, p) =(q1, . . . , qN ,

L

q1, . . . ,

L

qN

). :=B>;

( / :=B>; .

7!+

L # , JEBK !%+ %+ *

% DL() % , RN RN Legendre % ,,, % 9 * %DL() L %

JL =[IN 0NCN HN

], :=BB;

IN N N % 0N N N % CN N N % HN N N % D2qL(q, q) ' / Legendre L.

%+ * % JL % , % .

' 7!+

L

%+ % HN

!" # # %4

9 Legendre % ,,,N % + %

(L) = {(q, q) RN RN / D2qL(q, q) = 0}

/ !% +

% %

, JEC F> FBK #

Legendre %+ RN RN (L) % ! 2#

H

Legendre L

H : RN RN R R

H(q, p, t) =N

i=1

qipi L(q, q, t). :=B=;

6!%+ , 2#

H

dH(q, p, t) =N

i=1

H

qidqi +

Ni=1

H

pidpi +

H

tdt. :=BD;

, H !, %

dH(q, p, t) = d

(N

i=1

qipi L(q, q, t))

=N

i=1

(pi L

qi

)dqi +

Ni=1

qidpi N

i=1

(d

dt

(L

qi

))dqi L

tdt, :=BA;

# :=BD; :=B=; 0 % . :=BD; :=BA; %

% pi

pi =L

qi, i = 1 . . . N.

+% :=G; %

L(q, q, t) = T (q) U(q) = 12

Ni=1

miq2i U(q)

L

qi= miqi pi, i = 1 . . . N

qi 2# ,

%8

)' % # . :=BA; .#

:=BD; :=BA; % ! .# 2#

qi =H

pi, pi = H

qi:=B?;

! i = 1 . . . N # .%

L

t= H

t.

9 qi pi % ! ! % 9 .#

:=B?; %

# $ Poisson

# % Poisson % '%# RN %+

[u, v] =N

k=1

(u

qk

v

pk vqk

u

pk

)

u v %

! ! qk pk 9 .# %

!, , !# Poisson ! u !

v 2#

H

[qi,H] =N

k=1

(qiqk

H

pk Hqk

qipk

)=

H

pi, i = 1 . . . N.

! [pi,H] )' % ! 2# .# :=B?; ! !, ! # ,

qi = [qi,H], i = 1 . . . N,pi = [pi,H], i = 1 . . . N.

3 !# Poisson % Jacobi( # %

[u, v] = [v, u]

!" # # %;

%

[[u, v], w] + [[w, u], v] + [[u,w], u] = 0.

2# # 2# .# !! %

%

= (q, p, t) q, p !, ,

d

dt=

Ni=1

(

qi

H

pi Hqi

pi

)+

t= [,H] +

t. :=BE;

.

t = 0 ' ! Poisson [,H]%+

%

!" #

,

!, =B % .# % #

,

S = (S1, . . . , Sn) (

%+ + % % % = (1, . . . , n) n = 2N

$# !% M .# &3'

:=B; :=?; % # N

3# % #

% . # % !, %

!

: ;

)' ! % # !+

,! !/

9 ,+ % % % %

! % # # % N 1

2#

)9 ! + % %

.!! 7! . % .

# % , !#

. ! %

% ! ! ! ,

# , " !!

. , % . . , # #

%6 $ %&'# !"# ' ( &

# $ Poincare % &! !' %'

( , # Poincare % % !

2# 0 ! 2# #

% # , % A 1 ! # q1 q2 #

p1 p2 % '! % P # , % ! !% #

% : =B; " % , % !

% % P . : ,; % % % 2#

& =BI '.. % # , , #Poincare P

% ! , P % % ! % !

2#

H = H0 % % q2 = 0 : =B; 7 .% H(q1, 0, p1, p2) = H0 p2

p2 = p2(q1, p1;H0) :=BF;

! (q1, p1) ! % % ,! %

& %+ P # ! !! % 41 # ,

!" # # %7

( , # Poincare !! ! 2# N = 2 % 0 2# N # % % ! , H = H0 # , 2N 1 ' qN = !% ! # pN

.% H = H0 + % 2N 2 1 , # !q1, . . . , qN1, p1, . . . , pN1 % ! ( ,

% N = 2 % # % , # Poincare % % % 2N 2

' % 2# N = 2 # % %

I(q1, q2, p1, p2) = I0 = R :=BG;

.

H . [H, I] = 0 H(q1, q2, p1, p2) = H0 :==G; % .% %, p2

p1 = p1(q1, q2;H0, I0). :=BC;

' ! H0 I0 , #

Poincare % % % % :=BC; q2 = 0 . . # %

# # , # Poincare . ,% % %

% # " .%+ #

2# % , % F = 0 ! # , %

F = (H

p1, . . . ,

H

pN,H

q1, . . . , H

qN). :==>;

"

F =N

i=1

qi

(Hpi

)+

Ni=1

pi

(Hqi

)= 0 :==B;

# 2#

H % , # !!%

$ % & KAM

)9 , . N .

%. # %

! 2# ( %

, .# -

% %+I ( % !%

%@ ) && ! *+'# %&' &

! % % %

( ,# .# Hamilton %

/ .. !

% !% # % %

% # # !% % !

!# # ! # 2#

# % KAM .,%+

. % % %

J?AK

"% !% 2# #M)9 % ! % ,%+ # , ! %

% % # # !

% .

%

% % !# 9 ! %

% ,! % KAM .% %

!

' ! (% #

! ! # % %

2#

H H % # % ,%

1 !% J1, J2, 1, 2

H(J1, J2) = E,

E % ! J1, J2 % #

% # , %+ ! %

% ! ! % ! %

1 = 1t+ 10, 2 = 2t+ 20,

!# % .

!% % % 2

)9 # Liouville1 Arnold % % !% , # , % , Poincare p1, q1 : p2, q2; # # " J1 % # !% 1 :

2; !, % # % % , ,%

!" # # %2

==

0 J1 J2 , 1 = 1(J) 2 = 2(J) ! = 12 % % 0 = r/s r, s . # % % = 2/T % ,%+ % % T % % r , !% 1

s !% 2 ! %

%

( % % % !% % !%

% % #

# , % # ! N %

N % % 2N 1 # !% % N 1 , % !, N N !%

& ==I % # , % ! # % :a; ( % % J1 = J2 = :b; " % , # 2 = ! #

) *+ Lyapunov

( )' *+

9 Lyapunov %+ % 2# # ! % !%

): , # Lyapunov &

% , % ! Lyapunov % %+ !# # .. 9 !

# , !# # # JA?K

Zaslavski Chirikov :1972; JCGK Froeschle Scheidecker:1973; JDFK Ford :1975; JDDK

( % # Lyapunov :1907; # Oseledec :1968; JEDK # . # # # Benettin et al. :1976;JB>K Pesin :1977; JEFK 7%! ! % % % ! ! # Lyapunov Benettin et al. 1980 JG CK

&, Benettin et al. JG CK Lyapunov ! % x(t) ! &3' # .

dx(t)dt

= V (x(t)) :===;

M 1 '% # x(t) RM %+ .I % , :===; x(t) % ! x

(t) x0 x0+x0 %

=D # d(x0, 0) = ||x0|| x0 '%

& =DI 9 ! Lyapunov 3 ! % x0 x0 +x0%+ # . 9 . w(t)%+ ! Lyapunov

!" # # )%

9 .% .

d(x0, t) = x(x0, t) = ||x(t) x(t)||.

# # :===; x(t) x(t) + w w x = (x1,x2, . . . ,xM ) !# ! ||w(t)|| ||x(t)||

dw

dt= J(x(t))w, :==D;

J V/x % * % % V (x) ( :==D; + )variational equation* &3' :===;

9%+ # !# # # ,

(x0, w) = limt

(1t

)lnd(x0, t)d(x0, 0)

.

% % ! ! ! w : ||w|| x0; % JEEK ' % M 1 ei !

w # ! , ! w % % :!

;

i(x0) = (x0, ei), i = 1 . . .M

% + Lyapunov

9 % #

1 2 . . . M

% .

! # , '% , %

Lyapunov % # # ,

.%+ # w . !

x(t) % % .# # :==D; 0 ! ei ,% % % i = 00 % 2# x = (q, p) % V %

V (x) =(Hq

,H

p

)

)) , # Lyapunov &

# % Lyapunov %

i = 2Ni+1, i = 1 . . .M

M = 2N N % # % 2# '! % 2#

# Lyapunov %

Mi=1

i(x0) = 0,

# ! # ! # ,

%

# # Lyapunov # %

.I ! % % %

!# # % #

. " %+

Lyapunov !

# Lyapunov w % ! Lyapunov # . 9 Oseledec :1968; !% ! !/

. p 1 ! # , :p M; %+ p w1, . . . , wp 2#

Vp = w1 w2 . . . wp

! Vp p 1 % %

w1, . . . , wp # % !# .

p

% %+

(p)(x0, Vp) = limt

(1t

)ln

Vp(x0, t) Vp(x0, 0)

:==A;

! Lyapunov p 9 Oseledec :1968; JEDK Benettin et al. :1980; JG CK. (p) % p # Lyapunov . 1 )9 (x0, w) =1(x0) ! w # % (p) % p! # Lyapunov

(p) = (p)1 = 1 + 2 + . . .+ p :==?;

! Vp ( :==?; % ! !

!" # # )4

# Lyapunov 0 p = M . ! # ,

(M) =Mi=1

i(x0).

( $,

*+

# Benettin et al. :1980; JG CK !/ # # !% Lyapunov M 1 '%, !% 1, 2, . . . .

! , w e1, e2, . . . .%

% !% e1

% % . % p #

!% p 1 ! Vp(t) % %+ 9 Lyapunov (p)i . p % :==A; ! p = 1 . . .M :==?; ! Lyapunov 1, 2, . . . )9 % % ,%+ # , .% !% . !% !#

% ! ! ' ! ,!

t % #

! # ! %

!

( % % % ! !%+ .. p 1

! !%+ .. M

% Gram Schmidt ! .I )' wk1(0) , % x(0) wk1(t) , % x(kt)2# !%+ # ! t

d(1)k = w(1)k1(t)

w(1)k (0) =

w(1)k1(t)

d(1)k

.

)' !%+ ! j = 2 . . .M

u(j)k1(t) = w

(j)k1(t)

j1i=1

[w(i)k (0) w(j)k1(t)] w(i)k (0),

d(j)k = u(j)k1(t) ,

w(j)k (0) =

u(j)k1(t)

d(j)k

.

)8 - . # '+' & ' /SALI0

(k 1) 1 ! Vp . d(1)k d

(2)k . . . d

(p)k ( :==A; %

(p)1 = limt

(1Nt

) Ni=1

ln(d(1)i d

(2)i . . . d

(p)i

).

,# (p1)1

(p)1 # :==?; p 1

Lyapunov

p = limt

(1Nt

) Ni=1

ln(d(p)i

).

$ % # %+

/ ! .,% ! # Lyapunov' %+ % # #

Lyapunov ! ! ! , #Poincare 3% '! :SALI; : !,; # !, % # % , :

; # ,% !

%

, - .+ .% *+/ 0SALI1

( 3% '! :Smaller ALignment Index % SALI;% % %

# ! , % 2#

% # % (

! # , JFFK ,% % ! !

% 2# # % ,

J= D BE =F EA FG FC G> GB G= GE GFK

9 SALI % % % % % #

# # # ! ! ! 2

# N 1 # % N > 2 & %

= 1 % SALI % ! ! ! , # % %

. # JFFKI ' SALI % # ! # , #

: 1/t) ! '# # !, , SALI % 2

# # 0 ! 2N 1 # ,

!" # # );

2#

H H(q1(t), . . . , qN (t), p1(t), . . . , pN (t)) = E :==E;

qi(t),i = 1, . . . , N % ! pi(t),i = 1, . . . , N % +!% E % ! ( .. %

x(t) = (q(t), p(t))

x(t0) = (q1(t0), . . . , qN (t0), p1(t0), . . . , pN (t0))

%+ 2N # . ,# .# %

dqi(t)dt

=H

pi(t),dpi(t)dt

= Hqi(t)

, i = 1, . . . , N :==F;

%+ 2#

:==E;

0 % 3% '! SALI .# # ! % # Lyapunov 9 .# % % ! .# .# % ,#

.# :==E; :=BF; ! % , :==E; :=BF; %+

di(t)dt

= J (x(t)) i(t),i = 1, . . . , 2N :==G;

J (x(t)) % * % . :==F; !

x(t) = (q1(t), . . . , qN (t), p1(t), . . . , pN (t)) " i(t) = (i,1(t), . . . , i,2N (t)),i =1, . . . , 2N % ! & ! k(t) l(t) 1 k, l 2N, k = l %+ % '! :ALI; t t0

ALI(t) = k(t) k(t) l(t)l(t)

, ALI+(t) = k(t) k(t) + l(t)l(t) :==C;

# '% R2N 9 3% '! :SALI; %+

SALI(t) min{ALI(t),ALI+(t)},t t0 :=D>; t0 % ! .. JFFK JFGK

)9 % ALI(t) 0 % k(t) l(t)

% !% !! !% %

ALI+(t) 2 # ALI+(t) 0 % !%!! !% % # %

)6 1 . # 2 & ' /GALI0

ALI(t) 2' x(t) = (q1(t), . . . , qN (t), p1(t), . . . , pN (t)) %

!%+ , #

:=D>;

limt

SALI(t) = min{0, 2} = 0.

& % , x(t) = (q1(t), . . . , qN (t), p1(t), . . . , pN (t)) %!

limt

SALI(t) = 0

SALI # ! % (0,2] # ! % !% , % % #

% ! ! , # .%

%+ N .

% N #

% JFCK

# , , SALI ! % . . 2#

'% % , ! # k(t0) l(t0) % % + #

JFF FC G> G=K " SALI % ,% %! # % % % % .

% # # Poincare

#

2 - .+ 3 *+/ 0GALI1

( % %+ % SALI % %+

! # Lyapunov % % , $! # .

% . % ! ! # Lyapunov # t

0 !% SALI .. , +% % w1(0) w2(0) ,, !

% B

wi(t) =wi(t)

wi(t) , i = 1, 2 :=DB;

!" # # )7

# '% :; # % 1 " SALI %+

SALI(t) = min {w1(t) + w2(t) , w1(t) w2(t)} , :=D=;

, SALI(t) [0,2] ( SALI = 0 # !% % % % . %

!# .

'% # + #

:=D=; % !%+ !

P (t) = w1(t) + w2(t) w1(t) w2(t) , :=DD;

t $! # % :

% % ; % !

P (t)

P (t) % > % ! % ! % # : % ! #;

% !%+ SALI(t) :=D=; !%+ .! :wedge product; : $

EA; w1 w2 ! %

w1 w2 = w1 w2 w1 + w22 . :=DA;

% ! %+ ( . :=DA; % $

E? $! . !

% # ,% , ! #

! %+ w1, w2, . . . , wk 2 k 2N % % 2# N 1 # %

% 2N 1

)' % %

!% / # # # ,#

2# # # % $

!, !% ,# !%+ !

SALI ! % w1 w2I

B & % # # w1 w2 !%+

! Lyapunov SALI(t) % % .! , % !,

%.

SALI(t) e(12)t :=D?;

1 2 ! Lyapunov 1 > 2

)@ 1 . # 2 & ' /GALI0

= & % ! % %

w1 w2 % % , # #

t1 ! , # SALI # ! , ,% !! ! !

. !

' % % , % % ! !

% ! t1 w1 w2 , , #

.%+ # % % ,%+

% % % , # % '

% % !% !# .

SALI % # 1/t % , , # SALI % % % # % JFFK " !! ,

SALI ! ! .

, #

, '% , % %

!% # ,% % ! %

% !% % %

# % , , SALI !

:=D?; #I $# , %

% % ! % ! Lyapunov 1 2 % % %M $ % % 2#

# !,% :# ;N '

JDK % 2# B? #

% ! Lyapunov % % !! 1 2 0.0002 '% , , ! % % ,% !! SALI 2+ %

! # %

# SALI 0 ! % % % .

Lyapunov % $ % ,%

, , , #

% !, #

% % ! : ! sticky; # # 1 # :super integrable;

!" # # )2

- .+ " ,%

/GALI0

2# N # % 2#

H(q1, q2, . . . , qN , p1, p2, . . . , pN ) = h = :=DE;

qi pi i = 1, 2, . . . , N % ! +!% % %

%+ x(t) = (q1(t), q2(t), . . . , qN (t), p1(t), p2(t), . . . , pN (t)) xi = qi xi+N = pi i = 1, 2, . . . , N ( .. % .#% Hamilton

dx

dt= V(x) =

(H

p,H

q

), :=DF;

# .. w(0) = (dx1(0), dx2(0), . . . , dx2N (0))

! x(t) = (q(t), p(t)) % .# #

dw

dt= (x(t)) w , :=DG;

= V/x % * % V

)9 % SALI % % ! w1 w2 : % B; !% !# .

!+ %

( ! .

#

,# % !

0 .. k

w1 w2 . . . wk 2 k 2N % !% !# .

! !

% ! !%

. ! : wedge product; # : !, =EB $

EA ! . !;

)9 wi i = 1, 2, . . . , 2N 2N 1 , # 2# 2# #

#

e1 = (1, 0, 0, . . . , 0), e2 = (0, 1, 0, . . . , 0), . . . , e2N = (0, 0, 0, . . . , 1) :=DC;

wi % !,% ! #

wi =2Nj=1

wij ej , i = 1, 2, . . . , k :=A>;

4: 1 . # 2 & ' /GALI0

wij % !% %

2Nj=1

w2ij = 1. :=AB;

' !, :=A>; ,w1

w2

wk

=w11 w12 w1 2Nw21 w22 w2 2N

wk1 wk2 wk 2N

e1

e2

e2N

:=A=;

# .% :EAA;

EA . ! k

,+

w1 w2 wk =

1i1

!" # # 4%

wi %

# #

!% !# .

GALIk(t) # t !# % # ! .

% . %

% ! % , ! t

- )'% $+

-# , *! GALI + &1+

$ 2# N

% % ! # Lyapunov# , !, =E J?A =EK )9 .! %

(x(0), w) % , x(0) %

(x(0), w) = limt

1t

lnw(t)w(0) . :=AE;

' % 2N 1 {u1, u2, . . . , u2N} , # 2#

# (x(0), w) % % 2N :# ;

i (x(0)) = (x(0), ui) , i = 1, 2, . . . , 2N :=AF;

% % % Lyapunov . !

1 2 . . . 2N . :=AG;

9 !% , % %

% T & % % .# # :=DG; %

T 1 2N 2N % .# # % !,% ,

w(t) = (t)w(0) , :=AC;

(t) % ! # % JCB ED BD E?K # (0) = : 2N 2N %; &, % Floquet , w(t) %+ i % (T ) . |1| |2| |2N | % ui i = 1, 2, . . . , 2N % " w(0) = ui

w(nT ) = ni ui , i = 1, 2, . . . , 2N :=?>;

# :=AE; %

(x(0), ui) = limn

1nT

ln |ni | =ln |i|T

, i = 1, 2, . . . , 2N. :=?B;

4) 1 . # 2 & ' /GALI0

' !/

w(0) =2Ni=1

ci ui , :=?=;

:=?>; .. # w(0) t = nT%

w(nT ) =2Ni=1

ci einT ui , :=?D;

0 % %

$ Oseledec JEDK . . {u1, u2, . . . , u2N} # Lyapunov ! .!% % # , 2# % %

! % % )'

.. % !!% % .% , :=?D;

w(t) =2Ni=1

ci editui , :=?A;

ci di % !% % .

# % # , %

, ' di i = 1, 2, . . . , 2N % Lyapunov ! t i i = 1, 2, . . . , 2N % ui i = 1, 2, . . . , 2N % !

# w(0) = ui . % .. ! !

! Lyapunov 1 .%

L1(t) =1t

lnw(t)w(0) , # 1 = limtL1(t) , :=??;

w(0) w(t) % % t = 0 t > 0% JG CK

# % 2#

:=DE; N % )' L1 L2 L2N % Lyapunov % # ! 1 2 2N % % % Lyapunov '% !

2# % ! # Lyapunov % % !% % +! %

% % J?AK ' i = 2Ni+1 !i = 1, 2, . . . , N 1 2 N1 N = N+1 = 0 N+2 2N 6 Li i = 1, 2, . . . , 2N # !

!" # # 44

!! Li i .. wi % !,%

wi(t) =2Nj=1

cij ejt uj , :=?E;

cij % !% % uj !% %

j % ! !# .

"%+ # cij uj

.

# , ! # , %

6 # 1 > 2 % % % # . '%

! ! %

wi(t) |ci1|e1t . :=?F;

% C : $

:EAD;) # k

wi(t) = wi(t)/wi(t) i = 1, 2, . . . , k 2 k 2N # # {u1, u2, . . . , u2N} !%

C(t) = [cij ] =

s1

c12|c11|e(12)t c

13

|c11|e(13)t c12N|c11| e

(12N )t

s2c22|c21|e(12)t c

23

|c21|e(13)t c22N|c21| e

(12N )t

. . .

skck2|ck1 |e(12)t c

k3

|ck1 |e(13)t ck2N|ck1 | e

(12N )t

, :=?G;

si = sign(ci1) i = 1, 2, . . . , k j = 1, 2, . . . , 2N [w1 w2 . . . wk

]T= C

[u1 u2 . . . u2N

]T:=?C;

:T; # , % " . ! k

!%+ ! .% :=AD; ,

w1(t) w2(t) wk(t) =

1i1

48 1 . # 2 & ' /GALI0

.%+ # #

Sk =

1i1

!" # # 4;

# %

D1,j1,j2,...,jk1 =

s1c1j1|c11|e(1j1 )t c

1jk1|c11|

e(1jk1 )t

s2c2j1|c21|e(1j1 )t c

2jk1|c21|

e(1jk1 )t

skckj1|ck1 |e(1j1 )t c

kjk1|ck1 |

e(1jk1 )t

=

=

s1c1j1|c11|

c1jk1|c11|

s2c2j1|c21|

c2jk1|c21|

skckj1|ck1 |

ckjk1|ck1 |

e[(1j1 )+(1j2 )++(1jk1 )

]t

:=EA;

1 < j1 < j2 < . . . < jk1 2N ' .. D1,j1,j2,...,jk1 %+%

D1,j1,j2,...,jk1 e[(1j1 )+(1j2 )++(1jk1 )

]t. :=E?;

$ % %+ # % C :=?G;

,

Dj1,j2,...,jk =

c1j1|c11|e(1j1 )t

c1j2|c11|e(1j2 )t c

1jk

|c11|e(1jk )t

c2j1|c21|e(1j1 )t

c2j2|c21|e(1j2 )t c

2jk

|c21|e(1jk )t

ckj1|ck1 |e(1j1 )t

ckj2|ck1 |e(1j2 )t c

kjk

|ck1 |e(1jk )t

=

=

c1j1|c11|

c1j2|c11|

c1jk

|c11|c2j1|c21|

c2j2|c21|

c2jk

|c21|

ckj1|ck1 |

ckj1|ck1 |

ckjk1|ck1 |

e[(1j1 )+(1j2 )++(1jk1 )+(1jk )

]t

:=EE;

1 < j1 < j2 < . . . < jk1 < jk 2N ' # +# % % #

Dj1,j2,...,jk e[(1j1 )+(1j2 )++(1jk1 )+(1jk )

]t. :=EF;

$,# kk %+ ,%+ Sk :=EB; % # ! % # % C :=?G;

46 1 . # 2 & ' /GALI0

D1,2,3,...,k e[(12)+(13)++(1k)]t . :=EG;

)9 %+ :=E?; % !! %+ D1,2,3,...,k

, % ! % :=EG; '

% % Sk % :=EG; % !!

Sk(t) e[(12)+(13)++(1k)]t . :=EC;

' , :=AA; k 1 w1 w2 wk .% Sk % GALIk % %

GALIk(t) e[(12)+(13)++(1k)]t . :=F>;

$ ! k = 2 % GALIn(t) % SALI !/ =F & % :=F>; ! :=D?; . %

% % JG=K

! 1 > 2 #

% !!% .% :=?F; '

# m Lyapunov 1 < m < k % % . 1 2 m .% :=F>; !%

GALIk(t) e[(1m+1)+(1m+2)++(1k)]t , :=FB;

% !, % # $ ! k m < N GALIk % # % %+ % C % %+ & %

, . #

% # GALIk ( Lyapunov i = 0 # % % ! !,

-2 3 GALI '+ &1+

)9 % ! ! 2# N # %

:=DE; ! % N 1 &, % KAM ! % % %

!,% N :formal; # % . # ,%+ # %

1 !% # J1, J2, . . . , JN N

!" # # 47

# # .# Hamilton % ,

Ji = 0i = i(J1, J2, . . . , JN )

i = 1, 2, . . . , N. :=F=;

9 .# # %

Ji(t) = Ji0i(t) = i0 + i(J10, J20, . . . , JN0) t

i = 1, 2, . . . , N, :=FD;

Ji0 i0 i = 1, 2, . . . , N %

9%+ i i i = 1, 2, . . . , N % Ji i % !, .# # :=F=; ,

i = 0i =

Nj=1 ij j

i = 1, 2, . . . , N, :=FA;

ij =iJj

J0 i, j = 1, 2, . . . , N, :=F?; J0 = (J10, J20, . . . , JN0) = N 1 # ( .# # %

i(t) = i(0)

i(t) = i(0) +[N

j=1 ijj(0)]t

i = 1, 2, . . . , N. :=FE;

:=FE; w(0) ! i(0) i =1, 2, . . . , N i(0) i = 1, 2, . . . , N !% w(0) =(1(0), 2(0), . . . , N (0), 1(0), 2(0), . . . , N (0)) .% # ! # ! !% . !

, # % ! .

0 w(0)

Ni=1

i(0)2 +N

i=1

i(0)2 = 1 :=FF;

# .. %

w(t) =1 +

Ni=1

Nj=1

ijj(0)

t2 +2 N

i=1

i(0) Nj=1

ijj(0)

t

1/2

, :=FG;

4@ 1 . # 2 & ' /GALI0

# w(t) !,

w(t) =1

w(t)

1(0), . . . , N (0), 1(0) + N

j=1

1jj(0)

t, . . . , N (0) + N

j=1

Njj(0)

t .:=FC;

, :=FG; ! ! t . !

t :=FC; % , #

, ! % : !

; # t1 , % % ! 2# # %

JFCK

2# 2N 1 # , 2N % {v1, v2, . . . , v2N} # # N v1, v2, . . . , vN N

vN+1, vN+2, . . . , v2N N +!% !%

% wi i = 1, 2, . . . % !,%

wi(t) =1

w(t)

Nj=1

ij(0) vj +N

j=1

(ij(0) +

Nk=1

kjij(0)t

)vN+j

. :=G>;"%+ # ij i, j = 1, 2 . . . , N :=F?; .

# ! , ! .%+ %

# ui i = 1, 2, . . . , 2N .

! % .% %+ ei i = 1, 2, . . . , 2N .% :=DC; , .% :=E=;[

v1 v2 . . . v2N

]T= o

[e1 e2 . . . e2N

]T:=GB;

o %+ % )' {e1, e2, . . . , e2N} %! !/ .. ! qi

pi i = 1, 2, . . . , N 2# :=DE; # {v1, v2, . . . , v2N} %! !/ % .. 1

!% # .# % % :=F=;

& % .I %

, # % ! ..N $! %

i(0) = 0, i = 1, 2, . . . , N , :=G=;

!" # # 42

N

i=1

i(0)2 = 1, :=GD;

% .# :=FE;

i(t) = 0 , i(t) = i(0) :=GA;

% % '

,

w(t) = (0, 0, . . . , 0, 1(0), 2(0), . . . , N (0)) . :=G?;

# % k ! ! .

%

{w1, w2, . . . , wk} 2 k 2N 2# {v1, v2, . . . , v2N} %[

w1 w2 . . . wk

]T=

[v1 v2 . . . v2N

]T. :=GE;

% , # %

,

= [dij ] =1k

m=1 wm(t)

11(0) 1N (0) 11(0) +

Nm=1 1m

1m(0)t 1N (0) +

Nm=1 Nm

1m(0)t

21(0) 2N (0) 21(0) +N

m=1 1m2m(0)t 2N (0) +

Nm=1 Nm

2m(0)t

k1 (0) kN (0) k1 (0) +N

m=1 1mkm(0)t kN (0) +

Nm=1 Nm

km(0)t

, :=GF;

i = 1, 2, . . . , k j = 1, 2, . . . , 2N # ! (:=FF; 1 :=G>;) wi(t) ! ! .! t

Mi(t) = wi(t) t :=GG;

%+ 0,ki ki k 1

0,ki =

1i (0)2i (0)

ki (0)

, ki =1i (0)2i (0)

ki (0)

, :=GC;

8: 1 . # 2 & ' /GALI0

% :=GF; !

(t) =1k

i=1Mi(t)[

0,k1 . . . 0,kN

k1 +N

i=1 1i0,ki t . . .

kN +

Ni=1 Ni

0,ki t

]=

=1k

i=1Mi(t)0,k(t). :=C>;

# ! m !# .

m k , # % #m .% :=GE; &, % i :=C>; ,

m,ki =

00

0m+1i (0)m+2i (0)

ki (0)

:=CB;

# % m , # ! % %

' % :=C>; % %

(t) =1km

i=1 Mm+i(t)[

m,k1 . . . m,kN

k1 +N

i=1 1im,ki t . . .

kN +

Ni=1 Ni

m,ki t

]=

=1km

i=1 Mm+i(t)m,k(t), :=C=;

# % m,k(t) .# :=C>; :=C=; ! % m,ki

2# . .% :=E>; !, . ! k

w1(t) w2(t) wk(t) =

1i1

!" # # 8%

! !

Sk =

1i1-*!

# % , &

% % % % ,%+ Sk , % .% :=C>; )' %+ ,%+

Sk ! 1/k

i=1Mi(t) .% :=GG; # ,

1ki=1Mi(t)

1tk. :=C?;

8) 1 . # 2 & ' /GALI0

0 % # .. Sk + !! . %+ kk %+ %0,k :=C>; # . t

. k , #

2 k N 9 !! . % % N !/(k!(N k)!)%+ k % % % N %

0,k

0,kj1,j2,...,jk = kj1 +Ni=1 j1i0,ki t kj2 +Ni=1 j2i0,ki t kjk +Ni=1 jki0,ki t , :=CE;

1 j1 < j2 < . . . < jk N 2# +#

.. 0,kj1,j2,...,jk %+ % , +#

j1m10,km1t j2m20,km2t jkmk0,kmk t = tk ki=1

jimi 0,km1 0,km2 0,kmk tk, :=CF;

mi {1, 2, . . . , N} i = 1, 2, . . . , k mi = mj ! i = j ' :=C?; :=CF; % , +# %+ 0,kj1,j2,...,jk , Sk % :=C=; )9 %+ ,%+ Sk % ,

0,kj1,j2,...,jk

%

# N % 0,k ! :=C=; !#

tk % % % ,

GALIk(t) 0 % %+

0,km = 0,k1 0,km k1 +Ni=1 1i0,ki t kkm +Ni=1 kmi0,ki t , :=CG;

# m %0,k % %+ !

# k m # !% %+

% 0,k 1 < m k 9 # m 0,km % .

2# % +# ..

0,km %+ % .. +# 0,k1 0,k2 0,km 1i10,ki1 t 2i20,ki2 t kmikm0,kikmt tkm, :=CC; ij {m + 1,m + 2, . . . , N} j = 1, 2, . . . , k m ij = il ! j = l ' , +# 0,km % ! tkm/tk = 1/tm

:1 < m k; % # t . , k k %+ ,%+ Sk , :=CE; . t

k

, :=CG; . tkm , Sk %+ %+ . tk % :=C?; %

!" # # 84

GALIk(t) ! 2 k N. :=B>>;

& % k N < k 2N 9!! . %+ % % N

% 0,k

0,kj1,j2,...,jkN ,1,2,...,N = 0,kj1 0,kjkN k1 +Ni=1 1i0,ki t kN +Ni=1 Ni0,ki t ,

:=B>B;

1 j1 < j2 < . . . < jkN N 9 # k N 0,kj1,j2,...,jkN ,1,2,...,N ! # N 0,k % % .

)'

N !/((kN)!(2N k)!) %+ :=B>B; % !, kk +# % % % N #ki i = 1, 2, . . . , N 5 , ji

0,ki t i, j = 1, 2, . . . , N 3 % ) /

0,ki i = 1, 2, . . . , N ,%+ % , , % % %+ % % %+ !!

. % ! t

, t % 0,ki #

kN 0,kj1,j2,...,jkN ,1,2,...,N t ,%+ N (kN) = 2N k ,3, %+ % 0,ki , % %

9 k (2N k) (k N) = k N ki % ,%+ % , ' .. 0,kj1,j2,...,jkN ,1,2,...,N%+ % %+ , 0,kj1 0,kjkN 0,ki1 0,kikN ikN+1m10,km1t iN m2Nk0,ki2Nk t t2Nk, :=B>=; il {1, 2, . . . , N} l = 1, 2, . . . , N il = ij ! l = j ml {1, 2, . . . , N} l =1, 2, . . . , 2N k ml {j1, j2, . . . , jkN} ml = mj ! l = j )' %+ ,:=B>B; , .. Sk ! ! t

2Nk/tk =1/t2(kN) )9 %+ ,%+ Sk , 0,kj1,j2,...,jkN ,1,2,...,N ! % !! 1/t

2(kN) ,

kN , .

, 0,ki i = 1, 2, . . . , N ' Sk GALIk % #

,

GALIk(t) 1t2(kN)

! N k 2N. :=B>D;

B) 1'! ', m > 0 ",', -+, ,"5', 1>-*!

, GALIk m m k m N % , # & % %

88 3 4' +' &

% % ,%+ Sk , % :=C=; ' %+ ,%+ Sk !

1/

kmi=1 Mm+i(t) % # #

1kmi=1 Mm+i(t)

1tkm

. :=B>A;

'!# % % m = 0 % % 2 k N !! . kk %+ % m,k % , ki1 ki2 kim im+1n10,kn1 t im+2n20,kn2 t iknkm0,knkmt tkm, :=B>?; il {1, 2, . . . , N} l = 1, 2, . . . , k il = ij ! l = j nl {1, 2, . . . , N} l = 1, 2, . . . , km nl = nj ! l = j ' % , Sk ) GALIk%+ , +# , :=B>?; % +

:=B>A; %

GALIk(t) ! 2 k N. :=B>E;

( % N < k 2N m > 0 , % , # # !

% m = 0 % m < k N GALIk .% ! t2Nk/tkm = 1/t2(kN)m m k N %.

. %+ % ! tNm & % Sk GALIk .% # , tNm/tkm =1/tkN

$ GALIk ! ! k N # % ! k > N # % .

m :m N m k; % , & %.

GALIk(t)

2 k N

1t2(kN)m N < k 2N 0 m < k N

1tkN N < k 2N m k N

. :=B>F;

4 & .%

& !, , % !# #

%

, # % 9

! % %+ !

!" # # 8;

# !# 2#

% ! ! % !

%

" / !

.!! ! # #

! # # % Fermi, Pasta Ulam (FPU) 2# , Bose Einstein Condensation (BEC) . , A

%+ !, .+ !# ,# .#

J(t) = A(t)(t) :=B>G;

% # !#! (t) % 2n! 9 A(t) % T 1 2n 2n % # J %+

2n 2n %

J =

(0n InIn 0n

):=B>C;

0n n n % In % n n % " :=B>G; !% .# % 2# %

! % T .!!

# !

3 % %. :=B>G; % 2#

% .# ! % 2#

! ,

H2 =n

i,j=1

aij(t)ij :=BB>;

aij(t) % T 1

% %A(t)

9 % J :=B>C; %

J1 = J = JT :=BBB;J2 = I2n :=BB=;

JT % , % J )9 J

!% #

86 3 4' +' &

4 5 * %

" :=B>G; !,

(t) = (J1A)(t) :=BBD;(t) = P (t)(t), P (t) = P (t+ T ) :=BBA;

P (t) = J1A(t) = JA(t) ' % A(t) % % % Trace[P (t)] = 0 " :=BBA; % ! !

,+ ! n !# .

9 n

! nn % X(t) % +

:=BBA; ( # !

% X(t) % X(0) = In # % X(0) t = 0 ! Rn % X(t) ! (t) :=BBA; !, ,

(t) = X(t)(0), :=BB?;

(0) % ' % # % Z(t) :=BBA; !,

Z(t) = X(t)Z(0) :=BBE;

Z(0)

% % # % X(t) .% :=BBA;% J?=K

det[X(t)] = det[X(0)]e t0 Trace[P(t)]dt. :=BBF;

' Trace[P(t)] = 0 %+ # % X(t) % .

,

det[X(t)] = det[X(0)] :=BBG;

% 2# % .+ #

' # % X(t) ! :=BBA; %

det[X(t)] = det[X(0)] = 1,t 0 :=BBC;

! :=BBF; , X(0) = In Trace[P (t)] = 0 3 %+ # % :=BBA; % % % .

' det[X(T )] = 1 ! t = T T % % %P (t)

9 % X(T ) % ! % # !# + % ( % !!

!" # # 87

! (t) :=BBA; % (0) +

! t = T %

(T ) = [X(T )](0) :=B=>;

X(T ) = [X(T )] :=B=B;

T % % % P (t) 0 %. !+ .I , % P (t) % T 1 P (t) = P (t + T ) #% X(t) % , X(t) = P (t)X(t) ! :=BBA; :=BB?; )

X(t+ T ) = P (t+ T )X(t+ T ) :=B==;

% X(t+ T ) % # % :=BBA; %! :=BBE;

X(t+ T ) = X(t)X(T ). :=B=D;

' # :=B=D; t = T

X(2T ) = [X(T )]2 :=B=A;

!! ! t = 2T, . . . , T

X(T ) = [X(T )] , N. :=B=?;

4 + 5 *

& !, , % X(T ) # %+ # #

2# ! % ! .# %

, :=BBA; *% , ! , # # #

,

$ # # #

%

"> A )B2B% 2n 2n X

XTJX = J :=B=E;

8@ 3 4' +' &

J =

(0n InIn 0n

). :=B=F;

' X = J :=B=E; J :=BBB; :=BB=; 2n 2n % J % %

' :=B=E; J

J1XTJX = I J1XTJ = X1. :=B=G;

( % :=B=G; # % XT X1 % XT %! , ! # 0 )! , % J X

% / !! % XT X1 %

% XT % % X X1

% % % X i, i = 1, . . . , 2n # % X(t) :=BBA; I 0 i

# % X(t) % 2ni = 1i " % ! % X(T )

# % :=BBA; X(t) % ! % ! % X(T ) ' # .%

|X(T ) I| = 0 :=B=C;

! % % ! % .

+! ! +! ! ! +! !

&/%+ % ! 2#

:=BBA; ! % . 4 I

B ' % ! = 1 ! %, 1 = 1 " % ! = 1 % 1 1 % :;

= ' i % ! i > 1 %, 1i =

1i< 1 1i > 0 '% i % ! i < 1

1i =1i> 1 i < 0 : =A;

D ' % ! i C |i| = 1 :i = e, 2 = 1; +! ! i |i| = 1 :i = e; # ii = 1 % ! . % : =?;

!" # # 82

& =AI 9 ! |i| > 1 + . ! % +!I (a) 0 i > 1 0 < 1i < 1 (b) ! i < 1 1 < 1i < 0

& =?I 0 ! i C |i| = 1 :i = e, 2 = 1; +!

i |i| = 1 :i = e; % # ii = 1

A ' % ! i = e, R, = 1 : > 1; .3 I

:H; j = 1e

:H; k = e

:!H; l = 1e

% il = 1 jk = 1 :. %, +!;# % i = k j = l :. +! ! +!; : =E;

" % i, i i, j , k, l % % %

;: 3 4' +' &

& =EI 0 ! i = e, R = 1 : > 1; j = 1e

k = e l = 1e % % ! %,

+! ! +!

4# %, .

! !' %'

" # ! !# 2#

% % X(T ) ! .. :=BBA; P (t) ! % % T P (t) = P (t+ T )

)' % X(T ) n 1, . . . , n % , ! # 0 )! !% , nn / % W #

W1X(T )W = D :=BD>;

D =

1 0 . . . 00 2 . . . 0

0 0 . . . n

:=BDB; W % n n

1, . . . , n % X(T ) , n n % X !%! N

[X(T )] = WDW1. :=BD=;

:=BD>; :=BD=; :=B=>; !%

(T ) = [X(T )](0) W1(T ) = W1[X(T )](0) W1(T ) = D(W1(0)). :=BDD;

!" # # ;%

' % :=BDD; = W1

(T ) = D(0) :=BDA;

D =

1 0 . . . 00 2 . . . 0

0 0 . . . n

. :=BD?; :=BDA; i, i = 1, . . . , n %

X(T ) %+ :=B=>; , .. = W %+ # # ( (T ) !

i, i = 1, . . . , n % X(T ) % % % !

# %+ , % % i |i| > 1 ! ! % . (t) # . t

# % X(T ) % (1 = 2) !% :=BD>; %+ % X(T ) % D % , Jordan

D =

1 1 0 . . . 00 1 0 . . . 00 0 3 . . . 0

0 0 0 . . . n

. :=BDE;

& % !!

(D) =

1 11 0 . . . 0

0 1 0 . . . 00 0 3 . . . 0

0 0 0 . . . n

. :=BDF;

.

#

% % & % :=BDE; ,%

:=BDF; |i| = 1 ! i % D % %

& %. # JCFKI

;) 3 4' +' &

+!"A )B2B% P (t) = P (t+ T )! " "

(t+ T ) = (t) :=BDG;

# X(T )

1>0=!A # % :=BDG;

" t = 0

(T ) = (0). :=BDC;

)9 !%+ :=BB?; ! t = T

(T ) = X(T )(0). :=BA>;

&+ % %

(X(T ) I)(0) = 0 :=BAB;

% (0) : :=BDG;) %

! ! :=BAB; %

det (X(T ) I) = 0, :=BA=;

:=BDG; % % % X(T )

)' % # + ! %

% Floquet ! !# %

.I

+!"A )B2B) $ " # " %&!

" %' # t

1>0=!A :=BB?; ! % t

(t) = X(t)(0) (t+ T ) = X(t+ T )(0). :=BAD;

!" # # ;4

'.% :=B=D; :=BDC; :=BA>;

(t+ T ) = X(t+ T )(0) (t+ T ) = X(t)X(T )(0) (t+ T ) = X(t)(T ) (t+ T ) = X(t)(0) (t+ T ) = X(t)(0) (t+ T ) = (t). :=BAA;

% :=BDG;

% % % X(T )# :=BDG; %

.

:=BDG; % .

% X(T )

' , % Floquet

(t) = f(t)et :=BA?;

% ! .% :=BDG; = 1T ln f(t) = f(t+T ) :T 1

;( + Floquet ' :=BBA; , :=BA?; " % n i, i = 1, . . . , n% % , :=BA?; ! :=BBA; ,+

! n , :=BA?; # n !# .

% (0) = f(0) % ) ! !, !

, :=BA?;

(t) = c1f1(t)ea1t + . . . , cnfn(t)eant :=BAE;

i = 1T lni, i = 1, . . . , n ( % ,! Re(i) 0, i = 1, . . . , n

$,# = i = 0 :=BA?; = 2/T % # .

% # / = m1/m2 :=BA?; ,+ # % T T

= 2/ 6

% # % T M

0 . :=BBA; % (t) =(t+T ) (0) = 0 , (t) = (t+T ) t = 0 % (0) = (T )

;8 3 4' +' &

:=BB?; ! t = T

(T ) = X(T )(0) :=BAF;

(X(T ) I)(0) = 0 :=BAG;

) , =CB =C= :=BBA; T 1

% X(T ) %

# # %, ! %,

', # !# 2# %

! 2# % , :=B>G; (t)! % (0) + % i, i = 1, . . . , n % X(T ) , & % . ! (t) % ,

, fke

kt k > 0 ) ! : % ; ! % X(T ) % % % 6

% i = 1 i = 1 % %. % % , , f2(t)+ tf1(t) % % ( ! , % ,! .

% !

' ! % !

(t) ! 2# :=B>G; ! % ! (0) % ,! ! t % % ) % % %

' % % 1 1 % % % % % , % . !

! % ,!

(!!A )B2B% ( # # # )*

+ ,' " ( + - ." ."*

( # ." # ."

+ + - " " " /-

# 0 + ,' " )+ 0

# "

1 .

!" # # ;;

5 & Krein

6 "+ .

% !

& ! !, ! 2# ! % .!!

% !

# # !# 2# ( %

, # % X(T ) # X(T )+ ! % ! !

& !, % % % Krein % , # % # %

! 2#

% ,% ! # %

% !% ! % % 2#

# % % %

)' ! 2# :=B>G; A = A0 % T 1 2n 2n % 6 % % % % %

.

% 9+

& % # / %

2# ! ' ! 2#

,

J = (A0 + A1) :=BAC;

2n 2n % A1 % T 1 || 1 + ', :=BAC; % !

%X(T ) % , %! % X0(T ) % ! ! 2# :=B>G; L

% X(T ) ! :=BAC; %%! % ! % % X0(T ) ' ! !|| 1 % % % , !% " # % % %I 0 2#

, :=BAC; %

% ,!M

. % ! 2# ,

J = A(t, ), :=B?>;

;6 5 4' Krein

2n 2n % A(t, ) % T 1 % .

% ) % :=B?>; .

%

X(t, ) ' % ! ! ! 0 ! 2# :=B>G; %

" % # %I %

% ! 2# :=B>G;M

" % A #

% #

% X(T, 0) ! = 0 X(T, ) ! ! 1$ ! . I

I 0 I& % %

% !, =C= % % % :

=F;

& =FI %

II , I'# 1 = 2 3 = 4 % % % : =G:a); 3 % % % %

:; : =G:b); % ! . :; : =G:c); (

% % + !

III ( 1 1I'# % 1 = 2 = 1 1 : =C:a); ( % % % :; :

=C:b); % . :; : =C:c);

%

% ! ! 2# :%

!" # # ;7

& =GI 3 1 = 2 3 = 4 % % % :b; 3 % % % :; % :c; . :! ;

*; % % !% % !

% # :% ** ***; %

% % & # % ! # %

!% % Krein % .

"> A )B%:B% 2$ /-

/- # #

9 % ! ! 2# :=B?>; % %

: % %

.

% #; %

2# 1 - % % % % % !

% # ! % 2#

& % % # .,% !

2+ % ,

! # % !% "

! %+ %

;@ 5 4' Krein

& =CI :a; 3 1 = 2 = 1 :b; % % % :; % :c; . :; % ! % 1

! I

"> A )B%:B) + < x, y > + x y

< x, y >= i(Jx, y) :=B?B;

J 2n 2n (

0 InIn 0

)! x y 2n 1 !

i2 = 1 .# 3 (, ) 3 " + +

# % X(T ) ! ! 2# % m % h ! %

I

I < h,h >> 0I ( + # %

II < h,h >< 0I ( + %

III < h,h >= 0I ( + % ! h % # % #! % % .%

!" # # ;2

# Krein, Gelfand Lindskii JCFK I

(#A )B%:B% 2$ /- "

" 3 + "

"! # - # "

&, % Krein ! # # % % % !% % % % %

2# !% !

, % !% % % %

2# % !

.67 *! &% 8

' 51 7 8 $

19 # % !% ,

: ! ; !, #

. , !

)' % ! ! # %

% P ! V % T % !

PV = nRT n % :mole; % R = 8.31J/mol K % ! % : % % Kelvin 0 K %; 9 % %

% / !

. : ., ! V % % ;

! # !%

%

9 % !

%

% ,% % %

% . $# .%

! % M

9 S. Carnot 1824 !! ! ! , % .

: Carnot; R. Clausius 1865 !! % !, I $ %

.!% , %

%M

$ %

%M " !!

# % % , ! % % !

6: &67 '! 4'+ 8%'

% % ! ! ! %

,% % !% !

M

% # % # P. Duhem E. Mach .I'!! % ,!#

% + % # # ,

% !! % ! # % %

! # 4

%

# % %

%+ # % %

9 # !% % W. Herepath J. Waterston 1845! !% A. Kronig 1856 / % + @ 9 J. C. Maxwell % ! # . %, ! f(v) v

N % % %

%

f(v) = 4N(

m

2kBT

) 32

v2e mv2

2kBT :=B?=;

m % + % kB = 1.38066 1023JK1 % ! Boltz-mann T % % % Maxwell Boltzmann . % % % .%

S = kB log10W :=B?D;

% S % ! Boltzmann kB % ! W # ( %, .%

Boltzmann ! f(v)

% f/t = 0

:=B?=; Maxwell JCCK

( % Boltzmann !!

% # . t

&! # :recurrence theorem; Poincare ! ,! , % !!

% %

&, # % !

% %

, %

! , , %

% .%N 9 J. Loschmidt % ! / % !% ! !

%

% / % .#

% ! t t : vi vi;

!" # # 6%

9 Maxwell Boltzmann ! ..

% % % % % % %

&! Maxwell ! % % . %

9 Boltzmann % % .! .% ! !, !!

% % !! /

% % # . %

9 # , . # , &

! # % % % 2

20 # " 2 %. % # / % , ,

% % !# % %

%

% % N " ,%+

. % % ! 1 !

% % 9 Maxwell Boltzmann % % % ( ,%

% % ,%N

51 3,

& # % &

!

% ! Maxwell Boltzmann J. Gibbs !! :microcanonical ensemble; ( # + # % %

7! 2#

% ! %

% % ! t % ! , %+ 2#

'# %+ # %

. )' #

% ! ! ,% : microcanonical ensemble;

!, ! % 0% ! M $#

!%+

M

9 Boltzmann ! # #

! !

.! . ,# #

.. 0 !% .

6) &67 '! 4'+ 8%'

%

' A % = (x) % x = (q, p) % % :ensemble; % % !% A ,

< A >=1

N (E)N (E)i=1

(xi) :=B?A;

xi % i 1 ensemble N (E) # , ! E

% , A !%+

< A >=1V

(x)(x E)dx, :=B??;

V ! ! ! , % !% A ! x(t) !, #

A = limT

1T

T0(x(t))dt. :=B?E;

' % x(t) % % %

! , < A > A %N

( . ,# : microcanonical ensemble; # % ! ( < A >= A # % ! +

% 0 ! 2# %

% ! , ! ! %!

% % ! ! % %

:=B?E; ! 2# 0 !

+ '! 6 %+ % ! .

#

0 % %. %

! ! # %

.) . , # ( ' '!

% J. von Neumann ! G. Birkhoff . ! % ! : ;

% # , %

% % % %

!" # # 64

# . !

6

!!M % % % ! # # I (a) % (b)

# , ( . # !

! , ! !, %! +

!!

)' ! % Ya. Sinai # # ,#

% ! %

!! % #

! # Kolmogorov Arnold Moser (KAM) % : ! # ; ,%

!

% # ,

@ # # % N 1 # , % ! % % % 1

! % )9 N ,%

.,%+ !#

& '! 6 )'

% % #

!! # N

# ) 9

9 ! % # %

@ ! .! % ! , %

% # (

( , ,# ! %

% % , !! 0 ! % % !

! !

% ! . %

% % !

! % % & %

# % # %

)9 Rudolf Clausius !! % 1865 % # & % !

! % # , ! % % ( %

Clausius % % %+ . .%

68 &67 '! 4'+ 8%'

# # % !, %

% , ! !,

& L. Boltzmann % R. Clausius !

, :=B?D; ! % ( .% % % . # ! , ,

4 I dSdt 0 9 % !% % & # ,

, ! ! , # !#

!.#

7%! ! :BGFE 1 BGFG; @ $% Yale ($J. W. Gibbs & ! !/ , . .% % ( Gibbs !, , # % 0 !

% Gibbs !, %+ % % % % (

% % % %+ ! .%

% % , #

% ! % N $ I ( %

% % # %+ %

& % 1870 Gibbs !! % , ! !/ % : % Gibbs;

S = kB

i

pi log10 pi :=B?F;

%

pi ' %+ pi = 1/w .%:=B?F; %+ Boltzmann :=B?D; 0 ! Boltzmann :=B?D; :=B?F; , % Boltzmann 1 Gibbs

( % % % % !, %

% + .

! L %

! % N : . Avogadro NA 1023;# ! V % % ! % #

% % = N/V '

2# %

%+ E N # ! E/N = E % ! 2# N #

%

)9 % !, =E

, % %+

!" # # 6;

!# # # ! # #

Lyapunov ,% ! !% ! % Kolmogorov Sinai :Kolmogorov Sinai entropy; ! ,!

hKS =

i:i>0

i :=B?G;

i > 0 % Lyapunov JEEK & . :=B?G; , ,! % Kolmogorov Sinai hKS

i:i>0

i JD>K #

. Pesin % % # JEEK ! !! ,

9%! :

& # ,% + #

! # +

D , &! # ! #

! !% 2#

Legendre % 7!+ & ,

# # # , !

Poincare % % 2# '!/ 2# #

# , %, % KAM )9 / ! % # 2#

, ! %

& % ! % % % .

% ! # # %

0 # .%+ % %

% ! 2# # # %

% A , L #

, % ! # Lya-punov %+ !# # $# . % ,! % , ! 2#

%+ % SALI % !! . !

! Lyapunov &, SALI .. # !# .

!!

!!,% % # % % ,

!! ! !

, # Lyapunov !! 0 3%

66 9+!&&

'! . k :GALIk; !% SALI ( %% !% # #

% 2# N # %

, !! SALI % ! ! % % GALIk

! k > 2 !# .

% & . GALIk !%+ . !

k

0 ! # GALIk # , # .. k :!;

.# # , $

M ! # , 2#

% 2N 1 k % % 2N # GALIk % % )9 . %

!# .

% !% !# .. &

% ! GALIk %+N # % ! k > 2 , % % ! k > N # , # : ;

' . % % %+ ! GALIk

% # .

# #

Lyapunov : . :=F>;) % ! # GALIk 2 k N # ! # ! N < k 2N % # : . :=B>F;) 9 .

k N # m %

% , % %

& !/ # # 2#

+ % % Krein # % ,% % % &

,% ! &

' , %

% , % ,

!

"# $ % &

''()

$ !

-# $

SALI GALI

# #

! ,$ SALI

! /# $

+ ( " GALI 2 3 ! &

/ '!

- !: SALI / (6

)9 % , = =F % SALI % ! # # % # !

& !, .! !% SALI !%

EF

6@ $ . !'' SALI # # *'6#

% ! # # , )'

% % # 2#

% % ,! % %

! ! # , J?AK

)' # 2# 2 # % 2

# H % F .

H .

[H,F ] = 0, :DB;

[, ] # ! Poisson : ==D; & %%+

H = h, F = f , :D=;

h f % ' ! 2#

, # % ! %

fH = (Hpx ,Hpy ,Hx,Hy),

fF = (Fpx , Fpy ,Fx,Fy). :DD;

, H F % % .#

Hamilton JAK ' 4 1 # , 2 1 % %

H = (Hx,Hy,Hpx,Hpy),F = (Fx, Fy, Fpx , Fpy), :DA;

x y ! px py % +!% #

% # !#! % : Hx H/x;"

fH ,

fF :

H,F ; % !# .

.%

.

% H,F % )' %

fH =fH

fH, fF =

fF

fF+% H =

HH

, F =F

F:D?;

% 4 1 # , % .% % ! !# #

H, F = fH , fF =HxFy +HyFy +HpxFpx +HpyFpy

H F:DE;

% ! # H = fH F = fF , # !

)' # :D?; % ! ..

$ & $ :

'# 62

v1 !/

v1 = a1fH + a2fF + a3H + a4F :DF;

a1, a2, a3, a4 R

t & % 2 1 ! % % :1 1

; % JCAK :#

! fH H # ;

!% , % % ,

, n1 n % /

$ % % # 2# = # %

% = 1 v1 % , #

! fH fF # .% :DF; a3 0 a4 0 # a1 a2 % ! , L ! ,

# 2#

Van der Waals J=BK

H(x, y, px, py) =12(p2x + p

2y) E(x2 + y2) +A(x6 + y6) +B(x4y2 + x2y4), :DG;

E A B % ! 0 B = 3A E R 2# :DG; % # % %

, JDGK

F (x, y, px, py) = (xpy ypx)2. :DC;

& ! A = 0.25 B = 3A = 0.75 E = 1080 , !%+ .. #

a1 a2 a3 a4 .% :DF; -% # a1 a2 ,

# a3 a4 % )' ! ! +

DB $+ DB(b) % ! ! % |a3| |a4| t1 " .! =GA , k .%

2# N # %

DB % v1 % , # % .% , # ! fH fF

!% ! .

#

% ! , %

a1 a2 SALI : , = =?; # ! ! ,

% , D=:a; %+ .. SALI ! % x = 0.6 y = 0 px = 0 py = 1.99416 Poincare D=:b; "

7: $ . !'' SALI # # *'6#

& DBI ( .. # a1 a2 a3 a4 ! # a1 = 1 a2 = 1 a3 = 0 a4 = 1 : :DF;) (a) a1 : !; a2 :! !;(b) |a3| : !; |a4| :! !; ! % ( 2#:DG; x = 0.6 y = 0 px = 0 py = 1.99416 & (b) %+ % |a3| = 0.107 t0.962 |a4| = 1.067 t0.995

% v1 = fH + fF + F v2 = fH + fF + H ( .. v1% DB $ ! !

v2

- !: SALI / (6

)9 !/ =? ,% = 3% '! :SALI;% % ! !! !

, % 2# ( %+ , ,

SALI #I '% SALI # ! % ! # % SALI 0 # t ! ! ! Lyapunov %# # # % %. , =?

,% = ,+ SALI e(12)t 1, 2 ! Lyapunov (1 > 2) 2# % # ! % 2#

% : , ?;

! %!I

$ & $ :

'# 7%

& D=I ( .. SALI (a) ! % ! % Poincare y = 0 (b)

" 2# N = 2 %

H2 =12(p2x + p

2y) +

12(x2 + y2) + x2y 1

3y3 :DB>;

! H2 = 18 Henon & Heiles JA?K

2# N = 3 %

H3 =12(p2x + p

2y + p

2z) +

12(Ax2 +By2 + Cz2) xz2 yz2, :DBB;

A = 0.9 B = 0.4 C = 0.225 = 0.56 = 0.2 % J=B =DK 0 !#

)' 2# #

%! # % % # # , # Poincare :PSS; J?AK ( % ! 2# = # % , , Poincare % %2 & 3 # % # , % . , # Poincare % N , # % % % )' .%

% Poincare # )9 % %

! , % %

7) $ . !'' SALI # # *'6#

0 , SALI = D % . # ! & DD:a; % ! 2#

:DB>; % %+ x = 0 " % ! % % % :; % #

, # DD % % ,%+

% ,

( .. SALI # # %+ DD:b; &

% ! :! !; SALI , # % : !;

SALI ,% !% t 800 # , % ! :1016; % ! t 800 = !% % ,

% - ,

% % # %+

SALI , #, 108 t = 400 $ , + SALI ! ! 2# # #% :DBB; D?

& DDI :a; ( Poincare ! x = 0 y = 0.1 px 0.49058 py = 0 x = 0 y = 0.25 px 0.42081 py = 0 %! Henon & Heiles :DB>; ( ! % # % , % PSS :b;( .. SALI! :a; ! % ( !! % ! # !

" w = (dx, dy, dpx, dpy) !

$ & $ :

'# 74

DD % w1(0) = (1, 0, 0, 0) w2(0) = (0, 0, 1, 0) ! ! # !% , SALI (

%+ !I

'# % # #

DD # w1(0) = (0, 1, 0, 0) w2(0) 0 , +! w1(0)w2(0) !%+ T % SALI ! . ! #, 1012 ! T .

! ! # '! w2(0) , I $# % , Poincare :DD; .

!% # w1(0) = (0, 1, 0, 0) w2(0) = (0, cos , 0, sin )

& DA:a; %+ T

!% [0, ] )9% ! !% = 0 = T = 0 , % ! # , ( ! T

% % w2(0) %

& % ! w2(0) % % ! .

# !% : , !,

:DA); $ % ! (0, ) T + [400, 500] 9 ! !% T ! 1000 w2(0) % ! % % % : DA(b); DA % T .

! #

! ! Lyapunov 1 : =E; % ! % SALI !% , DE + .. L1(t) L1(t) 1 ! t : DE(a); SALI : DE(b); ! % :DB>; & t 1900 SALI , 1016 !% % ! L1(t) % :t 1900; % ! DE:a;@ % %! ! ,

t 1000 SALI 108 ' !% , SALI ! Lyapunov , ,% !% ! : % ! ; !! % ! ,

! Lyapunov !% : % %;

( , SALI # % #2# = # % !% JGK '% %. +!

% % , # %

.% % %+ = .

%

78 $ . !'' SALI # # *'6#

& DAI 9 T + SALI ! !% 1012 % DD # w1(0) =(0, 1, 0, 0) :a; w2(0) = (0, cos , 0, sin ) :b; w2(0) % ! % %% 9 T + :a;

!% :b;

% i % !

& D?I ( .. SALI ! % ! :! !; % : !; 2# D # % :DBB; x = 0.01054y = 0.05060 z = 0 px = 0 py = 0 pz 0.11906 x = 0.07310 y = 0 z = 0 px = 0.07695py = 0 pz 0.06760 %

= # % ' , w1(t) w2(t) %!# .

.

% # !

$ & $ :

'# 7;

& DEI ( .. :a; L1(t) L1(t) 1 t 1 ! Lyapunov :b; SALI ! x = 0 y = 0.01597 px 0.49974 py = 0 # % 2# :DB>;

& # # 2# N # % !

% , )9 .! # =GD ,% =

% .. w(t) Lyapunov # ! 1 2 2n % % Lyapunov ! )9 !/ , = Lyapunov # # 2# % !% % +! %

% % J?AK ' ..

w1(0) % !! % , ,

w1(t) =2ni=1

c(1)i e

Lit ei , :DB=;

c(1)i % ! !% % Li ei .

# !

x(t) # , % .%

0 # % .. SALI # % Lyapunov Li # ! !! Li i 3 , w1(t) ! :DB=;

w1(t) c(1)1 eL1t e1 + c(1)2 eL2t e2 c(1)1 e1t e1 + c(1)2 e2t e2. :DBD;

76 $ . !'' SALI # # *'6#

:DBD; ! !% !! %

w1(t)w1(t)

c(1)1 e

1t e1 + c(1)2 e

2t e2

|c(1)1 | e1t= s1 e1 +

c(1)2

|c(1)1 |e(12)t e2 , :DBA;

w2(t)

w2(t) c(2)1 e

1t e1 + c(2)2 e

2t e2

|c(2)1 | e1t= s2 e1 +

c(2)2

|c(2)1 |e(12)t e2 , :DB?;

si = sign(c(i)1 ) i = 1, 2

& DFI :a; ( .. L1(t) ! x = 0y = 0.25 px 0.42081 py = 0 # % 2# :DB>; :b; 9 SALI % :! !;

e1t : !; ! 1 = 0.047 9. % !

0 !% # SALI : =F; ,

.# :DBA; :DB?; %

SALI(t) = min w1(t)w1(t) w2(t)w2(t)

c(1)2|c(1)1 |

c(2)2

|c(2)1 |

e(12)te2. :DBE;9%+ c . .% #

!! e2 % %

SALI(t) c e(12)t . :DBF;

: . :=D?;) ( % SALI ! % !! # %+ , !

$ & $ :

'# 77

# Lyapunov

0 !. # 2

# = # % % % %

Lyapunov :1 > 0; # % % 2 = 0 )' .% :DBF;

SALI(t) c e1t . :DBG;

& DF(a) %+ ! % ! Lya-punov L1(t)

t ! 2# =# % :DB>; 9 ! Lyapunov , , . # .. !

! ! . !!% :t 10000; % 1 0.047 & DF(b) %+ SALI ! % :! !;# ! % ! t )9 ,%

:DBG; : !; !

% , SALI 1016 ! t 800 ) # SALI % , t 400 SALI !% 108

9 2# # # %

Lyapunov 1 2 )' ! !! , SALI :DBF; % !# 2 # # Lyapunov 0 # :DBF; !%+ # Lyapunov 1 2 ! %

2# # # % :DBB; ,+

Benettin et al. JG CK " + :DF;(a) ! %+ % Lyapunov L1 L2 # !% 1 2 )' ! t 106 ! ! H3 = 0.0141 % 1 0.0310 2 0.0064 2# DG(b) % SALI :! !; % % :DBF; : !; '% % , SALI !% !! % !! !

Lyapunov % ! ! !%

)' DF DG % %.

:DBF; ! , SALI

% 6 (;

9 SALI , % % ! % , 2# '#

!

!, # # ,

7@ $$ ''+# '6 *; #

& DGI :a; ( .. ! # # Lyapunov !

1 2 ! % x = 0 y = 0 z = 0 px = 0py = 0 pz 0.167929 2# # # % :DBB; :b; 9 .. SALI % :! !;

c e(12)t :!; ! 1 = 0.0310 2 = 0.0064 c = 101 9 . % !

% % !

& DC + , # Poincare x =0 2# = # % :DB>; )9 % , ! % ! % #

! % %. SALI % , % !I

% py = 0 ' 5000 , ! % !%+

SALI ! % " + DB> %+ SALI

y ! # !t = 1000 : DB>(a); t = 4000 : DB>(b); & % % ! # ! SALI % $ SALI ! : ! 104; % ! % % 6

SALI : 1012; % # % % &!% DB> DC %

$ % SALI ! 104 1012 % %

$ & $ :

'# 72

SALI :1012 < SALI 104; ! ! ! t = 1000 DB>(a) % ,

! ! ,%+

9 ! sticky % % ! ! # %

,% % # .# % '# % #

,% .% SALI % t = 1000 % 104 SALI 1012 # % . % sticky # %! % N

)' !% DB>(a) DB>(b) !% ! , % 1012 < SALI 104 ! t = 4000 DB>(b); # SALI ! % ! % %- % % #, %

, $

. . SALI ! #, % 104

& DCI ( Poincare x = 0 2# = # % :DB>; &+ % . py = 0

# DB> !

y 0.1 % % % ! %

SALI > 104 % % % % % Poincare DCN ' % y = 0.2088 ! SALI :> 0.1;

DB> # ! % SALI < 109 ! t =4000 % % % ! %

@: $$ ''+# '6 *; #

% ,% % ! !

Poincare - ,!

SALI % / % ,

',+ %

% Poincare + ! # ! SALI # ! , # & DBB:a; + %

! 2# = # % :DB>; # !

# # Poincare 9 # , %+ !% ! , SALI ! t = 1000 %+ , # .I SALI 1012 ! 1012 < SALI 108! 108 < SALI 104 ! SALI > 104

' DBB(a) % ) ! # % % ! % % & . #

# % % ! ! # !

sticky .%+ , # DBB(a) % % ! % % ,%

Poincare DC ! ! y 0.1 py 0 $ DBB !% ! t = 1000 % ! /# # ! ! !

:CPU; % !, # SALI % ! # # + DBB:a; % % % % %

G 1 C # ! Pentium 4 2GHz

' % + . %+

SALI % % # ,Poincare # ! DBB:a; = # CPU : ! Pentium 4 2GHz; ! ! # # 375750 , % py 0 t = 1000

0 2# D # % :DBB; % Poincare %A % = # %

% SALI % % % . 0 ! ! Poincare + , # :+% % ; !

SALI t = 5000 % . % % 2 1 4 1 Poincare

0 ! DBB(b) %+ ! ! H3 = 0.00765 y = 0 py = 0 :! z = 0 pz 0; 2# :DBB; #

$ & $ :

'# @%

& DB>I 9 SALI ! :a; t = 1000 :b; t = 4000 ! 2# = # % :DB>; % py = 0 Poincare DC

y !

% % % 2# = # % $ #

! : ! #; % : #; # %

sticky : ! ! #; % ' DBB ! ! ! #, ! ! # # # %

',+ % !% # ! 2# :DB>;

# # , # Poincare I :B; # % SALI < 108 ! SALI > 108 :=; #) % DBB(a) % JA?K ! % = # % :DB>; " ! #

+ DB= ' DB=(a) , # Poincare % SALI % % : #; % # # 55% %

Poincare %! ! %

, Poincare , : % SALI % %+ % ; + ! # & DB=(b) %+ % ! # # :DB>;

! H2 !% J?GK $ ! H2 = 1/8

% % % DB=:a;

@) $$ ''+# '6 *; #

& DBBI $ ,# # SALI (a) Poincare x = 0 2#= # % :DB>; ! t = 1000,H2 = 1/8 (b) y = 0 py = 0

Poincare z = 0 2# D # % :DBB; ! t = 5000,H3 = 0.00765 = %+ SALI 1012 ! 1012 < SALI 108 ! 108 < SALI 104 ! SALI > 104

& DB=I (a) " # # (SALI < 108)

Poincare 2# 2 # % :DB>; ! ! H2 = 18

t ( ! % % # ! ! 2 (b) 9 ! # # !% J?GK ! :DB>;

! E $ ! E = 1/8 % :a;

$ & $ :

'# @4

' !% Poincare % : DB=(a) % #;% ! ! # Poincare # CPU % % ! 0 !

!, Poincare DBB:b; % % %68300 CPU ! # # !! t = 10000 % ! % % % 2100 N

( SALI % % % ! . . Poincare # # L ! , :DBD; % %+

SALI % (y,H2) ! 2# = # % :DB>; '! x = 0 y [0.3, 0.4] py = 0 H2 = [0.08, 0.18] !%+

% SALI ! # t = 2000 " % % JDAK %

FLI JDEK ! !% . DBD # ! H2 2# . 0.08 0.18 y [0.3, 0.4] ! % KAM

:% %; ,% ,%

& ! ! H2 = 0.18 ! ! ! ,!

$ - !: SALI

@8 $< . !'' SALI =!# # '6#

& DBDI 9 SALI ! ! 500 500 , # # % (y,H2) ! 2# = # % :DB>; ! t = 2000 ' SALI < 1012 %+ SALI [1012, 108) %+ ! SALI [108, 104) %+ ! SALI [104,2] %+ ! #

= 0.25 % # % !% % , # Poincare

= 0.24 % Poincare (x, px), y = 0, py > 0

DBA % % % 7 % '! % ! x = 0, y = 0, px = 0.77, py 44.71473: ! :DBA;) # SALI , % !% % ! Lyapunov

)9 DB? ! SALI % I 0 ! DB?:a; t 100 % SALI % 104 # ! t 700 1012% .. ! Lyapunov % # % , % ! t > 100 !% % : DB?:a;)

+ # >=? % % SALI % DBA !

,% % !% JB?K

$ & $ :

'# @;

& DBAI % % Poincare (x, px), y = 0, py > 0 ! (x = 0, y = 0, px = 0.77, py 44.71473) %+ ! 9 :DBC; % = 0.24 E = 1000

& DB?I (a) ( .. SALI ! ! ! 2#:DBC; = 0.24 E = 1000 (x = 0, y = 0, px = 0.77, py 44.71473) (b)9 % ! Lyapunov

t ! %2# ! (a)

' "/ SALI 89 .+% ;"

% /

#: & 6 7 /5+ 7 90 & 9$ 5+,

" ! !, % SALI % %!! ! . !

@6 $) "'& SALI 89# +# ' ;" + %'

, # ( % SALI % ! # % .. % !%

+ , #

#

. :k > 2; , !% ! GALIk ! % , = ' ! # # % :N 1; ! GALIk % % SALI % ! % # !# # # #

SALI ! ,% ! , % !! ! .% ! ! Lyapunov Lmax % % ! ! Lmax % :

; % ! + SALI ! !% # DE DF DG $ ,%

!! ! SALI # .

# # Lyapunov % !!%+ % ,1 2 " % % . % !% SALI '% ,# # , ! > 1 B test !! Gottwald Melbourne JA>K

& JFFK SALI !% ! Fast Lyapunov Indicator(FLI) JD? DEK ,# !# ! # :stretchingnumbers; JCDK &! JCDK % SALI !% , FLI ! ! ! ! % % FLI )' SALI FLI % !! SALI % ! % % , #

FLI . ! ! ! ,

' % ,# ! !

SALI : DBB; .

# % FLI % FLI +# . )9 .! JFFK ! %

SALI % !! , # # ! ! % % " SALI % % ! # . % %

,# % ,%

& Sandor et al. JFAK !% ! . Relative LyapunovIndicator (RLI) SALI % % , # % ( RLI % , ! # Lyapunov !# # #

$ & $ :

'# @7

! % # ! , . Lmax

, JFAK "%+ ! RLI % ..

# : ! % ; # !

SALI % # , %

$ , Gottwald Melbourne JA>K !! .! , + > 1 B Test( % ! % ,% ) % !

% .# %

0 1 1 test % ! SALI ! % % !# .# % #

% .# # ( ,! 0 1 1 test %+

!# #

p(t) = t

0(s) cos((s)) ds , :D=>;

(s) % ! % %

(t) = kt+ t

0(s) ds , :D=B;

k % 9%+ !

p(t) %

M(t) = limT

1T

T0

(p(t+ ) p())2 d , :D==;

K(t) =log10(M(t) + 1)

log10 t, :D=D;

%+ %

K = limtK(t) . :D=A;

( %+ % ! ! K = 0 K = 1 ( ! # !, ,! # JA>K

0 !% 0 1 1 test SALI ,+ ! 2# Henon & Heiles :DB>; DD )9 %. % # # SALI !%+ , t 800 : SALI 1016; ! t 400 SALI 108 ,!

0 1 1 (t)

(t) = y(t) + py(t) , :D=?;

@@ $) "'& SALI 89# +# ' ;" + %'

: :D=B;) # ! k % k = 1.7 % JA>K ( ,! 0 1 1 test % % ! # (t) # . ! M(t) ! T :D==; % K K(t) ! t : :D=A;) & ! T = 90000 !%+ K(t) ! 10000 ! % ! # t = 105 $ % T t ! ! K(t) JA>K , K(t) ! ! ,% % % . # :

DBE;

& DBEI ',! > 1 B test ! DD ( K(t)

# # K = 0 ! ! :!!; # % K = 1 ! : !;

' K(t) . ! : ! DBE; % % , K = 1 # % % K = 0 ! ! :! ! DBE; + #

% , ! , # %+

> 1 B test ! % !

% / % SALI ( , ,% # ! ! K # + !

# :D=>; 1 :D==; ! .# +

! ! # !!% % :D==; :D=A;

% DBE

$ & $ :

'# @2

) *:/ .+% 3 *+/#

0GALI1

#( $,

% $+ 2 3 ;, ,

0 , GALI 2# # , %!I !

Henon & Heiles JAAK !, 2#

= #%

H2 =12(p2x + p

2y) +

12(x2 + y2) + x2y 1

3y3, :D=E;

# # %

H3 =3

i=1

i2

(q2i + p2i ) + q

21q2 + q

21q3, :D=F;

JC ==K & !, ! H2 = 0.125 H3 = 0.09 1 = 1 2 =

2 3 =

3 + , GALIk !

, k # # # ! #

A) "*+,# $!" 2 /"+, */

. # % 2 # % :DB>; x = 0 y = 0.25 px = 0.42 py = 0 & DBF(a) ..

! Lyapunov L1(t) L1(t) + !

# !!% % : , ;

L % !! ! Lyapunov :LCE; 1 ! t 105 % 1 0.047)9 !%+ 2# # %

LCE 1 , ! % 2 = 0 * % 3 = 2 4 = 1 :=F>; % !, .. GALIk ! %

GALI2(t) e1t, GALI3(t) e21t, GALI4(t) e41t. :D=G;

& DBF(b) %+ GALIk k = 2, 3, 4 !

t %+ t ! % #

:D=G; % GALI2 GALI3 GALI4 % % 1/ ln 1021/ ln 10 41/ ln 10 % DBF(b) % # % ! 1 = 0.047 !!%+ % !

GALI ( ! . # # #

2: $, '# + 2 & ' /GALI0

& DBFI (a) ( .. ! Lyapunov L1(t) :! ; L2(t) : ; L1(t) L2(t) :! ; ! %

x = 0 y = 0.25 px = 0.42 py = 0 2# 2# % :DB>; (b) ( .. GALI2 GALI3 GALI4 % 9 + !

! e1t :!!; e21t : !; e41t :! !; ! 1 = 0.047 9 .

% !%

,%+ % GALI4 / # % % , %

)9 , ! !, E? GALI2 % SALI !% , # DBG %+ , .

GALI2 SALI DBF

t 9 % %+ t 300 , , , %

! :1016; ! % . 105 : DBF(b);

# , GALIk ! % ! 2#

:DB>; :=B>F; % 2#

(N =)2 % GALI2 , # GALI3 GALI4 % # ! % .

m % ,

% "# ! % ! :DB>; % % !

# :=B>D; % GALI ,

GALI2(t) , GALI3(t) 1t2, GALI4(t) 1

t4. :D=C;

)9 .! =D 2#

$ & $ :

'# 2%

& DBGI ( , . GALI2 SALI DBF

t

& DBCI ( Poincare (PSS) %+ (a) x = 0 (b) y = 0 !

x = 0 y = 0 px = 0.5 py = 0 ! 2# Henon &Heiles :DB>;

% # % # % , Poincare (PSS) J?AK& = # % , % %

% # % ! : % = 1 ; ! %

.

2) $, '# + 2 & ' /GALI0

.I ! x = 0y = 0 px = 0.5 py = 0 & DBC %+ % % %+ x = 0, px 0 : (a)) y = 0 : (b))

,!% # # DBC ,%

e1 = (1, 0, 0, 0) e4 = (0, 0, 0, 1) % ,

& D=>I ( .. SALI :! ; GALI2 GALI3 GALI4 ! ! DBC ! % ! , m , I (a) m = 0 (b) m = 1 (c) m = 2 "%+ (a) SALI GALI2 % . & ! % %+! !

& D=> %+ .. SALI GALI2 GALI3 GALI4 ! DBC ! , ! .I

& D=>(a) % % ! #

$ & $ :

'# 24

% , & % SALI GALI2 #! # GALI3 GALI4 % # : :D=C;) & D=>(b) + m = 1 , : e1; & % % .% :=B>F;

SALI GALI2 # GALI3 1/t GALI4 1/t3 " D=>(c) % # m = 2 , : e1 e4; O / :=B>F;

% % , GALI3 1/t GALI4 1/t29 , , SALI : GALI2; ! !

% # ! . # ! ,

JBF ?C E> FA GBK 0 ! # %

% ! % % !

# ! SALI ! % % % . # , # 2#

, % JD GB GFK

" DBF(b) D=> % GALI3 GALI4 % !

! ! , # '

% % ! #

! ! ,

# ! . # ! # # %

, !! SALI GALI2I 0 ! :D=G; :D=C; GALI4 e41t ! GALI4 1/t4 ! ! % % % % ! % ! !%

% % ! ! ! ' % %

% : SALI

GALI2; tth + ! GALI4 , #, 1012 %+ % ! !

" ! % ! 2# = # %

Henon & Heiles :DB>; + D=B ! %

% A

# t = 500 GALI4 %! 1012 % % ! %+ ! # % ! 500 tth 400 % !+ ! ! %400 tth 300 %+ ! # # #

)' . % . , # , # % ! :

tth; 1 # ! % %+ !

! % :! tth; & . # ! !

28 $, '# + 2 & ' /GALI0

& D=BI $ , tth GALI4 ! !% 1012 Poincare %+ % x = 0 2#= # % Henon & Heiles :DB>; 0 % %

% %, % ! sticky &% !! !%