My Community

Dear All,
I am doing a staic calculation for SiO2 with HSE hybrid function. I follow the setting from manual: first, do the normal DFT calculation to get WAVECAR, then do the hybrid-funtion based on the WAVECAR. But error happen with indication "internal error in FOCK_ACC: number of k-points incorrect 14 100 2 2 2"

I use M-P mesh grid 5x5x4 for the calculation. (This is ok for the normal DFT one). So would someone can give me some advice for that. I am very grateful in advanced.
Below is my INCAR and KPOINTS setting.

..................................................
SYSTEM = SiO2-Primitive Cell

Start parameter for this run:
ISTART = 1
ICHARG = 2
NWRITE = 2
ISPIN = 2

Electronic Relaxation
ENCUT = 540
NELM = 100
EDIFF = 1E-05
LREAL = .FALSE.

Ionic Relaxation
EDIFFG = -0.03 Stopping-Criterion for IOM (HF-FORCE)
NSW = 0
IBRION = -1
ISIF = 2
POTIM = 0.5

LHFCALC = .TRUE.
HFSCREEN = 0.2
NBANDS = 24
ALGO = All
TIME = 0.4
ENCUTFOCK = 0
NKRED = 2

DOS related values:
ISMEAR = 0
SIGMA = 0.1

Write flags

LCHARG = .TRUE.
LWAVE = .TRUE.

.....................

Auto
0
M-K
5 5 4
0 0 0

The problem has been solved after removing the parameter "ENCUTFOCK = 0 NKRED = 2", maybe the problem comes from the setting for NKRED=2 which is not suitable for the present hexagonal cell.

thanks! I also faced this problem.