gch123 发表于 2021-8-16 10:27 是的 |
请问在这个近似的公式里,用vasp计算过振动能的计算之后,怎么从计算的结果中得到zpe和ts的值呢? |
您好,请问一下,这个公式算出来的自由能是标准吉布斯自由能吗? |
wzkchem5 发表于 2021-1-5 09:40 好的,谢谢 |
dali 发表于 2021-1-5 11:50 有的,你看一下那两项的表达式,把那两项的值对振动频率作图,就知道了。 如果你非要问是不是有现成的教材或者文献得出过这两项很小的结论,我也不知道 |
wzkchem5 发表于 2021-1-5 09:40 你好,请问忽略掉的项随反应的变化都很小是否有理论依据? |
笙歌落 发表于 2021-1-4 22:40 这个问题是这样的:忽略掉的项随反应的变化都很小,所以忽略掉问题也不大;但是一旦有声子频率了以后,计算这几项特别快,所以即便这些项特别小,也完全没有理由不算它们。也就是说,假如你现在要做计算,那么应该考虑这些项,但是如果你之前做过计算而且没考虑这些项,文章已经发出去了,那也没必要特意重新算一遍然后发个erratum。 |
dingniu2 发表于 2021-1-2 10:41 文献有错误,还是提出这个公式的人理解错了?记得有人说是Jens K. Nørskov提出的吧 |
许多vasp计算文献都根据文献J. Phys. Chem. B 2004, 108, 17886-17892,但按照G=Eele+ZPE-TS,相当于G(T)=U(0)-TS,如果这里是U(T)还好理解,就是凝聚态ΔG≈ΔH,但这里不是U(T)而是U(0),应还差一个∫CpdT,不知道这个∫CpdT是如何忽略的,忽略掉误差有多大? |
碧雪霜花剑 发表于 2021-1-1 11:34 分清楚U(0)和U(T)的差异 能算U(T)干嘛要用U(0) |
碧雪霜花剑 发表于 2021-1-1 11:34 电子的能量是通过自洽场计算得到的,ZPE和TS以及U(T)对声子频率进行分析就能得到。Phonopy官网是有公式的,如果还不明确或者进一步了解可以专门看统计力学方面的书籍,通常在理想晶体部分会提到。气体比固体多出的是平动和转动对S和U以及pV项(固体没有)的贡献,有些固体计算的程序并不会考虑这部分,对于吸附前的分子在处理过程要小心。 |
算固体的自由能,一般用简谐近似,也就是算声子 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-25 20:50 , Processed in 0.400296 second(s), 26 queries , Gzip On.