Barycentric Corrections for HST/STIS Data

技术摘要：仪器科学报告 STIS 2026-02

问题陈述
对于需要毫秒级精度的 HST/STIS 观测（如脉冲星计时），准确的质心时间修正至关重要。这些修正考虑了观测者（HST）与太阳系质心（SSB）之间的光传播时间差，即罗默延迟（Rømer delay），以及相对论效应（爱因斯坦延迟和夏皮罗延迟）。历史上，这些修正使用 IRAF 任务 odelaytime 进行计算。然而，随着 IRAF 的弃用，HST 社区缺乏专门用于此目的的工具。此外，遗留的 odelaytime 依赖于过时的历表数据（DE200，创建于 1981 年），并且需要通过已弃用的 FTP 协议手动下载 HST 轨道文件。过时的轨道数据、潜在的时间标准歧义（UTC 与 TDB）以及缺乏现代软件实用程序的结合，引入了数毫秒的计时偏差，这在光传播速度约为每毫秒 300 公里的情况下是非常显著的。

方法论
作者引入了 stistools.barycentric_correction，这是一个旨在取代 odelaytime 的新型 Python 实用程序。该方法涉及三个主要组成部分：

1. 位置确定：

   • 地球： 该工具不再使用静态 FITS 文件，而是通过 astroquery.jplhorizons 查询 JPL Horizons 系统，或利用内置的 astropy 历表（特别是 DE440）。这确保了相对于 SSB 的地球位置是最新的。作者强调了在查询过程中使用正确时间标准（TDB vs. UTC）的极端重要性，指出此处错误会导致数千公里的位置偏移。
   • HST： 该工具通过三种方法计算相对于地球地心的 HST 位置，以进行交叉验证：(1) 查询 JPL Horizons；(2) 从 HST 轨道历表（GCRS 坐标）进行插值；(3) 使用 SPICE 内核（via SpiceyPy）。该工具既能处理遗留的轨道文件，也能处理现代 JPL 查询，并针对每一个 TIME-TAG 事件（而非仅针对曝光开始/结束时间）更新位置。

2. 时间标准转换：
   该实用程序严格管理时间标准。它将输入时间（通常为 MJDUTC）转换为以 Barycentric Dynamical Time 标准表示的重心儒略日（BJDTDB）。这包括：

   • 使用观测者位置矢量与目标单位矢量的点积来计算几何罗默延迟 ($\Delta R_{\odot}$)。
   • 应用“有限距离”修正，以补偿太阳系内目标由于波前曲率产生的项，这一项在标准的 astropy 光传播时间计算中被忽略，但在原始的 odelaytime 中是存在的。
   • 应用时钟修正 ($\Delta C$)，以解释闰秒以及 TAI、TT 与 TDB 之间的偏移。
   • 在相关情况下包含夏皮罗（Shapiro）和爱因斯坦（Einstein）延迟。

3. 实现：
   该工具接受 FITS 请求头（tag, raw, flt, x1d, x2d）并输出修正后的时间。它通过以 1 分钟分辨率进行批处理并添加填充，来管理 JPL Horizons 的查询速率，以避免服务器过载。

核心贡献

• 现代化： 用基于 Python 的实用程序（stistools.barycentric_correction）取代了已弃用的 IRAF odelaytime，使其与当前的 HST 数据归约工作流兼容。
• 历表精度： 从 1981 年的 DE200 历表迁移到通过 JPL Horizons 和 Astropy 实现的现代 DE440 标准，将地球位置的不确定性从数千公里降低到几公里。
• 自动化 HST 定位： 通过集成直接查询 JPL Horizons 和 SPICE 内核，消除了手动下载 HST 轨道文件的需求，同时保持了对现有轨道文件的向后兼容性。
• 精度处理： 显式处理 UTC 与 TDB 之间的转换，管理闰秒，并包含针对太阳系目标的有限距离修正，旨在实现毫秒级的精度。

结果与验证
作者通过以下几次对比验证了该工具的准确性：

• 与遗留代码的比较： 当使用旧的 DE200 历表时，新工具与原 odelaytime 移植版本显示出 >3 ms 的差异。然而，当更新至 DE440 时，新工具与更新后的 odelaytime 移植版本在 1 ms 以内保持一致。
• 交叉验证： 通过计算 HST 位置（JPL Horizons 方法 vs. HST 轨道文件方法）进行的比较显示，两者在 1 ms 以内达成一致。
• 第三方检查： 该工具的结果在 BJDTDB 标准下与 barycorrpy.utc_tdb (Kanodia et al. 2018) 的结果相比，差异小于 1 ms。
• 绝对精度： 文中指出，虽然修正算法实现了毫秒级精度，但 STIS 时间戳的绝对精度受限于机载硬件。在计数率低于 32,000 ct/s 时，STIS 时间戳的精度约为 10.125 ms（结合了 $\le$10 ms 的 HST 时钟漂移和 125 $\mu$s 的精细定时器分辨率）。在高于此速率时，队列延迟可能会引入高达 82 ms 的额外不确定性。

意义
论文声称 stistools.barycentric_correction 恢复并现代化了对 HST/STIS 数据进行高精度重心修正的能力。通过解决基于 IRAF 的工作流过时问题，并使用当前标准（DE440）更新底层历表，该工具确保了用于脉冲星研究及其他毫秒级精度应用的计时数据不再受到过时轨道模型或时间标准歧义的影响。作者强调，该工具在处理河外目标时实现了与独立验证代码在 1 ms 级别的吻合，在处理太阳系目标时约为 5 ms 的吻合，满足了现代 HST 科学案例的精度要求。