VASP的FERWE报错求助

wendawei

自己试手固定电子的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
------------------分界线------------------

卡开发发

多k点还是只有Gamma点？

wendawei

卡开发发 发表于 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

卡开发发

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点的占据情况也要写在里面。

wendawei

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

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

原来如此，谢谢！

不明嚼栗子

楼主您好，我想请教下多k点的情况该怎么写FERWE和FERDO两个tag呀，万分感谢！

smalll_wu

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

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

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

卡开发发

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部分输出的占据信息一致。

smalll_wu

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

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

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

盈朕

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

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

smalll_wu

盈朕 发表于 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）

盈朕

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”是根据什么取得

smalll_wu

盈朕 发表于 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数目。

since-2010

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