24小时热门版块排行榜

>论坛更新日志 (3396)
>虫友互识 (294)
>导师招生 (258)
>文献求助 (175)
>休闲灌水 (142)
>教师之家 (106)
>考博 (91)
>基金申请 (88)
>论文投稿 (77)
>考研 (64)
>硕博家园 (59)
>论文道贺祈福 (54)
>找工作 (50)
>博后之家 (42)
>招聘信息布告栏 (40)
>公派出国 (24)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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计算资料

有机

自旋轨道耦合

» 猜你喜欢

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

关于VASP中MAGMON的磁矩设置问题已经有6人回复
vasp计算磁性的时候 OSICAR中有总磁矩的大小，但是OUTCAR中每个原子的磁矩均显示为0 已经有8人回复
怎么补算电声子耦合？已经有3人回复
VASP自旋轨道耦合已经有26人回复
NBANDS 的设置问题已经有6人回复
自旋轨道耦合第一性原理计算已经有11人回复
VASP 运行报错已经有11人回复
vasp自旋极化计算时，初始设置自旋能否不被优化已经有6人回复
VASP如何读取分子轨道？已经有18人回复
VASP+SOC后的能带计算已经有21人回复
求助vasp中的自旋轨道耦合已经有3人回复
ADF计算自旋轨道耦合已经有7人回复
VASP计算：如何处理波函数占用空间太大的问题已经有5人回复
请问MS里怎么考虑自旋轨道耦合作用？已经有9人回复
【求助】LAPACK: Routine ZPOTRF failed! 已经有5人回复
【求助】vasp 离子驰豫与电子自旋的关系已经有11人回复
【求助】VASP考虑自旋轨道耦合的话，如何编译已经有17人回复
【求助成功】charge-density read from file:unknown 是怎么回事呢? 已经有3人回复

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

查看全部散金贴

1楼 2010-12-01 23:20:53

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

lastfreedom

木虫 (小有名气)

应助: 4 (幼儿园)
金币: 4209.4
帖子: 99
在线: 658.9小时
虫号: 1087212

★
小木虫(金币+0.5):给个红包，谢谢回帖

引用回帖:

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

请教：这样计算和relax的时候就加soc，费的时间是不是会少点，对大体系呢