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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【问题已解决】vasp计算自旋轨道耦合读问题
111/1返回列表
查看: 4066  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 考研材料与化工，求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
关闭
请填处理意见
关闭