2025-11-23T04:58:16.808235

Assessing the star formation history of all-sky and part-sky 100pc white dwarf samples

Roberts, Tremblay, Bédard
Thanks to Gaia and large-scale spectroscopic follow-up surveys (4MOST, DESI, WEAVE, SDSS-V), it is now possible to build representative and minimally biased samples of the local white dwarf population. Here we analyse several volume-limited 100pc samples of white dwarfs, constructed from different surveys and studies, to evaluate their completeness and residual biases. We model the underlying star formation history and Galactic disc age via comparison with simulated populations of white dwarfs to quantitatively characterise completeness. We assess whether the benefit of Gaia XP spectra in datasets outweighs the reduction in sample size, and to what extent targeted, part-sky, and magnitude limited surveys can be used in comparison to all-sky volume limited surveys. Additionally, we simulate the 4MOST 100PC sub-survey and discuss its use to better understand the local star formation history.
academic

Assessing the star formation history of all-sky and part-sky 100pc white dwarf samples

基本信息

  • 论文ID: 2510.20624
  • 标题: Assessing the star formation history of all-sky and part-sky 100pc white dwarf samples
  • 作者: Emily K. Roberts, Pier-Emmanuel Tremblay, Antoine Bédard (University of Warwick)
  • 分类: astro-ph.SR (Solar and Stellar Astrophysics)
  • 发表时间: 2025年11月21日 (Preprint)
  • 期刊: MNRAS (Monthly Notices of the Royal Astronomical Society)
  • 论文链接: https://arxiv.org/abs/2510.20624

摘要

本文利用Gaia卫星数据及大规模光谱后续巡天（4MOST, DESI, WEAVE, SDSS-V）构建的多个100秒差距（pc）体积限制白矮星样本，系统评估其完备性和残余偏差。通过模拟白矮星种群并与观测样本对比，定量刻画星族历史（SFH）和银河系盘年龄。研究评估了Gaia XP光谱的价值、部分天区巡天与全天巡天的适用性，并预测了4MOST 100PC子巡天对理解本地星族历史的贡献。

研究背景与动机

核心问题

白矮星作为恒星演化终点产物，其冷却时间与温度、光度和质量存在明确关系，可作为"宇宙时钟"反推星族历史。但这一潜力的实现依赖于：

  1. 完备且无偏的白矮星样本：Gaia提供了前所未有的大样本，但不同选择标准和巡天覆盖导致样本差异
  2. 精确的物理模型：需要白矮星冷却模型、主序寿命、初始质量函数（IMF）、初始-终态质量关系（IFMR）等
  3. 观测偏差的定量评估：光谱完备性、视星等限制、天区覆盖等因素的影响

研究重要性

  • 银河系考古学：通过白矮星重建太阳邻域的星族历史，与主序星结果对比验证
  • 样本选择策略：为未来大规模巡天（如Vera Rubin）提供样本设计指导
  • 方法学基准：40pc样本已建立方法学框架（Roberts et al. 2025），100pc样本可检验其外推性

现有方法局限

  • 40pc样本：完备性高（~97%）但统计样本小（1073颗），泊松噪声大
  • 100pc样本：样本量大（1-1.5万颗）但完备性和偏差未充分量化
  • 部分天区巡天：光谱覆盖仅~30-50%天区，其代表性存疑
  • 不同选择标准：Gaia DR3衍生的多个白矮星目录（Gentile Fusillo, GCNS, GSPC-WD等）差异显著

核心贡献

  1. 系统比较多个100pc白矮星样本：定量评估6个主要样本（3个全天+3个部分天区）的完备性和偏差
  2. 验证星族历史重建方法的稳健性：证明光度函数法和绝对星等分布法在不同样本间结果一致
  3. 量化XP光谱的价值：GSPC-WD样本虽有光谱分类但存在红端不完备性（Teff<6000K）
  4. 评估部分天区巡天适用性：DESI DR1表现出与全天样本相当的无偏性，SDSS则存在明显偏差
  5. 预测4MOST 100PC子巡天性能：模拟显示其视星等限制（G<19.5）会导致暗端不完备，但足以区分不同星族历史模型
  6. 确认最佳星族历史形式：晚峰型星族历史（Roberts et al. 2025）或恒定星族率均能很好拟合数据

方法详解

任务定义

输入：100pc内白矮星样本的观测数据（Gaia视差、星等、光谱参数等）
输出：本地银河系盘的星族历史SFH(t)和年龄
约束：体积限制（100pc球形或锥形区域）、质量范围（0.54-1.34 M⊙，排除双星产物）

模型架构：单星种群合成

核心流程

  1. 初始恒星生成：
    • 根据Salpeter IMF（α=2.35±0.075）生成初始质量
    • 按假定SFH(t)分配形成时间
    • 赋予初始金属丰度（Z=0.0134±0.0104）
  2. 恒星演化：
    • 主序+巨星阶段寿命：tMSGP（Byrne et al. 2025模型，±3.5%不确定性）
    • 通过IFMR转换为白矮星质量
    • 白矮星冷却：Bédard et al. (2020)模型（±4.5%不确定性）
    • 晶化延迟：开/关两种情况
    • 并合延迟：tΔ（43%概率开启）
  3. 运动学演化：
    • 银河系盘标高：h(t) = 10.71 × t^0.65 pc
    • 年龄-速度弥散关系：离开样本概率 p_left = 1 - 65/h
    • 修正缺失的主序星（未演化为白矮星）
  4. 观测模拟：
    • 大气组成：25%氦主导（±6.5%），DC型<5000K假定为氢大气
    • 星等计算：考虑消光（平均AV=0.02±0.015 mag）
    • 部分天区：应用实际巡天足迹（图3）
    • 视星等限制：DESI G<20, 4MOST G<19.5

不确定性处理

运行50次蒙特卡洛模拟，每次在表2所列参数空间内随机采样：

  • 系统不确定性：种群年龄（10.6±0.5 Gyr）、IMF斜率、主序寿命、冷却模型等
  • 个体不确定性：金属丰度、Gaia测量误差

三种星族历史重建方法

1. 光度函数法（Luminosity Function）

  • 观测量：白矮星绝对热光度Mbol分布
  • 计算：需Teff和半径R → 依赖质量-半径关系和大气模型
  • 挑战：需已知大气组成（H或He主导）
  • 处理：无光谱样本假定纯氢大气；有光谱样本使用定制模型

2. 绝对G星等分布法（Absolute Magnitude Distribution）

  • 观测量：Gaia绝对G波段星等MG分布
  • 优势：仅需视差和视星等，无需光谱参数
  • 模拟：种群合成中编码25%氦大气比例
  • 消光修正：使用Gentile Fusillo et al. (2021)的AV值

3. 直接年龄法（Direct Age Method，附录A）

  • 计算：ttotal = tWD + tMSGP + tΔ
  • 应用：需完整的Teff、质量、光谱型
  • 修正：缺失主序星、年龄-速度弥散效应
  • 结果：归一化星族历史dN/dt

技术创新点

  1. 低质量问题修正：Teff<6000K时Lyman-α红翼不透明度误差导致光度质量低估，应用O'Brien et al. (2024)五阶多项式修正至中值0.6 M⊙
  2. 概率性体积成员：考虑视差不确定性的100pc成员概率（图1），但发现<2000颗边界天体且P100pc<0.5，对结果影响可忽略
  3. 双星排除策略：质量限制0.54-1.34 M⊙排除双简并星和非单星演化产物，虽降低完备性但减少系统偏差
  4. 部分天区模拟：解析化巡天足迹（图3）模拟银河标高和年龄-速度关系的方向依赖效应

实验设置

数据集

全天100pc样本（3个）

样本源目录初始数PWD≥0.75质量筛选后特点
GF21Gentile Fusillo et al. (2021)16,29816,28115,001Gaia EDR3，纯氢大气假定
GCNSGaia Collaboration (2021)301,92816,31914,699随机森林分类，98%重叠GF21
GSPC-WDVincent et al. (2024)14,90614,80213,235Gaia XP光谱，定制大气模型
  • GF21：基准样本，完备性~93%（与SDSS足迹对比）
  • GCNS：独立选择但与GF21高度重叠（14,687/14,699）
  • GSPC-WD：需XP光谱，G<20.5，红端不完备（Teff<5000K）

部分天区100pc样本（3个）

样本覆盖率白矮星数光谱完备性特点
SDSS (Kilic et al. 2025)~30%2,75075% (Teff>5000K达86%)可靠性优先，避冷星
DESI DR1~30%2,978随机采样无偏随机目标选择
4MOST (预测)~50% (南天)~8,000目标完备性G<19.5，Dec: 5°至-80°

评价指标

  1. χ²拟合优度：模拟与观测分布的逐bin比较
  2. 最佳种群年龄：使光度函数暗端下降匹配的年龄
  3. 星族历史形式：四种假定SFH的相对表现
    • 恒定星族率（Cukanovaite et al. 2023）
    • 晚峰型（Roberts et al. 2025）
    • 早峰型（Fantin et al. 2019）
    • 双峰型（Mor et al. 2019）

对比基准

  • 40pc样本（O'Brien et al. 2024）：960颗白矮星，光谱完备性100%
  • 不同SFH模型：图7引自Roberts et al. (2025)

实现细节

  • 模拟次数：每个配置50次蒙特卡洛运行
  • 误差计算：标准差（系统不确定性）与泊松误差（观测）的平方和
  • 消光处理：缺失AV时从N(0.02, 0.015²)采样
  • 并行策略：不同SFH和样本独立模拟

实验结果

主要结果：全天样本

GF21样本（图4）

光度函数：

  • 最佳拟合：Roberts et al. (2025)晚峰型SFH，其次Mor et al. (2019)
  • 峰值问题：所有模拟在Mbol≈15处预测小幅上翘，观测显示下凹
  • 年龄约束：最佳拟合11.0 Gyr（1σ内符合默认10.6 Gyr）

绝对G星等分布：

  • 趋势一致：晚峰型和双峰型优于恒定和早峰型
  • Fantin模型：明显预测过多暗星、过少亮星

GCNS样本（图5）

  • 高度相似：与GF21结果几乎完全一致
  • 最佳年龄：精确10.6 Gyr
  • 结论：两个独立Gaia选择方法的交叉验证

GSPC-WD样本（图6）

  • 不完备性明显：光度函数暗端过早下降
  • 年龄偏差：需降至8.5-9.3 Gyr才能匹配（偏离3σ）
  • 原因：Teff<6000K（红色）白矮星缺失，G≈20视星等限制

关键发现：XP光谱提供光谱分类价值，但样本不完备性引入显著年龄偏差

主要结果：部分天区样本

SDSS样本（图7）

  • 不完备性：需年龄降至7.5-9.3 Gyr
  • 梯度不匹配：下降坡度与任何模拟均不符
  • 原因：有意避开无特征DC型白矮星（Teff<5000K）

DESI DR1样本（图8）

  • 优异表现：10.6 Gyr完美拟合所有SFH
  • 随机性优势：tile和目标均随机选择，无颜色/亮度偏差
  • 视星等限制：经验G<20排除~5%模拟样本（图9）
  • 全天外推：部分天区结果可代表全天性质

突破性结论：DESI的随机采样策略实现了部分天区的无偏表示

4MOST预测（图10-11）

  • 天区影响：南天50%覆盖对光度函数无显著影响（顶图）
  • 星等限制：G<19.5在暗端造成不完备（底图）
  • SFH区分能力：误差棒内可区分四种星族历史形式（图11）
  • 预期样本量：~8,000颗白矮星

消融实验：直接年龄法（附录A）

归一化星族历史dN/dt（图A1-A2）：

  • GF21与GCNS：与40pc样本高度一致，低早期星族率→近5 Gyr高星族率
  • GSPC-WD：峰值时间提前，5-10 Gyr陡增（红端偏差导致）
  • DESI：介于全天和SDSS之间
  • 年龄上限：12-13 Gyr bin，但误差>2 Gyr

关键实验发现

  1. 样本间一致性：
    • GF21、GCNS和40pc三者星族历史完全一致
    • 不同Gaia选择标准（HR图位置、概率阈值）影响可忽略
  2. 完备性等级：
    DESI (无偏) > GF21 ≈ GCNS (轻微) > GSPC-WD (中度) > SDSS (严重)
  3. 光谱价值权衡：
    • 光谱分类改善大气模型（GSPC-WD, SDSS）
    • 但选择偏差抵消甚至超过模型改进收益
    • 光度法和绝对星等法对光谱依赖低
  4. 峰值之谜：
    • 所有模拟在Mbol≈15预测上翘
    • 所有观测显示平滑或下凹
    • 可能原因：冷却物理（晶化？）或选择效应
  5. 40pc→100pc外推：
    • 泊松噪声显著降低（图2中12.5 mag凸起在100pc样本消失）
    • 系统不确定性仍主导（表2）
    • 方法学稳健性验证

相关工作

白矮星星族历史研究

  • 早期工作：Winget et al. (1987)首次提出用白矮星光度函数测年
  • SDSS时代：Rowell (2013), Fantin et al. (2019)基于SDSS样本的早峰型SFH
  • Gaia革命：
    • Mor et al. (2019): 20pc样本，双峰SFH
    • Cukanovaite et al. (2023): 40pc样本，恒定SFH，年龄10.6±0.7 Gyr
    • Roberts et al. (2025): 三方法对比，晚峰型SFH

主序星对比研究

  • Gallart et al. (2024), Fernández-Alvar et al. (2025): Gaia主序星SFH
  • 目标：白矮星与主序星结果的交叉验证

Gaia白矮星目录

  • DR2: Gentile Fusillo et al. (2019) - 26万候选体
  • EDR3: Gentile Fusillo et al. (2021) - 35.9万，PWD概率
  • DR3新增：
    • GCNS (Gaia Collaboration 2021): 随机森林分类
    • GSPC-WD (Vincent et al. 2024): XP光谱+SDSS合成测光
    • Jiménez-Esteban et al. (2023): 独立XP分析

地面光谱巡天

  • SDSS I-IV: Kleinman et al. (2013), Kilic et al. (2025) - 北天~30%
  • LAMOST: Cui et al. (2012) - 中国巡天
  • DESI: Cooper et al. (2023), Manser et al. (2024) - 随机采样策略
  • 未来: 4MOST (de Jong et al. 2019), WEAVE (Jin et al. 2024), SDSS-V (Kollmeier et al. 2025)

冷却物理

  • Bédard et al. (2020): 质量-半径关系和冷却模型
  • Caron et al. (2023), O'Brien et al. (2024): Lyman-α红翼修正
  • 晶化延迟：仍有争议的物理过程

结论与讨论

主要结论

  1. 星族历史稳健性：
    • 晚峰型SFH（Roberts et al. 2025）或恒定星族率最佳拟合所有样本
    • 结果不依赖具体Gaia选择方法或天区覆盖
    • 40pc与100pc样本一致性验证方法外推性
  2. 样本完备性排序：
    DESI DR1 (最优) > GF21 ≈ GCNS > GSPC-WD > SDSS
    • DESI随机采样实现部分天区无偏表示
    • XP光谱价值被不完备性抵消
  3. 年龄约束：
    • 银河系盘年龄10.6±0.5 Gyr（与40pc一致）
    • 视星等限制阻碍暗端精确测量
    • 系统不确定性>2 Gyr主导误差
  4. 部分天区适用性：
    • 随机采样策略（DESI, 4MOST）可代表全天
    • 目标选择偏差（SDSS）引入严重系统误差
    • 天区覆盖50%足够（银河标高影响小）
  5. 未来巡天建议：
    • 4MOST 100PC将提供~8,000颗南天白矮星
    • G<19.5限制可接受（能区分SFH模型）
    • 深度巡天（Vera Rubin）对暗端更关键

局限性

  1. 物理模型不确定性：
    • 冷却模型±4.5%（晶化物理未完全理解）
    • 主序寿命±3.5%（金属丰度、旋转效应）
    • IFMR散射（~0.1 M⊙，未充分建模）
  2. 单星假设：
    • 质量限制排除双星但并合历史简化处理
    • 43%并合延迟比例基于40pc样本经验
  3. 大气组成：
    • 25%氦大气比例粗糙
    • DC型<5000K组成不确定性未传播
  4. 小尺度结构：
    • Mbol≈15峰值差异未解释
    • 无法约束SFH的短时标变化（<1 Gyr）
  5. 体积限制：
    • 100pc内晕族白矮星污染可忽略
    • 更大体积需多种群分解

未来方向

  1. 双星种群：
    • 4MOST完整光谱可测双星参数
    • 质量比、轨道周期分布
    • 质量转移历史约束
  2. 化学丰度：
    • 光谱金属线测量
    • 星族-化学演化联合建模
    • 与主序星Fe/H对比
  3. 更深样本：
    • Vera Rubin LSST: r<24.5，探测200pc
    • 约束银河系盘年龄至±0.5 Gyr
    • 晕族-盘族分离
  4. 物理改进：
    • 22Ne沉降对晶化的影响
    • 磁场对冷却的影响（10%白矮星）
    • 大气对流3D模型
  5. 机器学习：
    • 端到端SFH推断（跳过中间参数）
    • 模拟-观测域适应
    • 不确定性量化

深度评价

优点

  1. 系统性与全面性：
    • 首次定量比较6个主要100pc样本
    • 三种独立方法交叉验证（光度、星等、年龄）
    • 覆盖全天与部分天区、光度与光谱样本
  2. 方法学严谨：
    • 50次蒙特卡洛模拟量化系统不确定性
    • 详细误差预算（表2）
    • 低质量问题修正（关键技术细节）
  3. 实用价值高：
    • 为4MOST等未来巡天提供设计指导
    • 验证部分天区策略的可行性
    • 明确XP光谱的价值权衡
  4. 结果稳健：
    • 不同样本间一致性强
    • 与40pc基准吻合
    • 物理解释合理（DESI随机性优势）
  5. 写作清晰：
    • 表1完整总结样本特性
    • 图3直观展示巡天足迹
    • 附录A提供完整第三方法

不足

  1. 峰值差异未解决：
    • Mbol≈15系统性偏差缺乏物理解释
    • 可能暗示冷却模型缺陷或选择效应
    • 需独立数据（如光谱Teff）验证
  2. 4MOST预测简化：
    • 未模拟光纤分配算法（邻源排斥）
    • 时间压力导致的不完备性未量化
    • 实际完备性可能低于预期
  3. SFH形式受限：
    • 仅测试4种参数化形式
    • 短时标变化（<1 Gyr）不可分辨
    • 系统不确定性掩盖精细结构
  4. 双星处理粗糙：
    • 质量限制简单粗暴
    • 并合延迟二元化（开/关）
    • 未探索质量比分布影响
  5. 统计检验不足：
    • 缺乏定量χ²值或p值
    • "最佳拟合"基于视觉判断
    • 贝叶斯模型比较更严格

影响力

学术贡献：

  • 为Gaia时代白矮星星族学建立基准
  • 证明部分天区巡天的代表性（重要方法学结论）
  • 量化XP光谱在大样本中的边际效用

实用价值：

  • 直接指导4MOST、WEAVE观测策略
  • 为DESI后续数据释放提供分析框架
  • Vera Rubin等深度巡天的设计参考

可复现性：

  • 代码基于Roberts et al. (2025)公开框架
  • 数据来自公开目录（Gaia DR3, DESI DR1）
  • 模拟参数完整列于表2

潜在影响：

  • 引用Roberts et al. (2025)基础上的扩展
  • 可能成为100pc样本选择的标准参考
  • 推动随机采样策略在MOS巡天中的应用

适用场景

直接应用：

  1. 太阳邻域星族历史研究（<100pc）
  2. 白矮星样本设计与质量控制
  3. Gaia DR4/DR5白矮星目录验证

扩展应用：

  1. 银河系盘-晕分解（>100pc样本）
  2. 星团白矮星年龄测定
  3. 初始-终态质量关系校准

不适用场景：

  1. 银河系中心/旋臂等非盘族环境
  2. 短时标星爆事件（分辨率不足）
  3. 单个白矮星精确测年（系统误差>10%）

参考文献（关键引用）

方法学基础：

  • Roberts et al. (2025, MNRAS 538, 2548): 40pc样本三方法对比
  • Cukanovaite et al. (2023, MNRAS 522, 1643): 恒定SFH基准

Gaia白矮星目录：

  • Gentile Fusillo et al. (2021, MNRAS 508, 3877): GF21样本
  • Gaia Collaboration et al. (2021, A&A 649, A6): GCNS目录
  • Vincent et al. (2024, A&A 682, A5): GSPC-WD XP光谱样本

地面光谱巡天：

  • Kilic et al. (2025, ApJ 979, 157): SDSS 100pc样本
  • DESI Collaboration et al. (2025, arXiv:2503.14745): DESI DR1
  • de Jong et al. (2019, Messenger 175, 3): 4MOST巡天设计

物理模型：

  • Bédard et al. (2020, ApJ 901, 93): 白矮星冷却与质量-半径关系
  • O'Brien et al. (2024, MNRAS 527, 8687): 低质量问题修正
  • Byrne et al. (2025, MNRAS 537, 2433): 主序寿命模型

总体评价：这是一篇高质量的天体物理方法学论文，通过系统比较多个100pc白矮星样本，为Gaia时代的星族历史研究提供了重要基准。其核心贡献在于量化了不同样本选择策略的偏差，特别是证明了DESI随机采样的优越性和XP光谱在大样本中的有限价值。研究对未来巡天（4MOST, Vera Rubin）具有直接指导意义，但在物理模型细节和统计检验方面仍有改进空间。论文的稳健结论（晚峰型或恒定SFH）为银河系考古学提供了可靠约束。