[论文解读] A Bayesian Population Model for the Observed Dust Attenuation in Galaxies

Gautam Nagaraj, John C. Forbes+3 more|arXiv (Cornell University)|Feb 10, 2022

Galaxies: Formation, Evolution, Phenomena参考文献 97被引用 31

一句话总结

本文提出一种分层贝叶斯群体模型，利用恒星质量、SFR、金属丰度、红移和倾角作为预测变量，对来自3D-HST巡天的30,000个星系的尘埃消光特性（如消光斜率n和光学深度τ）进行联合推断。该模型采用灵活的5维线性插值，结合边缘化不确定性和基于投影的平均方法，以解决参数退化问题，揭示出τ随SFR增加而增大，随光学深度增大而趋于平坦，并且在红移和几何结构方面表现出显著依赖性，尤其在大质量、高SFR星系中更为明显。

ABSTRACT

Dust plays a pivotal role in determining the observed spectral energy distribution (SED) of galaxies. Yet our understanding of dust attenuation is limited and our observations suffer from the dust-metallicity-age degeneracy in SED fitting (single galaxies), large individual variances (ensemble measurements), and the difficulty in properly dealing with uncertainties (statistical considerations). In this study, we create a population Bayesian model to rigorously account for correlated variables and non-Gaussian error distributions and demonstrate the improvement over a simple Bayesian model. We employ a flexible 5-D linear interpolation model for the parameters that control dust attenuation curves as a function of stellar mass, star formation rate (SFR), metallicity, redshift, and inclination. Our setup allows us to determine the complex relationships between dust attenuation and these galaxy properties simultaneously. Using Prospector fits of nearly 30,000 3D-HST galaxies, we find that the attenuation slope ($n$) flattens with increasing optical depth ($ au$), though less so than in previous studies. $ au$ increases strongly with SFR, though when $\log~{ m SFR}\lesssim 0$, $ au$ remains roughly constant over a wide range of stellar masses. Edge-on galaxies tend to have larger $ au$ than face-on galaxies, but only for $\log~M_*\gtrsim 10$, reflecting the lack of triaxiality for low-mass galaxies. Redshift evolution of dust attenuation is strongest for low-mass, low-SFR galaxies, with higher optical depths but flatter curves at high redshift. Finally, $n$ has a complex relationship with stellar mass, highlighting the intricacies of the star-dust geometry. We have publicly released software (https://github.com/Astropianist/DustE) for users to access our population model.

研究动机与目标

• 通过将尘埃消光建模为多个星系属性的函数，解决Sed拟合中的尘埃-金属丰度-年龄退化问题。
• 通过考虑群体星系数据中相关参数和非高斯误差分布，改进对尘埃消光的推断。
• 开发一种灵活且公开可用的5维插值模型，用于描述依赖于恒星质量、SFR、金属丰度、红移和倾角的消光曲线。
• 证明群体水平建模优于简单贝叶斯拟合，通过整合数据密度和边缘化不确定性实现性能提升。

提出的方法

• 采用具有超参数θ的分层贝叶斯模型，θ代表在星系属性5维空间中控制点处的尘埃参数（n, τ）的插值值。
• 使用5维线性插值建模尘埃消光如何随M∗、SFR、Z∗、z和b/a（倾角）变化，实现平滑且灵活的预测。
• 通过对模型超参数后验分布和Prospector后验样本进行马尔可夫链蒙特卡洛采样，计算边缘化期望。
• 采用加权投影方法，利用数据密度作为权重函数，对被抑制的维度（如红移、金属丰度）进行平均，避免稀疏区域带来的偏差。
• 在聚合星系样本时，通过引入逆先验权重（1/p(tq, uq|α)）校正Prospector的中间先验。
• 通过在M个群体模型后验样本和Q个星系-Prospector样本对上对参数空间积分进行近似，实现可扩展的推断。

实验结果

研究问题

• RQ1 在统计上严谨的前提下，尘埃消光（n和τ）如何随恒星质量、SFR、金属丰度、红移和倾角变化？
• RQ2 年龄、金属丰度和尘埃之间的退化程度在多大程度上影响Sed拟合？在群体水平上如何缓解这一问题？
• RQ3 尘埃消光与星系属性之间的关系如何随红移演化，特别是在低质量、低SFR系统中？
• RQ4 灵活的5维插值模型能否在处理非高斯不确定性的同时，准确恢复合成数据中的已知趋势？

主要发现

• 消光斜率（n）随光学深度（τ）增加而变平，但程度弱于以往研究，表明尘埃几何结构响应具有非线性特征。
• 光学深度（τ）强烈依赖于SFR，但当log SFR ≤ 0时，τ在大范围内恒星质量下几乎保持不变。
• 仅当log M∗ ≥ 10时，侧向星系的τ才高于正面星系，反映出低质量星系缺乏三轴结构。
• 尘埃消光的红移演化在低质量、低SFR星系中最为显著，这些星系在高红移时表现出更高的τ和更平坦的曲线。
• n与恒星质量之间的关系复杂，揭示出非单调行为，凸显了星-尘几何结构的重要性。
• 该模型在测试模拟中成功恢复了合成趋势，系统偏差极小，且内在离散度估计准确（log σ恢复值为−4.95和−5.09，与模拟值−5.30接近）。