| 查看: 2678 | 回复: 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 |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有14人回复
-
为什么国自然不能直接公布
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有3人回复
-
系统查不到
已经有9人回复
-
科研孤儿太难了
已经有22人回复
-
怎么查啊
已经有4人回复
-
国合里面能看到了
已经有7人回复
-
国际合作可查了，中了面上
已经有18人回复
-
能否退出参与的面上项目解除限项
已经有25人回复
-
系统进不去
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算界面结构时原子弛豫的问题
已经有12人回复
- VASP弛豫收敛太慢
已经有18人回复
- 静态计算能收敛，同样的参数弛豫不能收敛，这是为什么？
已经有4人回复
- 做表面弛豫计算用多少时间呢？回帖就有金币哦。。。。
已经有21人回复
- vasp离子弛豫不收敛，求方法
已经有31人回复
- 弛豫时间的计算
已经有3人回复
- 欠松弛因子不同得到的收敛解也不同，且相差很大
已经有6人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助】NVT弛豫问题
已经有3人回复
- 【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算
已经有23人回复
amberkaka
木虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 3723.9
- 帖子: 42
- 在线: 51.8小时
- 虫号: 983640
- 注册: 2010-03-27
- 专业: 凝聚态物性 II ：电子结构
5楼2012-07-10 19:25:40
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
20