Page 1 of 1

nonlr.f bug error

Posted: Mon Dec 09, 2024 2:53 pm
by inioluwa_afolabi

Hello,

I am trying to run NEB calcs for 5 images for water on defective graphene surface. The first calcs, ran well without issues, but reruning the unfinished calculation after clearing my run directory with vfin.pl started giving this bug error.

Code: Select all

 -----------------------------------------------------------------------------
|                     _     ____    _    _    _____     _                     |
|                    | |   |  _ \  | |  | |  / ____|   | |                    |
|                    | |   | |_) | | |  | | | |  __    | |                    |
|                    |_|   |  _ <  | |  | | | | |_ |   |_|                    |
|                     _    | |_) | | |__| | | |__| |    _                     |
|                    (_)   |____/   \____/   \_____|   (_)                    |
|                                                                             |
|     internal error in: nonlr.F  at line: 596                                |
|                                                                             |
|     REAL_OPTLAY: internal error (1) 483402 492794                           |
|                                                                             |
|     If you are not a developer, you should not encounter this problem.      |
|     Please submit a bug report.                                             |
|                                                                             |
 -----------------------------------------------------------------------------

I have tried to increase NCORE in my incar and delete the vaspout.h5 in the folders before rerunning but it doesn't help.

Any advise on how to solve this bug will be appreciated. Thanks.


Re: nonlr.f bug error

Posted: Wed Dec 11, 2024 1:23 pm
by merzuk.kaltak

Dear inioluwa_afolabi,

it seems the reason for the failure is that one or more images have ions that are way too close to each other.
VASP is actually reporting this, for instance looking at the out.log file, I see following line:

Code: Select all

...
 06/POSCAR found :  3 types and      52 ions
 -----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|     The distance between some ions is very small. Please check the          |
|     nearest-neighbor list in the OUTCAR file.                               |
|     I HOPE YOU KNOW WHAT YOU ARE DOING!                                     |
|                                                                             |
 -----------------------------------------------------------------------------

 POSCAR, INCAR and KPOINTS ok, starting setup
...

I have looked at 01/POSCAR using vesta and found that the distance between the C-ion 32 and 35 is less than 0.42 A. Other ions seem also seem to be very close to each other. This is most likely the reason why your second NEB crashes.

Also, I am not convinced your first NEB calculation finished successfully.
For instance, inspecting the 01/OUTCAR of your first NEB shows already increasing forces after each ionic step.
The same holds true for other images.

I am not entirely sure what's the cause of this problem, though.

Maybe IVDW causes it? Maybe, you need to use an alternative ionic minimizer, such as IBRION=3 as described on SPRING wiki entry.
In any case, I would first try to converge with an simpler DFT functional, for instance plain PBE or Opt-PBE, and continue with the preconverged images using the preferred IVDW settings.