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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【问题已解决】vasp计算自旋轨道耦合读问题
111/1返回列表
查看: 4067  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +4 Liwangman 2026-03-15 4/200 2026-03-16 08:01 by wang_dand
[考研] 化学工程321分求调剂 +6 大米饭！ 2026-03-15 6/300 2026-03-16 07:58 by wang_dand
[考研] 本人考085602 化学工程 专硕 +8 不知道叫什么！ 2026-03-15 9/450 2026-03-16 07:45 by L135790
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
关闭
请填处理意见
关闭