长周期彗星动力学特性统计研究

基于JPL (Jet Propulsion Laboratory)和MPC (Minor Planet Center)的彗星观测数据, 对长周期彗星(Long Period Comet, LPC)的轨道进行逆向演化, 得到其原始轨道, 并分析了LPC的观测轨道和原始轨道的动力学特性. 研究表明: LPC观测轨道半长轴倒数(1/a)_obs和原始轨道半长轴倒数(1/a)_ori分布存在显著差异, 峰值分别位于(1/a)_obs～ 2 × 10^-5au^-1和(1/a)_ori～ 6 × 10^-5au^-1; 大近日点LPC(近日点距离q > 3.1 au)和小近日点LPC (近日点距离q 3.1 au)的原始轨道(1/a)_ori峰值位置相同, 均位于1 × 10^-4au^-1; JPL和MPC数据中LPC的原始轨道(1/a)_obs峰...     

小天体离散元仿真软件的搭建与测试

小天体探测是当今太阳系探测的一大热点, 对小行星演化的研究有助于人们了解太阳系的起源. 演化研究的一项重要内容是小天体结构的演化, 即小天体在多种力学机制作用下自身形状与结构的演化. 在小天体碎石堆结构的假设之下, 一种比较常见的模拟小天体结构演化的方法是离散元仿真算法(Discrete Element Method, DEM), 目前国内外有数个团队开发了相关算法软件. 本文介绍了团队开发的《基于DEM仿真算法的多粒子系统模拟软件》的基础理论、实现方法以及加速算法, 并使用二体接触模型、声速传播仿真、小行星内部压强、小行星旋转稳定性仿真算例验证了算法的可靠性.     

金属丰度对球状星团动力学演化的影响

动力学过程和恒星演化及二者的互相影响都会对球状星团的演化产生重要影响.由于金属丰度会影响恒星的演化轨迹,与之相伴随的恒星质量损失率的变化也会对球状星团的动力学过程造成影响.通过一系列N体模拟研究金属丰度对球状星团的质量损失率、半径等的影响,并分析其原因,同时研究了大质量恒星以及星团初始数密度分布的影响.模拟中采用的球状星团模型初始成员星数目N=50000,运行于类银河系的引力势中并考虑成员星的演化.结果显示,由于低金属丰度恒星拥有较快的演化时标,所以贫金属球状星团在早期会拥有较高的质量损失,但与此同时它们的核塌缩时间会比后者显著推迟,因此在核塌缩之后其质量损失会被富金属星团反超.另外由于大质量恒星演化导致的质量损失较大,所以大质量星的存在会使金属丰度更加显著地影响球状星团早期的扩张以及随后的核塌缩过程,同时星团的初始数密度分布也对该效应有着不可忽视的影响.     

彗星动力学研究综述

彗星是太阳系中天文学的重要研究对象．它涉及到天体物理、化学、天体力学等多种领域。近年的研究表明,彗星可能是揭开太阳系起源和演化之谜的突破口,还可能与生命的起源有关．本文旨在对彗星动力学的发展情况作一简单的回顾和评述．     

鹿林彗星的断尾事件研究

鹿林彗星在2009年2月4日(UT)发生了一次断尾事件,通过对2月4日前后的观测数据进行分析,发现在2月3日鹿林彗星就已发生过一次断尾事件.对2月3日发生的断尾事件进行研究可得出断尾移动的速度为68 km·s-1,断尾事件发生的时间为2009年2月3日(15.24土5.16)时.结合断尾事件发生前后鹿林彗星的轨道特征和STEREO-A飞船所测量到的太阳风数据,鹿林彗星2月3日发生的断尾事件可能是由于日冕物质抛射和彗星相互作用导致的磁重联引起的.     

近地天体望远镜的PSF分布研究

盱眙近地天体望远镜CCD视场为1.94°×1.94°,作为大视场望远镜由于存在场曲现象,因此根据单颗星像调焦会导致PSF(Point Spread Function)分布不均匀.研究所拍摄图像的PSF分布特征,并通过拟合半高全宽和离焦量的函数关系模拟出弯曲焦面,从而确定出最佳成像平面.制定出根据PSF的半高全宽直方图来判断最佳成像平面的调焦方法.     

修正的极大似然法拟合幂律频数分布

针对包含饱和样本数据的频数幂律分布拟合,提出一个新的幂律分布指数的极大似然估计方法的修正公式.对比研究显示,修正公式适用于剔除异常饱和值的幂律频数拟合.如果不剔除饱和值,幂律指数的估计只能使用修正前的公式,其误差随幂律指数变化,指数较小时误差较大.由此建议,对于包含饱和样本的频数分布拟合,首先剔除异常的饱和值,然后对剩余不含饱和值的子集使用修正公式进行参数估计.     

一种低表面亮度星系的自动搜索算法newline--YOLOX-CS

低表面亮度星系(Low Surface Brightness Galaxy, LSBG)的特征对于理解星系整体特征非常重要, 通过现代的机器学习特别是深度学习算法来搜寻扩充低表面亮度星系样本具有重要意义. LSBG因特征不明显而难以用传统方法进行自动和准确辨别, 但深度学习确具有自动找出复杂且有效特征的优势, 针对此问题提出了一种可用于在大样本巡天观测项目中搜寻LSBG的算法---YOLOX-CS (You Only Look Once version X-CS). 首先通过实验对比5种经典目标检测算法并选择较优的YOLOX算法作为基础算法, 然后结合不同注意力机制和不同优化器, 构建了YOLOX-CS的框架结构. 数据集使用的是斯隆数字化巡天(Sloan Digital Sky Survey, SDSS)中的图像, 其标签来自于$\alpha.40$-SDSS DR7 (40%中性氢苜蓿巡天与第7次数据发布的斯隆数字化巡天的交叉覆盖天区)巡天项目中的LSBG, 由于该数据集样本较少, 还采用了深度卷积生成对抗网络(Deep Convolutional Generative Adversarial Networks, DCGAN)模型扩充了实验测试数据. 通过与一系列目标检测算法对比后, YOLOX-CS在扩充前后两个数据集中搜索LSBG的召回率和AP (Average Precision)值都有较好的测试结果, 其在未扩充数据集的测试集中的召回率达到97.75%, AP值达到97.83%, 在DCGAN模型扩充的数据集中, 同样测试集下进行实验的召回率达到99.10%, AP值达到98.94%, 验证了该算法在LSBG搜索中具有优秀的性能. 最后, 将该算法应用到SDSS部分测光数据上, 搜寻得到了765个LSBG候选体.     

热海王星系统HD 106315的近$2:1$ 平运动共振捕获与轨道演化

通过结合理论分析和数值模拟方法,可以对热海王星系统HD 106315轨道迁移中的近2:1平运动共振捕获机制以及潮汐作用下的演化过程进行研究.在轨道迁移阶段,初始轨道半长径、初始偏心率以及行星c的偏心率衰减系数K会对系统轨道构型产生影响.数值模拟结果显示当初始轨道半长径分别为ab～0.4 au、ac～0.8 au,偏心率eb和ec均小于0.03时, HD 106315b和HD 106315c在中央恒星的引力作用以及原行星盘粘滞作用下向内迁移, 65000 yr左右两颗行星均可迁移至当前观测位置附近并形成近2:1平运动共振捕获.此外,中央恒星的潮汐效应也可能会对行星系统共振构型产生影响,理论分析表明当行星潮汐耗散系数Q=100时,潮汐效应造成的轨道半长径衰减使系统轨道周期比发生的变化可能是系统脱离共振构型的原因.数值模拟结果显示, HD 106315系统内两颗行星Q103时,来自中央恒星的潮汐效应并不会使行星系统产生明显的偏心率和轨道半长径衰减,不足以使HD 106315行星系统在剩余寿命内脱离2:1平运动共振轨道构型.     

On the origin of HNC in Comet Lee

We have modelled the chemistry occurring in the coma of Comet Lee and have critically evaluated the possible routes leading to HNC. We show that the observed levels of HNC cannot be produced by ion–molecule chemistry, or by reactions of energetic H atoms with HCN. Rather, it appears that HNC is injected into the coma following the photodestruction of an unknown precursor. We discuss the possible nature of the parent of HNC and conclude that photofragmentation of large HCN polymers, such as polyaminocyanomethylene (PACM), is responsible. The degradation of hydrogen cyanide polymers may constitute a common source of HNC in comets, accounting for HNC/HCN ratios in the range measured in Lee and Hyakutake (≈ 0.06–0.12). The high HNC/HCN ratio measured in Hale–Bopp (≈0.2) and its heliocentric variation may, however, require an additional source.     

Organic synthesis in the coma of Comet Hale–Bopp?

Several organic molecules have now been detected in the coma of Hale–Bopp. These species may either emanate from the nucleus, or, as has been suggested by Bockelée–Morvan et al., could be synthesized in the coma. We have modelled the gas phase chemistry which occurred in the coma of Hale–Bopp, concentrating on the observed organic molecules HCOOH, HCOOCH₃, HC₃N and CH₃CN. We find that gas phase chemical reactions are unable to synthesize the observed abundances of these molecules, so all these species are most probably present in the nuclear ice. We briefly discuss the implications of this result for the connection between cometary and interstellar ices.     

On the periodicity hypothesis of the ages of large impact craters

An analysis is made of the periodicity hypothesis of the ages of large craters, based on the compilation by Grieve with the addition of recently identified craters. A method earlier proposed by Broadbent is used to derive a period, and the significance of the derived period is tested by a Monte Carlo experiment. In accordance with the result of Stothers, the ages of large craters  ( D >30 km)  are shown to exhibit a period close to 37.5 Myr. Monte Carlo experiments show, however, that the derived period is far from being statistically significant. A subset of crater data earlier adopted by Napier for the purpose of similar investigation is also tested, and it is shown that they also exhibit a similar period at an almost identical level of confidence. A brief discussion is made of the relation between the derived period and that associated with faunal mass extinctions.     

Orbit of Comet C/1857 D1 (d'Arrest)

Comet C/ 1857 D1 (d'Arrest) is one of a large number of comets with parabolic orbits. Given that there are sufficient observations of the comet, 299 in right ascension and 279 in declination, it proves possible to calculate a better orbit. The calculations are based on a 12th order predictor‐corrector method. The comet's orbit is highly elliptical, e = 0.99982 and, from calculated mean errors, statistically different from a parabola. The comet will not return for at least 44000 years and thus represents no immediate NEO threat (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

New measures of the orientation of the linear structure in the sunward hemisphere of the coma of Comet 19P/Borrelly: apparitions of 1994, 2001, and 2008

The sky‐projected orientation (position angle) of the axis (line of maximum density or maximum brightness) of the long time‐known, linear structure (LS) in the sunward hemisphere of the coma of Comet 19P/Borrelly is measured on 45 photographs taken by different observers under different projection conditions and covering three consecutive apparitions (1994, 2001, and 2008) for a total time interval of 5174 days. The analysis of the results by a tomographic approach yields an LS axis constantly oriented towards a fixed point in the space, at Right Ascension 214°.4 ± 0°.5 and Declination –7°.0 ± 0°.5 (J2000), corresponding to an obliquity of 103°.5 ± 1° and an orbital longitude of 147°.2 ± 1°, throughout the relevant interval. Such coordinates are close to the ones found by other authors for the spatial orientation of the nucleus spin axis during the apparitions of 1994 and 2001. In the hypothesis of an LS orientation aligned with the nucleus spin axis, the new results confirm the previous ones and show that this orientation remained unchanged during the subsequent 2008 apparition (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)