← 返回 PaperDaily
视觉与图像
普林斯顿生成式框架让0.2m/s行星信号无处遁形
让AI看星星,居然能比人眼强26倍!普林斯顿团队用生成式模型从噪声中“捞”出系外行星信号,地球2.0的搜寻又近了一大步。
龙哥读论文
发布于 2026-08-14 09:10:36
阅读 3
查看原文
🐉 龙哥读论文知识星球来了! 公众号每日8篇拆解不够看?星球 无上限更AI领域论文、资讯、招聘、招博、开源代码, 一站式干货,每日2分钟刷完即赚!
👇扫码加入「龙哥读论文」知识星球,前沿干货、实用资源一站式拿捏~
龙哥推荐理由: 让AI看星星,居然能比人眼强26倍!普林斯顿团队用生成式模型从噪声中“捞”出系外行星信号,地球2.0的搜寻又近了一大步。
原论文信息如下:
龙哥敢说,99%的天文爱好者都低估了找一颗“地球2.0”的难度。别看望远镜一个比一个牛,仪器精度卷到0.1米/秒,可恒星自己就是个“话痨”——它表面那点活动(震荡、米粒组织、磁场斑块)引起的视向速度抖动,动辄就是1米/秒上下。地球这样的行星产生的引力信号才0.1米/秒,完全被“噪音”淹没了。再加上地球大气的水汽吸收、仪器自身的不稳定……想从这一团乱麻里把行星信号“拎”出来,简直像在暴风雨里听一声树叶落地。
最近,普林斯顿大学天体物理系Yan Liang团队在arXiv上挂出了一篇重量级论文——AESTRA II ,首次将生成式光谱建模框架应用到真实的太阳光谱数据上,用“去噪+解耦”的思路,一下子把系外行星的探测灵敏度从传统方法的0.5 m/s拉到了0.2 m/s以下,探测数量暴涨26倍!龙哥看完直呼:这才是AI+天文该有的样子。
新颖点:生成式光谱建模如何实现“去噪提纯”
传统方法处理恒星光谱,通常需要先做“减法”——用外部大气传输模型扣除地球大气吸收,用多项式拟合扣除连续谱,再用模板拟合恒星线形。问题在于:这些模型总有误差,尤其那些微弱的“微地球大气线”(深度<1%),很容易被漏掉或矫正过头。仪器本身的线型函数(LSF)也随时间、波长变化,固定模版根本搞不定。结果就是:恒星活动和行星信号搅在一起,传统的交叉相关函数(CFF)提取的视向速度序列,再用活动指标(如Ca II H&K发射线、CFF的宽度、深度等)去回归,丢掉了太多光谱细节,效果天花板很低。
AESTRA的核心思想是“不做减法,做分解”。它不是先扣除再分析,而是用一个生成式模型,把每一条观测光谱同时 分解成三个物理可解释的分量:
1) 地球大气吸收分量(Telluric) ——固定在地球参考系中的窄吸收线,强度随水汽含量变化;
2) 连续谱变化分量(Continuum) ——光滑的、低频的乘性变化,来自仪器响应和大气透射的慢变;
3) 恒星线形变化分量(Stellar) ——集中在恒星吸收线附近的线形扭曲,由恒星表面活动驱动。
这三个分量的分离,不需要任何外部大气模型、恒星模板或预先定义的线型函数——完全靠数据自己学出来。模型仅通过“弱结构约束”(比如地球大气分量必须窄且静止在地球参考系,连续谱必须平滑,恒星分量必须宽且在恒星参考系)来分离物理来源。这种设计,从根本上避免了传统方法中模型误差的引入。
图2:AESTRA工作流程概览。(a) 每个阶梯阶被分解为地球大气吸收、连续谱变化和恒星线形变化;移除地球大气和连续谱后,得到清洁光谱。(b) 清洁光谱用于测量视向速度,然后编码为恒星活动潜在向量,再映射到活动驱动的径向速度,与候选行星信号联合优化。
最妙的是,AESTRA在训练时只用了光谱重建作为主要损失,没有预先知道哪个分量代表什么。但学出来的结果却“自洽”地分开了——论文中做了严格的验证:提取出的地球大气分量与NEID数据管道(DRP)的可降水汽值(PWV) 相关系数超过0.9,说明它确实捕捉到了真实的大气变化;三个分量的形态也截然不同(图4),一个窄而稀疏,一个光滑,一个集中在谱线附近,没有混淆。
图1:NEID太阳观测中的原始光谱变化性。顶部:某阶的CF视向速度时间序列,高亮12个时段。中部:这些时段的平均光谱。底部:相对于全局中值光谱的残差,展示多尺度光谱变化性。
原理揭秘:从光谱到视向速度的分解与重建
我们具体看看AESTRA是怎么一步步把“脏”光谱变成干净信号的。整个流程分两大阶段:先是每阶光谱的逐阶分解与清洗 ,然后是全谱活动编码与视向速度推断 。
图3:地球大气分量示例。灰线为观测光谱,彩色线为推断的地球大气模型,颜色代表最深处线的强度。右侧小图显示提取的地球大气强度与NEID DRP可降水汽值(PWV)紧密相关,相关系数>0.9。
第一步:波长对齐与预处理
NEID光谱仪有121个阶梯阶,覆盖380-930 nm,但信号质量在两端较差。AESTRA只保留中间42个高质量阶(430-623 nm)。每阶9216个像素,但只有中间30-80%的像素在自由光谱范围(FSR)内是可靠的。数据预处理包括:用NEID DRP提供的blaze函数除噪,按中值归一化,移除负值像素和单像素坏点。然后用NEID的重心校正去除地球公转速度,将光谱对齐到恒星静止参考系中。火前的NEID数据有一个零点头偏,AESTRA额外用一个平滑的样条函数对火前数据做经验波长校正,让它与火后数据统一到同一个波长网格上。
第二步:生成式光谱分解
对每个阶、每次曝光,AESTRA将观测光谱建模为:
ymodel,i = (1 - yt,i) ⊙ (1 + yc,i) ⊙ (y ★ + Δy★,i) + bi
其中yt,i是地球大气吸收(乘性),yc,i是连续谱变化(加性乘性),y ★是训练的时不变恒星模板,Δy★,i是时变的恒星线形变化,bi是小偏置。所有变量都是一维光谱数组,⊙表示逐元素乘法。整个模型通过一个共享编码器+三个解码器分支实现,编码器将预处理后的光谱压缩成潜在向量,再分别解码出三个分量。
为了确保分量不混叠,模型加入了弱结构正则化:地球大气分量只能在波长点附近有非零值、且必须固定在地球参考系(即波长不随恒星运动移动);连续谱分量被惩罚为光滑的低频结构;恒星分量只能出现在恒星吸收线位置附近。这些先验非常弱,只是引导方向,不依赖任何外部模型。
训练完成后,只移除地球大气和连续谱分量——恒星分量被直接丢掉,因为它的作用是“兜底”吸收恒星变化,防止它们被误分到大气或连续谱里。这样得到的就是所谓的“活动保持光谱”(activity-preserving spectra),里面只留下了恒星线形变化。
第三步:合并与活动编码
清洗后的42个阶被合并成一个连续的一维光谱。然后,AESTRA再次使用编码器-解码器结构,将合并后的光谱压缩成一个低维潜在向量(latent vector)。这个潜在向量相当于对恒星活动状态的“紧凑描述”——它比传统的CF活动指标(宽度、深度、不对称度)携带了更多信息,因为它直接来自全谱、全阶的联合学习。
接下来,这个潜在向量通过一个小的神经网络映射到“活动驱动视向速度”(activity RV)。同时,模型还会联合拟合一个或多个开普勒信号作为候选行星——这样,活动信号和行星信号就在同一个概率框架下被分离了。注意,行星信号(正弦曲线)和活动信号(来自潜在向量)是同时优化的,可以自然地去除相关性。
图4:光谱分解的三个分量示例。从上到下:地球大气分量(窄且稀疏,仅在地球大气线位置有值)、恒星分量(集中在恒星吸收线附近,捕捉线形不对称变化)、连续谱分量(光滑低频变化)。三者形态完全不同,证明分离成功。
结果亮眼:探测能力相比传统方法提升一个数量级
为了公平评估探测能力,论文设计了一套极其严格的注入-恢复测试 :在真实NEID太阳观测数据中注入50个(实际重复多次总数为50次独立试验?原文说50个单行星注入测试,周期2.5-40天,半振幅0.1-0.7 m/s)不同参数的单行星信号。注意,这些信号是注入到真实太阳数据中,所以恒星活动、地球大气、仪器噪声都真实存在。然后,用AESTRA和传统CF+活动指标去趋势方法(基线方法)分别去探测,且探测过程是“周期盲”——分析者不知道注入的信号参数,必须自己搜索候选者。只有当候选者与注入信号在周期和相位上匹配时,才算正确恢复;其他被接受的信号都算假阳性。
为了保证比较公平,两种方法的假阳性率(FPR)被校准到完全相同(零假阳性)。在这个基准下,结果令人震惊:
传统CF方法: 只恢复了9颗行星,且没有一颗的半振幅低于0.5 m/s。在K=0.1-0.3 m/s的地球类比范围内,恢复数为零。
AESTRA: 恢复了238颗行星,其中包括13颗半振幅低于0.3 m/s的信号!
什么概念?传统方法在0.5 m/s以下完全失明,AESTRA却能在0.1-0.3 m/s区间恢复13颗。这相当于探测灵敏度直接提升了1个数量级(从0.5 m/s降到0.2 m/s)。论文还展示了恢复率与行星参数的详细关系,AESTRA在所有周期和半振幅区间都碾压传统方法。
表1:NEID太阳观测数据筛选和质量控制步骤。包括火后关闭期、硬件问题、日食、大气质量、云过境、低日照天数等筛选条件,最终保留约30,000幅高质量光谱用于分析。
局限与展望:迈向地球孪生探测的下一步
当然,AESTRA II还不是完美无缺的。目前所有的验证都是在太阳数据上做的——太阳是距离最近、信号最强的恒星,而且我们同时拥有高分辨率的全日面图像(如SDO/HMI)来验证活动指标。对于更暗的恒星,光谱信噪比会显著下降,AESTRA的分解能力和活动编码效果可能会打折扣。此外,AESTRA目前仅使用单个编码器处理合并后的全谱,没有显式地利用空间信息或时间相关性来约束活动演化。论文作者也提到,模型假设活动特征只在恒星吸收线位置出现,但某些活动现象(如斑点、斑黯)会同时影响连续谱和线形,目前的弱约束可能还不够强。
展望未来,AESTRA框架可以很自然地扩展到多波段联合建模、时间序列建模(如GP先验)以及更复杂的行星系统(多行星、含月卫)。更重要的是,它提供了一条清晰的路径:用生成式模型替代手工特征工程,直接在全谱级别分离信号源。如果能在更多恒星(如HARPS-N、EXPRES的数据)上验证成功,那“地球2.0”的发现真的可以提上日程了。
龙迷三问
Q1: 什么是视向速度法?为什么要在太阳上测试? A1: 视向速度法(Radial Velocity,RV)是通过测量恒星相对于地球的径向运动引起的谱线多普勒位移来间接发现行星的方法。太阳是唯一一颗我们可以同时获得高分辨率光谱和独立活动图(如黑子、光斑分布)的恒星,所以是验证新RV算法的理想“实验室”。
Q2: AESTRA中的“生成式”是什么意思?跟传统机器学习有什么区别? A2: 生成式模型(Generative Model)是指模型学习观测数据的联合概率分布,能够生成 与观测相似的样本。在AESTRA中,模型不仅学到如何重构光谱,还学到了不同物理过程如何生成 光谱中的特征(即三个分量的生成过程)。传统方法通常只是“判别式”——比如直接用活动指标回归RV,或者用模板匹配做去噪。AESTRA则是在“模拟”光谱的生成过程,因此能够实现更物理地分解。
Q3: 论文中说“零假阳性”,是怎么做到的? A3: 在注入-恢复测试中,作者通过调整行星信号检测的显著性阈值,使得在无注入信号的真实数据上探测到的候选者数量为零(即零假阳性率)。然后保持这个阈值不变,再去计算注入信号后的恢复率。这样确保了AESTRA和传统方法在相同的置信水平下比较,结果才公平。
如果你还有哪些想要了解的,欢迎在评论区留言或者讨论~
龙哥点评
论文创新性分数: ★★★★✰
将生成式光谱分解引入EPRV领域,实时数据上的首次应用,端到端无模板联合建模思路新颖,但分解层序结构已有先例。
实验合理度: ★★★★★
严格的零假阳性校准+周期盲恢复测试,对比基线选择合理(传统CF+活动指标),实验设计非常严谨,结果可信度高。
学术研究价值: ★★★★★
开辟了EPRV领域的新范式——从手工特征工程转向生成式全谱建模,对后续系外行星探测和恒星活动建模有重要启发意义。
稳定性: ★★★✰✰
在太阳高信噪比数据上表现优异,但尚未在暗弱恒星或低信噪比数据上验证,真实望远镜观测中的稳定性待考察。
适应性以及泛化能力: ★★★✰✰
目前仅针对太阳光谱(分辨率R>100,000,高信噪比),对不同光谱类型、不同仪器、不同活动强度的泛化能力尚未验证。
硬件需求及成本: ★★☆☆☆
需要训练生成式神经网络,GPU资源消耗较大;推理阶段也需要编码器-解码器,实时处理高采样率数据有挑战。
复现难度: ★★☆☆☆
代码未开源,需要自行实现生成式模型和训练流程;依赖NEID特定数据格式和预处理,复现门槛较高。
产品化成熟度: ★★✰✰✰
目前处于科学验证阶段,远未达到产品级;需要在更多望远镜和恒星上验证后,才可能集成到现有EPRV数据处理管线。
可能的问题: 生成式分离的三个分量可能存在退化(比如弱地球大气线和弱恒星线混叠),论文虽然验证了相关性,但未对极端情况(如极低水汽、强活动)做测试。另外,仅用光谱重建损失可能导致模型过度拟合噪声,但作者使用了弱正则化,这一问题有待进一步分析。
主要参考文献
Y. Liang et al., "AESTRA II: Generative Spectral Modeling of the Sun as a Star for Precise Radial Velocities," arXiv:2606.13574v1, 2026.
Y. Liang et al., "AESTRA: A Generative Spectral Modeling Framework for Exoplanet Radial Velocities," ApJ, 2024.(介绍AESTRA I的原始方法)
C. Schwab et al., "NEID: A High Precision Doppler Spectrometer for WIYN," SPIE, 2016.
E. B. Ford et al., "Solar Radial Velocity Variability with NEID," AJ, 2024.
*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击 "阅读原文", 查看更多原论文细节哦!
欢迎加入龙哥读论文粉丝群,
扫描下方二维码或者添加龙哥助手微信号加群 :kangjinlonghelper。
一定要备注:研究方向+地点+学校/公司+昵称(如 图像处理+上海+清华+龙哥) ,根据格式备注,可更快被通过且邀请进群。
『龙哥读论文』微信群目前包含:图像处理、大模型及智能体、自动驾驶及机器人、AI医疗及AI金融5个群