自己计算的自由能和文献差别很大

slash-youth · 发表于 Post on 2021-3-16 10:07:33

我想重复文献中的计算过程，文献中△G是10kcal/mol，我算出的是7.86个Hartree，
输入文件中的xyz坐标从Si文件中拷贝的，
请帮忙看一下，哪里出了问题

ms860309 · 发表于 Post on 2021-3-16 11:17:09

首先...
差到快8個hartree實在是太誇張了
我先猜你取值錯誤
多了重複的計算
你可以看看orca output的計算公式
U= E(el) + E(ZPE) + E(vib) + E(rot) + E(trans)
H = U + kB*T (kB is Boltzmann's constant)
T*S = T*(S(el)+S(vib)+S(rot)+S(trans))
G = H - T*S

slash-youth · 发表于 Post on 2021-3-16 12:47:52

ms860309 发表于 2021-3-16 11:17
首先...
差到快8個hartree實在是太誇張了
我先猜你取值錯誤

好的好的谢谢你了

sobereva · 发表于 Post on 2021-3-17 05:53:25

相差到Hartree数量级这种事，要么结构建错了（原子数不对应），要么犯了低级错误（基组不对应、数据读错了、计算表格写错了等），要么公式用错了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 自己计算的自由能和文献差别很大

浏览过的版块