Using Vasp as calculator in ASE

[karol_kawka]

Hello,

I'm trying to use VASP as a calculator in ASE to obtain minimum‑energy paths with the so‑called NEB method.
Previously I used another DFT code in ASE (Siesta) and it worked well. I’d like to compare the results from the two programs, but I’m having trouble doing the same with VASP.

It seems that during a single‑point calculation VASP produces very large forces in the first optimisation step, which show up as asterisks (*) for certain atoms. For example, the output from the second image (out of seven) has this issue—only img0 and img1 ran to completion. Is there a tag or setting that can prevent this error? This is how the NEB routine in ASE launches VASP; if there’s no workaround, I may have to switch to VASP’s own NEB.

[christopher_sheldon1]

Hello Karol,

Thank you for your question. It's not immediately clear from the image what is happening. What is the system that you are studying? I wonder if those particular lines are for a specific element or atoms that lie very closely together. This can sometimes cause problems with the images generated using NEB.

If you could upload the POSCAR, POTCAR, KPOINTS, INCAR, and OUTCAR, then it would help me to make more sense of the issue.

Best wishes,

Chris

[karol_kawka]

I uploaded files from image 2 of the procedure.

The issue is that ASE starts from the initial atomic positions and interpolates the path to the final positions. In my case, there are 7 images in total (from 0 to 6, so the uploaded one is the 3rd structure). As you can see, VASP completed the job, but the forces at the end of the OUTCAR file are shown as asterisks, so ASE couldn’t read them to determine the next step in the optimization path.

[christopher_sheldon1]

Hi Karol,

Thank you for the input and output files. The forces are far too large and meaningless. The energy is also positive (check your OSZICAR file to see that each elementary step is positive). Does this occur for all of your images?

I also notice that there seem to be some atoms missing at the edge of the cell (when viewing a 3x3x1 supercell). I suspect that this will cause big issues (though the positive energy is still unexplained).

Best wishes,

Chris