Large influence of SYMPREC on energies

[MBaeker]

I have a question regarding the SYMPREC-parameter: I have found a case where changing it from the standard 1.e-05 to 1.e-04 changes the energy quite significantly, from
E0= -.79707341E+03 (correct value, as far as I can tell) to
E0= -.79920458E+03

Looking at the OUTCAR (which I cannot upload because it is too large), the only differences I see in the setup (before iteration begins) are:
< SYMPREC = 1.00e-04
> SYMPREC = 1.00e-05
512c512
< Subroutine GETGRP returns: Found 16 space group operations
> Subroutine GETGRP returns: Found 8 space group operations
518d517
< The point group associated with its full space group is D_4h.
542c541
< Subroutine GETGRP returns: Found 16 space group operations
> Subroutine GETGRP returns: Found 8 space group operations
548d546
< The point group associated with its full space group is D_4h.
572c570
< Subroutine GETGRP returns: Found 16 space group operations
> Subroutine GETGRP returns: Found 8 space group operations
578d575
< The point group associated with its full space group is D_4h.
581c578
< Subroutine INISYM returns: Found 16 space group operations
> Subroutine INISYM returns: Found 8 space group operations

When the simulation starts (OSZICARs attached), energies are the same, but then they drift apart significantly.

I do understand that choosing a value of SYMPREC that is too large may force an incorrect symmetry on the charge distribution, but I am puzzled by the effect being so large here and by the fact that the energy with too large SYMPREC value is lower (shouldn't it be higher because there are less DOFs to minimize when we enforce a symmetry?).
I find it also a bit worrying considering that "change SYMPREC" is advice frequently given when a VASP run has trouble with the symmetry.

I've attached the INCAR and POSCAR and the 2 OSZICAR results; I have not attached the POTCAR because I think that's forbidden, the potentials used were
TITEL = PAW_PBE Ni_pv 06Sep2000
TITEL = PAW_PBE Nb_sv 25May2007

[fabien_tran1]

Hi,

As a test, could you please follow these steps (all with NSW=0, i.e. no geometry optimization)?

1) Do a calculation with "SYMPREC = 1.00e-04" and "LWAVE = .TRUE.".

2) In the same directory, do the calculation with "SYMPREC = 1.00e-05". This calculation should start using the WAVECAR file generated at step 1.

Is the energy obtained at step 2 the same as the one that you reported?

[MBaeker]

I tried with two different INCARs, one my standard one, one where I reduce AMIX and set NELMDL to get the initialization right.
In both cases, the job does not converge, the first ends with
| Error EDDDAV: Call to ZHEGV failed. Returncode = 23 1 24 |
| |
| ----> I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <----

the second (different mixing and NELMDL set) diverges and gives lots of
WARNING: Sub-Space-Matrix is not hermitian in DAV 60
6.535011242324960E-005
warnings, but with an energy of
DAV: 66 -0.145384783713E+10 -0.94014E+08 -0.32750E+08207888 0.604E+03 0.213E+03
that is probably not surprising.

I attach all slurm outputs:
slurm-9298808.out initial run with SYMPREC 1.e-4
slurm-9298931.out restart with SYMPREC 1e-5 and my usual INCAR parameters
slurm-9299562.out restart (from the first WAVECAR) with SYMPREC 1e-5 and my second try for INCAR parameters.

If there is anything else I can do to clear this up, please tell me.

BTW, I am very sure that from the two runs I posted yesterday, the result with SYMPREC=1.e-4 and the lower energy is incorrect.