Inquiry Regarding High Energy Barrier in CI-NEB Calculation

To share experience including discussions about scientific questions.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
emmanuel_adejumo
Newbie
Newbie
Posts: 4
Joined: Thu Feb 22, 2024 2:11 pm

Inquiry Regarding High Energy Barrier in CI-NEB Calculation

#1 Post by emmanuel_adejumo » Wed Feb 26, 2025 9:23 pm

Dear VASP Support Team,

I am performing CI-NEB calculations to investigate Li diffusion in amorphous Mo₃S₁₃. I inserted lithium into the system and created initial and final structures, both of which converged to the required accuracy during optimization before the NEB calculation. However, my NEB results indicate a significantly high diffusion barrier, which appears questionable despite the absence of experimental data for direct comparison. Given the results, I suspect that something may be incorrect in my setup or calculation procedure.

Here are the key energy values from my NEB calculation:

00.009413-1482.7993000.00000010.325824-1008.260400474.53890020.372460-1245.674300237.12500030.357279-1394.22290088.57640040.261281-1406.74070076.05860050.008624-1482.7350000.064300
Given these values, I am concerned that the computed barrier may not be physically accurate. I would appreciate any insights or possible explanations regarding this issue.

Please find attached my input files for your reference.

Thank you for your assistance.

Best regards,
Emmanuel

You do not have the required permissions to view the files attached to this post.

ferenc_karsai
Global Moderator
Global Moderator
Posts: 530
Joined: Mon Nov 04, 2019 12:44 pm

Re: Inquiry Regarding High Energy Barrier in CI-NEB Calculation

#2 Post by ferenc_karsai » Thu Feb 27, 2025 9:15 am

I've moved this topic to "from users to users".


ferenc_karsai
Global Moderator
Global Moderator
Posts: 530
Joined: Mon Nov 04, 2019 12:44 pm

Re: Inquiry Regarding High Energy Barrier in CI-NEB Calculation

#3 Post by ferenc_karsai » Thu Feb 27, 2025 9:25 am

I have no experience with NEB (and the most of us at VASP also have no experience), so I can only look at your input files and see if I see something there.

The file looks reasonable. Also the plotted reaction barrier you attached looks smooth.
How was your electronic relaxation, did it always converge? You used ALGO=VeryFast and Prec=Low, maybe you could try ALGO=normal and PREC=normal and see how the results change.

Otherwise I hope maybe someone in the forum here with NEB experience can give you an advice.


emmanuel_adejumo
Newbie
Newbie
Posts: 4
Joined: Thu Feb 22, 2024 2:11 pm

Re: Inquiry Regarding High Energy Barrier in CI-NEB Calculation

#4 Post by emmanuel_adejumo » Thu Feb 27, 2025 1:12 pm

Thank you for taking a look at my input files and reaction barrier plot. My electronic relaxation generally converged, but I will check the convergence behavior more carefully. I'll also try running with ALGO=Normal and PREC=Normal to see if that affects the results.

I appreciate your feedback, and I'll wait to see if someone with more NEB experience can provide further insights. Let me know if you notice anything else that might be worth adjusting!


Post Reply