24小时热门版块排行榜    

查看: 4069  |  回复: 10

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计算资料 有机 自旋轨道耦合

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)



zhang668(金币+1):谢谢参与
你好,我也遇到了相类似的问题,提示K点不一致。请问应该如何解决?

谢谢!
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计算也开启自诩轨道耦合和非共线,这样产生的电荷密度读入就没错了,说明书说的有问题。
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

自旋轨道耦合 开



直接这么算?
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

自旋轨道耦合 开



直接这么算?

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

sunyb83

新虫 (初入文坛)



小木虫(金币+0.5):给个红包,谢谢回帖
ISYM=0
关闭对称性操作。当启用自旋轨道耦合时,不同磁矩指向情况下体系的对称性存在差异,因此k点也会相应变化。
我用的办法是从共线计算开始就关掉对称性操作,然后用IBZKPT 代替KPOINTS,保证各个方向计算时使用相同的K点。
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,费的时间是不是会少点,对大体系呢
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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 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
信息提示
请填处理意见