24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2383  |  回复: 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的回帖

huxiao.gl

木虫 (小有名气)

引用回帖:
3楼: Originally posted by yeahgyz at 2014-07-19 09:01:42
我看黄铜矿文献上取的245,还有一篇取的是280。如果提高截短能,晶胞常数会不会更趋于实验值,误差更小些?别人为什么截短能200多就能收敛,我的不能呢?
...

245即使与实验值吻合,也是巧合,没有意义的。
不用把时间浪费在这个问题上了。
good luck
5楼2014-07-19 16:42:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
截断能至少350~400eV才合理. 300eV以下算出来的结果基本没有意义。

» 本帖已获得的红花(最新10朵)

good luck
2楼2014-07-19 02:48:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahgyz

铜虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by huxiao.gl at 2014-07-19 02:48:21
截断能至少350~400eV才合理. 300eV以下算出来的结果基本没有意义。

我看黄铜矿文献上取的245,还有一篇取的是280。如果提高截短能,晶胞常数会不会更趋于实验值,误差更小些?别人为什么截短能200多就能收敛,我的不能呢?

[ 发自小木虫客户端 ]
3楼2014-07-19 09:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

测试过程建议采用结构单点,一旦确定参数,做结构优化,并在优化结构的单点做参数的检验,这样可以节约时间。最后确定截断能和k不是随便设定一个,而是由测试确定,如果很小的截断能能够满足测试条件,直接使用也无妨(当然我觉得不太可能)。如果测试不能达到标准,或是标准较低,导致的结果就是计算“力”及“应力”会产生偏差,导致结构优化失败。除了这两个参数,空带数目是否充分也很重要。

原则上来说,截断能和k越充分结果是越好的,但是由于DFT的特性,于是不见得是这样。对于某些结构,PBEsol也许比PBE好得多。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-07-19 09:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见