求助2，6-NDA的晶格能计算的结果不合理

yangyangyang · 发表于 Post on 2023-5-9 09:26:26

老师好，我想要计算2，6-NDA的晶格能，E晶格=E超晶胞中单分子的平均能量-真空中单分子的能量，超晶胞结构是2*2，单胞结构是从COD数据库直接下的，没经过结构优化直接用castep计算的能量，真空中分子能量是用dmol3优化后计算的，但所得结果大的不合理，请教老师们，学生在哪一步是不合理或有问题的？

NDA.outmol (16.75 KB, 下载次数 Times of downloads: 0)

NDA-3.castep (80.31 KB, 下载次数 Times of downloads: 0)

sobereva · 发表于 Post on 2023-5-9 23:36:13

恰当设置k点
必须考虑色散校正

yangyangyang · 发表于 Post on 2023-5-15 22:19:53

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，学生想问我的更改是不是还有问题？

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

[CASTEP/Dmol3/MS] 求助2，6-NDA的晶格能计算的结果不合理