共振Kuiper带天体轨道稳定性的数值模拟 Numerical Simulation on Orbital Stability of RKBOs.pdfVIP

【文章编号】1001．246X(2009)03-0467．嘴 共振Kuiper带天体轨道稳定性的数值模拟 吴晓梅1， 聂清香2 (1．泰山学院物理与电子T程学院，山东泰安271000；2．山东师范大学物理与电子科学学院，山东济南250014) [摘要】 对当前观测到的具有可靠轨道根数的共振Kuiper带天体的轨道未来演化进行数值模拟发现，它们的 轨道在一亿年的时间演化中具有相对稳定性，即仍处于各自的共振区中，但稳定程度不同，且共振天体的轨道 稳定性和其初始轨道半长径、偏心率和倾角有关．根据其轨道半长径、偏心率和倾角随时间的演化行为，这些 共振天体可以分为具有“规则轨道”和“混沌轨道”两种类型． [关键词】 RKBOs；轨道半长径；轨道偏心率；轨道倾角；数值模拟 [中圈分类号】P136 [文献标识码] A O 引言 在黄道面附近，大多数在离太阳30AU一50AU之间的环带内(AU即天文单位，1AU等于地球到太阳的平均 距离)，少数分布在50 Au之外，远的达几百天文单位．Kuiper带是目前观测到的太阳系中最遥远的小天体 群，它可分成两个区域，主带和散射带．主带天体轨道半长径在30—50AU之间，这里的天体也被称为“Trans． Neptunian 构的残留成分存在形式，对Kuiper带天体的研究可以追溯到短周期彗星的来源．由于大行星长期摄动的影 响，Kuiper带的小天体可能穿越海王星的轨道，并与海王星发生交会，其中一部分会向外太阳系抛射；另外一 值积分初步证实了这一看法¨】．他们指出，由Kuiper带小天体演化形成的短周期彗星应比由长周期彗星转变 而成的更多． AU。48 带天体．CKBOs是指拥有小偏心率且远离海王星摄动的Kuiper带天体，通常分布在42AU之间，被认 带天体轨道特性的分析b1发现，在海王星平运动共振区有小天体聚集现象；位于该处的天体普遍具有较大轨 道偏心率．这些共振天体的轨道性质引起学者的普遍关注【6。1¨，有人猜测它们可能是短周期彗星的潜在来 源‘12，1引． 对海王星平运动共振带中天体的动力学结构的形成，讨论较多的是行星迁移理论模型，该模型通过假设 行星的迁移时间和距离以及小天体的初始轨道参数，成功地解释了冥王星的轨道特性1和大量KBOs在海 王星3：2共振带处的聚集[7’1¨，但是，该设想不能解释海王星2：1共振带处天体的分布，也不能解释为什么观 测到的KBOs的半数具有大于150的轨道倾角[1¨，即行星迁移模型当前仍面临着较大的困难．本文的目的不 是讨论RKBOs的过去演化历史，而是数值模拟它们的未来演化，试图由此及彼，从未来的演化过程中探索这 些天体在整个演化过程中的规律．这样做的优点是可利用观测的天体轨道数据作为初始条件，不必用假设参 数，结果较可信． [收穰日期]2007一lO一15；[修回日期】2008—09—02 【基金项目】国家自然科学基金资助项目 【作者简介】吴晓梅(1971一)，女，汉族，山东泰安，副教授．硕士，从事天体力学及动力天文方面的研究，山东泰安迎宾大道中段泰山学院 物理与电子工程学院271000． 万方数据万方数据 计 算 物 理 第26卷 1基本方程和数值方法 数值模拟，研究其轨道半长径、轨道偏心率和倾角的演化状况，及其轨道稳定性问题． 考虑到除海王星以外的其它大行星到Kuiper带的距 离较远，对其中小天体的引力摄动较小，而观测事实也说 明Kuiper带的共振天体都是对于海王星的口】，因此，本文 采用限制性三体问题的简化模型，只考虑太阳和海王星 对小天体的作用，其它大行星的作用和小天体间的相互 作用忽略不计．设太阳、海王星和小天体的质量分别为 M，m7和m，三个天体均看成质点，相对位置如图1所示． 天体m相对太阳的运动方程为 ；：=一笋r—Gm7(万A+；)， (1) 其中p=G(M+m)，G为万有引力常数，r，，7和厶的意 义见图1． 采用天体的六个轨道根数——轨道半长径口、轨道偏 图1天