24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2990  |  回复: 23

liang7035

金虫 (正式写手)

[交流] 【求助】镍110表面优化失败的原因在那?已有13人参与

我优化ni110 表面,优化3天,结果失败,请高手给分析一下原因!谢谢。。



晶胞切割110表面,depth:2.5 建了5层,真空层厚度11埃,然后优化计算
*Warning* max. SCF cycles performed but system has not reached the groundstate.
  
Current total energy, E     =  -27108.96036451     eV
Current free energy (E-TS)  =  -27109.09020830     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS) =  -27109.02528641     eV

  ****************************************************************************
  Warning: electronic minimisation did not converge when finding ground state.
  ****************************************************************************

Writing model to 3D_Atomistic__1_1_0_.check                                                                                                                                                                                                                                      
Error in geom_get_forces - electronic_minimisation of current_cell failed
[0] MPI Abort by user Aborting program !
[0] Aborting program!

[ Last edited by liang7035 on 2010-4-17 at 14:37 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一原理资料汇编 我的饕餮 ms计算模型及其方法

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一个一心向上的人整个世界都会给他让路!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

acridine

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
xiaohunhun(金币+10):非常好,非常仔细的解释,受教了,呵呵,这样说的话还是要计算对比一下 2010-04-20 09:07
引用回帖:
Originally posted by xiaohunhun at 2010-04-17 19:27:05:

请您解释一下这段话,谢谢
我觉得还是要优化,如果不优化,我怎么知道该界面能不能稳定存在
话说楼主连导入的模型都优化了。
Z轴变形,我觉得也是可以理解的,要是不是低能态,晶格畸变是正常的吧??

我的理解是这样,楼主处理的是slab,应该是用来模拟表面体系吧。模型在z方向加足够厚的真空,以忽略层间的作用,因此在优化的时候,至少在这个方向的晶格常数的优化是没有意义的。
至于x, y方向的优化,也是需要小心的。固然,表面的原子配位状态的改变会导致原子间结合的改变,进而导致晶格常数改变,但这种改变的幅度是否可以用优化slab的x,y方向来描述,这是个值得注意的问题。因为实际材料不会像slab模型那么薄,因此表面原子的影响力比较小,而体系晶格常数还是取决于体材的。而理论模拟的slab模型则不是这样,表面原子占很大比例,因此有可能对晶格常数造成较大影响,但这种情况下优化出来的晶格常数的变化,究竟是实际存在的还是由于slab模型过薄而引起的,是值得注意的。这是一个理论模型与实际模型之间的权衡
而“将晶格常数设置为体材值”,用于模拟实际的表面,还是有一定合理性的,就如同“slab中心的原子是要固定的”,也无非是为了要模拟“这里已经是体材区域”。既然中心的原子的固定是个很好的假设,那我觉得晶格常数设置为体材值也是个不错的近似
22楼2010-04-20 01:32:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★
liang7035(金币+5):谢谢参与
gavinliu7390(金币+1):谢谢交流! 2010-04-15 20:51
这是原因:
electronic minimisation did not converge when finding ground state.
寻找基态时,电子减小不收敛。
估计是模型建的有点问题~能把模型建立过程详细说说么?
思想重于技巧,内涵重于表象
4楼2010-04-15 18:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+3):谢谢交流! 2010-04-17 10:22
楼主的问题解决了没有,上一个回贴过于简单了,楼主还奖励金币, 惭愧啊惭愧~
首先说max. SCF cycles performed but system has not reached the groundstate.这条,楼上说的都有道理,但我想说的是从你的前几步优化都没有问题,表明各种参数的选择上应该是可以的。可能是scf的步数设置比较小,你可以看以下你的最后一步scf过程,在迭代的过程中能量是否有收敛的趋势。如果收敛趋势比较好,那就表明在增加几步scf次数还是可以收敛的,那么就不用改变其他参数了,增加max. SCF cycles就可以
这是计算停止的直接原因,另外就是我在上面的回贴说的,既然设置了真空层就不要优化晶格了,你可以看一下,在优化后模型的z轴应该已经改变了,这个方向上的改变是没有意义的,不仅浪费计算时间,而且会导致格子的形状发生改变,也有可能导致不收敛
12楼2010-04-16 23:41:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shrek826

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+2):谢谢解答! 2010-04-17 14:39
Ni优化不成功是很正常的,有两条建议:
第一,增加empty_bands
第二,增加smearing
第二条尤其重要
15楼2010-04-17 14:15:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzfish

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
不知道现在楼主的问题解决没有。
Ni应该有磁性的吧,不知道楼主是如何设置Ni的电子结构的?
也就是SPIN怎么设置,这对有磁性体系计算很重要的。
16楼2010-04-17 15:30:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

04nylxb

木虫 (正式写手)

★ ★
liang7035(金币+5):谢谢参与
gavinliu7390(金币+1):谢谢交流! 2010-04-15 20:50
把你构建的Ni(110)模型放上来,几层、supercell多少、束缚……
便于大家参考
集中精力发文章
2楼2010-04-15 17:29:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
liang7035(金币+5):谢谢参与
xiaohunhun(金币+2):请您解释仔细一点好么?说明一下原因,谢谢 2010-04-17 19:19
表面,slab模型,还要优化晶格,这是件值得商榷的事情
3楼2010-04-15 17:44:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)

建模过程

引用回帖:
Originally posted by souledge at 2010-04-15 18:37:29:
这是原因:
electronic minimisation did not converge when finding ground state.
寻找基态时,电子减小不收敛。
估计是模型建的有点问题~能把模型建立过程详细说说么?

谢谢你的回复,我先导入镍的晶胞对其优化,优化成功以后,再拿优化的晶胞切割110表面,depth:2.5 建了5层,真空层厚度11埃,然后优化计算。
一个一心向上的人整个世界都会给他让路!
5楼2010-04-15 19:00:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chensxxy

铁杆木虫 (正式写手)

★ ★
liang7035(金币+5):谢谢参与
gavinliu7390(金币+1):谢谢参与! 2010-04-15 20:51
估计有2个方面的原因:
1、模型存在不合理的地方
2、计算机内存不够或者其他配置不合理

你用的什么泛函???? 或者是 empty band 设置

[ Last edited by chensxxy on 2010-4-15 at 20:04 ]
尊崇自然,敬畏生命!
6楼2010-04-15 19:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)

内存应该还可以啊!empty band 我更改了用的是21
模型不合理的地方在哪里呢!
一个一心向上的人整个世界都会给他让路!
7楼2010-04-15 20:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzfish

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gavinliu7390(金币+1):谢谢参与! 2010-04-15 20:51
好象只驰豫原子坐标,不驰豫晶包参数。
8楼2010-04-15 20:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzfish

金虫 (正式写手)


gavinliu7390(金币+1):谢谢交流! 2010-04-16 09:53
引用回帖:
Originally posted by liang7035 at 2010-04-15 20:08:56:
内存应该还可以啊!empty band 我更改了用的是21
模型不合理的地方在哪里呢!

empty band 都增加到30%,可以多试试看。看看华南理工大学的教程吧。
9楼2010-04-16 09:23:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xjwang8621

金虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
华南理工大学的教程,在那儿下?
引用回帖:
Originally posted by hzfish at 2010-04-16 09:23:58:



empty band 都增加到30%,可以多试试看。看看华南理工大学的教程吧。

10楼2010-04-16 23:13:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liang7035 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见