| 查看: 2098 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
带电体系的弛豫计算： 难收敛问题
|
|||
|
最近在做一个磁性体系引入一个空穴的离子弛豫计算， 我试了不同的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 |
» 猜你喜欢
回收溶剂求助
已经有7人回复
-
职称评审没过，求安慰
已经有40人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- VASP弛豫收敛太慢
已经有18人回复
- 静态计算能收敛，同样的参数弛豫不能收敛，这是为什么？
已经有4人回复
- 做表面弛豫计算用多少时间呢？回帖就有金币哦。。。。
已经有21人回复
- vasp离子弛豫不收敛，求方法
已经有31人回复
- 弛豫时间的计算
已经有3人回复
- 欠松弛因子不同得到的收敛解也不同，且相差很大
已经有6人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助】NVT弛豫问题
已经有3人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3679.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 915.2小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
9楼2012-07-12 10:07:10
enola
捐助贵宾 (正式写手)
- 1ST强帖: 2
- 应助: 13 (小学生)
- 金币: 1485.9
- 散金: 1549
- 红花: 14
- 帖子: 847
- 在线: 102.5小时
- 虫号: 947564
- 注册: 2010-01-24
- 专业: 半导体光电子器件
2楼2012-06-26 11:34:18
3楼2012-06-26 22:59:36
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2012-06-27 11:16:46