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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【问题已解决】vasp计算自旋轨道耦合读问题
111/1返回列表
查看: 4069  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
信息提示
关闭
请填处理意见
关闭