小行星的搜寻和定轨

本文简要介绍了与小行星有关的一些基本知识和在小行星搜寻方面的国际进展情况，侧重与ASTROD项目有关的内容。并介绍了与小行星定轨有关的网络资源。     

小行星YORP效应的可探测候选体

YORP (Yarkovsky-O''Keefe-Radzievskii-Paddack)效应是小行星长期动力学演化的机制之一. 与碰撞、引力摄动等因素相比, YORP效应作用量级小, 短时标观测效应不明显, 这给直接测量YORP效应带来了很大的困难. 利用小行星光变数据库中已知的小行星数据, 统计了小行星的自转速率分布, 使用核密度估计以及Kolmogorov-Smirnov检验分别分析了近地小行星和主带小行星自转速率的分布特性, 分别给出了在近地小行星和主带小行星中寻找受YORP效应影响减速自转的最佳样本群; 基于7颗已被探测到YORP旋转加速度的近地小行星, 利用YORP强度估计方法和光变探测条件建立了筛选模型, 给出了未来可直接通过光变数据探测\lk YORP效应的10颗近地小行星.     

近地小行星发现数目与发现场景的统计分析

近地小行星是一类可能对地球安全造成潜在威胁的太阳系小天体, 目前绝大部分的近地小行星是由地基望远镜发现的, 且数目仍在不断增加. 为了对我国未来开展近地小行星发现监测提供参考和借鉴, 利用国际小行星中心公开的数据库对所有近地小行星首次发现时刻的观测资料开展了多维度统计分析. 发现望远镜探测能力的限制会对近地小行星的发现造成选择效应, 导致不同轨道类型近地小行星发现的相对比例逐年变化且与直径有关. 另外, 结合数值模拟获得的轨道数据, 对近地小行星首次发现时的观测场景进行了还原, 获得了发现时刻近地小行星位置在不同天球坐标系的分布, 分析了其分布特征与季节、测站纬度和小行星直径的依赖关系. 最后, 通过分析数据定量考察了太阳、月球和银道面对近地小行星发现的影响, 发现地基望远镜一般难以发现来自太阳方向90$^\circ$范围内直径140m以下的近地小行星, 并且随着小行星直径的减小该限制范围也将变大; 月光污染对近地小行星发现的影响也非常显著, 望月前后几天的观测限制可导致约29%的目标无法被发现, 而且分析表明农历上半月发现的目标一般比下半月发现的更难以被跟踪观测; 银道面特别是银心方向会对近地小行星发现产生影响, 使得黄道面附近存在与季节相关的观测``盲区''.     

小行星热物理及Yarkovsky效应和YORP效应的研究进展

小行星热物理的研究是近年来小行星研究中一个前沿方向。简单阐述了在小行星热物理这一研究领域中常用的热物理模型,利用这些热物理模型并结合空间或地面红外望远镜的观测数据,可获得小行星的一些热物理参数,如热惯量、几何反照率、有效直径、表面粗糙度、表面温度等。介绍了小行星的形状模型及红外观测,还介绍了利用这些模型得到的部分小行星热物理参数。通过这些热物理参数,可进一步研究小行星的Yarkovsky效应和YORP效应等,甚至可为人造航天器着陆小行星表面及小行星采样返回任务提供相关科学信息。     

Flora小行星族自转特性研究

小行星族作为灾变碰撞的残留物,其基础物理性质提供了其母体以及后续演化信息.其中轨道以及自转特性分别反映了Yarkovsky效应以及Yarkovsky-O’Keefe-Radzievskii-Paddack效应(YORP效应)对于小行星族演化的影响.基于小行星光变数据库(Asteroid Lightcurve Database),通过对Flora小行星族自转速率分布进行研究,发现随着直径减小,族成员自转速率倾向于主要集中在3–5 d~(-1)的范围内.同时,可以注意到Flora小行星族整体表现出更倾向于顺行自转状态的现象,但对于轨道半长轴小于2.2au的成员来说,其顺行自转与逆行自转状态成员数目比接近于近地小行星中顺逆行自转状态源1:3的比例;此外,对于轨道半长轴大于2.2 au且具有顺行自转状态的部分族成员,在轨道半长轴-绝对星等分布中表现出聚集现象,并在聚集区域中有9颗成员展现出类似Slivan状态特征.     

飞行器近小行星轨道动力学的特点及研究意义

该文在介绍小行星及对其探测意义的基础上,以NEAR、Muses-C(Hayabusa)小行星探测器为例,重点分析了飞行器与小行星距离较近时(near asteroid)的动力学环境及其对环绕小行星飞行的飞行器设计的影响.并以旋转的二阶二次重力场和小行星Castalia为例,进一步阐述近小行星轨道动力学的一些特点和结论.此外,还给出了小行星探测的一些建议和设想.     

小行星动态(2001.08～2001.09)

月份 获永久编号小行星数 小行星命名数2001.08 863(10 0) 123 (5 4)2001.09 1420(28 0) 52(1 0) (括号内数字分别指国家天文台施密特CCD小行星项目组和紫金山天文台的发现数) 在8月4日的小行星通报上,小行星中心改变了从年初开始的每两期公布一次新的小行星命名的惯例,而在这一期继7月份的小行星命名之后,继续给出了新的小行星命名。这使得9月2日的小行星通报上     

小行星与彗星—近地小行星,主带小行星和彗星之间相互关系探讨

小行星和彗星都是太阳系中的小天体,而近地小行星又是小行星中特殊的一类。近年来对近地小行星的观测和研究已全面开展,这已成为当今太阳系研究的前沿领域之一。本文初步探讨太阳系中不同小天体之间的关系:1.小行星与短周期彗星之间的关系;2.近地小行星与主带小行星和短周期彗星之间的关系。最近的观测研究表明某些小行星可能是由彗星演变来的,而某些近地小行星可能来源于主带小行星和死亡的彗核。     

在电子计算机上进行的小行星特别摄动计算——方法与精度问题

一、引言紫金山天文台于1954年开始从事小行星摄动的研究,八年多来进行了近10颗小行星的特别摄动计算和轨道改进.鉴于国内计算技术的迅速进展和小行星研究的迫切需要,1962年秋季,在经过较多的准备以后,我们利用中国科学院计算技术研究所的快速电子计算机,初步成功地建立了小行星特别摄动计算的工作程序.我们考虑了在太阳-内行星的质心综合引力下,小行星所受木星、土星的摄动,按照科     

用小行星观测检测星表的系统误差

本文提出了用小行星观测检测星表系统误差的方法。叙述了选取小行星的方法及小行星轨道真值的求解方法,介绍了怎样用小行星观测来体现和确定星表的系统误差。     

Physical properties of meteorites—Applications in space missions to asteroids

Abstract— Based on reflectance spectroscopy and chemical/mineralogical remote sensing methods, it is generally assumed that asteroids are parent bodies for most meteorites reaching the Earth. However, more detailed observations indicate that differences exist in composition between asteroids and meteorites resulting in difficulties when searching for meteorite‐asteroid match. We show that among other physical parameters the magnetic susceptibility of an asteroid can be determined remotely from the magnetic induction by solar wind using an orbiting spacecraft or directly using the AC coil on the lander, or it can be measured in samples returned to the laboratory. The shape corrected value of the true magnetic susceptibility of an asteroid can be compared to those of meteorites in the existing database, allowing closer match between asteroids and meteorites. The database of physical properties contains over 700 samples and was recently enlarged with measurements of meteorites in European museums using mobile laboratory facility.     

小行星(469219) Kamo`oalewa轨道的确定与误差分析

作为一颗与地球共轨道的小行星,(469219)Kamo'oalewa是一个具有很高研究价值的近地小天体,也是中国首次小行星探测计划的目标天体之一.针对其轨道特性,建立了兼顾太阳、地球和月球非球形引力作用的小行星动力学模型.并在该模型的基础上,利用国际小行星中心(Minor Planet Center,MPC)提供的2004|2018年间的光学观测数据对该小行星的轨道进行确定.拟合后观测残差的均方根误差约为0:2″(与美国喷气推进实验室的Horizons在线历表系统相当),其中2004年期间数据的观测残差有所改进.最后,对小行星(469219)Kamo'oalewa的轨道误差进行了详细分析,并预报了2020-2025年期间该小行星的轨道误差.     

Orbit Determination of the Asteroid (469219) Kamoòalewa and Its Error Analysis

As an Earth co-orbital asteroid, (469219) Kamoòalewa is a near earth object (NEO) with high value of research, and one of the targets explored by the first Chinese asteroid exploration mission. Given its orbit characteristics, we build a refined dynamical model for this asteroid, in which the effects induced by nonspherical gravitational fields of the Sun, the Earth, and the Moon are combined. On the basis of the dynamical model of the asteroid (469219) Kamoòalewa, its orbit is determined with optical data from 2004 to 2018 available on the Minor Planet Center (MPC) database. The root mean square error of post-fit residuals is about 0.2 arc second (comparable with that of the Jet Propulsion Laboratory (JPL)/Horizons), and the post-fit residuals of optical observations in 2004 are decreased. At the end, we implement error analysis on the asteroid (469219) Kamoòalewa's orbit in detail, and also predict its orbit error at the time interval between 2020 and 2025.     

Asteroid spin‐axis longitudes from the Lowell Observatory database

By analyzing brightness variation with ecliptic longitude and using the Lowell Observatory photometric database, we estimate spin‐axis longitudes for more than 350,000 asteroids. Hitherto, spin‐axis longitude estimates have been made for fewer than 200 asteroids. We investigate longitude distributions in different dynamical groups and asteroid families. We show that asteroid spin‐axis longitudes are not isotropically distributed as previously considered. We find that the spin‐axis longitude distribution for Main Belt asteroids is clearly nonrandom, with an excess of longitudes from the interval 30°–110° and a paucity between 120° and 180°. The explanation of the nonisotropic distribution is unknown at this point. Further studies have to be conducted to determine if the shape of the distribution can be explained by observational bias, selection effects, a real physical process, or other mechanism.     

Massive identification of asteroids in three-body resonances

An essential role in the asteroidal dynamics is played by the mean motion resonances. Two-body planet–asteroid resonances are widely known, due to the Kirkwood gaps. Besides, so-called three-body mean motion resonances exist, in which an asteroid and two planets participate. Identification of asteroids in three-body (namely, Jupiter–Saturn–asteroid) resonances was initially accomplished by Nesvorný and Morbidelli (Nesvorný D., Morbidelli, A. [1998]. Astron. J. 116, 3029–3037), who, by means of visual analysis of the time behaviour of resonant arguments, found 255 asteroids to reside in such resonances. We develop specialized algorithms and software for massive automatic identification of asteroids in the three-body, as well as two-body, resonances of arbitrary order, by means of automatic analysis of the time behaviour of resonant arguments. In the computation of orbits, all essential perturbations are taken into account. We integrate the asteroidal orbits on the time interval of 100,000 yr and identify main-belt asteroids in the three-body Jupiter–Saturn–asteroid resonances up to the 6th order inclusive, and in the two-body Jupiter–asteroid resonances up to the 9th order inclusive, in the set of ～250,000 objects from the “Asteroids – Dynamic Site” (AstDyS) database. The percentages of resonant objects, including extrapolations for higher-order resonances, are determined. In particular, the observed fraction of pure-resonant asteroids (those exhibiting resonant libration on the whole interval of integration) in the three-body resonances up to the 6th order inclusive is ≈0.9% of the whole set; and, using a higher-order extrapolation, the actual total fraction of pure-resonant asteroids in the three-body resonances of all orders is estimated as ≈1.1% of the whole set.     

Observations of asteroid occultations by the trailed-image method

The method of trailed CCD images for observations of asteroid occultations is described. This method was used to observe 9 asteroid occupations at the Lisnyky observational station in 2006 with the telescope AZT-8 (D= 0.7m and F= 2.8m) equipped by the CCD ST-8 XME. In the case of occultation of the star TYC 0587-00209-1U by the asteroid 76 Freia in November 4, 2006, the distance between the asteroid center and the star, as well as the time of asteroid occultation were determined. The size of asteroid 76 Freia is determined assuming that its shape is spherically symmetric. In other cases, the minimal distances between the asteroid center and the star are determined. The method makes it possible to observe asteroid occultations with high time resolution.     

Reliability of estimates for the sizes of spots on asteroid surfaces made by the spectral-frequency method

A new spectral-frequency method (SFM) for the study of solid body surfaces is briefly described. This method allows estimation of the sizes of various spots. Estimates for the sizes of spots on asteroid surfaces made by the SFM and other methods are compared and discussed. The sizes of spots on the surface of asteroid 1620 Geographos determined by the SFM are well consistent with those of the craters obtained from radar data. The sizes of hydrosilicate spots on the surface of asteroid 21 Lutetia found by the SFM agree with those of the craters determined by the Rosetta spacecraft. The size of a blue spot on the surface of asteroid 4 Vesta found by the SFM is consistent with the size of the well-known crater on the south pole of the asteroid. It is inferred that the SFM is a promising method for the estimation of the sizes of spots on asteroid surfaces.     

Catastrophic disruption of asteroids and family formation: a review of numerical simulations including both fragmentation and gravitational reaccumulations

In the last few years, thanks to the development of sophisticated numerical codes, a major breakthrough has been achieved in our understanding of the processes involved in small body collisions. In this review, we summarize the most recent results provided by numerical simulations, accounting for both the fragmentation of an asteroid and the gravitational interactions of the generated fragments. These studies have greatly improved our knowledge of the mechanisms that are at the origin of some observed features in the asteroid belt. In particular, the simulations have demonstrated that, for bodies larger than several kilometers, the collisional process not only involves the fragmentation of the asteroid but also the gravitational interactions between the ejected fragments. This latter mechanism can lead to the formation of large aggregates by gravitational reaccumulation of smaller fragments, and helps explain the presence of large members within asteroid families. Numerical simulations of the complete process have thus reproduced successfully for the first time the main properties of asteroid families, each formed by the disruption of a large parent body, and provided information on the possible internal structure of the parent bodies. A large amount of work remains necessary, however, to understand in deeper detail the physical process as a function of material properties and internal structures that are relevant to asteroids, and to determine in a more quantitative way the outcome properties such as fragment shapes and rotational states.     

Stability of Binary Asteroids

The stability and final outcome of a strongly interacting binary asteroid system is considered. We discuss the implications of the system transferring energy and angular momentum between rotational and translational motion while conserving the total system energy and angular momentum. Using these results we can develop a set of sufficient conditions for stability against escape and impact. These allow us to delineate several classes of final outcomes for a binary asteroid system, each of which may have implications for asteroid observations. The effects of energy dissipation on an asteroid binary system are also considered and are shown to be able to change the stability of the system against escape and impact. An example computation for the near-Earth asteroid binary 1996 FG3 is given along with a series of numerical explorations of an evolving binary system consisting of an ellipsoid and a sphere of equal mass.     

Effect of the finite size of an asteroid on its deflection using a tether–ballast system

Potentially hazardous near-Earth objects which can impose a significant threat on life on the planet have generated a lot of interest in the study of various asteroid deflection strategies. There are numerous asteroid deflection techniques suggested and discussed in the literature. This paper is focused on one of the non-destructive asteroid deflection strategies by attaching a long tether–ballast system to the asteroid. In the existing literature on this technique, very simplified models of the asteroid-tether–ballast system including a point mass model of the asteroid have been used. In this paper, the dynamical effect of using a finite size asteroid model on the asteroid deflection achieved is analyzed in detail. It has been shown that considering the finite size of the asteroid, instead of the point mass approximation, can have significant influence on the deflection predicted. Furthermore the effect of the tether-deployment stage, which is an essential part of any realistic asteroid deflection mission, on the predicted deflection is studied in this paper. Finally the effect of cutting the tether on the deflection achieved is analyzed and it has been shown that depending on the orbital properties of the asteroid as well as its size and rotational rate, cutting the tether at an appropriate time can increase the deflection achieved. Several numerical examples have been used in this paper to elaborate on the proposed technique and to quantitatively analyze the effect of different parameters on the asteroid deflection.