什么是 C/2023 A3（紫金山-阿特拉斯）彗星？让我们用 Python 来绘制它的路径

彗星的基本概念

  彗星（Comet），是指进入太阳系内亮度和形状会随日距变化而变化的绕日运动的天体，呈云雾状的独特外貌，也是中国神话传说的扫帚星（星官名）。彗星分为彗核、彗发、彗尾三部分。彗核由冰物质构成，当彗星接近恒星时，彗星物质升华，在冰核周围形成朦胧的彗发和一条稀薄物质流构成的彗尾。彗星的质量、密度很小，当远离太阳时只是一个由水、氨、甲烷等冻结的冰块和夹杂许多固体尘埃粒子的“脏雪球”。当接近太阳时，彗星在太阳辐射作用下分解成彗头和彗尾，状如扫帚。彗星以椭圆形、抛物线和双曲线这三种轨道围绕太阳运行。

C/2023 A3（紫金山-阿特拉斯）彗星

  2023年3月1日，国际小行星中心发布一颗新彗星C/2023 A3（Tsuchinshan-ATLAS）的确认公告。该目标由中国科学院紫金山天文台近地天体望远镜于1月9日首次发现，后由南非ATLAS计划在2月22日报告存在彗星特征，进一步通过美国帕洛玛ZTF望远镜观测资料的回溯检测，确认为一颗已经开始活动的彗星。彗星的临时编号为A10SVYR，根据彗星命名系统，这颗彗星获得了两个天文台的名称，并被正式命名为C/2023 A3（Tsuchinshan-ATLAS）。

观测

• 10月14日至20日：0.5到1.7等，亮度减弱，但距离太阳更远，因此更容易被看到。日落后约一小时，肉眼可见。您可以看到反彗尾，即一条明亮的条纹，似乎指向太阳，与其他彗尾相反。10月15日，彗星将经过M5球状星团1.4°，提供良好的拍照机会。彗星穿过巨蛇座。
• 10月20日至31日：亮度从1.9到3.9等，穿过蛇夫座，位于西南天空相对较高的位置，晚上可见。用双筒望远镜和望远镜可见。在这几天里，彗尾将开始快速增长。
• 11月：亮度从4到8等，晚上可见。日落后在北半球升得更高。
• 12月：亮度从8到10等，在天空中逐渐靠近太阳，升到地平线以下。在南半球不可见。
  

绘制

import numpy as np
from matplotlib import pyplot as plt
from matplotlib.collections import LineCollection

from skyfield.api import Star, load
from skyfield.constants import GM_SUN_Pitjeva_2005_km3_s2 as GM_SUN
from skyfield.data import hipparcos, mpc, stellarium
from skyfield.projections import build_stereographic_projection

# 加载天文历
ts = load.timescale()
t_comet = ts.utc(2024, 10, range(1, 32))
t = t_comet[len(t_comet) // 2]  # middle date

# 加载太阳系天体数据（行星及以上）
eph = load('de421.bsp')
sun = eph['sun']
earth = eph['earth']

# 加载小行星中心数据，创建彗星对象
with load.open(mpc.COMET_URL) as f:
    comets = mpc.load_comets_dataframe(f).set_index("designation", drop=False)
n = 'C/2023 A3 (Tsuchinshan-ATLAS)'
row = comets.loc[n]  # 提取 Tsuchinshan-ATLAS 彗星参数
comet = sun + mpc.comet_orbit(row, ts, GM_SUN)

# 加载恒星目录（Hipparcos）
with load.open(hipparcos.URL) as f: