天文动力学和天体力学

叙述了与Astrod工程有关的天体力学和天文动力学的基本结果和学科概况 ,内容包括二体问题、摄动理论、人造地球卫星运动、限制性三体问题等     

HD 12661行星系统构形的形成与稳定性

最近的多普勒观测表明恒星HD 12661周围存在两颗中等偏心率轨道上运行的行星,内行星的最小质量为2．3木星质量,轨道周期为263．6天;外行星的最小质量为1．57木星质量,轨道周期为1444．5天．该系统的稳定性要求两颗行星处在平运动轨道共振．用数值方法研究了该系统形成初期在恒星气体盘作用下的轨道迁移与稳定性,计算了行星在迁移中被平运动共振俘获的概率．发现这两颗行星目前很可能正处在11:2平运动共振边缘,且运动是混沌的,从而澄清了关于系统目前构形的不同说法,并且很可能在系统形成后行星迁移到目前构形时,气体盘几乎消失了．     

HD 82943行星系统的2:1共振与稳定性

用直接数值积分方法通过模拟不同的行星构型探讨了HD82943行星系统(由两颗共振的巨行星组成)的长期动力演化,同时,还研究了在相空间的轨道运动.在对系统长达107年数值积分中,发现所有的稳定轨道均与2:1共振相关联.典型地,在相同的时标内,两个共振幅角θ1和θ2同时(或其中之一)存在秤动.由于共振幅角在一定范围内的秤动,因而使两颗行星轨道半长径被约束而表现为规则运动模式.另外,利用分析模型(包含了外行星偏心率e2的因素),还讨论了对于不同取值的e2和相对近星点经度θ时,内行星在相空间的运动,并发现2:1轨道共振对于相对较小的e2以及当θ=0°时易于保持.此外,适中的e2将导致系统的两颗行星进入深度共振状态.再者,分析模型和数值计算的结果吻合得很好,两者都揭示了行星系统的2:1共振结构.     

天文动力学和天体力学

叙述了与Astrod工程有关的天体力学和天文动力学的基本结果和学科概况，内容包括二体问题，摄动理论，人造地球卫星运动，限制性三体问题等。     

空间碎片轨道预报中动力学模型基准的选取

在不同的轨道预报场景中, 使用的动力学模型也不同. 例如, 在低轨空间碎片的预报中大气阻力是十分重要的摄动力, 而到了中高轨, 大气阻力就可以忽略不计. 如何为不同轨道类型的空间碎片选择最优(满足精度要求下的最简)动力学模型还没有系统、详尽的研究. 将对不同精度需求、不同轨道类型下的大批量轨道进行预报, 通过比较不同动力学模型下的预报结果, 给出各种预报场景的最优动力学模型建议. 可以为不同轨道类型的空间碎片在轨道预报时选择基准动力学模型提供参考或标准.     

小行星主带上单个行星的质量上限

介绍了N体模拟的Hermite算法,并利用该算法研究了不同质量行星在小行星主带上轨道的演化情况．采用的演化模型是太阳系N体模型(N=7),即把水星、金星、地球的质量加到太阳上,忽略冥王星,同时在小行星主带附近增加一个假想行星,系统演化时间为1亿年．数值模拟显示能够稳定存在于小行星主带上的单个天体的质量上限其量级为10~(25)kg．模拟同时还显示在某些情况下,假想行星与木星之间的低阶共振可以增强系统的稳定性．     

星系振动的动力学研究进展

概括地介绍了盘星系的运动学边缘弯曲现象、盘星系中心的HI发射线空洞和棒旋星系中核球的长轴和棒的长轴不互相平行也不互相垂直等星系振动的问接观测证据;综述了星系振动动力学的理论研究情况,主要介绍了利用非静态维里方程和N体数值模拟得到的有关结果;展望了该领域今后的发展趋势,指出了建立引力自治的星系振动模型的必要性。     

PMOE精密行星历表研究进展

本文扼要介绍了PMOE行星历表研究当前的进展和工作 ,给出了PMOE2 0 0 2与DE40 5历表的比较。     

窄吸积环动力学不稳定性的数值解

本文在细环近似和共转半么位移不为零的条件下，考虑吸积环中非轴对称动力学不稳定性的线性扰动过程，采用数值计算方法求得不稳定性的线性增长率和共转半径位移随波数的变化关系，发现线性增长率受共转半径不为零的影响较小，而共转半径位移项随波数的色散关系与线性Ｋｄｖ方程的色散关系相同，说明窄吸收环在动力学不稳定性数值模拟中出现的“行星状解”很可能是类似于Ｋｄｖ方程中的孤子解。     

太阳系动力学中当前的重要课题

以当前太阳系动力学中的重要课题以及研究方法进行讨论，并提出一些看法，课题中主要讨论动力学模型，轨道共振，行星环，混沌和长期演化，近地天体运动，Kupiper带，太阳系中的引力理论，以及其他有关问题。     

内日冕的绿线强度分布研究

日冕是太阳大气活动的关键区域,是日地空间天气的源头.受观测限制,对日冕低层大气等离子体结构和磁场状态的研究非常欠缺,国际上对于可见光波段日冕低层大气的亮度分层研究很少.利用丽江日冕仪YOGIS(Yunnan Green-line Imaging System)的日冕绿线(Fe_ⅩⅣ5303?)观测资料,对内日冕区域(1.03R_☉-1.25R_☉,R_☉表示太阳半径)亮结构及其中冕环进行了有效的强度衰减分析.对亮结构的强度在太阳径向高度上进行了指数衰减拟合,比较这些拟合结果发现所得到的静态内冕环的衰减指数在一固定值附近.然后将比较明显的冕环提取出来,通过对不同高度的绿线强度进行指数拟合,得出的衰减指数与亮结构中也比较相近,这对进一步研究日冕中的各项物理参数演化提供了参考.     

基于无量纲的轨道参数开展卫星性质统计研究

目前已发现了285颗围绕太阳系八大行星公转的卫星, 它们的轨道和物理性质呈现了丰富多样性. 目前为止, 几乎所有的卫星研究工作都基于单个卫星系统或者卫星群, 似乎缺少统一的研究. 提出了一个新的与行星性质无关、只与恒星半径有关的轨道参数n, 定义为以太阳半径为单位的轨道半长轴的自然对数. 不同行星的卫星的n值都存在双极分布, 绝大部分卫星在$n\gtrsim2$区间, 其次在$n\lesssim-1$区间, 位于中间区域的行星则很少. 从卫星物理参数和轨道参数与n的关系中发现, 分属六大行星的卫星有明显的共同特征. 首先, 轨道偏心率和轨道倾角偏大的卫星的n值都在3.5左右, 它们都是巨行星的不规则卫星. 其次, n值在-1和1之间的卫星绝大部分体积大、质量大、反照率高、自转速度慢. 从文献中找到11颗系外卫星候选体, 获得了它们轨道n值和卫星质量, 发现后者也是在-1< n< 1区间最大,其他区间偏小.这些统一的 规律暗示,太阳系内不同行星的卫星形成机制以及太阳系外卫星的形成机制可能一样或类似.     

Planetary migration to large radii

There is evidence for the existence of massive planets at orbital radii of several hundred au from their parent stars where the time-scale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.     

Analysis of Error Sources in STEP Astrometrytwo

The space telescope Search for Terrestrial Exo-Planets (STEP) employed a method of sub-pixel technology which ensures that the astrometric accuracy of the telescope on the focal plane is at the order of 1 μas. This kind of astrometric precision is promising to detect the earth-like planets beyond the solar system. In this paper, we analyze the influence of some key factors, including the errors in the stellar proper motion, parallax, the optical center of the system, and the velocity and position of the satellite, on the detection of exoplanets. We propose a relative angular distance method to evaluate the non-linear terms in the variation of star-pair's angular distance caused by the possibly existing exoplanet. This method could avoid the direct influence of measuring errors of the position and proper motion of the reference stars. Supposing that there are eight reference stars and a target star with a planet system in the same field of view, we simulate their five-year observational data, and use the least square method to get the parameters of the planet orbit. Our results show that the method is robust to detect terrestrial planets based on the 1 μas precision of STEP.     

The Rise of New Planets: Super-Earths and Sub-Neptunes

With the ceaseless progress of detecting technology, over 3500 exo-planets have been discovered. It is interesting but unexpected that the majority of the detected exoplanets are unlike any planets in our solar system. They have a size and mass between the Earth and the Neptune, and thus they are called as Super-Earths or Sub-Neptunes. In this article, I introduce these newly rising planets and review our current knowledge of their physical properties, orbital properties, and origins. Finally, I discuss the promising and exciting prospects.     

On the stability of hypothetical satellites coorbital to Mimas or Enceladus

In the present work, we study the stability of hypothetical satellites that are coorbital with Enceladus and Mimas. We performed numerical simulations of 50 particles around the triangular Lagrangian equilibrium points of Enceladus and Mimas taking into account the perturbation of Mimas, Enceladus, Tethys, Dione, Titan and the oblateness of Saturn. All particles remain on tadpole orbits after 10 000 yr of integration. Since in the past the orbit of Enceladus and Mimas expanded due to the tidal perturbation, we also simulated the system with Enceladus and Mimas at several different values of semimajor axes. The results show that in general the particles remain on tadpole orbits. The exceptions occur when Enceladus is at semimajor axes that correspond to 6:7, 5:6 and 4:5 resonances with Mimas. Therefore, if Enceladus and Mimas had satellites librating around their Lagrangian triangular points in the past, they would have been removed if Enceladus crossed one of these first-order resonances with Mimas.     

Lunar surface potential and electric field

The Moon has no significant atmosphere, thus its surface is exposed to solar ultraviolet radiation and the solar wind. Photoemission and collection of the solar wind electrons and ions may result in lunar surface charging. On the dayside, the surface potential is mainly determined by photoelectrons, modulated by the solar wind;while the nightside surface potential is a function of the plasma distribution in the lunar wake. Taking the plasma observations in the lunar environment as inputs, the global potential distribution is calculated according to the plasma sheath theory, assuming Maxwellian distributions for the surface emitted photoelectrons and the solar wind electrons. Results show that the lunar surface potential and sheath scale length change versus the solar zenith angle, which implies that the electric field has a horizontal component in addition to the vertical one. By differentiating the potential vertically and horizontally, we obtain the global electric field. It is found that the vertical electric field component is strongest at the subsolar point,which has a magnitude of 1 V m-1. The horizontal component is much weaker, and mainly appears near the terminator and on the nightside, with a magnitude of several mV m-1. The horizontal electric field component on the nightside is rotationally symmetric around the wake axis and is strongly determined by the plasma parameters in the lunar wake.     

Astrometric observations of the Uranian satellites with the Faulkes Telescope North in 2007 September

The results of astrometric observations of the main Uranian satellites taken with the Faulkes Telescope North are presented. A median filter algorithm was applied to subtract a scattered-light halo caused by Uranus. The Two-Micron All-Sky Survey (2MASS) and USNO-B1.0 were used as reference catalogues. The mean value of the differences between the equatorial coordinates of the satellites determined with 2MASS and USNO-B1.0 is close to 200 mas. A comparison of the observed equatorial coordinates of the satellites and their relative positions with ephemerides based on different combinations of theories of motion of Uranus and its satellites (DE405+GUST86, DE405+GUST06, INPOP+GUST86, INPOP+GUST06) was performed. The satellites' positions obtained with respect to 2MASS are in better agreement with theories. The values of (O−C) of the equatorial coordinates determined with the 2MASS are mainly less than 100 mas. The majority of (O−C) of relative positions are within ±50 mas. The mean values of the standard errors of (O−C) are within 20 to 60 mas.     

Dynamical behaviour of planetesimals temporarily captured by a planet from heliocentric orbits: basic formulation and the case of low random velocity

Planetesimals encountering with a planet cannot be captured permanently unless energy dissipation is taken into account, but some of them can be temporarily captured in the vicinity of the planet for an extended period of time. Such a process would be important for the origin and dynamical evolution of irregular satellites, short-period comets, and Kuiper-belt binaries. In this paper, we describe the basic formulation for the study of temporary capture of planetesimals from heliocentric orbits using three-body orbital integration, such as the definition of the duration and rate of temporary capture, and present results in the case of low random velocity of planetesimals. In the case of planetesimals initially on circular orbits, we find that planetesimals undergo a close encounter with the planet before they become temporarily captured. When planetesimals are scattered by the planet into the vicinity of one of periodic orbits around the planet, the duration of temporary capture tends to be extended. Typically, these capture orbits are in the retrograde direction around the planet. We evaluate the rate of temporary capture of planetesimals, and find that the ratio of this rate to their collision rate on to the planet increases with increasing semimajor axis of the planet. Similar results are obtained for planetesimals with non-zero but small random velocities, as long as Kepler shear dominates the relative velocity between the planet and planetesimals. For larger initial random velocities of planetesimals, temporary capture in both prograde and retrograde directions with much longer duration becomes possible.