应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于spin orbit coupling不收敛的问题 求助各位高手啦
91/1返回列表
查看: 1856  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

liaoren5

木虫 (正式写手)

[求助] 关于spin orbit coupling不收敛的问题 求助各位高手啦已有1人参与

大家好，我最近在算WS2单层体系，需要用到spin orbit coupling，于是求助上海超算上的管理员帮我编译了SOC版的并行VASP5.2
但是在计算WS2单层的时候，一直不收敛，非常郁闷，求助各位大家啦！
之前尝试过不用SOC版的VASP，可以收敛，但是用了SOC版的VASP之后，不管加不加SOC的相关设置都不收敛。
也不知道是啥问题，是设置问题还是编译问题，因为我自己在上海超算上总是编译不成功，所以求助他们的管理员帮我编译的。
谢谢大家啦
这是我的INCAR，在SOC版不收敛
SYSTEM = WS2

   ISTART = 0
   ICHARG = 2

   IBRION = 2
   ISIF   = 2
   NSW   = 100
   
   PREC   = Medium
   ENCUT  = 400
   NELM   = 1000
   NELMIN = 2
   EDIFF  = 1.0e-4
   ISMEAR =  0
   SIGMA  = 0.1
   LWAVE = .FALSE.
  LCHARG = .FALSE.
   IALGO  = 48

   LSORBIT = .TRUE.
   LMAXMIX = 4
   NBANDS = 32
   MAGMOM = 0 0 0 0 0 0 0 0 0
   GGA_COMPAT = .FALSE.
   ISYM = 0
   SAXIS = 0 0 1

   LREAL = .FALSE.

不用SOC版的VASP能很好的收敛，INCAR如下（也拿这个INCAR在SOC版试了，也不收敛）：
SYSTEM = WS2

   ISTART = 0
   ICHARG = 2

   IBRION = 2
   ISIF   = 2
   NSW   = 100
   
   PREC   = Medium
   ENCUT  = 400
   NELM   = 1000
   NELMIN = 2
   EDIFF  = 1.0e-4
   ISMEAR =  0
   SIGMA  = 0.1
   LWAVE = .FALSE.
  LCHARG = .FALSE.
   IALGO  = 48

   LREAL = .FALSE.
回复此楼

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

自旋轨道耦合

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-12-21 13:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)
能否說一下是怎樣的不收斂?
能量上下震盪? 還是出現 error message? warning message?
赞一下 回复此楼
2楼2013-12-21 18:33:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liaoren5

木虫 (正式写手)
引用回帖:
2楼: Originally posted by jpchou at 2013-12-21 18:33:59
能否說一下是怎樣的不收斂?
能量上下震盪? 還是出現 error message? warning message?

就是能量一直震荡不收敛，这个是截取到三百多步，其实已经算了2000多步就是不收敛啊

RMM: 372     0.232601230704E+03    0.39695E+02    0.16783E+03  1398   0.173E+02    0.305E+00
RMM: 373     0.200120444466E+03   -0.32481E+02    0.41439E+02  1385   0.198E+02    0.204E+00
RMM: 374     0.224944037297E+03    0.24824E+02    0.13004E+03  1397   0.176E+02    0.309E+00
RMM: 375     0.203419852701E+03   -0.21524E+02    0.54658E+02  1394   0.192E+02    0.252E+00
RMM: 376     0.211581013622E+03    0.81612E+01    0.10757E+03  1372   0.181E+02    0.252E+00
RMM: 377     0.202071962315E+03   -0.95091E+01    0.76228E+02  1387   0.181E+02    0.248E+00
RMM: 378     0.184480670422E+03   -0.17591E+02    0.80134E+02  1366   0.181E+02    0.232E+00
RMM: 379     0.209147679050E+03    0.24667E+02    0.12268E+03  1403   0.172E+02    0.284E+00
RMM: 380     0.186569067088E+03   -0.22579E+02    0.80967E+02  1384   0.183E+02    0.207E+00
RMM: 381     0.197445569875E+03    0.10877E+02    0.96235E+02  1386   0.178E+02    0.234E+00
RMM: 382     0.207623184187E+03    0.10178E+02    0.11965E+03  1396   0.179E+02    0.288E+00
RMM: 383     0.202025366728E+03   -0.55978E+01    0.74417E+02  1382   0.187E+02    0.246E+00
RMM: 384     0.208244957474E+03    0.62196E+01    0.10661E+03  1405   0.180E+02    0.262E+00
RMM: 385     0.179394007863E+03   -0.28851E+02    0.45520E+02  1378   0.186E+02    0.162E+00
RMM: 386     0.200009314610E+03    0.20615E+02    0.13549E+03  1388   0.171E+02    0.205E+00
RMM: 387     0.176457202955E+03   -0.23552E+02    0.54682E+02  1372   0.185E+02    0.175E+00
RMM: 388     0.213454065941E+03    0.36997E+02    0.14114E+03  1397   0.169E+02    0.273E+00
RMM: 389     0.201833939226E+03   -0.11620E+02    0.55763E+02  1367   0.193E+02    0.189E+00
RMM: 390     0.225960922607E+03    0.24127E+02    0.15142E+03  1395   0.165E+02    0.298E+00
RMM: 391     0.190925928204E+03   -0.35035E+02    0.24449E+02  1391   0.191E+02    0.195E+00
RMM: 392     0.202571566992E+03    0.11646E+02    0.11356E+03  1366   0.174E+02    0.226E+00
RMM: 393     0.189771788659E+03   -0.12800E+02    0.59809E+02  1389   0.183E+02    0.221E+00
赞一下 回复此楼
3楼2013-12-22 19:51:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流 2013-12-23 21:33:51
這震盪的幅度有夠誇?�
先懷疑一下也許管理員沒有編好

測試一下
1)  沒 SOC 的 INCAR 用一般 VASP 算一次
2)  沒 SOC 的 INCAR 再用能計算 SOC 的 VASP 算一次
理論上這兩個結果要一模一樣
如果沒 SOC 的 2) INCAR 計算後還是不收斂~或結果與 1) 不一致
就可以說明了這個有 SOC 的 VASP 是有問題的
赞一下 回复此楼
4楼2013-12-23 01:56:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liaoren5

木虫 (正式写手)
引用回帖:
4楼: Originally posted by jpchou at 2013-12-23 01:56:09
這震盪的幅度有夠誇?�
先懷疑一下也許管理員沒有編好

測試一下
1)  沒 SOC 的 INCAR 用一般 VASP 算一次
2)  沒 SOC 的 INCAR 再用能計算 SOC 的 VASP 算一次
理論上這兩個結果要一模一樣
如果沒 SOC 的 2) ...

多谢啦，这问题很麻烦啊！
赞一下 回复此楼
5楼2013-12-23 09:24:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SplendorZhao

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by jpchou at 2013-12-23 01:56:09
這震盪的幅度有夠誇?�
先懷疑一下也許管理員沒有編好

測試一下
1)  沒 SOC 的 INCAR 用一般 VASP 算一次
2)  沒 SOC 的 INCAR 再用能計算 SOC 的 VASP 算一次
理論上這兩個結果要一模一樣
如果沒 SOC 的 2) ...

你好，我也遇到了这种问题，SOC的程序没问题，因为我是调应力，到10%之后就不收敛了，怎么办？新人求指点
赞一下 回复此楼
6楼2015-04-02 21:40:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liaoren5

木虫 (正式写手)
引用回帖:
6楼: Originally posted by SplendorZhao at 2015-04-02 21:40:31
你好，我也遇到了这种问题，SOC的程序没问题，因为我是调应力，到10%之后就不收敛了，怎么办？新人求指点...

这我也不知道啦
回复此楼
7楼2015-04-02 22:18:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)
突然發現
你的能量是"正"的耶
那表示你這計算根本有問題
如果你確定一般版本計算這系統(同樣 POSCAR,POTCAR,INCAR,KPOINTS) 沒有問題
但 SOC 卻是這樣結果
那很名顯可以看出
這個 SOC 版本 VASP 應該是沒有正確編譯
赞一下 回复此楼
8楼2015-04-03 03:06:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liaoren5

木虫 (正式写手)
引用回帖:
8楼: Originally posted by jpchou at 2015-04-03 03:06:30
突然發現
你的能量是"正"的耶
那表示你這計算根本有問題
如果你確定一般版本計算這系統(同樣 POSCAR,POTCAR,INCAR,KPOINTS) 沒有問題
但 SOC 卻是這樣結果
那很名顯可以看出
這個 SOC 版本 VASP 應 ...

高手，确实没有正确编译
后来编译个正确的就行了
赞一下 回复此楼
9楼2015-04-03 09:18:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liaoren5 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭