°´¼üÅÌÉÏ·½Ïò¼ü ¡û »ò ¡ú ¿É¿ìËÙÉÏÏ·ҳ£¬°´¼üÅÌÉ쵀 Enter ¼ü¿É»Øµ½±¾ÊéĿ¼ҳ£¬°´¼üÅÌÉÏ·½Ïò¼ü ¡ü ¿É»Øµ½±¾Ò³¶¥²¿£¡
¡ª¡ª¡ª¡ªÎ´ÔĶÁÍꣿ¼ÓÈëÊéÇ©ÒѱãÏ´μÌÐøÔĶÁ£¡
¡¡is¡¡directly¡¡related¡¡to¡¡the¡¡semimajor¡¡axis¡¡a£©¡¡might¡¡be¡¡much¡¡less¡¡affected¡¡by¡¡perturbing¡¡plas¡¡than¡¡is¡¡the¡¡angular¡¡momentum¡¡exchange¡¡£¨which¡¡relates¡¡to¡¡e£©¡£¡¡Hence£¬¡¡the¡¡eccentricities¡¡of¡¡Venus¡¡and¡¡Earth¡¡can¡¡be¡¡disturbed¡¡easily¡¡by¡¡Jupiter¡¡and¡¡Saturn£¬¡¡which¡¡results¡¡in¡¡a¡¡positive¡¡rrelation¡¡in¡¡the¡¡angular¡¡momentum¡¡exchange¡£¡¡On¡¡the¡¡other¡¡hand£¬¡¡the¡¡semimajor¡¡axes¡¡of¡¡Venus¡¡and¡¡Earth¡¡are¡¡less¡¡likely¡¡to¡¡be¡¡disturbed¡¡by¡¡the¡¡jovian¡¡plas¡£¡¡Thus¡¡the¡¡energy¡¡exchange¡¡may¡¡be¡¡limited¡¡only¡¡within¡¡the¡¡Venus¨CEarth¡¡pair£¬¡¡which¡¡results¡¡in¡¡a¡¡negative¡¡rrelation¡¡in¡¡the¡¡exchange¡¡of¡¡orbital¡¡energy¡¡in¡¡the¡¡pair¡£
¡¡¡¡As¡¡for¡¡the¡¡outer¡¡jovian¡¡plaary¡¡subsystem£¬¡¡Jupiter¨CSaturn¡¡and¡¡Uranus¨CNeptune¡¡seem¡¡to¡¡make¡¡dynamical¡¡pairs¡£¡¡However£¬¡¡the¡¡strength¡¡of¡¡their¡¡upling¡¡is¡¡not¡¡as¡¡strong¡¡pared¡¡with¡¡that¡¡of¡¡the¡¡Venus¨CEarth¡¡pair¡£
¡¡¡¡5¡¡¡À¡¡5¡¡¡Á¡¡1010¡yr¡¡integrations¡¡of¡¡outer¡¡plaary¡¡orbits
¡¡¡¡Since¡¡the¡¡jovian¡¡plaary¡¡masses¡¡are¡¡much¡¡larger¡¡than¡¡the¡¡terrestrial¡¡plaary¡¡masses£¬¡¡we¡¡treat¡¡the¡¡jovian¡¡plaary¡¡system¡¡as¡¡an¡¡independent¡¡plaary¡¡system¡¡in¡¡terms¡¡of¡¡the¡¡study¡¡of¡¡its¡¡dynamical¡¡stability¡£¡¡Hence£¬¡¡we¡¡added¡¡a¡¡uple¡¡of¡¡trial¡¡integrations¡¡that¡¡span¡¡¡À¡¡5¡¡¡Á¡¡1010¡¡yr£¬¡¡including¡¡only¡¡the¡¡outer¡¡five¡¡plas¡¡£¨the¡¡four¡¡jovian¡¡plas¡¡plus¡¡Pluto£©¡£¡¡The¡¡results¡¡exhibit¡¡the¡¡rigorous¡¡stability¡¡of¡¡the¡¡outer¡¡plaary¡¡system¡¡over¡¡this¡¡long¡¡time¡span¡£¡¡Orbital¡¡nfigurations¡¡£¨Fig¡£¡¡12£©£¬¡¡and¡¡variation¡¡of¡¡eccentricities¡¡and¡¡inclinations¡¡£¨Fig¡£¡¡13£©¡¡show¡¡this¡¡very¡¡long¡term¡¡stability¡¡of¡¡the¡¡outer¡¡five¡¡plas¡¡in¡¡both¡¡the¡¡time¡¡and¡¡the¡¡frequency¡¡domains¡£¡¡Although¡¡we¡¡do¡¡not¡¡show¡¡maps¡¡here£¬¡¡the¡¡typical¡¡frequency¡¡of¡¡the¡¡orbital¡¡oscillation¡¡of¡¡Pluto¡¡and¡¡the¡¡other¡¡outer¡¡plas¡¡is¡¡almost¡¡nstant¡¡during¡¡these¡¡very¡¡long¡term¡¡integration¡¡periods£¬¡¡which¡¡is¡¡demonstrated¡¡in¡¡the¡¡time¨Cfrequency¡¡maps¡¡on¡¡our¡¡webpage¡£
¡¡¡¡In¡¡these¡¡two¡¡integrations£¬¡¡the¡¡relative¡¡numerical¡¡error¡¡in¡¡the¡¡total¡¡energy¡¡was¡¡¡«10£¿6¡¡and¡¡that¡¡of¡¡the¡¡total¡¡angular¡¡momentum¡¡was¡¡¡«10£¿10¡£
¡¡¡¡5¡£1¡¡Resonances¡¡in¡¡the¡¡Neptune¨CPluto¡¡system
¡¡¡¡Kinoshita¡¡&£»amp£»¡¡Nakai¡¡£¨1996£©¡¡integrated¡¡the¡¡outer¡¡five¡¡plaary¡¡orbits¡¡over¡¡¡À¡¡5¡£5¡¡¡Á¡¡109¡¡yr¡¡¡£¡¡They¡¡found¡¡that¡¡four¡¡major¡¡resonances¡¡between¡¡Neptune¡¡and¡¡Pluto¡¡are¡¡maintained¡¡during¡¡the¡¡whole¡¡integration¡¡period£¬¡¡and¡¡that¡¡the¡¡resonances¡¡may¡¡be¡¡the¡¡main¡¡causes¡¡of¡¡the¡¡stability¡¡of¡¡the¡¡orbit¡¡of¡¡Pluto¡£¡¡The¡¡major¡¡four¡¡resonances¡¡found¡¡in¡¡previous¡¡research¡¡are¡¡as¡¡follows¡£¡¡In¡¡the¡¡following¡¡description£¬¦Ë¡¡denotes¡¡the¡¡mean¡¡longitude£¬¦¸¡¡is¡¡the¡¡longitude¡¡of¡¡the¡¡ascending¡¡node¡¡and¡¡£¿¡¡is¡¡the¡¡longitude¡¡of¡¡perihelion¡£¡¡Subscripts¡¡P¡¡and¡¡N¡¡denote¡¡Pluto¡¡and¡¡Neptune¡£
¡¡¡¡Mean¡¡motion¡¡resonance¡¡between¡¡Neptune¡¡and¡¡Pluto¡¡£¨3£º2£©¡£¡¡The¡¡critical¡¡argument¡¡¦È1=¡¡3¡¡¦ËP£¿¡¡2¡¡¦ËN£¿£¿P¡¡librates¡¡around¡¡180¡ã¡¡with¡¡an¡¡amplitude¡¡of¡¡about¡¡80¡ã¡¡and¡¡a¡¡libration¡¡period¡¡of¡¡about¡¡2¡¡¡Á¡¡104¡¡yr¡£
¡¡¡¡The¡¡argument¡¡of¡¡perihelion¡¡of¡¡Pluto¡¡¦ØP=¦È2=£¿P£¿¦¸P¡¡librates¡¡around¡¡90¡ã¡¡with¡¡a¡¡period¡¡of¡¡about¡¡3¡£8¡¡¡Á¡¡106¡¡yr¡£¡¡The¡¡dominant¡¡periodic¡¡variations¡¡of¡¡the¡¡eccentricity¡¡and¡¡inclination¡¡of¡¡Pluto¡¡are¡¡synchronized¡¡with¡¡the¡¡libration¡¡of¡¡its¡¡argument¡¡of¡¡perihelion¡£¡¡This¡¡is¡¡anticipated¡¡in¡¡the¡¡secular¡¡perturbation¡¡theory¡¡nstructed¡¡by¡¡Kozai¡¡£¨1962£©¡£
¡¡¡¡The¡¡longitude¡¡of¡¡the¡¡node¡¡of¡¡Pluto¡¡referred¡¡to¡¡the¡¡longitude¡¡of¡¡the¡¡node¡¡of¡¡Neptune£¬¦È3=¦¸P£¿¦¸N£¬¡¡circulates¡¡and¡¡the¡¡period¡¡of¡¡this¡¡circulation¡¡is¡¡equal¡¡to¡¡the¡¡period¡¡of¡¡¦È2¡¡libration¡£¡¡When¡¡¦È3¡¡bees¡¡zero£¬¡¡i¡£e¡£¡¡the¡¡longitudes¡¡of¡¡ascending¡¡nodes¡¡of¡¡Neptune¡¡and¡¡Pluto¡¡overlap£¬¡¡the¡¡inclination¡¡of¡¡Pluto¡¡bees¡¡maximum£¬¡¡the¡¡eccentricity¡¡bees¡¡minimum¡¡and¡¡the¡¡argument¡¡of¡¡perihelion¡¡bees¡¡90¡ã¡£¡¡When¡¡¦È3¡¡bees¡¡180¡ã£¬¡¡the¡¡inclination¡¡of¡¡Pluto¡¡bees¡¡minimum£¬¡¡the¡¡eccentricity¡¡bees¡¡maximum¡¡and¡¡the¡¡argument¡¡of¡¡perihelion¡¡bees¡¡90¡ã¡¡again¡£¡¡Williams¡¡&£»amp£»¡¡Benson¡¡£¨1971£©¡¡anticipated¡¡this¡¡type¡¡of¡¡resonance£¬¡¡later¡¡nfirmed¡¡by¡¡Milani£¬¡¡Nobili¡¡&£»amp£»¡¡Carpino¡¡£¨1989£©¡£
¡¡¡¡An¡¡argument¡¡¦È4=£¿P£¿£¿N£«¡¡3¡¡£¨¦¸P£¿¦¸N£©¡¡librates¡¡around¡¡180¡ã¡¡with¡¡a¡¡long¡¡period£¬¡«¡¡5¡£7¡¡¡Á¡¡108¡¡yr¡£
¡¡¡¡In¡¡our¡¡numerical¡¡integrations£¬¡¡the¡¡resonances¡¡£¨i£©¨C£¨iii£©¡¡are¡¡well¡¡maintained£¬¡¡and¡¡variation¡¡of¡¡the¡¡critical¡¡arguments¡¡¦È1£¬¦È2£¬¦È3¡¡remain¡¡similar¡¡during¡¡the¡¡whole¡¡integration¡¡period¡¡£¨Figs¡¡14¨C16¡¡£©¡£¡¡However£¬¡¡the¡¡fourth¡¡resonance¡¡£¨iv£©¡¡appears¡¡to¡¡be¡¡different£º¡¡the¡¡critical¡¡argument¡¡¦È4¡¡alternates¡¡libration¡¡and¡¡circulation¡¡over¡¡a¡¡1010¡yr¡¡time¡scale¡¡£¨Fig¡£¡¡17£©¡£¡¡This¡¡is¡¡an¡¡interesting¡¡fact¡¡that¡¡Kinoshita¡¡&£»amp£»¡¡Nakai's¡¡£¨1995£¬¡¡1996£©¡¡shorter¡¡integrations¡¡were¡¡not¡¡able¡¡to¡¡disclose¡£
¡¡¡¡6¡¡Discussion
¡¡¡¡What¡¡kind¡¡of¡¡dynamical¡¡mechanism¡¡maintains¡¡this¡¡long¡term¡¡stability¡¡of¡¡the¡¡plaary¡¡system£¿¡¡We¡¡can¡¡immediately¡¡think¡¡of¡¡two¡¡major¡¡features¡¡that¡¡may¡¡be¡¡responsible¡¡for¡¡the¡¡long¡term¡¡stability¡£¡¡First£¬¡¡there¡¡seem¡¡to¡¡be¡¡no¡¡significant¡¡lower¡order¡¡resonances¡¡£¨mean¡¡motion¡¡and¡¡secular£©¡¡between¡¡any¡¡pair¡¡among¡¡the¡¡nine¡¡plas¡£¡¡Jupiter¡¡and¡¡Saturn¡¡are¡¡close¡¡to¡¡a¡¡5£º2¡¡mean¡¡motion¡¡resonance¡¡£¨the¡¡famous¡¡¡®great¡¡inequality¡¯£©£¬¡¡but¡¡not¡¡just¡¡in¡¡the¡¡resonance¡¡zone¡£¡¡Higher¡order¡¡resonances¡¡may¡¡cause¡¡the¡¡chaotic¡¡nature¡¡of¡¡the¡¡plaary¡¡dynamical¡¡motion£¬¡¡but¡¡they¡¡are¡¡not¡¡so¡¡strong¡¡as¡¡to¡¡destroy¡¡the¡¡stable¡¡plaary¡¡motion¡¡within¡¡the¡¡lifetime¡¡of¡¡the¡¡real¡¡Solar¡¡system¡£¡¡The¡¡send¡¡feature£¬¡¡which¡¡we¡¡think¡¡is¡¡more¡¡important¡¡for¡¡the¡¡long¡term¡¡stability¡¡of¡¡our¡¡plaary¡¡system£¬¡¡is¡¡the¡¡difference¡¡in¡¡dynamical¡¡distance¡¡between¡¡terrestrial¡¡and¡¡jovian¡¡plaary¡¡subsystems¡¡£¨Ito¡¡&£»amp£»¡¡Tanikawa¡¡1999£¬¡¡2001£©¡£¡¡When¡¡we¡¡measure¡¡plaary¡¡separations¡¡by¡¡the¡¡mutual¡¡Hill¡¡radii¡¡£¨R_£©£¬¡¡separations¡¡among¡¡terrestrial¡¡plas¡¡are¡¡greater¡¡than¡¡26RH£¬¡¡whereas¡¡those¡¡among¡¡jovian¡¡plas¡¡are¡¡less¡¡than¡¡14RH¡£¡¡This¡¡difference¡¡is¡¡directly¡¡related¡¡to¡¡the¡¡difference¡¡between¡¡dynamical¡¡features¡¡of¡¡terrestrial¡¡and¡¡jovian¡¡plas¡£¡¡Terrestrial¡¡plas¡¡have¡¡smaller¡¡masses£¬¡¡shorter¡¡orbital¡¡periods¡¡and¡¡wider¡¡dynamical¡¡separation¡£¡¡They¡¡are¡¡strongly¡¡perturbed¡¡by¡¡jovian¡¡plas¡¡that¡¡have¡¡larger¡¡masses£¬¡¡longer¡¡orbital¡¡periods¡¡and¡¡narrower¡¡dynamical¡¡separation¡£¡¡Jovian¡¡plas¡¡are¡¡not¡¡perturbed¡¡by¡¡any¡¡other¡¡massive¡¡bodies¡£
¡¡¡¡The¡¡present¡¡terrestrial¡¡plaary¡¡system¡¡is¡¡still¡¡being¡¡disturbed¡¡by¡¡the¡¡massive¡¡jovian¡¡plas¡£¡¡However£¬¡¡the¡¡wide¡¡separation¡¡and¡¡mutual¡¡interaction¡¡among¡¡the¡¡terrestrial¡¡plas¡¡renders¡¡the¡¡disturbance¡¡ineffective£»¡¡the¡¡degree¡¡of¡¡disturbance¡¡by¡¡jovian¡¡plas¡¡is¡¡O£¨eJ£©£¨order¡¡of¡¡magnitude¡¡of¡¡the¡¡eccentricity¡¡of¡¡Jupiter£©£¬¡¡since¡¡the¡¡disturbance¡¡caused¡¡by¡¡jovian¡¡plas¡¡is¡¡a¡¡forced¡¡oscillation¡¡having¡¡an¡¡amplitude¡¡of¡¡O£¨eJ£©¡£¡¡Heightening¡¡of¡¡eccentricity£¬¡¡for¡¡example¡¡O£¨eJ£©¡«0¡£05£¬¡¡is¡¡far¡¡from¡¡sufficient¡¡to¡¡provoke¡¡instability¡¡in¡¡the¡¡terrestrial¡¡plas¡¡having¡¡such¡¡a¡¡wide¡¡separation¡¡as¡¡26RH¡£¡¡Thus¡¡we¡¡assume¡¡that¡¡the¡¡present¡¡wide¡¡dynamical¡¡separation¡¡among¡¡terrestrial¡¡plas¡¡£¨&£»gt£»¡¡26RH£©¡¡is¡¡probably¡¡one¡¡of¡¡the¡¡most¡¡significant¡¡nditions¡¡for¡¡maintaining¡¡the¡¡stability¡¡of¡¡the¡¡plaary¡¡system¡¡over¡¡a¡¡109¡yr¡¡time¡span¡£¡¡Our¡¡detailed¡¡analysis¡¡of¡¡the¡¡relationship¡¡between¡¡dynamical¡¡distance¡¡between¡¡plas¡¡and¡¡the¡¡instability¡¡time¡scale¡¡of¡¡Solar¡¡system¡¡plaary¡¡motion¡¡is¡¡now¡¡on¡going¡£
¡¡¡¡Although¡¡our¡¡numerical¡¡integrations¡¡span¡¡the¡¡lifetime¡¡of¡¡the¡¡Solar¡¡system£¬¡¡the¡¡number¡¡of¡¡integrations¡¡is¡¡far¡¡from¡¡sufficient¡¡to¡¡fill¡¡the¡¡initial¡¡phase¡¡space¡£¡¡It¡¡is¡¡necessary¡¡to¡¡perform¡¡more¡¡and¡¡more¡¡numerical¡¡integrations¡¡to¡¡nfirm¡¡and¡¡examine¡¡in¡¡detail¡¡the¡¡long¡term¡¡stability¡¡of¡¡our¡¡plaary¡¡dynamics¡£
¡¡¡¡¡ª¡ªÒÔÉÏÎĶÎÒý×Ô¡¡Ito£¬¡¡T¡£&£»¡¡Tanikawa£¬¡¡K¡£¡¡Long¡term¡¡integrations¡¡and¡¡stability¡¡of¡¡plaary¡¡orbits¡¡in¡¡our¡¡Solar¡¡System¡£¡¡Mon¡£¡¡Not¡£¡¡R¡£¡¡Astron¡£¡¡Soc¡£¡¡336£¬¡¡483¨C500¡¡£¨2002£©
¡¡¡¡ÕâÖ»ÊÇ×÷Õß¾ý²Î¿¼µÄһƪÎÄÕ£¬¹ØÓÚÌ«ÑôϵµÄÎȶ¨ÐÔ¡£
¡¡¡¡»¹ÓÐÆäËûÂÛÎÄ£¬²»¹ýÒ²¶¼ÊÇÓ¢Îĵģ¬Ïà¹Ø¿ÎÌâµÄÖÐÎÄÎÄÏ׺ÜÉÙ£¬ÄÇЩÂÛÎÄÏÂÔØһƪҪ¾ÅÃÀÔª£¨¡¶Nature¡·ÕæÊDZ©Àû£©£¬×÷Õß¾ýдÕâƪÎÄÕµÄʱºòÒѾ»Ø¼Ò£¬²»ÔÚ¼ì²âÖÐÐÄ£¬ËùÒÔûÓÐÊý¾Ý¿âµÄʹÓÃȨ£¬Ï²»Æ𣬾Ͳ»ÌùÉÏÀ´ÁË¡£
¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡
¡¡¡¡µÚÒ»°ÙÈÕ£¨2£©¶Ú»ý·à±ã·¢¼ÒÖ¸»
¡¡¡¡ÌÆÔ¾×øÔÚ·à¶ÑÀÓüôµ¶¼ô¿ª´ó±ãµÄÕæ¿Õ°ü×°£¬ÓÉÓÚÍÑË®ºÍµÍΣ¬´üÖеķà±ãÈ«²¿¸ÉÔï¶øËɴ࣬¿´ÉÏÈ¥¼áÓ²£¬Êµ¼ÊÉÏÇáÇáÒ»µ·¾Í»áÏñɳ×ÓÒ»ÑùËéÁÑ¡£
¡¡¡¡Ëû³¤Õâô´óҲûÏë¹ý×Ô¼º´ËÉú»áÓÐÐÒÓë´ó±ãÇ×ÃܽӴ¥£¬ÒÔÍùÌÆÔ¾¶¼ÊÇÀÍêÌá¿ã³åË®×ßÈË£¬Ï൱ÎÞÇ飬´ÓÀ´²»»á»ØÍ·¿´Ò»ÑÛ¡£
¡¡¡¡ÌÆÔ¾°ÑÕæ¿Õ´üÖеĸÉÔï·à±ãµ¹ÔÚÌúÅÌ×ÓÀËûÉõÖÁ¿ÉÒÔ¸ù¾Ý·à±ãÖеÄʳÎï²ÐÔü·Ö±æ³ö×Ô¼ºÄÇÌì³ÔµÄÊÇʲô£¬Ðí¶à¸»º¬´ÖÏËάµÄʳÎïÔÚÈ˵ij¦Î¸Öеò»µ½³ä·ÖÏû»¯£¬»áËæ×Å·à±ãÒ»ÆðÅųöÀ´£¬ÌÆÔ¾ÓÐÒ»Õó×Ó¾³£³ÔÊ߲˹ÞÍ·£¬ÄǸöʱÆڵĴó±ãÖо;³£ÄÜ¿´µ½Î´ÍêÈ«Ïû»¯µÄÖ²ÎïÏËά¡£
¡¡¡¡ÌÆÔ¾´÷×ÅÁ½²ã¿ÚÕÖ£¬ºÃÔÚÕâЩ´ó±ã¶¼ÒѾ¸ÉÔïÁ˺ܳ¤Ê±¼ä£¬Ã»Ê²Ã´Æøζ¡£
¡¡¡¡¡°·à±ãÖеÄζµÀÖ÷ÒªÀ´×Ô3¡¼×»ùßÅßᣬҲ¾ÍÊÇ·à³ôËØ£¬ÕâÖÖ¶«Î÷Èç¹ûÕ´ÔÚÒ·þÉÏ£¬ÄÇÑùÄãÎÅÉÏÈ¥Ò²»áÏñÒ»Ûç´ó±ã¡£¡±ÀÏè×øÔÚ±ßÉÏ°ï棬ËüÒ»±ßµ·Ëé·à±ãÒ»±ßËéËéÄî¡£
¡¡¡¡ÌÆÔ¾¾õµÃËü¸ü´óµÄÐËȤÔÚÓÚ·à±ã±¾Éí¡£
¡¡¡¡¡