[论文解读] Partition functions 1: Improved partition functions and thermodynamic quantities for normal, equilibrium, and ortho and para molecular hydrogen

A. Popovas, Uffe Jørgensen|arXiv (Cornell University)|Jul 15, 2016

Advanced Chemical Physics Studies参考文献 31被引用 15

一句话总结

本论文针对分子氢在1–20,000 K温度范围内的高精度配分函数和热力学量进行了研究，通过严格求和实验与理论计算的基态和激发电子态的能级数据，实现了高精度建模。研究结果表明，常见的简化方法（如简谐振子/刚性转子近似或人为设定截断）可引入高达40%的误差，因此本文主张采用杜纳姆系数（Dunham coefficient）建模方法，并完整包含束缚态与准束缚态。

ABSTRACT

Aims. In this work we rigorously show the shortcomings of various simplifications that are used to calculate the total internal partition function. These shortcomings can lead to errors of up to 40 percent or more in the estimated partition function. These errors carry on to calculations of thermodynamic quantities. Therefore a more complicated approach has to be taken. Methods. Seven possible simplifications of various complexity are described, together with advantages and disadvantages of direct summation of experimental values. These were compared to what we consider the most accurate and most complete treatment (case 8). Dunham coefficients were determined from experimental and theoretical energy levels of a number of electronically excited states of H$_2$ . Both equilibrium and normal hydrogen was taken into consideration. Results. Various shortcomings in existing calculations are demonstrated, and the reasons for them are explained. New partition functions for equilibrium, normal, and ortho and para hydrogen are calculated and thermodynamic quantities are reported for the temperature range 1 - 20000 K. Our results are compared to previous estimates in the literature. The calculations are not limited to the ground electronic state, but include all bound and quasi-bound levels of excited electronic states. Dunham coefficients of these states of H$_2$ are also reported. Conclusions. For most of the relevant astrophysical cases it is strongly advised to avoid using simplifications, such as a harmonic oscillator and rigid rotor or ad hoc summation limits of the eigenstates to estimate accurate partition functions and to be particularly careful when using polynomial fits to the computed values. Reported internal partition functions and thermodynamic quantities in the present work are shown to be more accurate than previously available data.

研究动机与目标

• 解决现有分子氢配分函数计算中存在的显著不准确性，这些误差会传播至热力学量及天体物理学模型中。
• 识别并量化广泛使用的简化方法（如简谐振子与刚性转子近似）以及人为设定的求和截断所存在的缺陷。
• 通过完整包含H₂基态与激发电子态的所有束缚态与准束缚态，提供更精确、更完整的内部配分函数处理方法。
• 报告H₂基态与激发电子态的杜纳姆系数，以实现对振转能级的精确表征。
• 在宽温度范围内，为平衡态、正常态、正氢与仲氢提供改进的、高精度的配分函数与热力学量（内能、熵、比热等）。

提出的方法

• 对H₂基态及多个激发电子态的实验测量与理论计算能级（包括束缚态与准束缚态）进行直接求和。
• 从高精度能级数据中确定杜纳姆系数（Yₖₗ），以建模振转能级，避免依赖经典振子-转子类比。
• 分别计算平衡态、正常态、正氢与仲氢的配分函数，考虑核自旋统计与简并度。
• 在配分函数求和中包含所有相关电子态，而不仅限于基态，以确保完整性。
• 采用所有可访问能级的完整玻尔兹曼求和：Q(T) = Σ gₙ exp(–Eₙ / kT)，其中gₙ为简并度，Eₙ为能级n的能量。
• 通过标准统计力学关系从配分函数推导热力学量（U, Cᵥ, S, F）：U = kT² ∂lnQ/∂T, Cᵥ = ∂U/∂T, 等等。

实验结果

研究问题

• RQ1 常见简化方法（如简谐振子与刚性转子近似）在计算分子氢配分函数时引入的定量误差有多大？
• RQ2 对转动与振动量子数人为设定的截断如何影响配分函数及其导出的热力学量的准确性？
• RQ3 对计算配分函数数据进行多项式拟合在多大程度上无法捕捉H₂在高温下的真实行为？
• RQ4 包含激发电子态（束缚态与准束缚态）对H₂总内部配分函数有何影响？
• RQ5 与分子常数相比，基于实验与理论能级数据导出的杜纳姆系数在多大程度上提升了配分函数建模的准确性？

主要发现

• 常见简化方法（如简谐振子与刚性转子近似）在配分函数计算中可引入高达40%或更多的误差，尤其在高温区域更为显著。
• 对量子数人为设定的求和截断会导致配分函数显著低估，尤其在1000 K以上温度区域。
• 对配分函数数据进行多项式拟合不足以满足高精度应用需求，因其无法捕捉真实配分函数的非多项式行为。
• 与仅基于基态的计算相比，包含H₂所有激发电子态的束缚态与准束缚态能级，可使20,000 K时的配分函数提高最多达15%。
• 本文报告的新配分函数与热力学量比以往文献值更为精确，且在天体物理温度区间内具有更高的内在一致性。
• 本研究以高精度报告了H₂基态与激发电子态的杜纳姆系数，为未来H₂能级与配分函数的高精度建模提供了基础。