计算化学公社

标题: 将氢替换为氘研究反应的同位素效应必须要重新进行freq计算吗 [打印本页]

作者
Author: seenbefore 时间: 5 day ago
标题: 将氢替换为氘研究反应的同位素效应必须要重新进行freq计算吗
各位老师好，现在我有一个氢转移的过渡态的频率计算文件，我希望将氢替换为氘研究反应的同位素效应，我必须要重新进行freq计算吗？是否有降低耗时的方法，比如读取fchk等？因为体系比较庞大。

作者
Author: UW_0728. 时间: 5 day ago
看http://sobereva.com/423
热力学数据计算结果和同位素设定密切相关，因此需要重新进行freq计算
可以读取之前同级别任务的fchk当波函数初猜节省SCF阶段时间，因为根据BO近似原子核没有波函数，不同同位素对应的波函数应该没啥差别；后续的Hessian矩阵得重新计算，这个时间节省不了

作者
Author: beefly 时间: 5 day ago
本帖最后由 beefly 于 2026-1-18 22:29 编辑

UW_0728. 发表于 2026-1-18 21:54
看http://sobereva.com/423
热力学数据计算结果和同位素设定密切相关，因此需要重新进行freq计算
可以读 ...

BO近似下的Hessian矩阵没有同位素效应，只有在对角化算频率的时候才代入不同同位素的原子量。频率要重算，但是Hessian矩阵不必重算，尤其是分子较大计算比较耗时的时候

作者
Author: wal 时间: 5 day ago

beefly 发表于 2026-1-18 22:27
BO近似下的Hessian矩阵没有同位素效应，只有在对角化算频率的时候才代入不同同位素的原子量。频率要重算 ...

所以这种情况是可以freq=readfc的对嘛

作者
Author: wxyhgk 时间: 4 day ago

wal 发表于 2026-1-18 23:13
所以这种情况是可以freq=readfc的对嘛

是的，我在公司经常算氘代的

作者
Author: sobereva 时间: 4 day ago
用guess=read freq=readfc结合新的同位素设定重新过一遍就完了，非常快。readfc读的是Hessian，不同同位素下Hessian、势能面极小点、电子波函数都是相同的

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