24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

umklapp

银虫 (小有名气)

应助: 9 (幼儿园)
金币: 463.6
散金: 58
红花: 2
帖子: 123
在线: 59.1小时
虫号: 1181441
注册: 2011-01-01
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 带电体系的弛豫计算：难收敛问题

最近在做一个磁性体系引入一个空穴的离子弛豫计算，我试了不同的SMEARING method(gaussian,fermi), 不同的弛豫方法（quasi-newton, CG）, 不同的Hamiltonian-subspace-diagonal( RMM, Davison), 不同的K-mesh, ENCUT, 甚至是不同的初始结构，总是得不到收敛的结果。能量看上去不是一直减小，而是有些震荡；（如下）急求大神指点
19 F= -.26223183E+03 E0= -.26223183E+03  d E =-.121863E-02  mag=    6.0000
  20 F= -.26223168E+03 E0= -.26223168E+03  d E =0.155725E-03  mag=    6.0000
  21 F= -.26223277E+03 E0= -.26223277E+03  d E =-.109502E-02  mag=    6.0000
  22 F= -.26223208E+03 E0= -.26223208E+03  d E =0.689328E-03  mag=    6.0000
  23 F= -.26223209E+03 E0= -.26223209E+03  d E =-.807981E-05  mag=    6.0000
  24 F= -.26223231E+03 E0= -.26223231E+03  d E =-.224072E-03  mag=    6.0000
  25 F= -.26223291E+03 E0= -.26223291E+03  d E =-.595735E-03  mag=    6.0000
  26 F= -.26223369E+03 E0= -.26223369E+03  d E =-.779615E-03  mag=    6.0000
  27 F= -.26223048E+03 E0= -.26223048E+03  d E =0.321320E-02  mag=    6.0000
  28 F= -.26222845E+03 E0= -.26222845E+03  d E =0.202155E-02  mag=    6.0000
  29 F= -.26222844E+03 E0= -.26222844E+03  d E =0.953689E-05  mag=    6.0000
  30 F= -.26222861E+03 E0= -.26222861E+03  d E =-.170270E-03  mag=    6.0000
以下是我的incar 设置
PREC = accurate
ISTART = 0
ICHARG = 2
ISPIN = 2
NELECT = 176 #177 for nuetral
IALGO =48
LREAL = Auto
NBAND =150
EDIFF = 1E-6
ISYM =-1
NELMDL = -5
NELMIN = 8
NELM = 40
LMAXMIX = 4
IBRION = 1 # I also tried 2
ENCUT = 530
NSW = 60
EDIFFG = -0.005
ISMEAR = 0
SIGMA = 0.05
ISIF = 0
VOSKOWN =1
MAGMOM =5 -5 -5 -5 5 -5 5 -5 5 -5 5 -5 18*0
LDAU = .TRUE.
LDAUTYPE = 1
LDAUL = 2 0 0
LDAUU = 4.0 0 0
LDAUJ = 0.58 0 0

回复此楼

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

VASP

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

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

1楼 2012-06-26 02:27:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

amberkaka

木虫 (初入文坛)

应助: 2 (幼儿园)
金币: 3723.9
帖子: 42
在线: 51.8小时
虫号: 983640
注册: 2010-03-27
专业: 凝聚态物性 II ：电子结构

★
xueht987: 金币+1, 鼓励交流！ 2012-07-11 17:28:39

引用回帖:

6楼: Originally posted by umklapp at 2012-07-10 23:59:24
NBANDS 是我自己设的，比valence band多40%；
IALGO 38我也试过了，没有多大的帮助；
奇怪的是
1）这个收敛问题我在bulk state(纯反铁磁) 没有遇到，而是在single-spin-flip
下难收敛。
2）而且如果我选取 ...

抱歉，不知道什么是single-spin-flip。。求科普。
这样看很有可能是初始构型不合理？
可以考虑先进行非磁计算，得出CHG*和WAVECAR，再读入它们进行磁性计算。