24小时热门版块排行榜    

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

yeahgyz

铜虫 (初入文坛)

[求助] 黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决已有1人参与

小弟按照文献所述参数对黄铜矿进行了几何结构的优化,但是却一直不能收敛是怎么回事,而且算得的晶胞常数与实验值相差超过了2%,甚至达到了4%或5%,是什么地方出现了问题吗?下面吧具体模拟过程及文献提供给各位大神分析,谢谢!

文献中的设置如下:黄铜矿,理想黄铜矿的晶胞参数为:a=b=5.289 Å,c=10.423 Å,α=β=γ=90°。采用GG-PBE梯度修正函数,在进行几何优化时,选取合适的截断能和k点,收敛标准为:原子位移为0.005 Å,原子间作用力为0.01 eV/nm,原子间的内应力为0.2 GPa,体系总能量的变化不大于5.0×10-4 eV/atom。
截断能和k点选择过程及数据如图一和图二:

我按照245计算,计算到98步进行不下去了,受不了就停止了。
手头上有一个270ev的数据,收敛过程如图三图四,设置参数如图五图六。

另外我算了E点为310,350和400的,k点位5x5x2,很快收敛,但是优化后的晶胞常数误差超过2%,最大的达到4%,应该是不能用的。
请问1.为什么我按照文献设置不能收敛?
       2.为什么提高E点可以收敛,但是晶胞常数误差过大?是不是K点设置问题,还是 其他问题?

谢谢指教!

黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决
360截图20140718223046593.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-1
360截图20140718223114573.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-2
360截图20140718223437784.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-3
360截图20140718223452372.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-4
360截图20140718223536488.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-5
360截图20140718223610451.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahgyz

铜虫 (初入文坛)

送红花一朵

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

huxiao.gl

木虫 (小有名气)

信息提示
请填处理意见