24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1107  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

pgi1706

木虫 (小有名气)

[交流] 晶格常数和能量是这种形式?已有2人参与

晕了,算出来的晶格常数和能量是这种形式?
a= 10.891 1 F= -.22064566E+03 E0= -.22064566E+03 d E =-.220646E+03
a= 10.892 1 F= -.22064519E+03 E0= -.22064519E+03 d E =-.220645E+03
a= 10.893 1 F= -.22064626E+03 E0= -.22064626E+03 d E =-.220646E+03
a= 10.894 1 F= -.22064570E+03 E0= -.22064570E+03 d E =-.220646E+03
a= 10.895 1 F= -.22064509E+03 E0= -.22064509E+03 d E =-.220645E+03
a= 10.896 1 F= -.22064517E+03 E0= -.22064517E+03 d E =-.220645E+03
a= 10.897 1 F= -.22064448E+03 E0= -.22064448E+03 d E =-.220644E+03
a= 10.898 1 F= -.22064530E+03 E0= -.22064530E+03 d E =-.220645E+03
a= 10.899 1 F= -.22064452E+03 E0= -.22064452E+03 d E =-.220645E+03
a= 10.900 1 F= -.22064591E+03 E0= -.22064591E+03 d E =-.220646E+03
所谓能量此起彼伏,何为最优结构?
这是采用4*4*4K点优化的结果,我想得到比较精确的能量。

原来采用单G点优化时,a=10.894时是最低能量,但是采用4*4*4后,能量却是现在这个模样?
可有人遇到这种情况,能否交流一下?谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
贺仪(金币+5): 多谢指教! 2011-09-09 17:03:41
说两点:
1,k点可能不够,要做收敛测试,也就是不断的增加k点数目进行晶格常数的计算,得到一系列的晶格常数值,当k点足够多时,直到得到的晶格常数趋近稳定,这个时候的k数目才是足够的。同样,也应该对Ecut做类似的测试。
2,也是更重要的,要对自己计算精度有一定的认识,这也是和上面说的收敛测试密切相关的。一般来说,能量收敛到meV的量级,晶格常数收敛到0.01埃以下,就已经是比较不错的精度了。从k点数目来看,你使用的参数并非精度极高的计算吧,如果a是原子单位的话,你所追究的这些0.001的差别可能已经超出了你的计算参数所能达到的精度极限了,能量也是同理。因此你所说的这些”能量此起彼伏“可能都是在误差范围之内的。
3楼2011-09-09 16:19:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pgi1706 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见