24小时热门版块排行榜    

查看: 5263  |  回复: 6

dx0620

木虫 (正式写手)

[求助] vasp 考虑自旋轨道耦合计算能带的过程? 已有2人参与

现在算的体系需要考虑自旋轨道耦合(SOC)对能带、态密度的影响,那请问
1)具体的步骤及INCAR、KPOINTS的设定是如何的?
2)下面是我认为的顺序:
      1、正常结构优化
      2、正常自洽计算得到电荷密度(1、2两步与一般计算能带过程一样吧?!)
      3、加上SOC计算能带、态密度
3)每一步(尤其是自洽和计算能带两步)的KPOINTS如何选取?

我的计算经常出现如下问题:
====================================================
magnetization density of overlapping atoms calculated
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
  WARNING: DENMP: can't reach specified precision
  Number of Electrons is NELECT =   44.9843750000000
DAV:   1    -0.607409496817+115   -0.60741+115   -0.41430+115  4680   0.151+114
  WARNING: DENMP: can't reach specified precision
  Number of Electrons is NELECT =   45.0312500000000
DAV:   2    -0.762479432673+116   -0.70174+116   -0.75190+116  7368   0.225+113
WARNING: Sub-Space-Matrix is not hermitian in DAV            6
WARNING: Sub-Space-Matrix is not hermitian in DAV            6
  9.401988633190671E+099
WARNING: Sub-Space-Matrix is not hermitian in DAV            6
  9.401988633190671E+099
WARNING: Sub-Space-Matrix is not hermitian in DAV            6
====================================================
回复此楼

» 猜你喜欢

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

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

whywhy01

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
dx0620(ljw4010代发): 金币+2, 谢谢指导解答! 2014-11-17 18:54:56
第2步的时候应该设定NBANDS和LMAXMIX,我看后面有NELECT的提示,我猜可能是NBANDS的问题。第2步的NBANDS应该是第1步的两倍
2楼2014-11-17 17:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
2楼: Originally posted by whywhy01 at 2014-11-17 17:17:16
第2步的时候应该设定NBANDS和LMAXMIX,我看后面有NELECT的提示,我猜可能是NBANDS的问题。第2步的NBANDS应该是第1步的两倍

手册上写的不是第3步是第2步的两倍吗?
第2步设定LMAXMIX了,这个就是一个数吧?如果两种元素如何设置?
3楼2014-11-17 17:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

引用回帖:
3楼: Originally posted by dx0620 at 2014-11-17 17:27:00
手册上写的不是第3步是第2步的两倍吗?
第2步设定LMAXMIX了,这个就是一个数吧?如果两种元素如何设置?...

对的我记错了....选LMAXMIX比较大的那个数
4楼2014-11-17 19:26:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zh11111

银虫 (正式写手)

---

我在进行第三步soc计算时INCAR为SYSTEM
PREC=High
ENCUT=520
EDIFF=1e-5
ISIF=3
ISYM=0
NSW=0
ISMEAR=-5
ISTART=1
ICHARG=11
LWAVE=TRUE
LCHARG=TRUE
NPAR=4
EDIFFG=-1e-2
ISPIN=2
LORBIT=11
LMAXMIX=4
LSORBIT=T
LORBMOM=T
SAXIS=1 0 0
NBANDS=192
NGXF=120
NGYF=54
NGZF=126
结果显示错误:LDA part: xc-table for Ceperly-Alder, standard interpolation
found WAVECAR, reading the header
  number of k-points has changed, file:    39 present:   125
  trying to continue reading WAVECAR, but it might fail
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(           1 )
reading WAVECAR
reading wavefunctions of collinear run, up
ERROR: while reading WAVECAR, plane wave coefficients changed       15319
       15285
这个问题您知道如何改进吗?
---
5楼2015-05-26 07:37:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jlmxtyy

木虫 (正式写手)

引用回帖:
5楼: Originally posted by zh11111 at 2015-05-26 07:37:39
我在进行第三步soc计算时INCAR为SYSTEM
PREC=High
ENCUT=520
EDIFF=1e-5
ISIF=3
ISYM=0
NSW=0
ISMEAR=-5
ISTART=1
ICHARG=11
LWAVE=TRUE
LCHARG=TRUE
NPAR=4
EDIFFG=-1e-2
ISPIN=2
LORBIT=11
LMAXM ...

ISIF=3?你这个是要优化还是干啥啊?
另外算能带的时候,可以不用之前的WAVECAR,只用CHGCAR,就不会有错误提示
6楼2016-06-13 16:56:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wobure

新虫 (小有名气)

【答案】应助回帖

不想麻烦的话,第二步就开自旋轨道耦合计算,K点和第一步一样,第三步用第二步的Charge计算能带,k点line mode,自定义都可以。 读入标准dft的共线计算,个人认为只是为了增加收敛速度而已。个人经验,仅供参考!
7楼2016-06-15 13:49:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dx0620 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:50 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:23 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:22 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:22 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:20 by Kddn1pB6pAzX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:17 by Kddn1pB6pAzX
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:11 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:02 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 10/500 2026-08-01 17:51 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:51 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:42 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:40 by Kddn1pB6pAzX
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +16 低垂的野花 2026-07-31 22/1100 2026-08-01 16:39 by Cathyau
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见