使用广义分数阶导数改进双原子分子振动能量的建模

为何分子微小的振动至关重要

你每一次呼吸、每一簇燃烧的火焰以及太空中的每一朵云中，都充满了不断振动的简单双原子分子。这些微观振动决定了分子如何吸收光、储存能量并参与化学反应。为理解和预测这些行为，科学家建立了描述分子运动的数学模型。本文介绍了一种更精确地描述这些振动的新方法，它将对经典而强大的分子模型应用一种现代化的微积分变体。

科学家通常如何描绘振动的分子

双原子分子——由仅两个原子组成——可以想象成由一根弹簧连接的两个球体。当原子间距离变近或变远时，它们在能量势阱上上下移动。近一个世纪以来，一种称为 Morse 势的简单公式一直是描述该能量景观的主力。它能捕捉到化学键的伸长、变硬以及最终断裂的行为。将其与量子力学中的薛定谔方程结合，Morse 势使研究者能够计算允许的振动能级，这些能级以分子光谱中的锐线形式出现在实验室测量或来自遥远恒星的观测中。

经典描绘的不足之处

标准的 Morse 模型在分子不自旋且激发不太高时效果最好。然而真实分子在振动的同时也会旋转，高能级下的行为更加复杂。为处理旋转效应，通常要加入一个额外的能量势——称为离心项——但这会使方程难以精确求解。多年来，人们发展了许多巧妙的近似和数值技巧，但在理论与精确实验数据之间仍存在小的偏差，尤其是在较高振动态和具有不同电子态的分子中更为明显。

对分子运动的新分数阶处理

作者采用了来自“分数阶”微积分的最新思想，将普通导数推广为非整数阶。在这一框架下，系统的响应可以带有对过去的记忆，或等效地编码微妙的非局域相互作用。基于一种保留了标准微积分关键规则的广义分数阶导数，研究团队将一种著名的解析解法（Nikiforov–Uvarov 方法）扩展为他们所称的广义分数阶 NU 方法。他们利用该方法从 Morse 势出发，并对旋转势采用精炼近似，推导出任意空间维度下分子的振动能级的解析公式。

将新模型付诸检验

为检验分数阶方法的有效性，作者将其应用于二十二种不同的双原子分子，包括在化学、材料科学和天体物理学中重要的物种，如 CO、Na₂、AlH、SiO⁺、TaO 和 TaS。对每一种情形，他们根据已知的分子常数构建势能曲线，并将预测的振动能级与来自 Rydberg–Klein–Rees（RKR）分析的高质量基准数据进行比较。他们将分数阶阶数视为一个可微调的小参数，用以有效概括简单模型未捕捉到的复杂相互作用。通过对每个分子系统地最小化平均百分比误差，结果显示分数阶公式以显著精度重现观测到的振动能量，通常相比经典极限将差异缩小两倍或更多，并且在与最先进的替代势和数值方法竞争时表现良好。

为什么更小的误差意义重大

除了平均误差度量外，研究还检视了理论与实验在各振动能级上的差异随级数的变化。对于许多分子，传统模型随着振动量子数的增大而逐渐偏离实验数据。相反，分数阶模型将这些偏差保持得更小且更均匀，尤其是在光谱上部——对准确预测要求最高的区域。作者还将其解析解用于 CO 分子的旋转态，并与高度精确的数值基准进行比较，发现良好的一致性，表明其近似在相当大的旋转运动下仍然可靠。

这对理解分子意味着什么

通俗地说，这项工作表明，通过“温和地”让微积分的常规规则发生弯曲——允许导数为非整数阶——可以得到对简单分子振动和旋转更灵活、更忠实的描述。通过调整单一的分数阶参数，该模型能够捕捉先前需要更复杂势或繁重数值计算才能表达的微妙效应。这使得它成为解释实验室分子光谱、设计新材料以及探测富含双原子分子的天体环境成分的一个高效且强大的工具。

引用: Khokha, E.M., Abu-Shady, M., Omugbe, E. et al. Improved modelling for vibrational energies of diatomic molecules using the generalized fractional derivative. Sci Rep 16, 12037 (2026). https://doi.org/10.1038/s41598-026-39091-5

关键词: 双原子分子, 振动光谱, Morse 势, 分数阶微积分, 分子光谱学