计算化学公社

标题: Seminario方法中键角力常数单位转换求助 [打印本页]

作者
Author: lucky1999 时间: 2024-7-23 08:46
标题: Seminario方法中键角力常数单位转换求助
各位老师，目前我在学习基于量化得到的Hessian矩阵计算用于分子动力学模拟力常数的Seminario方法(Int. J. Quantum Chem.: Quantum Chemistry Symposium, 30, 1271 (1996))的实现方式。请问对于键角力常数，如何将其单位转换为gmx中的kJ/mol/rad^2？下图中是该方法计算键角力常数最后一步对应的公式，观察该公式并没有发现关于角度的量纲，请问如何将该公式中计算得到的k_theta的单位转换为有角度的单位，谢谢各位老师。

作者
Author: sobereva 时间: 2024-7-26 04:30
之前我在写sobtop（http://sobereva.com/soft/Sobtop）程序中的Seminario代码的时候也对其键角力常数的公式产生过疑惑，后来想通了，当时写了以下笔记
(, 下载次数 Times of downloads: 13)

作者
Author: lucky1999 时间: 2024-7-26 08:11

sobereva 发表于 2024-7-26 04:30
之前我在写sobtop（http://sobereva.com/soft/Sobtop）程序中的Seminario代码的时候也对其键角力常数的公式 ...

谢谢sob老师的解答！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)