I am running some calculations with a charged adsorbate above a slab, and I am wondering about how to use the DIPOL flags correctly (...or if I need to). The slab is in the a-b direction, so I think that IDIPOL=3 is the flag that I need. The part that I am not sure about is how to set the DIPOL=X X X flag. The center of mass of the system is positioned at 0.5 0.5 0.5.
I have tried both setting DIPOL to be the COM of the adsorbate, and to be the COM of the whole system of slab+adsorbate. In both cases there is a very big difference between corrected and uncorrected energies ( >20 eV with a vacuum slab of 20 \AA ); looking at the OUTCAR file the trace of the quadrupole is very large: ~900, and this appears to be the reason for the big correction. Since my system is charged, I assume the DIPOL corrections are necessary. However, the un-corrected energies in fact decrease linearly with the inverse of the vacuum slab (over a range of 10-45 \AA), so is it unreasonable to simply extrapolate these values to 1/(vacuum slab) --> 0? Presumably the corrected and uncorrected values ultimately converge to the same number (if DIPOL is set correctly), so if the uncorrected value is converged it is OK to use this number for the total energy?
IDIPOL=3 and slab calculations
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 12
- Joined: Wed Jul 07, 2010 4:12 pm
IDIPOL=3 and slab calculations
Last edited by CStar on Fri Jul 23, 2010 2:51 pm, edited 1 time in total.
-
- Newbie
- Posts: 12
- Joined: Wed Jul 07, 2010 4:12 pm
IDIPOL=3 and slab calculations
OK so the problem was that calculations with charged adsorbates cannot be done, as it mentions in the last point in http://cms.mpi.univie.ac.at/vasp/vasp/FAQ.html
Last edited by CStar on Wed Aug 11, 2010 11:36 pm, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
IDIPOL=3 and slab calculations
you are correct, the answer to your question is given in the manual's FAQs
Last edited by admin on Thu Aug 26, 2010 1:44 pm, edited 1 time in total.