24小时热门版块排行榜    

查看: 2149  |  回复: 0

bjwang

铁杆木虫 (正式写手)

[求助] 二维体系加电场的ISIF取2还是3???

尝试计算电场对二维体系电子结构等的影响,借鉴了小木虫上的一些帖子(比如:http://muchong.com/bbs/viewthread.php?tid=8223517&fpage=1:)
主要设置下面三个参数:
  EFIELD = 0.05
  IDIPOL = 3   #z direction
  LDIPOL = .TRUE.
我的计算步骤如下:
(1)先不加电场进行结构优化,此时ISIF=3
(2)利用上面的结构,加上电场(在INCAR里加入上面三个参数),此时ISIF=3时,优化出错(ZBRENT: fatal error in bracketing      please rerun with smaller EDIFF, or copy CONTCAR   to POSCAR and continue),然而,当我把ISIF取为2时,结果很快就收敛了!
我的问题是:
(1)对于二维体系,ISIF是不是只能取2?
(注:大神valenhou001在回复关于“金属薄膜体系”加电场的一个帖子————http://muchong.com/bbs/viewthread.php?tid=8223517&fpage=1时说”对于表面体系的超原胞的计算,不要采用ISIF=3. 正确或合理的步骤:先优化和计算好块体材料的结构参数(晶格常数,原子坐标),再切表面构造出表面体系的超原胞,对这样得到的表面体系的超原胞进行只优化原子位置(也就是ISIF=2),平行于表面的晶格常数固定,表面法向方向上的晶格常数根据真空层厚度以及原子层厚度来取。“),

侯老师上面说的”表面体系的超原胞“和二维体系的意思类似吗?也就是说对于二维体系,ISIF也只能取2???

(2)如果(1)不合适,ISIF=3出错,应该是我的结构在cell改变时,结构发生不合理的变化造成的,如何解决?

(3)在计算中,采用ISIF=2ISIF=3在考虑电场的影响上来说,差别大不?

  谢谢各位大神!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

It is a shame for the soul to be first to give way in this life,when thy body does not give way.
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bjwang 的主题更新
信息提示
请填处理意见