My Community

Dear VASP Community,

I have tried HSE06 calculations with W.cI2, Cu.cF4, WO3 (Pearson symbols aP32, mP16, mP32, tP8), WO3-aP32 example enclosed

, starting from PAW PBE-relaxed structures from WAVECAR (ISTART = 1).

In the cases of WO3 structures (nonmagnetic), HSE06 did not even started. As for W.cI2 and Cu.cF4, these elemental HSE06 relaxations converged well with identical INCAR file.

What am I doing wrong? Are you able to successfully HSE06-rerelax the file I enclosed?

Eva

Dear Eva,

I was checking the calculation you sent. When executing the simulation I receive the following output message

 -----------------------------------------------------------------------------
|                     _     ____    _    _    _____     _                     |
|                    | |   |  _ \  | |  | |  / ____|   | |                    |
|                    | |   | |_) | | |  | | | |  __    | |                    |
|                    |_|   |  _ <  | |  | | | | |_ |   |_|                    |
|                     _    | |_) | | |__| | | |__| |    _                     |
|                    (_)   |____/   \____/   \_____|   (_)                    |
|                                                                             |
|     internal error in: fock_dbl.F  at line: 1440                            |
|                                                                             |
|     internal error in FOCK_FCC: number of k-points incorrect 27 125 2 2     |
|     2                                                                       |
|                                                                             |
|     If you are not a developer, you should not encounter this problem.      |
|     Please submit a bug report.                                             |
|                                                                             |
 -----------------------------------------------------------------------------

This problem can be overcome by removing the tag NKRED from your INCAR file.

All the best Jonathan

I still encounter problems with HSE06 runs with binary W-O oxides after NKRED tag has been removed.

I tried to relax several dozens of oxides within HSE06 setup and here I choose sample with the smallest 36-atom periodic unit cell: O2W-oP36, where the calculation failed at the very beginning just after the warning (in output):

"found WAVECAR, reading the header
number of bands has changed, file: 256 present: 136
trying to continue reading WAVECAR, but it might fail
WAVECAR: different cutoff or change in lattice found"

I do not understand, why "different cutoff or change in lattice found", sice I did not change the ENCUT = 550eV and I copied PBE CONTCAR (*.PBE files enclosed) to HSE06 POSCAR (*.HSE files enclosed) and I copied and read also WAVECAR (ISTART = 1), I used the same KPOINTS and POTCAR. I could not enclose the WAVECAR, because it is too large a file for the VASP forum.

I also do not understand, how it happened that "number of bands has changed".

I arrived at exactly the same failure (before the first electronic step) in dozens of cases with distinct W-O stoichiometries.

What else (except for placing NKRED tag in INCAR) am I doing wrong?

Hi Eva,

in your PBE calculation you do not set NPAR, which results in an automatic setting of NPAR to the number of processes. So far no harm done.
For performance reasons, VASP sets NBANDS to the number of processes (256?!), so you end up with 256 bands in your PBE simulation.

NPAR is not set (I'm not even sure if allowed) in HSE. You end up with 136 bands. The lowest number whoch you can divide by 8 and fulfills the requirements of your system.

Hence all the mess.

Cheers,

alex

Hi Alex,

could you see any way out? I would like to run HSE from well relaxed PBE sample, but HSE fails at the very beginning. Is there anything I should try?

Eva

Hi Eva,

Alex has a point that the number of bands changed between PBE and HSE06 because of your parallelization settings changed: in both calculations you use 128 MPI ranks but in HSE06 you have NPAR=8 while in PBE you don't.
In principle, if you set NPAR=8 in both calculations it should work because you will get the same NBANDS.

If this does not work still, then I noticed that you have a slightly different POSCAR for the PBE and HSE run. Is this intentional?
If you want to compute HSE with the HSE lattice parameter then I would suggest calculating the WAVECAR with PBE using the HSE lattice parameters.

Hope this helps!

Hello Eva,

please use the same NPAR in both INCARs, or set NBANDS to the highest observed value for starters and report back if it does not work, please.

Best regards,

alex