跳到主要内容
返回时间线
arXiv来源发表:

标准汽笛宇宙学不再每次从先验重跑:ASPIRE复用旧后验,暗汽笛分析提速约29倍

核心概要

标准汽笛宇宙学要联合推断宇宙学与致密双星族群参数,每换一个星系表、质量模型或事件目录,通常都要从先验重新跑一遍完整分析。ASPIRE把已有后验样本经归一化流拟合后,用序贯蒙特卡洛退火到新假设下的目标分布,在GWTC-5.0的暗汽笛算例中给出与嵌套采样一致的后验,并保留贝叶斯证据估计。以谱汽笛结果为起点加入GLADE+星系表信息时,采样时间从134.4小时降到4.6小时。

AI-generated editorial illustration: Efficient Standard Siren Cosmological Inference with Posterior Reuse

深度剖析

标准汽笛宇宙学分析可以不再从先验重新开始,而是把上一次分析留下的后验样本直接改造成新假设下的后验。 此前每探索一个族群模型、星系表或加权方案,都要独立跑一次完整推断,先前算出的信息在目录更新后被丢弃。 在GWTC-5.0的暗汽笛算例中,ASPIRE给出的后验与嵌套采样dynesty的结果一致,全部23个参数都落在设定的Jensen-Shannon散度阈值之内。

做法是先用归一化流把旧后验拟合成一个可采样、可求值的提议分布,再沿一串中间分布把它推到新目标。 归一化流是一种神经网络,学习简单基分布与目标分布之间的可逆映射;序贯蒙特卡洛则按逆温度逐步退火,每一步用重要性权重重采样并做MCMC移动。 退火步长按目标有效样本量自动确定,到最后一个温度步时提议分布已不留在后验中;在谱汽笛结果上加入GLADE+星系表信息的算例里,采样耗时4.6小时,而dynesty为134.4小时。

启示与展望

演示范围限于GWTC-5.0数据、GLADE+星系表与FullPop-4.0质量模型下的若干暗汽笛与谱汽笛算例,且报告的加速比假定提议分布已经存在,不计生成提议的代价。

把GWTC-4.0的谱汽笛结果更新到GWTC-5.0时,后验内部相关性发生变化,需要把四个与哈勃常数相关的参数改为从先验抽取;从较窄的固定宇宙学后验扩展到更宽目标时,重加权与重采样无法在提议没有样本的区域生成样本。