Efficient Orbit Integration by Kustaanheimo-stiefel ...

Æ

!

"

!

#$#

% & ! ! & '

( & )'(* + ,

+ & + & - ! + !. ! & ! ! ! & +. ( & + & ++ ! , & + + '(+ , .

#/#

!

! Æ +

& ! ! + & )0 . $111*

2 3 &

4567 + 4(8. )* + + ! ! '

( & )'(* + , ) $119&* )$* - ! + ! & '( + , ) $119+*

! )/*

!

& + ! ! ( !

:

& ) $11*. + ! &

&! ! &

! ! ) $11/ $119 !* &

& + & ! ! ; - + + !

! 3 - & ! & & ! , & + +

+

! !

+ 85?3 7@A 3 03 3. + &

& & ! ! '( + , ( !

:

& + + ! & . !+ B 1C. 0 & 8 7D. 0 & 8 D. 0 & 8 D7. 0 & 8 7 + 4777 .

#9# !

+ ! + + E & & ++ . 7 ! + + )* ! ! ! ! 83?3 ! )! * + )* 2E! , + + ! ! /1 )* !! + + @ ++: E ! )* !

& ,

! + & ! ! +

! ! & ,

! )* ! & ! ! ! &

& + ! ! ( !

:

& / & - ! + !

! & . 7 2+ & ! ! F( +: & + + ! ! '( + , + FG !: & - ! + ! ! '( + , + FH(: & + ! ! ( !

:

&

! F & H(: &

& + ! ! ( !

:

& . ? &

&

! ! '( + , + & F'(:. 7 + 2+ ! !

E ! + , I/1 ! +. & + + & + B 1C & + !. 3

! ! '( + , ! Æ + + & ! ) * . + !

! - ! & + ! +

! ! '( + ,

## + ,

+ & 1 .

! +

& + ! & ! + . 7 & - ! + ! + &

2 + ! !! & + 8 7D + D7

+ & ! ) *

! & Æ & & & + ,. 7 ! ! !J + !

! + & + + ! E 2 !

& & & & - ! & ! & & . ! ! Æ ! . ! !J + ! + !J + ! + . A & ! & ! & + '( & " E & ! &

++ + ! !

! & + &

+ & ! & + '( &

! & ! & - ! + ! 2 . E & !

K & + ! !. ( +

!

+ . 7 ! 2 ! & ! E ! ! + & 2 . % & ! + & & + . 6 + &

7 + ! +

!

.

#C#

. ?Æ & 3E 3 & + ,! ' > ? & 7 + /1 ( 8! & & B 1C

½¼ ¼

!

½¼ ½

!%

)

*

! " #$ & '(

& '(

& ('

+% !

!

! +%

½¼ ¾

! " #$

& '

! " #$

, " Æ " -% " . /&-

0 - " 1% 2 - 0 " & .

3

# 4/ " & "5 6 . & ! -% " - 7¼ 3

¬ ¬

¬

¼ ½ ¾ ¾ ¬ #

. &

" & 0" - -"% %

#L# & ! 2 + ! & + 2 !+ . 7 ! &

! !! + & + ! ! . + ! ( &

! (& )ML* ! !

!

! !! ! " ! & .

! &

& !2 ! ( !

:

E & &

! '( . ( 3-. )* . E ! ! & E + + ! '( + , !

! & ! ! .

+ + & + '( & & - ! + &

+

. ( +

! + . ! !! ! ! !! FH:. + !

+ & + + + & ( &

! (& )ML* & - & & + ! & & '( & + ' + .

! !

&

&! ,! + + ! ! 8 47 + 4777 ! ! ! + + !. K 2 ! + ! '( + , - & & & +. +

! ++! !! FH8: & + & + '( & + .

## 7

! ! & 2

! +

! . 0 & ' . % !

!J + ! &. 7 ! & + 2

! & ! ! '( + , + - ! + ! . 7 & + , E + ! + '( + , $

! E J + E + ! /.

& + !

&!

! ) + * ' +

! '

( & + , & ! ! )( &

! (& ML*. 7 E !2 !

& &

!

H )*

!J 2 .

5 - & & '( ) * E!

H B $

( 3-.)* & 8 7D.

B 9

)$*

#M#

!

B

)/*

+ & '( & &

! + . A + & - & & 3-. )/* &

!

B H

)9*

$ ) *$H

)*

)C* ) * % + & 3-. )* & 3-.)9* & & - & & . H9

( &

! (& )ML* ! E & + $ ) * H B

)L*

) H * ! +

& ! . ! ! )( 8 7D*. & *

B )* $

% + & & & 3-.)9*

! ! & & - & & .

E & + K E & B $ + ' + E! B )$*. *

# 1 # 7 ! B 1 +

! ! & 3-. )9* !

2 .

+

+ - ( &

! (& )ML* !. & ! + & ! Æ

& . ! ! + '( ! B

)M*

& - ++ & !! & ! & + . + & + ! . 7 8 7D + D7 ! + ! E &

+

! ! ) * E! B H & + ! & & !. = + & + & ! + E ! +

, + &

! ! E & 3-. )* ! & ! . !!! E + !

! & ! ) * & +

! B H

! + ! ! ! ! . ( !E 6 & )$11* & ! .

# #

! "# $ 5 E J & ! & + '( &

!

+ + - ! +. 6& + + & E + .

& 8 7 + D7 + !J & & ! + ! ! ! & ,. 7 & + ,! '( + , ( !

:

& !J &

+& !J . ! ! + + . ( !E 6 & 8 7D & ! .

% # & $ ' $ 7 ! ! & + & '(+ ,! +

! !. - & & !! H B 1 $

)1*

&

!

B

)*

& + !

! . 7 ) +* + ! ) * ) *

)$*

# $ # + & ! ! B $ ) * H

)/*

! +

! ! + !. ( 8 7D & E & &! . 7 !2 !

H

)9*

7 - & &

) * B H $

!

! ! & !

B $

)* )C*

! ! + & E !J " ) * ! + & 3- )1*

! )* & + )* ! + & 3- )1*

! )* + )* ! + & 3- )1*

! )C* + )$ * ? ) * + 3-. )$* ! + )$* ? )* +

! )$*

? )* +. 3=7 > 85?3 7@A 3 $ 03 3. + $ & ! ) * )*

! )$ *. % & ! ,

+ & ! + ! & . 7 2+ ! ! ! &

! & & ) *

! )$ *. ! !

# / # !2 !

)N* H $)N * H $) *

B

)N* H $)N *

)L*

% ! & & 2 ! & ! ! & & & + . ? ) * ++! FA !:

! )* ! FA !H8:. ! & & ! 2 . ! & + !. 0 ! & ) + * ! & ! & 2 + . & ! & ! & ! + & .

! + & & + .. ? )* !

+ - ! + . ! ! &" )* + 2 2+ )$* + & & ! & & - & &

! )/* + + - ! 2 & + .

! ! & 3-. )C* +

& & + ! & !! + 2+.

E !J

& + '(+ ,! ! +

! . 7 ! &

! ! & & - & + & ! .

# 9 # 3=7 > 85?3 7@A 3 / 03 3. 6& ! + !!

& +"

& & + . + / E ! & ! ! 83?3 ! )! *

& & & ! & + ! '( + , . + !J + ! !

! ! . !J & ! ! !J '(+ ,! - ! + !

! ! ! ! + C & 8 D.

% #

! . ! ! & . ! ! . + ! + ! + & ! '(+ , " ) * + '( + , ! & )* '( + , - & & )* '( + , - &

! ! & )$ * '( + , + + ! & )$* '( + , + +

! + !

! ! & )/ * '( + , - ! + ! &

! )/* '( + , - ! +

! + !

! ! & . + & 2 ) *

! )* + - !

# # . & - ! + & + & & & ! E + & ! . 5 !J 2 & + + . !&!+ !J + ! + + .. )* )$ *

! )$*

! )/ *

! )/*. !J E ! !J + ! & + " & 2 +

! / & - ! ! . +! & I/1 ! & ! ! IC1 ! & ! . , 2

! 1. !

. !J & ! & + ! !J & / !+ + !. E ! ! + + ! .. ) * )$ *

! )/ * ! ! & )* )* )$*

! )/*. !

! & + ! & ! . 3=7 > 85?3 7@A 3 9 03 3. 3=7 > 85?3 7@A 3 03 3. 3=7 > 85?3 7@A 3 C 03 3.

# C #

! !J

& + + . +

! + E + ! & . 0 . ( + 9 + C &

! E !2! . ? ! & + / & 8 7D & ) * + C & 8 7D & )$ *

! + / & 8 D & )/ *

+ , ! ! IM1 ! +.

& + &

+ ! !J ! !. !J ! ! & ! & +

! !J & +. 7 ! .. + 9 K E + ! & & ! & + ! &

+ ! ! 2 . + . 0

! + ! & ! ) * 2 + - & & 2 &

! ! 2

! . 7 & + + .. + K E

& + ! & . J & + + ' + & '( . 0 &

! ! + - ! . % ! & !

! & + + ! & - ! + ! + 2

+ & ! + &

+ ! .

# L # & + &

! ) * &

! . ! ! & + ! )* + 2

& ' +

! )$* & + !

! & &

! )/* + 2

+ & &

. & - ! + .. + K

E . + !

! & + .

! + ! & & & & & &

! ! 2 + . 7

& ! + ! !. % 2 ! ! + 2

+ + ! & & , . & ! - ! + ! + + &

& ! .

(

% ! . 3=7 > 85?3 7@A 3 L 03 3. & 2 ! ! & ! Æ . ( + L + & ! ! : !! .

# #

E & . % ! ! + + 2

!J

& ! & - ! + + ! ! : ! . 3=7 > 85?3 7@A 3 03 3. E E & & ! + . + + & ! ! & ! +

! ! & . + ! + & & & '(+ ,! + !

!! & & - & . & ! '( + , + & 2 + !

& &

! !. & + + ! & !. 7

& - ! + + &

+!. 3=7 > 85?3 7@A 3 M 03 3. 5 ! & + . + M + & ! ! & 3 + & & 3 : 2 . ! , &

! + +

+ &

! ! . ! +

!J + ! E ! . 7 & + & ! + &

# M # &

! !. & + ! + - ! . 3=7 > 85?3 7@A 3 1 03 3. E J & . + 1 + ! " ! 7 & + B 1$L ! . (

! !! & ! + ! + . 4 &- & ! ; & ! & ! + & !!

! ! & . 7 + 2+ ! & 7

! $111.1. & 2E! + $111.1. +! , + $1 ! . 7 2 ! ! &

2+. 3=7 > 85?3 7@A 3 03 3. !! + + + J . + & ! 7 ! ( :

. 2! & ! + & + ! . + 2 ! + 2

!J + ! & . 3=7 > 85?3 7@A 3 $ 03 3.

# $1 # E !J . + $ + L ! & , !! ! . 7 ! & + 9 & - ! + !

!

& + ! ! ( !

:

&

! & & . /I$ & & + '( + ,

!

+ & + ! & & - & . + & + ! !.

& + + - & & ! + ! & . &

! + + ! ! & . 8 & !.

) $

6 + + ! & & &

! '( !

! + '(+ ,

! + + - ! +. !! & ++ . ! + ! & '( + , " & !

!

+ & + & + - ! & ! '( + , . 7 & ! ! ! &

! + &

!

! +. 7 ! ! & + &

! + + 2

&

# $ # &

. & ! , &

+ ! '(+ , + ! .

# $$ #

*! . $11 $ LC . $11/ $C 1ML )8 7* . $11/ $C $CL )8 77* . $11/ $C // )8 777* . $119 $L /C/ )8 74* . $119 $ M$1 )8 4* . $119 $ 99C )8 47* . $119! ! ?. .

! = . . )8 477* . $119 $ 9 )8 4777* . $119& $ .C )8 7D* . $119+ $ .C )8 D* . $11 $M .$ )8 D7* 0 0. . $111 8( $ M '

8.

! ( & 3.5. MC . +. . $ $19 ( & 3.5.

! (& @. ML 5

! + ?

) O> ( +*

! ( 5)3D .$.

# $/ #

log10 |∆L0|

+ ,- #

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

30 Steps/Period Kepler, e=0.6

KS 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S

10−38 Single

10−44 Quad 10−52 KS+T 10−60 Single+T

10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 log10 (Physical Time/Period)

+. .P 3 + & + ) B 1C* '

+ ! + , !. 7 !

+ ! & & ++ . F'(:

'

( & + , ! F( +:

& + + ! + FG !:

& - ! + ! + F'(:

'

( & + , ! FH(:

( !

:

& ! FH8:

+ & &

! FH:

+ & .

# $9 #

log10 (Error)

+ ,- ' $ 2 0 Harmonic Oscillator ∆t Unscaled -2 60 Steps/Period -4 10−6∆t Unscaled+PT -6 -8 -10 -12 -14 -16 -18 -20 10−12∆t Scaled -22 -24 10−20∆ρ Unscaled -26 -28 -30 -32 10−24∆ρ Scaled -34 -36 -38 -40 -42 -2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 log10 (Fictitious Time/Period)

+. $.P 3 + & !. ( + !

! ! &

! & ! N )N* H $)N *

! & ! N. % ! !

! +. & !!! & ! + ! & + ! & & & !+ !

!! + FH8:. , 2 2E! IC1 ! ! &

! ! /1 . ! & ! ! + + ! " & !

! / & ! .

# $ #

. - # $ - $/ ' $ 22 Highest Stable Order

20

Harmonic Oscillator Implicit Adams, PECE

18 16 14 12 10

Scaled

8 6

Unscaled 8

16 32 64 128 256 512 1024 Steps/Period in Fictitious Time

+. /.P ( & . 7 ! + ! & 83?3 ! & ! ! + +

& & 2E! ,. % + !! + E + & & 2 1$9 !.

# $C #

log10 |∆E|

# / ! $ -2 -4 AM10, 30 Steps/Period o |∆x|/r -6 Kepler, e=0.6, I=31 KS+T -8 10−4∆L0 -10 -12 10−8e∆ϖ -14 -16 10−12∆e -18 -20 -22 -24 -26 10−16(sinI)∆Ω -28 -30 10−20∆I -32 -34 10−30(∆a)/a -36 -38 -40 -42 -44 -46 -48 -2 -1 0 1 2 3 4 5 6 7 8 9 10 log10 (Physical Time/Period)

+. 9.P 3 & '

+ ,! '(

& +

! +. ( + !

! !2! .

# $L #

log10 |∆E|

# / $ -2 -4 AM10, 30 Steps/Period o |∆x|/r -6 Kepler, e=0.6, I=31 -8 Single+T 10−4∆L0 -10 -12 -14 -16 10−8e∆ϖ -18 -20 10−12∆e -22 -24 10−16(sinI)∆Ω -26 -28 -30 10−20∆I -32 -34 -36 -38 -40 10−30(∆a)/a -42 -44 -46 -48 -2 -1 0 1 2 3 4 5 6 7 8 9 10 log10 (Physical Time/Period)

+. .P ( + 9 + + ! + .

# $ #

log10 |∆E|

# / 0 # $ -2 -4 AM13, 30 Steps/Period o -6 Kepler, e=0.6, I=31 -8 Quad+T |∆x|/r -10 -12 10−4∆L0 -14 -16 -18 10−8e∆ϖ -20 -22 10−12∆e -24 -26 10−16(sinI)∆Ω -28 -30 10−20∆I -32 -34 -36 -38 -40 -42 -44 10−30(∆a)/a -46 -48 -2 -1 0 1 2 3 4 5 6 7 8 9 10 log10 (Physical Time/Period)

+. C.P ( + 9 - ! + ! + .

# $M #

log10 |∆L0|

1 / 23 - . (

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

Moon’s Perturbation KS 30 Steps/Period HALCA 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S

10−38 Single 10−44 Quad 10−52 KS+T 10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Days

+. L.P ( + & ! ! : . ! & '( + , & . & + 2E! '

.

# /1 #

log10 |∆L0|

1 / & (

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

Air Drag 30 Steps/Period HALCA

KS 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S

10−38 Single 10−44 Quad 10−52 KS+T

10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Days

+. .P ( +. L ! ! ! ! +.

# / #

log10 |∆L0|

1 / $$ (

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

J2 Perturbation 60 Steps/Period HALCA

KS 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S 10−38 Single 10−44 Quad 10−52 KS+T 10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Days

+. M.P ( +. L ! & 3 . & ! B / . , ! ! . ! & ! ! !

+! & + ! ! ( !

:

& / & - ! + !

!

& + ! ! ( !

:

&

! & & . Æ

# /$ #

log10 |∆L0|

1 / - . (

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

Jupiter’s Perturbation KS 30 Steps/Period Icarus 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S

10−38 Single 10−44 Quad 10−52 KS+T 10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Years

+. 1.P ( +. L + ! ! 7 ! : !! . ! & , +.

# // #

log10 |∆L0|

1 / + 4$ (

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

General Relativity 30 Steps/Period Icarus

KS 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S 10−38 Single 10−44 Quad 10−52 KS+T 10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Years

+. .P ( +. 1 ! ( : + .

# /9 #

log10 |∆L0|

1 / 1

2 0 -2 -4 -6 -8 -10 -12 -14 -16 -18 -20 -22 -24 -26 -28 -30 -32 -34 -36 -38 -40 -42 -44 -46 -48 -50 -52 -54 -56 -58 -60 -62 -64 -66 -68 -70 -72 -74 -76 -78 -80 -82 -84 -86

Moon’s Perturbation 60 Steps/Period HALCA

KS 10−6 KS+PT

10−20 True+S

10−30 Antifocal+S

10−38 Single 10−44 Quad

10−52 KS+T 10−60 Single+T 10−66 Quad+T

-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 log10 Days

+. $.P ( +. L , ! & !J

!

. ! & ! ! 83?3 ! 9 & - ! + !

!

& + ! ! ( !

:

&

! & & .