论文标题：PyRAMD Scheme: A Protocol for Computing the Infrared Spectra of Polyatomic Molecules Using ab Initio Molecular Dynamics

论文链接：https://www.mdpi.com/2813-446X/2/3/12

期刊名：Spectroscopy Journal

期刊主页：https://www.mdpi.com/journal/spectroscj

研究背景

分子光谱是探索物质结构、揭示微观世界运动规律的核心工具之一。其中，红外光谱在化学、材料科学、环境监测和生命科学等领域具有广泛应用。理论模拟是解读和预测红外光谱的重要手段，特别是从头算分子动力学方法，能够处理复杂体系并包含非谐性效应。然而，该方法在追求高精度时往往伴随着巨大的计算成本，包括对核量子效应的准确描述、为获得足够频率分辨率所需的长时程模拟，以及为确保数值稳定而设置的微小积分步长。因此，如何在保证模拟可靠性的前提下，发展高效的计算方案来降低这些成本，是计算光谱学领域一个亟待解决的实际问题。本文提出并验证了一个完整的、基于PyRAMD软件的“方案”，旨在系统性地解决上述问题，为复杂分子的红外光谱计算提供一个平衡精度与效率的实用工具。

研究内容

该研究方案的核心在于整合了三个关键组成部分，旨在系统性地优化红外光谱的从头算分子动力学模拟流程。首先，为在经典分子动力学中有效且低耗地考虑核量子效应，研究引入了简化的维格纳采样方法。传统维格纳采样依赖于参考平衡结构，限制了其在模拟过程中的应用，而SWS方法通过对坐标和动量引入符合一定分布的位移，可以模拟零点振动效应，并且与热浴兼容。研究还基于分子的谐波振动频率平均值，推导并验证了一种更高效、无需多次试算的参数选择准则，从而降低了确定SWS关键参数的计算开销。其次，为了在保证精度的同时提高计算效率，作者明确提出了通过采用大积分步长来降低模拟成本的可能性。这通常会导致高频振动峰的数值蓝移，但研究通过分析推导给出了针对Verlet积分等常用算法的频率校正公式，可将模拟得到的频率修正回真实值。这有效扩展了可行模拟步长的上限。最后，在频谱后处理阶段，研究提出了两套提升光谱质量的方案。一方面，针对传统快速傅里叶变换频率分辨率受限于总模拟时长的局限性，引入了正则化最小二乘频谱分析算法，该算法能在无需显著延长模拟时间的前提下，获得具有更高表观频率分辨率的平滑光谱。另一方面，借鉴计算谐波频率的成功经验，研究建议对MD计算得到的光谱也应用尺度因子来校正系统误差。作者基于一组小分子训练集，给出了BLYP-D3(BJ)/6-31G、PBE-D3(BJ)/6-31G和PBEh-3c这三个常用泛函水平的预参数化尺度因子，为后续研究者直接应用提供了便利。

图 4. 在 PBEh-3c 理论水平下计算得到的甲烷（CH4）、乙烷（C2H6）、甲胺（CH3NH2）和甲醇（CH3OH）的振动光谱。

该PyRAMD方案的有效性在多个分子体系中得到了验证。研究表明，结合上述三个部分的协同作用，可以显著改善模拟光谱的质量。例如，在计算质子化甲烷CH5+的红外光谱时，未经校正的MD光谱与实验结果存在差距。但在应用了频率校正、强度校正（针对特定实验方法）以及预计算尺度因子后，理论谱图与实验谱图的吻合度得到大幅提升。此外，通过对乙酸和茚的案例分析，研究对比了经过尺度因子校正的谐波频率光谱、MD光谱与实验光谱。结果显示，尽管未校正的MD光谱因其内在包含非谐性效应，已比未校正的谐波光谱更接近实验；而经各自的尺度因子校正后，两者在谱峰位置上都与实验更为接近，展现了该方案在提供理论光谱参考方面的实用性。尽管目前提出的MD尺度因子训练集较小，但其思路为后续发展更精确的校正参数指明了方向。

研究总结

总之，本文系统性地提出了一个名为“PyRAMD方案”的计算框架，用于高效模拟多原子分子的红外光谱。该方案整合了基于SWS的核量子效应处理方法、适用于大积分步长的频率校正技术、以及可提升频谱质量的后处理策略，在精度与效率之间实现了良好的平衡。虽然从理论严格性上，此方案并非最优，例如在考虑核量子效应方面不及路径积分分子动力学，但它提供了一个在经典MD高效率和路径积分MD高精度之间的有效折衷方案。研究表明，该方案能够为CH5+、乙酸、茚等体系生成与实验数据吻合良好的理论光谱。因此，该方法有望成为连接实验测量与理论计算的可靠桥梁，为研究者提供一个在常规计算资源下即可获得高质量红外光谱预测结果的实用工具，在物理化学、催化、光谱学等多个研究领域具有广泛的应用潜力。

Spectroscopy Journal 期刊介绍

主编： Prof. Dr. Clemens Burda, Department of Chemistry, Case Western Reserve University, Millis Science Center, 10900 Euclid Ave., Cleveland, OH 44106, USA

我们鼓励科学家在光谱技术的各个方面、特性表征、理论以及其他光谱发展领域发表他们的实验和理论研究成果。光谱学涉及物质与电磁频谱任何部分之间的相互作用，并应用于所有学科，包括物理学、化学、生物化学、生物学、空间科学、材料科学和工程学等领域。来自非光子实验（如电子、中子和质子实验）的贡献同样受欢迎。