计算化学公社

标题: 自己计算的自由能和文献差别很大 [打印本页]

作者
Author: slash-youth 时间: 2021-3-16 10:07
标题: 自己计算的自由能和文献差别很大
我想重复文献中的计算过程，文献中△G是10kcal/mol，我算出的是7.86个Hartree，
输入文件中的xyz坐标从Si文件中拷贝的，
请帮忙看一下，哪里出了问题

作者
Author: ms860309 时间: 2021-3-16 11:17
首先...
差到快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

作者
Author: slash-youth 时间: 2021-3-16 12:47

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

好的好的谢谢你了

作者
Author: sobereva 时间: 2021-3-17 05:53
相差到Hartree数量级这种事，要么结构建错了（原子数不对应），要么犯了低级错误（基组不对应、数据读错了、计算表格写错了等），要么公式用错了

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