关键词:
双星
轨道拟合
数据处理
恒星基本参数
数值积分
摘要:
双(多)星系统是普遍存在的,并且双星的轨道拟合能够导出独立于模型的恒星参数,为探索恒星理论演化模型提供了新的基准恒星样本。近些年得益于地面和太空望远镜的发展,高分辨率光谱和高精度测光观测的实施,高精度的观测数据正在不断积累。与此同时,双星系统的证认及其轨道拟合依然是天文研究中的一个焦点。
在本论文中,首先我们利用郭守敬望远镜低分辨率光谱数据(LAMOST-LRS)对单谱线分光双星(SB1)通过轨道拟合的方式进行证认。这里先使用χ2检验对LAMOST-LRS的一个子样本进行筛选,得到视向速度变星。随后通过与变星星表的交叉匹配,丢弃其中的本质变星,得到SB1候选源。为了可靠地证认出SB1s,我们结合 Thejoker、rvfit 和 Levenberg-Marquardt 算法来给出双星轨道的最优解,如果其满足可接受轨道解的判别标准,则被证认为SB1s。最终获得255个SB1s,其中168个是新发现的。交叉验证的结果表明,我们给出的轨道周期与其他星表给出的周期在1σ不确定度内相一致。相应周期-偏心率分布图表明了大多数短周期双星偏心率较小,这与轨道圆化效应预期的结果一致。此外,与母样本相比,SB1s星表表现出相对较高的矮星比例和略低的金属丰度,这一统计特征与其他研究者得出的结论一致。进一步地,我们从上述SB1s星表中发现了一个有趣的系统,该系统为双谱线不接食双星(EPIC 211391083)。其同时拥有APOGEE多次观测的高分辨率光谱数据和K2两次观测的光变曲线数据。我们利用BF方法从连续归一化光谱中提取子星的视向速度数据,由rvfit程序求解了系统的分光轨道。然后我们同时对光变曲线和视向速度曲线进行建模,联合拟合这两种数据导出的主星和伴星的质量分别为M1=1.46647 M⊙和M2=1.37578 M⊙;半径分别为R1=1.61768 R⊙ 和R2=1.46257 R⊙。我们也利用已有的信息对子星的有效温度和表面重力进行了估计,并且将PARSEC恒星演化模型计算的数据与我们的结果进行比较,分析了子星的演化阶段和年龄。最后考虑到大约10%的双星系统本身是三星或四星系统的成员。研究这些多星系统的长期动力学稳定性,需要用到效果良好的数值积分方法。这里我们基于Nacozy的流形改正算法发展了一种新的流形改正算法,其可应用于多星系统。在一个纯二体问题中,存在七个积分守恒量。我们使用七个比例因子来改正积分后的位置和速度,以确保调整后的解严格满足七个守恒量。我们采用了结合奇异值分解的牛顿迭代法来计算这些因子。这种改正算法可以应用于受摄二体和多体问题中。这时与每个系统相关的七个量会随时间缓慢变化。通过同时积分这些量的积分不变关系和运动方程,可以得到这七个缓变量某个时刻的值。数值测试表明,新方法可以显著减小所有轨道参数数值误差的增长速率。