[论文解读] Equilibrium chemistry down to 100 K - Impact of silicates and phyllosilicates on carbon/oxygen ratio

P. Woitke, Ch. Helling+6 more|arXiv (Cornell University)|Dec 4, 2017

Astro and Planetary Science参考文献 104被引用 58

一句话总结

本文介绍 GGchem，一种快速的热化学平衡计算代码（可降至 100 K），涵盖气相平衡化学与凝结，并分析硅酸盐和层硅酸盐凝结如何改变气相 C/O 比以及尘埃成核（尤以钨为例）的形成种子作用。

ABSTRACT

We introduce a fast and versatile computer code, GGchem, to determine the chemical composition of gases in thermo-chemical equilibrium down to 100 K, with or without equilibrium condensation. We review the data for molecular equilibrium constants, kp(T), from several sources and discuss which functional fits are most suitable for low temperatures. We benchmark our results against another chemical equilibrium code. We collect Gibbs free energies, dG(T), for about 200 solid and liquid species from the NIST-JANAF database and the geophysical database SUPCRTBL. We discuss the condensation sequence of the elements with solar abundances in phase equilibrium down to 100 K. Once the major magnesium silicates Mg2SiO4[s] and MgSiO3[s] have formed, the dust/gas mass ratio jumps to a value of about 0.0045 which is significantly lower than the often assumed value of 0.01. Silicate condensation is found to increase the carbon/oxygen ratio (C/O) in the gas from its solar value of ~0.55 up to ~0.71, and, by the additional intake of water and hydroxyl into the solid matrix, the formation of phyllosilicates at temperatures below ~400 K increases the gaseous C/O further to about 0.83. Metallic tungsten (W) is the first condensate found to become thermodynamically stable around 1600 - 2200 K (depending on pressure), several hundreds of Kelvin before subsequent materials like zirconium dioxide (ZrO2) or corundum (Al2O3) can condense. We briefly discuss whether tungsten, despite its low abundance of ~2.E-7 times the silicon abundance, could provide the first seed particles for astrophysical dust formation. The GGchem code is publicly available at https://github.com/pw31/GGchem.

研究动机与目标

• 开发并验证一个快速、通用的代码（GGchem），以在热化学平衡下（可含或不含平衡凝结）确定气相丰度，温度降至 100 K。
• 回顾并基准化来自多个来源（NIST-JANAF、SUPCRTBL、Barklem & Collet）的分子和凝聚态热化学数据，以评估一致性和低温可靠性。
• 分析凝结序列（尤其是 Mg-硅酸盐）如何改变元素预算以及气相中的 C/O 比，并讨论在 T < 400 K 时层硅酸盐的作用。
• 研究第一批凝结物（如钨）及其作为天体尘埃形成种子的潜在作用。
• 向社区公开 GGchem 代码及数据。

提出的方法

• 通过最小化 Gibbs 自由能（可含或不含凝结）来求解气相化学丰度。
• 使用 Gibbs 自由能数据推导分子形成的平衡常数 k_p(T)（k_p 来自 ΔG_f°），并采用温度相关拟合（选取 Stock 等式形式）。
• 通过超饱和比 S_j 和凝结蒸汽压 p_vap(T) 纳入相平衡；通过凝结迭代地降低气相中的元素丰度。
• 采用牛顿-拉夫森迭代，嵌套调用平衡化学与四倍精度算术，在低温（降至 100 K）下确保数值稳定性。
• 使用 NIST-JANAF 和 SUPCRTBL 的 ΔG_f°(T) 数据来获取凝结物的 ΔG_f°(T)，并用鲁棒的函数对多数据集进行拟合，以实现低温外推。
• 提供公开的 GGchem 实现（FORTRAN-90），包含更新的数据和新增的预迭代方案。

实验结果

研究问题

• RQ1 平衡凝结如何在温度降至 100 K 时改变气相元素丰度及总体化学组分？
• RQ2 对于太阳组成，元素的凝结序列是什么，硅酸盐和层硅酸盐的形成如何影响气相的 C/O 比？
• RQ3 哪些凝结物最先形成（如 Mg2SiO4、MgSiO3、钨），在何种条件下它们变得热力学稳定？
• RQ4 尽管钨的丰度较低，钨是否可以作为早期的尘埃形成种子？
• RQ5 在 100–6000 K 范围内，不同的热化学数据源在预测分子和凝聚相平衡方面有何差异？

主要发现

• 凝结的主要镁硅酸盐（Mg2SiO4 和 MgSiO3）引发尘埃/气体质量比跃升至约 0.0045。
• 硅酸盐凝结将气相 C/O 从太阳值约 0.55 提升至约 0.71，而在 ~400 K 以下的层硅酸盐形成可进一步将 C/O 提高到 ~0.83。
• 钨（W）在约 1600–2200 K（压力相关）时成为首个在热力学上稳定的凝结物，领先于 ZrO2 与 Al2O3，且可能成为天体尘埃形成的种子。
• GGchem 代码，采用四倍精度和更新的预迭代策略，能在 100 K 处给出准确结果，且可公开获取于 GGchem GitHub 存储库。
• 对来自多源的 k_p(T) 数据进行全面比较，发现某些分子在 ~10 kJ/mol 量级存在显著偏差，强调数据选择和拟合形式对低温外推的重要性。
• 凝结过程的耦合显著影响气相丰度和光谱预测，强调在冷却/系外行星环境中建模平衡凝结的必要性。