计算化学公社

标题: VASP的FERWE报错求助 [打印本页]

作者
Author: wendawei 时间: 2020-10-12 17:25
标题: VASP的FERWE报错求助
自己试手固定电子的occupation，尝试使用FERWE和FERDO参数。已知体系电子数65个（LaCeB2O6,1*11+1*12+2*3+6*6=65，使用了Ce赝势而不是Ce_3赝势），up-spin和down-dpin打算顾定成33和32个电子，但是报错了，故求问。
另外我试过其他体系也是同样失败，但是我试着把Ce3+离子（改ELECT，加了背景电荷）放在10×10×10埃的格子中心然后使用FERWE和FERDO参数，却能跑起来。。。好迷
需要提供更多信息的话请告诉我，谢谢了！

------------------分界线------------------
错误信息如下：
Error reading item FERWE from file INCAR.                            |
|    Error code was IERR= 7 ... . Found N= 40 data.
------------------分界线------------------

我的INCAR文件
ISTART =  1          (Read existing wavefunction; if there)
ISPIN =  2          (Spin polarised DFT)
LREAL  = Auto       (Projection operators: automatic)
PREC =  Normal    (Precision level)
LWAVE  = .TRUE.       (Write WAVECAR or not)
LCHARG = .TRUE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid; helps GGA convergence)
ISMEAR =  -2          (gaussian smearing method)
SIGMA  =  0.05       (please check the width of the smearing)
LORBIT =  11          (PAW radii for projected DOS)
NELM =  100          (Max electronic SCF steps)
EDIFF  =  1E-06       (SCF energy convergence; in eV)
BMIX = 0.01
FERWE = 33*1.0 7*0.0
FERDO = 32*1.0 8*0.0
NCORE = 4
------------------分界线------------------

作者
Author: 卡开发发 时间: 2020-10-12 19:18
多k点还是只有Gamma点？

作者
Author: wendawei 时间: 2020-10-12 19:21

卡开发发发表于 2020-10-12 19:18
多k点还是只有Gamma点？

KPOINTS文件用vaspkit生成的，内容如下：

KPT-Resolved Value to Generate K-Mesh: 0.040
0
Gamma
4 5 6
0.0 0.0 0.0

作者
Author: 卡开发发 时间: 2020-10-12 19:50

wendawei 发表于 2020-10-12 19:21
KPOINTS文件用vaspkit生成的，内容如下：

KPT-Resolved Value to Generate K-Mesh: 0.040

手册：

If the occupancies are fixed by you, there should be a tag
FERWE = f1 f2 f3 ... f(NBANDS)
and for spin-polarized calculations
FERDO = f1 f2 f3 ... f(NBANDS)
in the INCAR file supplying the partial occupancies for all bands and k-points.

所以其他k点的占据情况也要写在里面。

作者
Author: wendawei 时间: 2020-10-12 20:49

卡开发发发表于 2020-10-12 19:50
手册：

所以其他k点的占据情况也要写在里面。

原来如此，谢谢！

作者
Author: 不明嚼栗子 时间: 2020-12-24 23:21
楼主您好，我想请教下多k点的情况该怎么写FERWE和FERDO两个tag呀，万分感谢！

作者
Author: smalll_wu 时间: 2021-8-31 17:01

卡开发发发表于 2020-10-12 19:50
手册：

所以其他k点的占据情况也要写在里面。

老师您好，请问如果有多个K点，该怎么设置相应的FERWE和FERDO的参数

作者
Author: 卡开发发 时间: 2021-9-2 05:03

smalll_wu 发表于 2021-8-31 17:01
老师您好，请问如果有多个K点，该怎么设置相应的FERWE和FERDO的参数

手册这段话要好好看：

in the INCAR file supplying the partial occupancies for all bands and k-points. The band-index runs fastest. The partial occupancies must be between 0 and 1 (for spin-polarized and non-spin-polarized calculations).
Mind: Partial occupancies are also written to the OUTCAR file, but in this case they are multiplied by 2, i.e. they are between 0 and 2.

The band-index runs fastest.的意思是说先循环每个能带的占据，再是不同的k，大概是这么个意思：
k1-b1 k1-b2 ... k1bn k2-b1 k2-b2 ...其中k1-b1指的是第一个k点的第一个能带。k应该和OUTCAR部分输出的占据信息一致。

作者
Author: smalll_wu 时间: 2022-10-19 16:19

卡开发发发表于 2021-9-2 05:03
手册这段话要好好看：

The band-index runs fastest.的意思是说先循环每个能带的占据，再是不同的k， ...

好久没登录了，多个k点确实是循环写这个占据情况，谢谢老师。

作者
Author: 盈朕 时间: 2022-10-20 19:44

smalll_wu 发表于 2022-10-19 16:19
好久没登录了，多个k点确实是循环写这个占据情况，谢谢老师。

您好，这个手册在哪里可以看到

作者
Author: smalll_wu 时间: 2022-10-20 20:37

盈朕发表于 2022-10-20 19:44
您好，这个手册在哪里可以看到

手册中可能不明显，我之前没找到。不过你可以试试在k点为2*1*1时，把1*1*1时的参数循环一遍，这样应该就能顺利计算了。同理，很多k点时，就循环写很多遍这些参数。k点能利用对称性简并的情况下我还没测试，你可以测试一下。（比如1*1*1时，FERWE = 5*1 4*0，2*1*1时就要写成FERWE = 5*1 4*0 5*1 4*0）

作者
Author: 盈朕 时间: 2022-10-21 11:30

smalll_wu 发表于 2022-10-20 20:37
手册中可能不明显，我之前没找到。不过你可以试试在k点为2*1*1时，把1*1*1时的参数循环一遍，这样应该就 ...

谢谢您的指点，如果说K点是2*2*1时，FERWE=5*1 4*0 5*1 4*0 5*1 4*0 5*1 4*0就要写成这样吗，“5”和“4”是根据什么取得

作者
Author: smalll_wu 时间: 2022-10-27 15:05

盈朕发表于 2022-10-21 11:30
谢谢您的指点，如果说K点是2*2*1时，FERWE=5*1 4*0 5*1 4*0 5*1 4*0 5*1 4*0就要写成这样吗，“5”和“4 ...

FERWE应该是那样写的，不过我没有测试。5和4只是我随便举的例子，具体看体系的总电子数，和设置的NBANDS数目。

作者
Author: since-2010 时间: 2024-12-30 22:04
我想请问下，我如果直接列出每个K点每个能带的占据书，总是报错Error reading item 'FERWE' from file INCAR。如果我写成n*1，m*0...... 这种形式的话是OK的，为什么呢？我的占据数是其他软件算出来了，不是各项同性的，所以只能人工输入，有什么建议嘛？
我在vasp论坛上问了，但是没人回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)