VASP激发态优化报错

maoxinxina

为了计算发射能，这里计算了晶体的激发态，通过人为控制电子的排列。在用PBE泛函的时候，可以很快地优化出来结果，在用PBE0泛函的时候，就会报错停止。
INCAR文件如下：
System = triplet_opt
NPAR = 4
#ISTART = 1
#ISYM = 2
ENCUT = 200
EDIFF = 1E-05
EDIFFG = -0.02
LREAL = Auto
PREC = Normal
ISPIN = 2
ISMEAR = -2
NBANDS= 216
FERWE = 103*1 1*0 1*1 111*0   #The syntax is BANDS*OCC
FERDO = 104*1   112*0  
AMIX = 0.4
BMIX = 1.0
NSW = 200 #sets the maximum number of ionic steps
NELM = 350 #sets the maximum number of electronic sc
ISIF = 2
ALGO = VeryFast
#LORBIT = 11
SIGMA = 0.1
IBRION = 3 #determines how the ions are updated and moved
NGX = 86
NGY = 42
NGZ = 30
LHFCALC = .TRUE.
GGA = PE
#ISYM = 0
#LORBIT = 11
POTIM = 0.1
#VOSKOWN = 1
#LMAXMIX = 4
#SAXIS = 0 0 1
#LNONCOLLINEAR = .TRUE.
#LPLANE = .FALSE.


KPOINTS如下：
KPT-Resolved Value to Generate K-Mesh: 3
0
Gamma
1   1   1
0.0  0.0  0.0


POSCAR如下
Cs2 Cu4 Cl6
   1.00000000000000
    15.8497257233000006    0.0000000000000000    0.0000000000000000
    -1.9435727601999999    7.6828364750000002    0.0000000000000000
     0.0000000000000000    0.0000000000000000    5.3450770378000003
   Cs   Cu   Cl
     4     8    12
Direct
  0.3466712934236256  0.3049552001793430  0.2500000000000000
  0.8466713234236209  0.3049552001793430  0.2500000000000000
  0.1533287215763792  0.6950448298206524  0.7500000000000000
  0.6533287365763840  0.6950448298206524  0.7500000000000000
  0.0862212967527469  0.1742094154131664  0.0002576158215035
  0.5862213037527439  0.1742094154131664  0.0002576158215035
  0.4137786962472489  0.8257905695868430  0.9997423841784965
  0.9137786962472561  0.8257905695868430  0.9997423841784965
  0.4137786962472489  0.8257905695868430  0.5002576158215035
  0.9137786962472561  0.8257905695868359  0.5002576158215035
  0.0862212967527469  0.1742094154131664  0.4997423841784965
  0.5862213037527439  0.1742094154131664  0.4997423841784965
  0.0512364027911190  0.8988172285028995  0.2500000000000000
  0.5512363917911216  0.8988172285028995  0.2500000000000000
  0.4487636082088784  0.1011827944970918  0.7500000000000000
  0.9487636082088784  0.1011827944970918  0.7500000000000000
  0.0792033381111708  0.4037651834189973  0.2500000000000000
  0.5792033611111691  0.4037651834189973  0.2500000000000000
  0.4207966688888263  0.5962347865810003  0.7500000000000000
  0.9207966388888309  0.5962347865810003  0.7500000000000000
  0.2990865795871827  0.8399924508728631  0.2500000000000000
  0.7990866095871922  0.8399924508728631  0.2500000000000000
  0.2009134054128126  0.1600075641271346  0.7500000000000000
  0.7009133904128078  0.1600075641271346  0.7500000000000000


在电子步优化到第15步的时候，就会自动停止，然后就不动了。
想知道这是什么问题。

zisheng

qq870624836 发表于 2022-3-3 13:29
谢谢老师

你好，请问你会计算激发态了吗

qq870624836

卡开发发 发表于 2022-3-3 13:17
对。

谢谢老师

卡开发发

qq870624836 发表于 2022-3-3 13:12
那这个能量差是不是可以在第三步的结果用grep energy OUTCAR搜索得到能量，再用同样的方法在第四步里得到 ...

对。

qq870624836

卡开发发 发表于 2022-3-3 13:06
对。

那这个能量差是不是可以在第三步的结果用grep energy OUTCAR搜索得到能量，再用同样的方法在第四步里得到D的能量然后相减就可以了把？

卡开发发

qq870624836 发表于 2022-3-3 13:02
我在https://zhuanlan.zhihu.com/p/454376223看到了优化激发的方法，里面C态和D态的能量差是否就是发射谱 ...

对。

qq870624836

卡开发发 发表于 2021-7-11 16:29
建议读取基态下的SCF的波函数，250的动能截断太低了。

我在https://zhuanlan.zhihu.com/p/454376223看到了优化激发的方法，里面C态和D态的能量差是否就是发射谱需要的激发能？

maoxinxina

山里面的和尚 发表于 2021-8-7 15:39
楼主这个体系我也尝试过，我算出来的和文献有区别

我用PBE0泛函算的话感觉跑不动，你那里能跑动吗。

山里面的和尚

楼主这个体系我也尝试过，我算出来的和文献有区别

山里面的和尚

请问一下楼主，你的FERWE和FERDO怎么设置的？还有你的K点设置的伽马点，但是你的体系VBM和CBM就在伽马点吗？

successjiang

请教一下，发射能是怎么定义的？是指一个电子从基态跃迁到费米能级以上空轨道的能量吗？

maoxinxina

卡开发发 发表于 2021-7-11 16:29
建议读取基态下的SCF的波函数，250的动能截断太低了。

嗯嗯，谢谢了。

卡开发发

maoxinxina 发表于 2021-7-9 10:15
把截断能改为250的时候，vasp.out 显示的结果如下：

POSCAR found type information on POSCAR  Cs Cu  ...

建议读取基态下的SCF的波函数，250的动能截断太低了。

maoxinxina

把截断能改为250的时候，vasp.out 显示的结果如下：

POSCAR found type information on POSCAR  Cs Cu Cl
POSCAR found :  3 types and      24 ions
scaLAPACK will be used
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
WAVECAR not read
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.165722427423E+04    0.16572E+04   -0.37310E+04   432   0.397E+02
RMM:   2     0.986186551870E+03   -0.67104E+03   -0.78497E+03   432   0.116E+02
RMM:   3     0.888158585797E+03   -0.98028E+02   -0.16101E+03   432   0.648E+01
RMM:   4     0.926437039734E+03    0.38278E+02   -0.51284E+02   432   0.560E+01
RMM:   5     0.946873408267E+03    0.20436E+02   -0.20777E+02   432   0.254E+01
RMM:   6     0.947110829056E+03    0.23742E+00   -0.15127E+02   432   0.270E+01
RMM:   7     0.955333499953E+03    0.82227E+01   -0.13755E+02   432   0.182E+01
RMM:   8     0.950508432298E+03   -0.48251E+01   -0.13857E+02   432   0.241E+01
RMM:   9     0.951373338494E+03    0.86491E+00   -0.17450E+02  1109   0.192E+01
RMM:  10     0.947741687115E+03   -0.36317E+01   -0.14042E+02  1103   0.117E+01
RMM:  11     0.949073351592E+03    0.13317E+01   -0.14217E+02  1107   0.120E+01
RMM:  12     0.946443590385E+03   -0.26298E+01   -0.13925E+02  1110   0.115E+01    0.330E+01
RMM:  13     0.956546733447E+03    0.10103E+02   -0.23449E+02  1017   0.252E+01    0.177E+01
RMM:  14     0.982609270047E+03    0.26063E+02   -0.22624E+02  1003   0.328E+01    0.846E+00
RMM:  15     0.960674816614E+03   -0.21934E+02   -0.27600E+02  1140   0.152E+01    0.488E+00
-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

卡开发发

暂时没办法判断，不行把文件发上来供大家一块讨论。另外200eV的动能截断有些离谱。