24小时热门版块排行榜    

查看: 3014  |  回复: 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 ]
回复此楼
一个一心向上的人整个世界都会给他让路!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖
查看全部 24 个回答

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的回帖

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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见