计算化学公社

标题: 求助2，6-NDA的晶格能计算的结果不合理 [打印本页]

作者
Author: yangyangyang 时间: 2023-5-9 09:26
标题: 求助2，6-NDA的晶格能计算的结果不合理
老师好，我想要计算2，6-NDA的晶格能，E晶格=E超晶胞中单分子的平均能量-真空中单分子的能量，超晶胞结构是2*2，单胞结构是从COD数据库直接下的，没经过结构优化直接用castep计算的能量，真空中分子能量是用dmol3优化后计算的，但所得结果大的不合理，请教老师们，学生在哪一步是不合理或有问题的？ (, 下载次数 Times of downloads: 23) (, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 18) (, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 15)

作者
Author: sobereva 时间: 2023-5-9 23:36
恰当设置k点
必须考虑色散校正

作者
Author: yangyangyang 时间: 2023-5-15 22:19

sobereva 发表于 2023-5-9 23:36
恰当设置k点
必须考虑色散校正

谢谢sob老师的解答，抱歉我前几天去干别的事了，没有及时回复您，我重新设置了一下，结果还是要比预估的大的离谱，（电脑传图片和文件失败了，抱歉），
我在dmol3里将Functional更改为GGA-PBE,点选了Grimme计算DFT-D矫正，最后读取 Ef -762.434411Ha -4.9567064Ha 6.76E-06 0.4m 14
在castep里同样将Functional更改为GGA-PBE,点选了Grimme计算DFT-D矫正，此外在electronic/more/k-points/separation从0.08更改到0.04，0.05与0.04所得能量得差值小于0.01ev，在0.05读取 Dispersion corrected final energy* = -14816.48811528 eV，在0.04读取 Dispersion corrected final energy* = -14816.48605583 eV
计算晶格能为-14816.48605583/4-（-726.434411*27.2114）=16063.1758175 ev，学生想问我的更改是不是还有问题？

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