将氢替换为氘研究反应的同位素效应必须要重新进行freq计算吗

seenbefore

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

UW_0728.

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

beefly

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

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

wal

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

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

wxyhgk

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

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

sobereva

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