. . 0
ÈÊõ ÇFùñÇ
Ï Ï
göö24/÷ 27ú û
´üQø ËVø Ë
{&þ:ÿsþ¡ £¢¥¤§¦©¨ ¨
ÿ ´7YyT7T7aÝæ_DÝ`Þ[\^_I·lô·ÄÞ{aX&öGßU`Ü\±ZÙY[´uaÞ{aD´<´7YËÝOU`X_Kðz_D\ ïW_Í]WYË_bÜ`´7ayïW\^X}_D]_DX2YyÝ ZÙY˶ ê¯P7ß:ÿ>]`Y
O
°^_b´7ÚUW´<_ ²xÞ{YyÝ ZÙ´<_D]_KYyX¿xê µ ß:ÿàê °^_?ßYËÞ[_b´<_DÞ{ZÙY[´\$y_D]`_ÖÜãaD´uU`X ´7Y[Ú\^X2Yé]WY aD´7V\±Z7_DT÷Ú\±´<\~ÝWaß
_ba ´7Yy]WaD´u]WaToÜãaÝ ZÙaTS¡Û_bÚD´<_DÝ`ÚDYy_DÝWaT¡ YV¡ ]`YóöGßU`Ü\±ZÙY[´y²¥Z7_D°xíUWY»a _DÝWÚU`°^aæ ´7YyT7T7aÝ`_DÝ ZÙY
î ï
M
?:CADCæ EkH=E+Fá°^\±V´<_2YyXZÙaD´<ÝWau]WY¢kb¶ ²C´7YyT7ÜrYyÞ{Z7\^äj_DX2YyÝOZÙYDêxÿ:T7TÙa Þ[\^_D]`_÷_}YyT7Z7_o´ ® Y{ï \^T7ZÙYUX}_ÜraDÜCU`°^_råÞ _ba2]WYX}_D\^TG]WY¢kb¶D¶}_DTÙZÙY[´Cßa\^]WYyTSíUWY´7YyÞ{Y[VãYa)Ý`aX2Ye]WYþÛ´ íOUWYxÞ{Y[´<Þ[_)]WYlkb¶ì ]WYy°±YyTaD´7V\±Z7_eYyX ZÙaD´<Ý`a]WYl¡ ²YyÝ`íOU_DÝ ZÙaíOUWYEP¶ì aáðz_LxYyX ZÙaD´<ÝWa]WYE¡ ê î ï ÿ6YyTÙZ7_bV\^°~\^]`_D]WY_ó°±aÝWÚDaóÜ`´<_L[aé]WaT)þ´7a\^_DÝWaTàZÙYyX T7\^]Wa»YyTÙZ7U`]`_D]`_ìÜãaD´d¡ÛY[ä\^T7aÝ9² úWaY{î ® vzv]w X}_L DY[´Y úWaYyX}_L DY[´oøg· ü²MíOU`Y)X}aTÙZÙ´<_b´<_DX4íUWY)YyTÙZÙYyTá_DT7ZÙY[´ßa\^]WYyTT _baõYyT7Z`ß_bäDYy\^Te_baì°±aÝ`ÚDa ]`_Ä\^]`_D]`Y»]Wa \^TÙZÙYyX}_ a°^_b´yê ³T÷þ´7a\^_DÝWaToTÙY[´<\^_DXÒaDV çÙY[ZÙaToÜ´<\^X2aD´<]\^_D\^T)Ý`_Ä´ ` ® ·lô·b²rÜ´7aläb_yäDYy°^X}YyÝ ZÙY}Þ[_bÜ`Z7U`´<_D]WaT_D\^Ý`]`_o]`UW´<_DÝOZÙYà_DT:Y[Z7_bÜ_DT Ý`_D\^T¼]_ðñaD´ ® vzv ø»Uì_b´ty_b´\Y ÞúWa°^°¥· yüê \^XU`°^_råÞ aDYyT:ÝOUXKßY[´<\^Þ[_DT:Z7_DXáVGßYyX X2aTÙZÙ´<_DXíUWYàYyT7ZÙYyT_DTÙZÙY[´ßa\^]WYyT ® vzv]w ZÙY[´<\^_DX4TÙaDðñ´<\^]WaõU`X_uY[äDa°^UMåÞ _baìÞ{a°^\~T7\±aÝ`_D°M\^X2ÜãaD´7Z7_DÝOZÙYìø»U»_b´ty_b´<\Y[Zá_D°mêu· ü²íUWY}]WY[äDY[´<\~_ vzv TÙYeX}_DÝ`\±ðñYyTÙZ7_b´¼Ý_)Y{ï \^T7Z{æYyÝÞ[\^_}]WYðz_DX ß °^\~_DTeø»Uì\^°~_DÝ`\· ýüê Ø YyT7ZÙY)Þ[_bÜß Z7U`°±au_bÜ`´7YyTÙYyÝOZ7_DX2aTaT´ Ø þÛ´7a\~_DÝWaT[ê aFÜ´<\^X2Yy\±´ ]`\^Ýæ_DX}\~Þ[_Ë_Ä°±aÝWÚDaKÜ`´<_L[aK]_Ë´7YyT7TÙaÝQæ_DÝ`Þ[\^_ù·lô·bê©ÿTÙY[ÚU\±´y²GYyTÙZÙY»X}a]WYy°^a Y»_bÜ°~\^Þ[_D]WaÄÜ_b´_ß ]WY[ZÙY[´ Ø íOUWYuÜãY[´ ® TÙY[ÚU`Ý`]WaõZÙ´_bV_D°^úWaW²M_ ]`\^ÝQæ_DX}\^Þ[_õ]`_õ´ ·sÓDÔ Ó Ó Ô Y X Y X ¹ [Zt]× ºKZ_¸ ¸¡Ú Öh ÖÚ(Z"!$# ¬¹»º Õ ¼ ¼ ÜC°^_DÝWY[Z`ß_b´<\^_ ²x_DÝ`_D°^\^T7_DÝ]Wa Þ{aX2aÖYyTÙZÙYéðñYyÝæaX2YyÝWaÆÜãa]`Y»_bðñY[Z7_b´ _KYyT7Z7_bV\^°^\^]_D]WYì]WaTuþ´7a\^_DÝWaT ® ´ {èþ¬ýÅþZ\æ¨ ) ¨'e)+æ. ¢ , ¤§. ¨ 9 :9 Ø aT7TÙaËX2a ]WYy°±aóTÙYyX}\ î_DÝ_D°mß Z7\^Þ{aóÜ_b´<_ó]WYyT7Þ{´7Y[äDY[´}_ó]`\^ÝQæ_DX}\^Þ[_ì_é°±aÝ`ÚDa»Ü`´_L[aé]`aT_DT7ZÙY[´ßa\^]WYyT ® ® þ´7a\^_DÝ`aTxVC_DTÙYy\^_jîTÙY)YyXZÙ´yæYyTe_bÜ`´ ® ® X\±ZÙYyX4ZÙ´<_DÝ`T7ðaD´ CADC ® v _DT ® ® _ìÞ{aÝ`TÙZÙ´ ® T7YyÝ Z7_b´+aX2a ä\~X2YyÝ ZÙae]_DT2ßaD´7V\±Z7_DTMÚ\±´\^ÝWaxÝ_x´ ÜãY[´7Z7U`´7V_bZ7\±äDa)Þ{aX _)ZÙY[aD´<\~_2]WaTS\^ÝOäj_b´<\~_DÝ ZÙYyTS_D]`\~_bV¼ß_bZ7\^Þ{aT[²`UWZ7\~°^\$y_DÝ`]Waàa2ðñ_bZÙao]WYíOUWYaT¼]`\±ðñY{î ® ´ Ü´ßaDÜ`´<\^_DTyêMÿÇTÙY[ÚU`\^´y²`_DÝ`_D°^\^T<_DX2aTSÞ{aX X}_D\~T/]WY[Z7_D°^úWYeÞ[_D]_}U`X}_2]WYy°~_DT[ê "!$#&% &')(*#,+¢-/.01 2#,+ ³>Ü´<\^Ý`Þ[\±ÜC_D°¼Ü`´ ÜãY[´7Z7U`´7V_bZ7\±äDaTußY)íUWYáaÄæ_DÝWÚU`°±a´7YyT7TÙaÝ`_DÝOZÙYõøzÝWa÷Þ[_DT7aß? üGÜãaT7T7U`\ÛU`X}_}ðñ´7Yyí¦¥UãæYyÝ`Þ[\^_oÜ`´ßaDÜ`´\^_ CADC ® XU`\±ZÙa»Ü`´ßabïW\^X}_ ]`Yp[Y[´7aWê TÙZÙa»ðz_Lu_bÜ_b´7YyÞ{Y[´)aTÜrYyíUWYyÝWaTà]`\±ä \^TÙaD´7YyTàÝ`_ìT7a°^UMåÞ _baóÜraD´TOßY[´<\±Y ® ]`aËÜ`´7aDV°±YyX_ ²+aÄíOU`YõZÙaD´<Ý`_éaÄXÆßY[ZÙa ]Wa ÜãY[´7Z7UW´ ® Ü´7aDV°±YyX}_ ²eßY÷Ý`YyÞ{YyT7T`ß_b´<\±a»YyÝ Z _baóZÙ´<_DÝ`T7ðaD´ ðñ´7Yyí¦¥UãæYyÝ`Þ[\^_àTÙYçÙ_]`\±ðñY[´7YyÝOZÙY¼]WYE[Y[´7aWêÿ X2a ]Waà]`Y/Y{ïYyX}Ü°±aW²äj_DX2aTGÞ{aÝ`T<\^]WY[´<_b´aáÜ`´7aDVC°±YyX}_e]WY ZÙ´yæYyTÞ{aD´<ÜraT:´7YyT7ZÙ´<\±ZÙa÷Yy°mß Ü`Z7\^Þ{aoÜ°^_DÝ`_b´XÆßYy]`\±aWêá³ Ü`´7aDVC°±YyX}_}ZÙYyX4]`a\^T¼ÚD´_DU`T]WY°~\±VãY[´<]`_D]WYYD² Ý`a2ZÙ´<_bZ7_DX2YyÝOZÙaoÜãY[´7Z7UW´ C @ C @ C @ C C C @ C @ ìÆø43 f*3 f65 f65 üæ_ì87Dø43 f*3 f65 ü3¯:9»ì ø43 f*3 f65 f65 ü øGPWê^·lü T7YyÝ`]Wa;9:< ·óU`XÒÜrYyíUWYyÝWaKÜ_b´jæ_DX2Y[ZÙ´ O Z7U`´7V_D]WaD´y²QY=3)>8f65?>Säj_b´\yß_bäDYy\^TÞ{aÝjçÙU`Ú_D]`_DT[ê X aUWZÙ´_DTàÜC_D°^_yä´<_DT[² ì87}Þ{aD´<´7YyTÙÜãaÝ`]WYna_»Ü_b´7ZÙY C Þ[\^´<Þ[U`°^_b´2]WaéÜ`´7aDVC°±YyX}_»Yðì _éÜC_b´7ZÙYõYy°mßa Ü`Z7\~Þ[_ êKÿ û:_DX}\^°^ZÙaÝ`\^_DÝ`_äì87÷ÜãaT7T7U`\¥U`X ÚD´<_DUÖ]WY ` O C °~\±VãY[´<]`_D]WYD²ÜãaD´7Z7_DÝ ZÙaÖßYá\^ÝOZÙY[ÚD´ß_bäDYy°ê °^_2]WY Ý`Ye_2ZÙaDÜãa°±aDÚ\^_2V¼ß_DT<\^Þ[_2]`_2´7YyT 7 7 7 ® YU`X4æß _DÝWÚU`°±aíOUWY:°^\±V`´_ êÿ:T7T7\^Xõ²]`_D]`_DT¥Þ{aÝ`]`\gåÞ aDYyT¥\^Ý`\~Þ[\^_D\^TA3 f*3 f65 ²OÜ´7aÞ{Yy]`Y{îTÙYx__DÞ<ú`_b´S_DT C @ C 7 ® C C @ C 7 äb_b´<\yß_bäDYy\~T/_råÞ _baY÷æ_DÝ`ÚU`°±a 3 f 5 f 3 æB3 ]WYïì ²]`YeðaD´X}_2íOU`Yá_bÚDaD´<_ 5 YàU`X>æß _DÝWÚU`°^a}íOUWY @ ýDµ Ó Ó Y tz× ]×àSà!FZ ºK[ZKº!mDC ¸ [ ¹w¬Ú(Z ½ Õ ® O Þ[\±´<Þ[U`°^_2Þ{aX ðñ´7Yyí¦¥UãæYyÝ`Þ[\^_Þ{aÝTÙZ7_DÝ ZÙY2]`\^ðY[´7YyÝOZÙYá]WYÑ[Y[´7aWê TÙZ7_DT:äj_b´<\lß_bäDYy\^T¼T _bau°±aDÚDau\^ÝOZÙ´7a]Uy\ î C C ]`_DT:YyX ì ²9ðñ_L[YyÝ`]`abîTÙY2_÷T7Y[ÚU`\±´U`X}_uXÆßYy]`\^_uTÙaDV´7Y 5 ²-]`YàðñaD´ ` vzv C @ @ °^\±VãY[´<]`_D]WY¥Y Þ[_b´Þ{aX¬U`X}_eû:_DX}\^°±ZÙaÝ\^_DÝ`_ ìÍø 3 f 3 f 5 ü+íOUWYußY\^ÝOZÙY[ÚD´ß_bäDYy°øzûxYyÝW´<_b´]ì· ¶ üê ® ÿ ]WYyTÙäb_DÝOZ7_bÚDYyX ]WYyTÙZÙYXKßY[ZÙa]Wa2YyT7Z`ß_àÝWaðz_bZÙa)]`YxZÙY[´SíUWY:]`Y[ZÙY[´ ® _DÝWÚU`°±aá]`Yæ mì87/Ü_b´<_áÞ[_D]_Þ{aÝ]`\gåÞ _baà\~Ý`\^Þ[\^_D°² aíUWYxÞ{aÝTÙaX2YXàU`\^ZÙaáZÙYyX2Üãaà]`Y:Þ9ß_D°^Þ[U°±aWêÿx°qßYyX ® CGF C C ]`\^T7TÙaW²WÝ`_}Ü_DT7T7_bÚDYyXE5 5 ÜrY[´]WY{îTÙYeíU`_D°^íUWY[´/\^ÝWðñaD´ ® ® O C C ZÙ´7a]U¥Þ{aX}Ü°^\^Þ[_råÞ aDYyTÝ`_úWaD´_:]`Y¥\^ÝOäDY[´7ZÙY[´+_TÙa°~UMåÞ _baÜ_b´<_:aDV`ZÙY[´A3 ø,I7üDf65 ø,I7üê X ÝWaT7TÙaZÙ´<_jî ® V_D°^úWaW²b_bÜ`´7YyTÙYyÝOZ7_DX2aTMU`X}_x_bÜ`´7alï \^X_råÞ _ba:]`\^ðY[´7YyÝOZÙY_ba:Ü`´7aDVC°±YyX}_ ²jaÝ`]WYYyT7ZÙYyTQ\~Ý`Þ{aÝ äDYyÝ\±YyÝ ZÙYyT ® ® ® T _baõ´7YyTÙa°±ä \^]WaTT7U`VTÙZ7\±Z7U`\~Ý`]Waa Þ9ß_D°~Þ[U`°±aõ]WY2äb_b´<\yß_bäDYy\^T_råÞ _ba Y»æ_DÝ`ÚU`°±aõÜãYy°^_oZÙ´<_DÝ`TÙðñaD´ äj_b´<\lß_bäDYy\^TS°±a Þ[_D\^T[²`Þ{aX2aoY{ïÜ°~\^Þ[_DX2aTS_}TÙY[ÚU\±´yê ® C 7 ÿ]`\^ÝQæ_DX}\^Þ[_ó]WYßì Y ß Þ[_b´<_bZÙY[´<\y_D]`_éÜãYy°^_»°~\±V`´<_råÞ _baé]WYJ5 YyX ZÙaD´<ÝWaË]WYJ5?K{êU»_D\^T)Ü`´ ® C C C C Þ[\^T7_DX2YyÝOZÙYD²C\^ÝOZÙ´7a ]`Uy\^Ý]Wa2_DT¼äb_b´<\yß_bäDYy\^TML ¯ON4P æRQ µ"3 Y{ïÜÛø,N/5 ü²ã_}°^\±V`´_råÞ _ba_DÞ{aÝ ZÙYyÞ{Y2YyX ® ZÙaD´<ÝWaì]WaìÜraÝOZÙaSL ¯TNUP êuÿxT K K Þ{aÝ`T7\^TÙZÙYYyXðñ_L[Y[´:U`X}_)ZÙ´<_DÝ`T7ðY[´yæYyÝ`Þ[\~_ ò C C C C L ¯;N P æ?øVL HWL ü3¯ONøVP HWP ü øGPWê µü K K ® C ]WY2ðaD´X}_ íUWY}aìÝWa äDaÖæ_DÝ`ÚU`°±a 5 Ý _ba»X}_D\^Te°^\^V`´<_÷YyX ZÙaD´<Ý`a»]WYX5 ²X}_DTàÞ[\±´<Þ[U`°~_uYyX ZÙaD´<ÝWa K ® ]WYuø L f P üæøz¶gf7¶ üê T7ZÙaFßY)ðñYy\±ZÙaõX2Yy]`\^_DÝOZÙY2U`X}_oZÙ´<_DÝ`TÙðñaD´ K K C @ C @ C @ ìÆø 3 f 3 f 5 f 5 ü²WaÝ]WYaTáæ_DÝWÚU`°±aT 5 f 5 Þ[\±´<Þ[U°^_DX YZ[æYyX ð´7Yyí¥U9æYyÝ`Þ[\^_DT ÝWaläb_2ûx_DX\^°±ZÙaÝ`\^_DÝ_ Ü`´ßaDÜ`´<\~_DT]\±ðY[´ ® ÜrY[´ ® T7\^X2Ü°^YyT/YZÙ´<_DÝ`TÙÜ_b´ ® ® ® íOUWYeT7Yá]WY[ZÙY[´ K K 7 aTSÜãaÝ ZÙaTx]WYYyíOU\^°mß V`´\±aà]`_2û:_DX}\^°^ZÙaÝ`\^_DÝ`_cì ê Y&Z&¥[\(*+1(*]^'_.)-`'_#ab(*]^cd.J "+$+e#,af%()c4!e#/0?. ® ßxX}_äDY_\~Ý ZÙ´7a ]`Uy\~]`_DTS_DT¼äb_b´<\yß_bäDYy\~Tx°±a Þ[_D\^T[²CaÜ_DT ® X2a]`Yy°±a Y ß _DÞú`_b´U`X}_óY{ï Ü`´7YyT7T _baÄ_DÝ`_D°mß Z7\^Þ[_»Ü_b´<_é_Ëû:_DX}\^°±ZÙaÝ\^_DÝ`_äì ²¥íOUWY÷ÜãaT7T7_éTÙY[´}X}_jî Ø Ý`\±ÜU`°~_D]`_o_D°^ÚDY[V`´<\^Þ[_DX2YyÝOZÙYDê _ìX}_D\±aD´\^_õ]WaTeYyT7Z7U`]WaT_DÝ`_D°mß Z7\^Þ{aTáTÙaDV`´7Y_ì´7YyT7TÙaÝQæ_DÝ`Þ[\^_Ä·lô·b² ® ë ® UWZ7\^°^\$y_DX2îTÙY)Y{ïÜ_DÝT aDYyTáYyX>Þ{aOaD´]WYyÝ`_D]`_DT _b´7ZÙYyT7\^_DÝ_DTáíUWYÝ _ba ´ ® _ ]WY[ZÙY[´ ® ® äj_b´<\lß_bäDYy\^T÷æ_DÝWÚU°±a»X2aX}YyÝ ZÙaW²+]Wa»Z7\±Üãa³¡Û_bÜC°^_DÞ[\^_DÝWaW²Ý _baóT _baóÞ{aÝOäDY[´7ÚDYyÝOZÙYyT}ÝWaóÞ[_DTÙaó]_ì´7YyTÙî ýDý ·sÓDÔ Ó Ó Ô Y X Y X ¹ [Zt]× ºKZ_¸ ¸¡Ú Öh ÖÚ(Z"!$# ¬¹»º Õ ¼ ¼ T7aÝæ_DÝ`Þ[\^_õ·lô·bê ® ¯Q_b´_Y[ä\±Z7_b´YyTÙZÙYyT+\^Ý`Þ{aÝOäDYyÝ`\±YyÝOZÙYyT[²YyX ÝWaT7TÙaX2a ]WYy°±a]WYyTÙYyÝOäDa°±äDYyX2aTU`X}_:Y{ï Ü_DÝ`T _baYyX TßY[´\±Yo]WY÷þ_+£°^aD´Ùî:cWaUW´<\^Y[´]`_ìûx_DX}\~°±ZÙaÝ`\^_DÝ`_ ²QV_DTÙYy_D]_ìÝWaó]WYyTÙYyÝOäDa°±ä\~X2YyÝ ZÙaó_DT7T7\^XÆßY[ZÙ´\^Þ{a»]`_ ® v v ðzU`ÝMåÞ _baoÜãY[´7Z7UW´7V_D]`aD´<_ øucWY[´7´<_L{îãUìYy°^°±a÷Y _bZÙa»· y üêxÿ?\~]ßYy\^_]WYyT7ZÙYà]WYyTÙYyÝOäDa°±ä \^X2YyÝOZÙaùßYíOUWY ® aóX2alä \^X2YyÝOZÙaé]WYo°^\^V`´<_råÞ _baìYyX>ZÙaD´<ÝWaé]Wa»ÜraÝOZÙagL fhP Üãa]WYoTÙY[´à´ K K ® O T7YoU`X}_ Y{ïÜ_DÝT _baó]_ ûx_DX\^°±ZÙaÝ`\^_DÝ_õYyX>TßY[´\±Yo]WY÷þ_+£°^aD´àYyX ZÙaD´<ÝWa»]`YiL fhP ê X>ÝWaT7TÙa K K Þ[_DT7aW²/_ û:_DX}\^°±ZÙaÝ`\~_DÝ`_ Y»Y{ï Ü_DÝ`]\^]`_éÜ`´<\~X2Yy\±´7aéYyXß TOßY[´<\±Y»]`Yìþ_+£°^aD´÷ÝWaToX2aX2YyÝOZÙaT 3¥>YD² ` ® ® ® ÜãaTÙZÙY[´\±aD´ ® 5j>²T _baY{ïÜ_DÝ]`\^]WaTGYyXTßY[´\±Yá]WY¢cWaU`´<\±Y[´yê Dkld]^'1 "!e#2m "] 0n#V o qpY#V ¥rs% &c4# 01 ® ÿÍZÙ´_DÝ`TÙðñaD´ ÜãY[´7Z7U`´7V_bZ7\±äb_DT[êR:Yeðñ_bZÙaW²r_2ûx_DX\^°±ZÙaÝ`\^_DÝ_uøGPWê^·lü¼_D]WaDZ7_}_}ðaD´ C @ C @ C @ C C @ C @ ìÆø 3 f 3 f 5 f 5 ü&æ ìg7Dø 3 f 3 ü/¯ut ì ø 3 f 3 f 5 f 5 ü øGPWê ýü ® aÝ]WYvtS<$·b² TÙYyÝ`]WaàÜãaT7Tyß äDYy°`_DÞú`_b´GU`X}_á_bÜ´7ayïW\^X}_råÞ _ba)\^ÝOZÙY[ÚD´ß_bäDYy°C]`Y ì4_bZÙ´<_yä+ßYyTS]WaàXÆßY[ZÙa ]Wa ® Ø ìg7}Þ{aÝ`TÙZ7\±Z7U\a8fK DY[´<Ý`Yy°hì]WaóXKßY[ZÙa]`aWê a»YyÝ Z7_DÝOZÙaW² Þ[°lß_DT7T7\^Þ{aì]`YoûxaD´<\²ÛaÝ]WYo_»TÙa°^UMåÞ _baó]WY Ý`a)Þ[_DT7a}]`_´7YyT7TÙaÝæ_DÝÞ[\^_u·lô·»ßYÜraT ÜãaT7T ÿ:]`\^_bVxß_bZ7\^Þ{aTSÜ_b´<_)´ C @ C @ YçÙ_DX Ü fhÜ _DTðñ´7Yyí¦¥UãæYyÝ`Þ[\^_DT)Þ[_b´<_DÞ{ZÙY[´lß T7Z7\^Þ[_DT2]WY 5 f 5 ²´ ` @1x C @ @ íUWYvw æöÜ Ü <$·b²ßY/ÜãaT7Tyß äDYy° ï _b´ 3 f 5 YS´7YyTÙa°±äDY[´_ûx_DX}\~°±ZÙaÝ`\^_DÝ`_ ìÞ{aX2aáTÙYSðñaT7TÙY/U`X C C T<\^TÙZÙYyX}_]WYxT`ßaáU`X?ÚD´_DU]WY/°^\^VrY[´]`_D]WYD² 3 f 5 ê+ÿ¼Ü°~\^Þ[_DÝ`]WaaáXÆßY[ZÙa ]Waà]`Y¼û¼aD´\`_YyTÙZÙY:T7\^TÙZÙYyX}_ ² y y y f ²OZ7_D\^TíUWY/_eÝ`aläb_áû:_DX}\^°±ZÙaÝ`\~_DÝ`_ìßY ü²OY _DÞú`_DX)îTÙY:ÝWaläb_DT+äb_b´<\yß_bäDYy\~T 3 y f 5 3^y f 5 ìgybø 3^y f 3^y f 5 C @ @ C @ C @ Üãa ]WYT7Y[´Q´ C Y¼UX}_Þ{aÝ`TÙZ7_DÝOZÙY¼_bZOßß Y/aD´<]WYyXzwy²DYS´7YyÞ{Y[VãYSaeÝWaX2YS]WYS\^ÝOäj_b´\^_DÝ ZÙYS_D]`\~_bV¼ß_bZ7\^Þ{aWê&Uì_DT]WYSðñ_bZÙaW² 3^y C @ ë ® y 3 y f 5 _bZÙ´<_yä+ßYyT¥]`YSZÙY[´ @ @ y Üãa ]WYyXÇTÙY[´Yy°^\^X}\~Ý`_D]WaTðñ_L[YyÝ`]`abîTÙY/U`X}_:XKßYy]`\^_TÙaDV´7YaÜãY[´lß a]Wa:]WY 5 ²aDV`ZÙYyÝ`]WabîT7YGUX äb_D°±aD´ @ 3{y íUWY}Üãa ]WY2TÙY[´àÞ{aÝT7\^]WY[´<_D]`auÞ{aX}aìa Yy°±YyX}YyÝ ZÙa»Ü`´ßaDÜ`´<\^a÷_DT7T7aÞ[\^_D]`aì_baìÜ`´\^X2Yy\±´7a÷ÚD´<_DU ]`Y C ® ]`Yu°~\±VãY[´<]`_D]WYDêVcM\^Ý`_D°^X2YyÝOZÙYD²_ËTÙa°~UMåÞ _baË]WY ìgyðaD´ÝWYyÞ{YõaéYy°±YyX2YyÝOZÙaËÜ`´ßaDÜ`´<\^a»_DT7TÙa Þ[\^_D]WaË_ba T7Y[ÚU`Ý`]Wa2ÚD´<_DUì]WYe°^\±VãY[´<]`_D]`YDê O TÙZÙY÷ðaD´X}_D°^\^T7X}a Üãa]WYoTÙY[´àYyT7ZÙYyÝ`]`\^]Waó_ìíU`_D°^íUWY[´àT<\^TÙZÙYyX}_óÞ{aXm%&ÚD´<_DU`T]WYo°^\^VrY[´_D]WY ý Ô Ô Ó Ó XsY X Y tz×RQz×S|¤[º8!FºK¸{¹uZ Z $Z º} ! [ Õ ß O @ C Z7_D\^TeíUWYLÜ_~z< Þ_Þ_Þ^< Ü < Ü ê aìÞ[_DTÙa»]`aTeþÛ´ ?:CADCËßYu]`_õaD´<]WYyX>]WYì·,xb¶ó_DÝWaT[²QaìÜãY[´lß a ]Wa ]WYLJ Y÷]`_õaD´<]`YyX>ýß zxb¶D¶ó_DÝWaT[²MY}aìÜãY[´lß a]Wa ]`Y ï ® O Y]`_õaD´]WYyX ]WY»·[¶ß _DÝWaT[êuÿ¬äj_DÝOZ7_bÚDYyX>]`YyTÙZÙYðñaD´ íOUWY¼anfK DY[´<Ý`Yy°há]WaXKßY[ZÙa]`a]WY:û¼aD´\SßYxX_D\^TÞ{aX}Ü°±Y[ZÙaà]`aàíUWYì87j²yçWß_íOU`Y¼Þ{aÝ ZßYyX YyXàVUWZ7\~]Wa O C ZÙY[´ ® ® Ø aT ® ]WY¥öGßUWÜ\±ZÙY[´y²b_DT7T7\~X Þ{aX2aaTQZÙY[´ ® öDalä \^_DÝWaTyê T7ZÙaÄÜãa]WY[´\^_óÚDY[´<_b´o_K_bÜ_b´<\gåÞ _ba ]`Yõ´7YyT7TÙaÝæ_DÝÞ[\^_DTTÙYyÞ[U°^_b´7YyT2YyÝOZÙ´7Y aTÚD´<_DU`T]`Y ® ® °^\±VãY[´<]`_D]WY/]`aÜ´7aDV°±YyX}_ ²a)íOUWY\^ÝOäb_D°^\^]`_b´<\~___bÜ´7ayïW\^X}_råÞ _ba}_D]`\~_bV¼ß_bZ7\^Þ[_]`_)TÙa°^UMåÞ _baWê+¯QaD´ ßYyXõ² T7_bVãYyX2aT/_}ÜC_b´7Z7\±´¼]WYáaUWZÙ´7aT/YyT7Z7U`]WaTàø»UìaD´<_D\^T:µb¶D¶W·lüxíOUWYà_IßU`Ý`\^Þ[_)´ C/ _DÞ{aÝ ZÙYyÞ{Y Ý_ ´7YyT ® O TÙZ7_à´7YyT7TÙaÝQæ_DÝ`Þ[\^_)X}_DÝ`\±ðñYyTÙZ7_jîTÙYÝ`_DT¥ÚD´<_DÝ]WYyTG_DX2Ü°^\±Z7U]WYyT]WY°^\^V`´<_råÞ _baW²aÝ`]WYíOU`_DT7YÝ`YyÝ`úOUX þÛ´7a\~_DÝWa2´7Yy_D°:ßYeaDVT7Y[´7äj_D]`aWê . ' {&þYXÅþG )+æ)+¢¥¤§¨/'1^ ¦,¨3^ ¦x.0¨/'z) æ , 9 9 68 ßxT7_DÝ]WaKaÖÝWaT þÛ´7a\~_DÝWaT+]WYxöGßUWÜ\±ZÙY[´l²\^Ý`Þ[°^U`\~Ý`]Wa:aT+aDV çÙY[ZÙaTÝU`X2Y[´_D]WaT+Y¼XU`°±Z7\^aDÜraT<\^Þ[\±aÝ`_D\^TMÞ{aÝ`ú`YyÞ[\^]WaT+YyX ® ® ]WY_[YyXáV´7a2]WYµb¶D¶D¶ ê¥ÿ Ü`´ aD´<]WYyX ]`Yචê^·+ì øñ´yê Xõê T[ê ü²raoíUWYà´7Y[Ü`´7YyTÙYyÝOZ7_÷U`XY[´7´ v w íOUWYauaDV`Z7\^]Wau_÷Ü_b´7Z7\±´]WYYyTÙZ7U`]`aTÜUW´<_DX2YyÝOZÙY)ÝOU`XÆßY[´\^Þ{aT}øuTG\±YyÝéY Þ<úUWV_b´7Z)· y âUì\^°~_DÝ`\ vzv · ýüê ë ` aX>YyTÙZÙYyT)Yy°±YyX2YyÝOZÙaT)Ü`´ßaDÜ`´\±aT[²M\^]WYyÝOZ7\ Þ[_DX2aT_DTðñ_DX ß °^\^_DTá]`Yo_DTÙZÙY[´ßa\~]WYyT[²+_bZÙ´<_lä+ßYyT]`a ® vzv XKßY[ZÙa]WaË]WYuÿxÚ°±aX2Y[´<_råÞ _baÄû:\±Y[´ß_b´<íU`\^Þ[_Æø'®ã_bÜ`Ü_D°Ja_ Y[Z2_D°êÖ· ¶ üê ÿ:Þ<ú`_DX2aT2Þ[\~Ý`Þ{a»ðz_DX ß °~\^_DT ® YyX¨ZÙaD´<ÝWaK]WY³¡ Y»_bÜrYyÝ_DT}]`U`_DTðñ_DX ß °^\^_DT)YyX¨ZÙaD´ÝWa ]WY³¡ êÀÿ:TX}_D\^T}´7Yy°±Y[äb_DÝ ZÙYyT÷T _ba _DT î ï O w v ðñ_DX ß °^\^_DT]WYË·,k î ® O X2YyXáV´7aT´ ® Ø þÛ´7a\~_DÝWaT]`Yl¡ _ba2°±aÝWÚDa2]`_\^]`_D]WY:]Wa \^TÙZÙYyX}_ a°^_b´yê a)YyÝ Z7_DÝOZÙaW²Ý _baða\rÜãaT7Tyß äDYy°]WY[ZÙYyÞqî î ® Z7_b´e_bÚ°±aX2Y[´<_råÞ aDYyTT7\^X}\~°^_b´7YyTxYyXZÙaD´<ÝWa ]WY5¡ ê¯M_b´<_õY{ï Ü°^\^Þ[_b´\^TÙZÙaW²-UX}_u_D°±ZÙY[´<Ý`_bZ7\±äb_õTÙY[´<\~_ ² ï T7\^X2Ü°^YyT7X2YyÝOZÙYD²íOUWY/Þ[_D]_eÚD´ ðñ_DX ß °^\^_DT9YyX ¡¦ï+Z7_DXáVGßYyX Üãa ]WY[´<\^_/T7Y[´QÞ{aÝ`T7Yyí¦¥U9æYyÝ`Þ[\~_:]`YY[äDYyÝOZ7U`_D\^TM]`\±ðñY[´7YyÝMåÞ[_DTQÝ`_xYyTÙZ7_bV\^°^\~]`_D]WY _à°±aÝ`ÚDaáÜ`´<_L[a)]WYÞ[_D]`_ÜãaÝ ZÙad¡Û_bÚD´<_DÝWÚDYy_DÝ`aWê R:Yxðz_bZÙaW²WYyTÙZ7_DTG]`\±ðñY[´7YyÝMåÞ[_DTY{ïW\^TÙZÙYyX6ÝWauæ_DXàV\±ZÙa ® vzv vzv ]`_á]`\~Ýæ_DX}\^Þ[_àÝ _babîÞ{aÝ`TÙY[´7äb_bZ7\±äb_õø»UìUW´<´<_,£ì· PâBo:aX2YyT· yü²`Y:Üãa]WY[´\^_DX Þ{aÝOZÙ´<\±VU\±´Ü_b´_ ýzx ·sÓDÔ Ó Ó Ô Y X Y X ¹ [Zt]× ºKZ_¸ ¸¡Ú Öh ÖÚ(Z"!$# ¬¹»º Õ ¼ ¼ Y[äb_bÜãaD´<_b´2_DT)ðñ_DXìß °^\^_DTàðaD´X}_D]`_DTYyX ¡ ²X}_D\~Tà´_bÜ\^]`_DX2YyÝOZÙYu]WaËíOUWYu_DT)ðñaD´ ï î ÿ ]`\±ðñY[´7YyÝMåÞ[_ ÝWaéÝßU`X}Y[´7a»]WYoðñ_DXìß °^\^_DTáYyÝ ZÙ´7Y¡ YS¡ Üãa ]WY[´<\^_õYyT7Z7_b´ä \^Ý`Þ[U°^_D]`_ _baóðñ_bZÙaË]`_ î ï ® ÜãaDÜU°^_råÞ _ba»aDVCTÙY[´7äb_D]`_Ë]WYQ¡¤îuTÙY[´}X}_D\±aD´}]`aéíUWYõ_é]WYQ¡¦ïb²+aÄíUWYõÞ{aÝ`TÙZ7\^Z7U`\GU`X Ü`´7aDV°±YyX_ _D\~Ý`]`_YyX _bVãY[´7ZÙaWê {èþ|{&þ .n[/¦ ¨^ ¢)I¤, ¦x. )T)'»¤ 2. . )T ¨/'S¦x¨/. ¢¨¥' 90 9 ¯M_b´<_óYyÝ`Þ{Y[´<´<_b´Ý`aT7TÙaóYyTÙZ7U`]Waé]`_»´7YyT7TÙaÝæ_DÝÞ[\^_K·lô·b²ðñ_L[YyX2aT2UX}_»_Dݼß_D°^\^TÙY÷]`_ìYyT7Z7_bV\^°^\^]_D]WY ® ]`YyTÙZ7_ ´7YyT ® Y÷íOUWYD²Þ{aX2a»Þ{aÝ`T7Yyíß ¦¥U9æYyÝ`Þ[\~_ ]`_ X}\^ÚD´<_råÞ _baW²ÛaTáÜ°^_DÝWY[Z7_DTZÙY[´\^_DX _bZÙ´<_läDYyT7_D]WaéÝWaõÜ_DT7_D]WaìÜraD´ äGß_b´<\^_DT+´7YyT7TÙaÝæ_DÝÞ[\^_DTX UWZ7U`_DT]WY/X2alä \^X2YyÝOZÙaTXÆßß Yy]\±aT[êÿIÜ_DT7T7_bÚDYyX ÜraD´YyT7Z7_DT+´7YyT7T7aÝæ_DÝ`Þ[\^_DT[² aU_Y[äDYyÝ Z7U_D°CÞ[_bÜ`Z7U`´<_ZÙYyX2ÜãaD´ß_b´<\~_Ý`_DT+X2YyT7X_DT[² ðñY_/aT+ÜC°^_DÝWY[Z7_DTY[äDa°^U\±´+ZÙYyX2ÜãaD´<_b´<\~_DX2YyÝ ZÙY ` YyXßaD´7V\±Z7_DTÞ[_ßaDZ7\^Þ[_DTy² ÜU`]`YyÝ`]Wa:ÚDY[´<_b´\^Ý`T7Z7_bV\^°^\^]_D]WYyTÛT7\^ÚÝ`\ Þ[_bZ7\±äb_DTMÝWaT_DTÙZÙY[´ßa\~]WYyT+þ´7a\^_DÝWaT[ê ® ` ÿ ÜãY[´7ÚU`ÝOZ7_KßYDòíU _baT7\±ÚÝ\ Þ[_bZ7\±äb_DTDê ® ¯Q_b´_ì´7YyT7ÜraÝ]WY[´á_ó\^TÙZÙaìðz_L[YyX2aT)U`X}_»TOßY[´<\±Yu]WYoT<\^XàU°^_råÞ aDYyTàÝU`XÆßY[´<\~Þ[_DTáaÝ`]`Y÷öGßU`Ü\±ZÙY[´áY ® _bZ7UW´<ÝWaìT _baóÞ{a°±a Þ[_D]WaTnfÙ_D]Kú`aÞ_h»YyX>]\±ðY[´ ® \~TÙZÙa÷_bðñY[Z7_õ_÷Y[äDa°^UMåÞ _baõ]`\^ÝQæ_DX}\^Þ[_]`YàÜ_b´ ® ® ß O \^X}ÜraD´ ® Þ{a°^aÞ[_DX2aT)Yy°^YyTàÝ`_DTÜãaT7Tyß äDYy\^Te´7YyT ® °^Y[äj_DX YyXÞ{aÝOZ7_Ë_DT)ÜãY[´7Z7UW´7V_råÞ aDYyT2]`aT2íOU`_bZÙ´ ® T _baùÞ{a°±aÞ[_D]`aTõYyX$´7YyT7TÙaÝQæ_DÝ`Þ[\^_ÍX UWZ7U`_ êß cQ\~Ý`_D°^X2YyÝOZÙYD²:]`Y[äDYyX2aT °±Y[äb_b´õYyXèÞ{aÝ Z7_FíOUWYéaT ® Ý`aT7TÙaTQ´7YyT7U`°^Z7_D]WaTMT _bae\^ÝOZÙY[´7Ü`´7Y[Z7_D]WaTM_DT7T ® íU`_DÝ`]`a}_2X}\±ÚD´<_råÞ _ba}a Þ{aD´7´7YyU-ê ® Ø w _DT÷ÝWaT7T7_DTuT7\~XàU`°~_råÞ aDYyT}YyTÙZ7U`]_DX2aTZÙ´yæYyTu´ ® v vzvzv Þ{aX aTX2a]`Yy°±aTÞ[°yß_DT7T7\^Þ{aT]`Y)X\±ÚD´<_råÞ _ba øucWY[´<ݼß_DÝ`]WY_2Y Üù· y Ø YyT7Z7_DTG´7YyT7TÙaÝæ_DÝÞ[\^_DT¥Üra ]WY[´<\~_DX ZÙY[´ST7\^]`a_bZÙ´<_läDYyT7_D]`_DTSÜãYy°±a)T7\^T7ZÙYyX}_öGßUWÜ\±ZÙY[´7î _bZ7UW´<ÝWaWê aT7TÙaT v ´ ® ]`Y·,xb¶D¶÷Ü_b´7Zyß Þ[U`°~_DT/]WYáZÙYyT7ZÙY\^Ý\^Þ[\^_D°^X2YyÝOZÙYeÝ`_2´ ë ® xb¶D¶D¶_DÝ`aT[ê aXÇaTMÜ°^_DÝWY[Z7_DTÝ`_x´7YyT7T7aÝæ_DÝ`Þ[\^_xôbµ² _:´7Y[Ú\ _bae]`aTQþ´7a\^_DÝ`aTáßYGYyTÙäb_Ly\^_D]`_íU`_DTÙY ZÙaDZ7_D°~X2YyÝ ZÙY¼YyX?_bÜ`´7ayïW\^X}_D]_DX2YyÝ ZÙY·[¶ _DÝWaTyêÅcM\^Ý_D°^X2YyÝOZÙYD²Þ{aX aTÜ°^_DÝWY[Z7_DTMÝ`_:´7YyT ® w ôbý²_x]WY[Ü°±Y åÞ _baõßYVãYyX©X}_D\~T-ðñ´<_DÞ[_ ²bYíOU_DTÙYÒkb¶ì©]`_DT-ÜC_b´7Zyß Þ[U`°^_DT]WY+ZÙYyTÙZÙYG\^Ý`\~Þ[\^_D\^TãTÙaDV`´7Y[ä \±äDYyX ® ® Ý_)´7Y[Ú\ _ba]WaT¼þÛ´7a\~_DÝWaT/_D\^Ý`]_2_bܼßaTx2X}\~°^ú aDYyTS]WYe_DÝWaTyê ® ® Ø aT7T<_÷Þ{aÝÞ[°^U`T _baIßY2íOU`Y)öGßUWÜ\±ZÙY[´:Y _bZÙ´ ýzk Ó Ó Ô Ó Ó Ó Ó) Ó Ó X Y X X tz×Rtz×SÍ\ ÚzÛ Z [ ¸{º*¹ ººD¹ Í[ º Z v Z" ¸NZ Ù ½ ½ ½ ® ]`UW´<_DÝOZÙY_÷X}\±ÚD´_råÞ _baWêá¯QaD´ßYyXõ²Û_÷Y{ï \~TÙZ{æYyÝ`Þ[\^_u]WaTþÛ´ w Y[äDYyÝ Z7U`_D°MÜC_DT7T7_bÚDYyX ]WaT:Ü°~_DÝWY[Z7_DTÜãYy°^_´ ´<\^_DX TÙaDV`´ w ý ·sÓDÔ Ó Ó Ô Y X Y X ¹ [Zt]× ºKZ_¸ ¸¡Ú Öh ÖÚ(Z"!$# ¬¹»º Õ ¼ ¼ ýzy ,)(à1 r( ¢>57./,)./; ïã5B./;$¢ 1á02(;¥y1¼*õ)(Ù¢ 1&9571x, 1¥×:)( ì*91G¡./, .&I§ r(+*91:570 *91÷(+* (;7(¢ (,m?./,)0 ¥y.-¢ 5;<5<( v ý ·sÓDÔ Ó Ó Ô Y X Y X ¹ [Zt]× ºKZ_¸ ¸¡Ú Öh ÖÚ(Z"!$# ¬¹»º Õ ¼ ¼ P¶ Icarus 153, 391–415 (2001) doi:10.1006/icar.2001.6699, available online at http://www.idealibrary.com on A Semianalytical Model for the Motion of the Trojan Asteroids: Proper Elements and Families C. Beauge´1 Observatorio Astronomico,´ Universidad Nacional de Cor´ doba, Laprida 854, 5000 Cor´ doba, Argentina and F. Roig Instituto Astronoˆmico e Geof´ısico, Universidade de Sao˜ Paulo, Avenida Miguel Stefano 4200, 04301-904 Sao˜ Paulo, SP, Brazil Received January 29, 2001; revised June 21, 2001 1. INTRODUCTION In this paper we develop a semianalytical model to describe the long-term motion of Trojan asteroids located in tadpole orbits The Trojans (or Jupiter Trojans) are a group of asteroids loca- around the L4 and L5 jovian Lagrangian points. The dynamical ted in the vicinity of Jupiter’s orbit when observed in the semima- model is based on the spatial elliptic three-body problem, including jor axis domain. Their orbits are characterized by an oscillation the main secular variations of Jupiter’s orbit and the direct pertur- of the angle 1 around one of the equilateral Lagrangian bations of the remaining outer planets. Based on ideas introduced points L4,L5( and 1 are the mean longitudes of the asteroid by A. H. Jupp (1969, Astron. J. 74, 35–43), we develop a canoni- and Jupiter, respectively). These orbits are usually referred to cal transformation which allows the transformation of the tadpole as “tadpole” orbits due to the shape of the zero-velocity curves librating orbits into circulating orbits. The disturbing function is of the three-body problem in the synodic reference frame. The then explicitly expanded around each libration point by means of a Trojan swarms can be further divided into two groups: the Taylor–Fourier asymmetric expansion. Making use of the property in which the different degrees of free- Greeks, orbiting the L4 point, and the genuine Trojans, orbit- dom in the Trojan problem are well separated with regard to their ing L5. periods of oscillation, we are able to find approximate action-angle Even though theoretical studies of equilateral equilibrium variables combining Hori’s method with the theory of adiabatic in- configurations of the three-body problem date back to Lagrange variants. This procedure is applied to estimate proper elements for in the late eighteenth century, the first asteroid in such a loca- the sample of 533 Trojans with well determined orbits at December tion was observed only in 1906. It was later designated as (588) 2000. The errors of our semianalytical estimates are about 2–3 times Achilles and was found orbiting L4. The same year a second larger than those previously obtained with numerical approaches body, (617) Patroclus, was found in L5. Since then, an ever in- by other authors. creasing number of asteroids have been discovered. The present Finally, we use these results to search for asteroidal families number of Trojans (as of December 2000) with well determi- among the Trojan swarms. We are able to identify and confirm ned orbits contained in the Asteroids Database of Lowell Ob- the existence of most of the families previously detected by Milani (ftp.lowell.edu/pub/elgb/astorb.dat) (1993, Celest. Mech. Dynam. Astron. 57, 59–94). The families of servatory amounts to 533. If we include the bodies with poorly determined Menelaus and Epeios, both around L4, are the most robust candi- dates to be the by-product of catastrophic disruption of larger as- orbits, this number grows to more than 800. Even this number teroids. On the other hand, no significant family is detected around may be only the tip of the iceberg. Levison et al. (1997) predict L . c 2001 Academic Press that as many as two million asteroids (with size larger than one 5 Key Words: asteroids, dynamics; asteroids, Trojans; celestial kilometer) may in fact lie around these points, thus rivaling the mechanics; resonances. population found in the main belt. The basic idea of the present study is to develop a semiana- lytical model for the dynamics in the 1 : 1 mean motion reso- nance. We focus on the long-term behavior of bodies in tadpole orbits in the vicinity of the equilateral Lagrangian points, and we present some results concerning the long-term dynamical evolu- 1 Present address: Instituto Nacional de Pesquisas Espaciais, Avenida dos tion of these asteroids. Our aim is to construct a model that can Astronautas 1758, (12227-010) Sao˜ Jose´ dos Campos, SP, Brazil. be applied to a qualitative study of the resonant structure in the 391 0019-1035/01 $35.00 Copyright c 2001 by Academic Press All rights of reproduction in any form reserved. 392 BEAUGE´ AND ROIG tadpole regime (e.g., secular resonances inside the libration zone asteroids, found several cases of what he refers to as “stable of each Lagrange point), but the main goal of the present paper chaos,” i.e., orbits with positive Lyapunov exponents (in some is the determination of proper elements for all known Trojans. cases even quite large), but which do not exhibit any gross insta- Although we center our study on Jupiter Trojans, it is worth not- bility over very long timescales. Only very few asteroids were ing that the scope of this work is not restricted to this subsystem. found with significant instabilities on timescales of the order of Specifically, one of the main advantages of analytical studies is 106 yr, all of them lying in the vicinity of secular resonances of the universality of the model: It may be applied to the case of the node. Similar results were also found by Pilat-Lohinger et al. any other perturber. (1999) and by Marzari and Scholl (2000). Giorgilli and Skokos We note that, up to now, all the determinations of Trojan (1997) used Nekhoroshev theory to study the stability of the proper elements have been carried out numerically. Perhaps the tadpole libration region in the Sun–Jupiter–asteroid model. Al- most complete study to date is due to Milani (1993). He per- though they found a zone around the Lagrange points which is formed a numerical integration over timescales of 106 yr, and effectively stable over the age of the universe, this region is too applied a Fourier analysis to the output to determine the funda- small and only includes a few of the real Trojans. Most of the mental frequencies of the free oscillations and their amplitudes present Trojan population seems to lie outside this stable region, (i.e., synthetic theory). These latter constitute his proper ele- which implies that chaotic diffusion and global instability can- ments. A similar approach was recently employed by Burger not be ruled out for these asteroids. A different type of study was et al. (1998) and Pilat-Lohinger et al. (1999), although their undertaken by Levison et al. (1997). They performed numeri- study was more concerned with chaotic orbits. On the other cal integrations of real and fictitious bodies over timescales of hand, analytical and semianalytical approaches have only been 109 yr and for various initial conditions. Although a number of made in the case of main belt asteroids (e.g., Williams 1969, particles showed lifetimes much shorter than the age of the So- Milani and Knezeˇ vic´ 1990, 1994, Lema´ıtre and Morbidelli 1994, lar System, about 90% of the initial conditions compatible with Knezeˇ vic´ et al. 1995, Knezeˇ vic´ and Milani 2001). the Trojan swarm survived the complete simulation, exhibiting From a theoretical point of view, proper elements are inte- (apparently) stable behavior. Thus, although at present we are grals of motion of a given dynamical system. They are values of not able to provide a precise quantification of the stability of the certain actions of the system that remain constant in time (see real Trojans, there is a certain confidence that whatever instabil- Lema´ıtre 1993 and Knezeˇ vic´ 1994 for discussions). However, in ity or chaotic diffusion exist must be extremely slow. For most practice almost all real dynamical systems are nonintegrable, so of the Trojan population, the present dynamical behavior should these integrals of motion do not really exist. In the best case we remain essentially invariant for timescales of at least 108–109 yr. can find only quasi-integrals, which are only approximately con- Over these timescales, proper elements are certainly meaningful stant in time, provided the chaos is sufficiently localized and slow and can be considered good indicators of the original dynamical so that the calculated proper elements still make sense and con- structure. tain meaningful information about the dynamics. In some sense, The present paper is organized as follows: Section 2 presents studies of proper elements constitute a kind of “archaeology” of an application of the method of adiabatic invariance to Hori’s the Solar System. Within the present observed distribution of averaging method in canonical systems (Hori 1966). This is the bodies we search for relics of their original dynamical structure: approach that will be used to determine the solution of the sets parameters that have remained almost unchanged for hundreds of canonical transformations, eventually leading to the determi- of millions of years. And it is through these relics that we hope nation of the proper elements. In Section 3 we introduce the to deepen our understanding of the origin of these bodies. general semianalytic expansion of the Hamiltonian of the three- When determining proper elements for a given dynamical sys- body problem in the vicinity of the 1 : 1 resonance. Section 4 tem, the main sources of approximation for the results are: (i) discusses the hierarchical separation of the different degrees of nonintegrability (chaos) of the system, which causes the real mo- freedom of the problem and their successive elimination. De- tion of the system to be nonquasiperiodic in nature; (ii) very long termination of proper elements follows in Section 5, as well as period perturbations, resulting from quasi-commensurabilites their comparison with previous studies. Identification of possi- between the different frequencies of the system; and (iii) limi- ble asteroidal families is treated in Section 6. Finally, Section 7 tations of the dynamical model, which include approximations is devoted to conclusions. both in the analytical model and in the model of the Solar Sys- tem. While the two latter approximations generate periodic vari- 2. AVERAGING METHODS WITH ations of the calculated proper elements (thus influencing only ADIABATIC INVARIANCE their precision), the first one concerns their very existence and can make them meaningless. Many characteristics of the Trojan asteroids complicate the The question that arises is: What is the magnitude of the chaos elaboration of an analytical model for their long-term motion. in the Trojan belt? Although this question has been addressed They have moderate-to-low eccentricities (e 0.2), but they several times in past years, we still do not have a complete an- can reach very high inclinations (up to 40 ) with respect to swer. Milani (1993), in his study of proper elements of Trojan Jupiter’s orbit. The resonant angle 1 may show large = SEMIANALYTICAL MODEL FOR MOTION OF TROJAN ASTEROIDS 393 amplitudes of libration: D 40 (we define D as half the dif- where the first-order generating function has the form ference between the maximum and minimum values of , i.e., D ( max min)/2). Moreover, in several cases the longitude Ai, j,k,l i j √ 1(k l ) B √ 1 (J ) (J ) E 1 2 , (4) of perihelion= $ can show what is usually referred to as “kine- 1 1 2 + = i, j k,l 0 k 1 l 2 matic” or “paradoxal” libration. In other words, $ does not take X X6= + all the values from zero to 2 , although from the topological and the new Hamiltonian is F (J , J ) F (J , J ). The new point of view the motion is related to a circulation. Fortunately, 1 2 0 1 2 frequencies , are given, in terms of=the new actions, as there is one feature that counteracts these difficulties. It refers to 1 2 the fact that the different degrees of freedom of the system are ∂ F (i 1) j well separated with respect to their periods. In other words, while iAi j00(J ) (J ) 1 ∂ J = , , , 1 2 the period of libration associated with the resonant angle is 1 i, j X (5) typically about 150 yr, the period of oscillation of the longitude ∂ F i (j 1) of perihelion $ is of the order of 3500 yr, while the period of 2 jAi,j,0,0(J1 )(J2 ) . 5 6 ∂ J2 = i, j the longitude of the ascending node is even longer: 10 –10 yr. X It is this property we exploit in the modeling of our problem. It is worth noting that B1 is a function of order 1, so the In this section we concentrate on the development of a gen- eral procedure to analyze multidimensional Hamiltonian sys- difference between (J, )and (J , )isalso of order . More- tems having this kind of “hierarchical” separation of the differ- over, the difference between the old unperturbed frequencies i ent degrees of freedom. In the next section the results obtained in and the new ones i is of order , and whenever 1 and 2 are this manner will be applied to the particular case of the Trojans. not commensurable, so should 1 and 2 be. In other words, we can assure that up to order there are no small integers k, l 2.1. System with Two Degrees of Freedom such that k 1 l 2 0. Then, the solution of system (3) can be found by an iterati+ ve'procedure of successive approximations. We begin by supposing a generic two-degree-of-freedom sys- For the discussion that follows, it is important to mention that tem defined by a Hamiltonian function F, it is not necessary to restrict Hori’s method to the first order in . Nevertheless, all the calculations are easier and the procedure F F(J, ) F0(J1,J2) F1(J1,J2, 1, 2), (1) = + developed in our model becomes much clearer. where is a small parameter and (J, ) are action-angle vari- 2.2. Hierarchical Structure and Adiabatic Invariance ables of F0. We will assume that the unperturbed frequencies i ∂ F0/∂ Ji are finite and large, i.e., neither 1 nor 2 are Let us concentrate on the solution of system (3) and the search = resonant angles. Furthermore, we will suppose that there are for the new action-angle variables (J , ). Let us forget, for the no significant commensurabilities between these frequencies. time being, that this system corresponds to a canonical transfor- In other words, there are no small integer values k, l such that mation and think about it as a system of four algebraic equa- k 1 l 2 0. Then, we can solve this system by using a clas- tions corresponding to two different sets of variables (degrees + ' sical averaging process such as Hori’s method. We search for of freedom). Knowing (J1, J2, 1, 2), we wish to determine a generating function B(J , )ofthe transformation (J, ) (J1 , J2 , 1 , 2 ) so as to satisfy (3). Instead of taking all equa- → (J , )tonewcanonical variables (J , )such that the trans- tions simultaneously, we will divide the system into two parts, formed Hamiltonian is F F (J ). In order to perform all adopting a hypothesis of adiabatic invariance as follows. Let calculations explicitly, imagine= that we have F written as a trun- us suppose that the unperturbed frequencies of each degree of cated Fourier–Taylor series of the type freedom satisfy the condition i j √ 1(k l ) F(J, ) A (J ) (J ) E 1 2 , (2) 1 2, (6) i,j,k,l 1 2 + = i,j,k,l X and let us introduce the small parameter 2/ 1.Weknow, x where E exp(x) and the coefficients Ai j k l (constant with = = , , , for example, that this kind of relation holds in the Trojan case. respect to the variables) are, for the time being, undetermined. Since the first degree of freedom is much faster than the second, The transformation equations between both systems of variables, we can solve it separately, assuming fixed values for (J2 , 2 ), up to first order, are given by and writing the solution in terms of these values. Let us rewrite the subsystem corresponding to the first degree J1 J ∂ B1/∂ = 1 + 1 of freedom explicitly as 1 ∂ B1/∂ J = 1 1 (3) J J ∂ B /∂ ∂ √ 1Ai, j,k,l √ ¯ 2 2 1 2 ¯ ¯ i j 1(k 1 l 2 ) = + J1 J 1 (J 1 ) (J2 ) E + = ∂ ¯ k ¯ l 2 2 ∂ B1/∂ J2 , 1 i, j k,l 0 1 2 ! = X X6= + 394 BEAUGE´ AND ROIG ˆ ¯ ¯ ∂ √ 1Ai, j,k,l √ ¯ where B1(J 1 , 1 ; J2 , 2 ) is the generating function with fixed ¯ ¯ i j 1(k 1 l 2 ) 1 1 (J 1 ) (J2 ) E + , ¯ (J2 , 2 ), corresponding to the single-degree-of-freedom = + ∂ J 1 i, j k,l 0 k ¯1 l 2 ! X X6= + Hamiltonian defined by (7) j √ 1l 2 i √ 1k 1 where we suppose (J2 , 2 )tobefixedparameters (with respect Fˆ(J1, 1;J2, 2) Ai,j,k,l(J2) E (J1) E to time). Since this is not true in the real system (3), we will = i,k ( j,l ) X X designate the results obtained by this approximation as (J¯ , ¯ ). 1 1 i √ 1k 1 Aˆi,k (J1) E , (12) In other words, (J1 , 1 ) will be the real (averaged to the first = ¯ ¯ i,k order) action-angle variables, and (J 1 , 1 ) will be the approx- X imations obtained from (7). Solving this system by iterations, 0 0 with (J2, 2)fixed. and denoting (J1 , 1)asthe initial values at t 0, we obtain = In other words, if we consider (J2, 2)asslowly changing ex- ¯ ¯ 0 0 ternal parameters, then the calculation of the action-angle vari- J 1 J 1 J1 , 1;J2 , 2 = (8) ables (J¯ , ¯ ) of the reduced Hamiltonian (12), expanded only ¯ ¯ 0 0 1 1 1 1 J1 , 1;J2 , 2 . in the first degree of freedom (i.e., with fixed values of J ), = 2, 2 is equivalent up to order to the determination of (J , )from ¯ ¯ 1 1 Now, we need to relate (J 1 , 1 ) to the original (J1 , 1 ). From the the complete Hamiltonian (2). theory of adiabatic invariants (see Henrard, 1993) we know that ¯ The action J 1 is an invariant of the “frozen” Hamiltonian (12), the difference between the solution of the system where (J2 , 2 ) but not of the full Hamiltonian (2). Since (J2, 2) vary slowly are fixed and the one where these quantities are slowly varying ¯ with time, so does J 1 , and according to the adiabatic theory, this with time is of the order of , i.e., variation is such that ¯ K ¯ ¯ J 1 J1 (J 1 , 1 ; J2 , 2 ) d J¯ ∝ (9) 1 2. (13) ¯ L(J¯ , ¯ ; J , ), dt 1 1 ∝ 1 1 2 2 where K and L are functions of order unity (see Henrard 1970 Then, for very small values of this second order variation can ¯ and Henrard and Roels 1974 for a detailed explanation). Thus, be neglected, and the resulting “constant” value of J 1 is called since 1, we can suppose that both sets of solutions are an adiabatic invariant of Hamiltonian (2). It is worth noting approximately the same. In other words, using the adiabatic that these second order corrections to the adiabatic invariant approximation, the new action-angle variables (J , )are de- are periodic with the same period of (J2, 2). Thus, they could 1 1 ¯ termined up to order , and this parameter defines the degree of be eliminated by a suitable averaging of J 1 over a period of precision of the method. (J2, 2). As we see below, averaging the corrections provides a better approach to the adiabatic invariant than neglecting them. 2.2.1. Solution for the first degree of freedom. Let us recall the subsystem (7). Since 1 2 (and so ), and the 2.2.2. Solution for the second degree of freedom. We now 1 2 indices k, l are bounded (recall that Eq. (2) is a truncated series), have expressions for the action-angle variables J1 J1 (J2 , 2 ) and (J )(determined up to order ) corresponding= we can approximate k 1 l 2 k 1 . Therefore, 1 1 2 , 2 + to the first= degree of freedom. Let us pass to the second one. We recall the second half of system (3) as ∂ √ 1Ai, j,k,l i j √ 1(k ¯ l ) J J¯ (J¯ ) (J ) E 1 2 1 1 ¯ 1 2 + ' ∂ 1 i, j k 0 k 1 ! X X6= J2 J2 ∂ B1/∂ 2 = + (14) ∂ √ 1Ai, j,k,l i j √ 1(k ¯ l ) ∂ B /∂ J . ¯ (J¯ ) (J ) E 1 2 2 2 1 2 1 1 ¯ 1 2 + . = ' + ∂ J 1 i, j k 0 k 1 ! X X6= (10) Introducing the solution of (J1 , 1 )into the generating function (4), we have The difference between this expression and (7) is also of the order . We can easily see that (10) can be rewritten as Ai, j,k,l B10 √ 1 = i, j k,l 0 k 1 l 2 ∂ X X6= + ¯ ˆ ¯ ¯ i j √ 1(k 1 (J1 ,J2 ) l 2 ) J1 J 1 B1(J 1 , 1 ; J2 , 2 ) (J1 (J2 , 2 )) (J2 ) E + . (15) ' + ∂ ¯ 1 (11) ¯ ∂ ˆ ¯ ¯ We write B instead of B because (once again) there is a differ- 1 1 B1(J 1 , 1 ; J2 , 2 ), 10 1 ' ∂ J1 ence of order from the hypothesis of adiabatic approximation. SEMIANALYTICAL MODEL FOR MOTION OF TROJAN ASTEROIDS 395 Explicitly, system (14) can be written as where 2 / . Comparing this with system (18), we can 1 = 1 see that the action-angle variables (J2 , 2 ) can be thought of as action-angle variables of the one-degree-of-freedom system ∂ √ 1Ai, j,k,l J2 J2 defined by the Hamiltonian F˜ , if calculated in terms of the av- = ∂ k l 2 i,j k,l 0 1 + 2 eraged variables: F˜ ( J , ). X X6= 2 1 2 1 In other words, wehcani avheragei the original Hamiltonian (19) i j √ 1(k (J ,J ) l ) (J (J , )) (J ) E 1 1 2 + 2 over a reference orbit of the first degree of freedom (which is 1 2 2 2 ! obtained by adiabatic approximation assuming that the second (16) degree is fixed), and then we can use this averaged Hamiltonian ∂ √ 1Ai, j,k,l 2 2 to solve the second degree of freedom. = + ∂ J k l 2 i, j k,l 0 1 + 2 X X6= 2.3. Extension to Many Degrees of Freedom i j √ 1(k 1 (J1 ,J2 ) l 2 ) (J1 (J2 , 2 )) (J2 ) E + , The procedure described above does not need to be restricted ! to two degrees of freedom. Imagine a general case with N de- grees of freedom (J , J ,...,J , , ,..., ), such that the which corresponds to a one-degree-of-freedom nonautonomous 1 2 N 1 2 N unperturbed frequencies of each angular variable are finite system, since is a linear function of time. i i 1 and large, and satisfy the condition of adiabatic invariance with Let us now define the averaging as respect to the previous one. In other words, let us call i j j / i , = and let us suppose that 1 2 N. Once again, the 1 2 X Xd (17) first degree is much faster than the second, the second much h i 1 = 2 1 Z0 faster than the third, an so on. Then, we can repeat the above procedure. First, we solve the Hamiltonian F for (J1 , 1 ) as- and apply it to both sides of (16). To first order in , this yields suming (J2 , 2 ,...,JN , N ) fixed (adiabatic approximation). ¯ We use this solution to average F over 1 obtaining F, we solve ¯ ∂ √ 1Ai, j,0,l F for (J2 , 2 )assuming (J3 , 3 ,...,JN , N )fixed, and we use J2 J 1 2 this solution to average F¯ over obtaining F¯, etc. h i = ∂ 2 i,j l 0 l 2 2 XX6= We can conclude that a general N-degree-of-freedom Hamiltonian system in which the different degrees of freedom i j √ 1l 2 (J1 (J2 , 2 )) (J2 ) E ! are well separated in periods can be approached and solved one (18) degree of freedom at a time. We solve a cascade of single-degree- ∂ √ 1Ai, j,0,l of-freedom Hamiltonians, each of them averaged over the faster 2 1 2 degrees of freedom and with fixed values of the slower degrees h i = + ∂ J2 i, j l 0 l 2 X X6= of freedom. The consequence of this is twofold. First, as in usual perturbation techniques, the modeling does not need to be done i j √ 1l (J (J , )) (J ) E 2 , 1 2 2 2 at once over all the dimensions of the problem, with the re- ! sulting simplification. Second, for each degree of freedom, the “unperturbed” Hamiltonian used to determine the action-angle which corresponds to a one-degree-of-freedom autonomous sys- variables is much more complete than in the usual perturbation tem. theories. On the other hand, there is a price to be paid: The 2.2.3. Reduction to the averaged Hamiltonian. Let us now accuracy of the approximate action-angle variables is directly return to our initial expansion of the Hamiltonian proportional to . i j √ 1(k 1 l 2) 3. APPLICATION TO THE TROJAN ASTEROIDS F(J, ) Ai,j,k,l(J1) (J2) E + , (19) = i,j,k,l X Let us now apply this method to the case at hand. Although and let us introduce the solution to the first degree of freedom the present work is semianalytical in nature, the procedure we adopted for the evaluation of the proper elements requires an J1 J1(t, J2, 2), 1 1(t, J2, 2). Furthermore, let us average = = explicit expression for the Hamiltonian (averaged over short- the resulting expression with respect to 1 and call this new function F˜ . Thus, period terms) for a massless body located in tadpole orbits in the vicinity of equilateral Lagrange points. This, in turn, implies 1 1 finding an expansion of the averaged disturbing function in terms i j √ 1(k 1 l 2) F˜ (J2, 2) Ai, j,k,l (J1) (J2) E + dt, of an appropriate set of variables. = 1 0 i, j,k,l Z X It is worthwhile mentioning that, although significant efforts (20) have been undertaken to find expansions for mean-motion 396 BEAUGE´ AND ROIG resonances p/q with p q, the same is not true in the case of In terms of variables (21), the Hamiltonian of the restricted the Trojans. Practically all6= the analytical studies of the 1 : 1 reso- three-body problem, in the extended phase space, can be written nance use either variational equations or expansions in Cartesian as coordinates (x, y, z) centered at the Lagrange points. This oc- curs because the equilibrium solutions are easier to represent as 2 F n1 L n13 R, (22) fixed points in the (x, y, z) space, and also because the disturb- = 2L2 + ing function in Cartesian coordinates is extremely simple and does not have the limitations of the Laplacian-type expansions, where R is the disturbing function, n1 is the perturber’s mean motion, and k2, this last denoting Gauss’s constant. which are not convergent in the Trojan case. = Although models based on Cartesian coordinates can yield precise results, they are not suitable for the determination of 3.2. Local Variables and Jupp’s Transformation proper elements. For this purpose, it is better to use an ex- In order to apply the averaging method described in Section 2, pansion of the disturbing function in terms of orbital elements we need to have an explicit expression for the Hamiltonian F in (or their canonical counterparts). An example of this kind of terms of nonresonant angular variables. Unfortunately, this is not expansion was recently given by Morais (1999), and it is based the case with Eq. (22), because is in fact a resonant angle and its on a local expansion around the resonant semimajor axis. Nev- unperturbed frequency is close to zero. So, our first step should ertheless, local expansions are not new. Woltjer (1924) devised be to find a canonical transformation from (L, W, T, ,$, ) an asymmetric expansion for the Trojans, although it seems that to new variables (J, W¯ , T¯ , ,$,¯ ¯) where all the angles are this work has been forgotten for almost 80 years. In the next nonresonant. section we present an alternative expansion which is, in many We begin with the first degree of freedom, associated with ways, similar to these. the subset (L, ). First, we split the Hamiltonian function in the form 3.1. The Hamiltonian for the 1 : 1 Resonance The first step in the expansion of the disturbing function con- F F0(L, W, ) F1(L,W,T, ,$, ). (23) = + sists in the choice of an adequate set of variables for the non- averaged system. We adopt the following set Here, F0 is simply the Hamiltonian of the circular–planar case (for which T 0), and F1 contains the remaining terms (includ- ing the dependence= on the perturber’s orbital elements, which is 1; L = implicit). Let us solve this “unperturbed” Hamiltonian F and $ ; W G L 0 = (21) search for its action-angle variables. Usually this could be done ; T H G = via an averaging process such as Hori’s method. But the fact that 1; 3, has an unperturbed frequency close to zero makes this impos- sible, even though F0 is a single-degree-of-freedom system. Let where is the mean longitude, $ the longitude of perihelion, us recall that classical averaging methods are not valid when a and the longitude of the ascending node. Elements belonging resonant angle exists, and in such cases one usually uses nu- to the perturber will be designated with a subscript 1. The mean merical algorithms (e.g., Henrard 1990, Morbidelli and Moons longitude of the planet, 1, is a short-period angle, and will 1993, 1995) to obtain the action-angle variables of F0.How- therefore be eliminated during the first averaging process. The ever, this has the drawback of being very CPU-time consuming canonical conjugates W, T are written above in terms of the (especially when this result has to be further introduced into usual Delaunay variables, and 3 is the conjugate to 1. the remaining degrees of freedom at F1) and does not explicitly We believe variables (21) are the best choice for the elabo- yield the proper element associated with this degree of freedom. ration of a model of proper elements, mainly because the main Here we choose a different approach which, to our under- frequencies of each angular variable, namely , $, , are standing, has several advantages. This approach is based on well separated from each other. An example of this property, a canonical transformation originally devised by Jupp (1969, which is the key piece of this work, is shown in Fig. 1. In this 1970) for the case of the ideal resonance problem. The idea is as way, it is possible to introduce the adiabatic approach to the follows. Let us think about the libration region of a resonance problem, defining the parameters 12 / , 23 / , as a set of invariant curves around a center (i.e., the libration = $ = $ and 13 / , which are all very small. point). If we only concentrate on this region and disregard the Another= advantage of variables (21) is due to the fact that, in structure of the phase space outside the separatrix, we can think the planar–circular problem, $ and practically do not show of these orbits as (distorted) circulations around a center which periodic variations associated with the libration period of . is displaced from the origin of the coordinate system. Now, if we This becomes very important when we perform the canonical find a canonical transformation (L, ) (J, ) that is simply transformation to obtain the action-angle variables of the planar– a translation of the origin to the libration→center, we will obtain circular problem, as we show below. a new angle having a frequency different from zero, and the SEMIANALYTICAL MODEL FOR MOTION OF TROJAN ASTEROIDS 397 FIG. 1. Numerical simulation of (659) Nestor over 5000 years. The different periods of each degree of freedom are noticeable. integral of J along any orbit will be the action of that trajectory. by the initial conditions (remember that W becomes an integral 2 In other words, we will have an angle that librates transformed of motion in the F0 approximation). The values of (Lc, c) into another angle that circulates. These new variables will are nothing but the equilibrium points of F0 and can be easily have properties of being “nonresonant” (even though they are obtained numerically. a simple translation), and we can use an averaging method to For the present work we choose a transformation different determine the action-angle variables. from (24) which, although a bit more complicated (because it Although there are many ways of determining (J, ), possibly can no longer be thought of as a simple translation) is still based the simplest consists of a series of transformations, on the same idea. The change of variables is represented by the relationship (L, ) (K,H) √2L(cos , sin ) → = 1/2 ˆ 2 ˆ 2 1/2 (K, H) (X, Y ) (K Kc, H Hc) (24) X 0 K Kc H → = = (25) (X Y ) √2J(cos sin ) (J ) 1 2 , , , , Y 0 / Hˆ , = → = where (Kc, Hc) √2Lc(cos c, sin c) marks the center of = libration corresponding to a predetermined value W W 0 given 2 Note that H here is not the classical Delaunay variable. = 398 BEAUGE´ AND ROIG of angular variables (Eq. 21) is that 1 and 2 are practically zero. This is because the frequency of does not have any significant contribution to the power spectra of $ and . Thus, we can con- sider the equalities $ $¯ and ¯ without introducing any important error in the =transformation.= Nevertheless, in order to guarantee the internal consistency of the transformation, we will maintain the canonical corrections 1 and 2 in our calculations. In order to simplify the notation, we use the set (W, T,$, ) instead of (W¯ , T¯ , $,¯ ¯). Now, using transformation (26), it is possible to write (23) in FIG. 2. Transformation to local variables as defined by Eq. (25). (a) Group the form of five invariant curves of the planar circular problem, corresponding to librations around L4. (b) Same invariant curves, in terms of the new (X, Y) variables. F F0(J, W, ) F1(J,W,T, ,$, ). (28) = + where (Kˆ , Hˆ ) √2L(cos( c), sin( c)) and 0 0(Kc) = = It is worth recalling that the set (J, W, T, ,$, )iscanonical is a scaling factor which modifies the shape of the trajectories. by construction, and that is a circulating angle with frequency This transformation is canonical and is valid as long as max | . c <