A generalized q-deformed Eckart model for the spectroscopic and thermodynamic properties of diatomic molecules

技术摘要：一种用于双原子分子的广义 q-变形埃卡特模型

问题陈述
准确的相互作用能模型对于解释双原子分子的光谱学、化学键合以及温度依赖的热力学性质至关重要。虽然量子化学方法（如密度泛函理论、耦合集群理论）能够提供极高的精度，但在生成完整的势能曲线或处理多种物种及温度点时，其计算量非常庞大。解析势能函数提供了一种补充途径，能够直接以核间距离表示相互作用能。然而，现有的解析势能函数往往缺乏完整的光谱参数化，无法将解离极限、平衡条件、势阱曲率以及高阶转动-振动关系联系起来。具体而言，广义 q-变形埃卡特（GqDE）势能尚未通过这种全面的参数化被确立为一种双原子相互作用能函数，也未通过独立的实验光谱或电子结构参考数据对其预测的曲线和能量进行验证。

研究方法
本研究开发了一种参数化的 GqDE 势能，并应用该模型来计算光谱学和热力学性质。

1. 势能参数化： GqDE 势能由受强度、范围和变形参数控制的双曲 q-变形函数定义。这些参数通过满足四个物理约束来推导：

   • 解离极限： 势能在无穷远处趋于解离能 ($D_e$)。
   • 平衡条件： 一阶导数在平衡键长 ($r_e$) 处为零。
   • 曲率： $r_e$ 处的二阶导数与由谐振振动波数 ($\omega_e$) 推导出的力常数相匹配。
   • 光谱形状： 三阶导数通过 Varshni 关系进行约束，该关系将势能形状与转动常数 ($B_e$) 和振动-转动相互作用常数 ($\alpha_e$) 联系起来。
     这些条件使得所有势能参数都可以用实验可获得的分子常数来表达。

2. 能量谱计算： 求解 GqDE 势能下的径向薛定谔方程。对于非零转动量子数 ($\ell \neq 0$)，离心项使用在与势能相同的变量下的二次 Pekeris 展开进行近似。这种变换使得方程可以简化为高斯超几何形式，从而得到振动-转动能量谱的解析表达式。

3. 热力学计算： 利用推导出的能量谱构建分子配分函数，包括振动、转动和平动贡献。

   • 振动配分函数通过对离散束缚态进行求和来计算，直至达到最大量子数 ($n_{max}$)，该量子数由能量等于解离极限的条件确定。
   • 引入了一个与温度无关的配分函数偏移量 ($\eta$)，并针对每个分子根据 NIST-JANAF 吉布斯自由能数据进行校准，以修正模型未明确建模的统计贡献。
   • 从配分函数导出标准摩尔熵 ($S$)、焓增 ($\Delta H$)、吉布斯自由能函数 ($\Phi$) 以及基于压力的平衡常数 ($K_p$)。

4. 验证： 该模型针对几种双原子分子对参考数据进行了测试：

   • 光谱学验证： 将 GqDE 产生的势能曲线和振动能量与光谱重建的转折点数据 (RKR) 及参考振动项值进行对比（涉及 $Br_2$、$Cl_2$、$ICl$ 和 $RbCs$）。
   • 热力学验证： 将计算得到的 $ICl$、$BrCl$ 和 $SO$ 的摩尔性质与 NIST-JANAF 热化学数据在 300–6000 K 温度范围内进行对比。

主要结果

• 光谱学精度： 该模型能够高精度地重现势能曲线和振动能量，尽管性能因分子而异。

  • 势能： 百分比平均绝对偏差 (PAAD) 分别为 1.2266% ($Br_2$)、0.6461% ($Cl_2$)、0.2457% ($ICl$) 和 2.5430% ($RbCs$)。$ICl$ 的符合度最好，$RbCs$ 的偏差最大，主要体现在内侧（短程）分支。
  • 振动能量： PAAD 值分别为 1.4933% ($Br_2$)、1.0259% ($Cl_2$)、0.5743% ($ICl$) 和 1.6687% ($RbCs$)。偏差通常随振动量子数的增加而增大，这表明模型很好地重现了势能最小值的曲率，但在波函数向远离平衡位置延伸的高能区域略有发散。
  • 离心近似： Pekeris 近似能够准确重现 $r_e$ 附近的离心因子，但在较大的核间距离处显示出增加的相对偏差，特别是在 $RbCs$ 中。

• 热力学性质：

  • 在校准了配分函数偏移量后，该模型在熵和吉布斯自由能函数方面与 NIST-JANAF 数据表现出高度一致。平均偏差普遍较低（例如，$BrCl$ 的熵和吉布斯函数的平均偏差分别为 0.0998% 和 0.1443%）。
  • 焓： 焓增在高温度下显示出较大的差异。对于 $SO$，在 6000 K 时，其焓增与参考值的差异为 $-9.027 \text{ kJ mol}^{-1}$。
  • 平衡常数： 对 $ICl$、$BrCl$ 和 $SO$ 生成过程的计算在低温（300–900 K）下显示出显著差异，$\log_{10} K_p$ 的百分比偏差在某些情况下超过了 100%。然而，在高温下（对于 $ICl$ 和 $BrCl$ 高于 3000 K，对于 $SO$ 高于 900 K），一致性显著提高。

意义与主张
本文声称提供了第一个针对双原子分子的广义 q-变形埃卡特势能的完整光谱参数化。其主要贡献在于将势能参数直接与解离能、平衡键长、曲率以及 Dunham 振动-转动关系联系起来，从而实现通过单一解析框架预测振动光谱和热力学性质。

作者强调，本研究成功地将光谱参数化与分子热力学计算联系了起来。然而，作者对模型的预测能力保持了审慎的态度：

• 热力学性质的一致性是在校准了温度无关偏移量之后实现的；因此，结果评估的是能量谱的一致性，而非对绝对配分函数归一化的独立验证。
• 平衡常数在低温下的巨大偏差以及 $SO$ 特有的高温焓误差，突显了当前模型的局限性，特别是针对高激发态的离心近似以及势能在长程处的特定形状。
• 作者总结道，虽然该模型在特定范围内有效地重现了参考数据，但在对反应定义、元素基准态的验证以及针对离散求和的配分函数近似有效性进行更深入的验证之前，尚不能对其热力学预测准确性做出更强的结论。