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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【问题已解决】vasp计算自旋轨道耦合读问题
北京石油化工学院2026年研究生招生接收调剂公告
111/1返回列表
查看: 4154  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhang668

木虫 (正式写手)

[交流] 【问题已解决】vasp计算自旋轨道耦合读问题

如题，先进行关闭自旋轨道耦合的共线的scf计算后，得到CHGCAR和WAVCAR，然后设置ICHARGE=11，计算非自洽时候，总出现错误，提示为：
running on   64 nodes
distr:  one band on    1 nodes,   64 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  O  Fe Bi
POSCAR found :  3 types and   10 ions
LDA part: xc-table for Ceperly-Alder, standard interpolation
found WAVECAR, reading the header
  number of k-points has changed, file:    19 present:   125
  trying to continue reading WAVECAR, but it might fail
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(           1 )
reading WAVECAR
reading wavefunctions of collinear run, up
reading wavefunctions of collinear run, down
the WAVECAR file was read sucessfully
charge-density read from file: BiFeO3                                 
magnetization density read from file 1
LAPACK: Routine ZPOTRF failed!           1
LAPACK: Routine ZPOTRF failed!           1
LAPACK: Routine ZPOTRF failed!           1
LAPACK: Routine ZPOTRF failed!           1

就是K点数变化，出现问题。不管第一步有没设置ISYM=0,第二步的计算的K点数都和第一步不一样，所以总是出现问题。只有第二步计算不读取波函数和电荷密度，计算才可以进行下去。

[ Last edited by zhang668 on 2010-12-2 at 22:21 ]
回复此楼

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

vasp计算资料 有机 自旋轨道耦合

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-12-01 23:20:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
★
zhang668(金币+1):谢谢参与
你好，我也遇到了相类似的问题，提示K点不一致。请问应该如何解决？

谢谢！
赞一下(1人) 回复此楼
2楼2011-07-19 11:34:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

deligent

禁虫 (初入文坛)
★
zhang668(金币+1):谢谢参与
本帖内容被屏蔽
3楼2011-07-19 15:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
引用回帖:
Originally posted by deligent at 2011-07-19 15:50:30:
可能是NBANDS设置的问题，NBANDS = 2 * number of bands of collinear run。

设置了，但是还不行。
赞一下 回复此楼
4楼2011-07-24 22:17:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)
★ ★
zhang668(金币+1):谢谢参与
bingmou(金币+1): 感谢参与 2011-07-31 17:55:16
引用回帖:
Originally posted by byin at 2011-07-24 22:17:19:
设置了，但是还不行。

第一步scf计算也开启自诩轨道耦合和非共线，这样产生的电荷密度读入就没错了，说明书说的有问题。
赞一下(2人) 回复此楼
5楼2011-07-30 23:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by whw19850730 at 2011-07-30 23:43:00:
第一步scf计算也开启自诩轨道耦合和非共线，这样产生的电荷密度读入就没错了，说明书说的有问题。

你的意思是直接一步计算，对吗？

istart=0
icharge=2

自旋轨道耦合 开



直接这么算？
赞一下(1人) 回复此楼
6楼2011-08-01 00:20:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by byin at 2011-08-01 00:20:05:
你的意思是直接一步计算，对吗？

istart=0
icharge=2

自旋轨道耦合 开



直接这么算？

是的读入非共线的电荷密度就不会有错了。你可以测试一下
赞一下(1人) 回复此楼
7楼2011-08-01 10:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyb83

新虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
ISYM=0
关闭对称性操作。当启用自旋轨道耦合时，不同磁矩指向情况下体系的对称性存在差异，因此k点也会相应变化。
我用的办法是从共线计算开始就关掉对称性操作，然后用IBZKPT 代替KPOINTS，保证各个方向计算时使用相同的K点。
赞一下(1人) 回复此楼
8楼2011-12-06 17:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastfreedom

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by sunyb83 at 2011-12-06 17:23:47:
ISYM=0
关闭对称性操作。当启用自旋轨道耦合时，不同磁矩指向情况下体系的对称性存在差异，因此k点也会相应变化。
我用的办法是从共线计算开始就关掉对称性操作，然后用IBZKPT 代替KPOINTS，保证各个方向计算时 ...

请教：这样计算和relax的时候就加soc，费的时间是不是会少点，对大体系呢
赞一下(1人) 回复此楼
9楼2011-12-06 17:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ffzx

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
1943759楼: Originally posted by lastfreedom at 2011-12-06 17:42:07
请教：这样计算和relax的时候就加soc，费的时间是不是会少点，对大体系呢...

大体系本来就慢，直接加SOC relax岂不更慢？
赞一下 回复此楼
10楼2012-06-02 20:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

国花1989

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主你好，你的问题解决了吗？？我的自旋轨道耦合没有加进去，是不是参数设置错啦？？
赞一下 回复此楼
11楼2013-06-09 08:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhang668 的主题更新
111/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 282求调剂 +6 aaa车辆 2026-04-02 10/500 2026-04-05 12:10 by arrow8852
[考研] 材料专硕283求调剂 +11 试试看呗 2026-04-04 12/600 2026-04-05 11:36 by 猪会飞
[考研] 化学357分，考研调剂 +10 .Starry. 2026-04-04 11/550 2026-04-05 10:57 by cql1109
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:57 by 544594351
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-31 8/400 2026-04-04 22:08 by hemengdong
[考研] 材料调剂 +12 一样YWY 2026-04-03 12/600 2026-04-04 21:46 by hemengdong
[考研] 考研调剂 +4 zybz冲冲冲 2026-04-03 6/300 2026-04-04 13:08 by zybz冲冲冲
[考研] 343求调剂085601 +6 要努力学习x 2026-03-29 7/350 2026-04-03 19:49 by 百灵童888
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦！ 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 化学工程专硕324分，一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 本科211生物医学工程085409求调剂339分 +7 里子木yy 2026-03-29 7/350 2026-03-31 14:35 by fmesaito
信息提示
关闭
请填处理意见
关闭