ENCUT Convergence
Moderators: Global Moderator, Moderator
-
Mausumi
- Newbie
- Posts: 3
- Joined: Wed Feb 15, 2023 4:56 am
ENCUT Convergence
I am a beginner in VASP, please help me in clarifying the error .We are trying to do dft calculations on N,S dual doped graphene in vasp. The ENCUT values are not getting converged for set of values ranging from 200 to 650 eV. We tried with Monkhorst pack grid 2*2*1 and Gamma centered pack grid 13*13*1 (as per the literature). Could you please look into this and give some appropriate suggestions?
You do not have the required permissions to view the files attached to this post.
-
alex
- Hero Member
- Posts: 585
- Joined: Tue Nov 16, 2004 2:21 pm
- License Nr.: 5-67
- Location: Germany
Re: ENCUT Convergence
Hello Mausumi,
if I got your question correctly you are after basis set convergence wrt. basis set size?!
In that case you need to know that VASP extrapolates to the basis set limit by adding a cut-off dependend energy per atom to correct for basis set incompleteness. You'll find this in your(!) OUTCAR.
PAW_PBE C 08Apr2002 :
energy of atom 1 EATOM= -147.1560
kinetic energy error for atom= 0.0854 (will be added to EATOM!!)
PAW_PBE S 06Sep2000 :
energy of atom 2 EATOM= -276.8230
kinetic energy error for atom= 0.0051 (will be added to EATOM!!)
PAW_PBE N 08Apr2002 :
energy of atom 3 EATOM= -264.5486
kinetic energy error for atom= 0.2274 (will be added to EATOM!!)
If you remove these energies you'll find a nice behaviour E_tot = f(ENCUT).
Hth,
alex
if I got your question correctly you are after basis set convergence wrt. basis set size?!
In that case you need to know that VASP extrapolates to the basis set limit by adding a cut-off dependend energy per atom to correct for basis set incompleteness. You'll find this in your(!) OUTCAR.
PAW_PBE C 08Apr2002 :
energy of atom 1 EATOM= -147.1560
kinetic energy error for atom= 0.0854 (will be added to EATOM!!)
PAW_PBE S 06Sep2000 :
energy of atom 2 EATOM= -276.8230
kinetic energy error for atom= 0.0051 (will be added to EATOM!!)
PAW_PBE N 08Apr2002 :
energy of atom 3 EATOM= -264.5486
kinetic energy error for atom= 0.2274 (will be added to EATOM!!)
If you remove these energies you'll find a nice behaviour E_tot = f(ENCUT).
Hth,
alex