IDIPOL=3 and slab calculations

Problems running VASP: crashes, internal errors, "wrong" results.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
CStar
Newbie
Newbie
Posts: 12
Joined: Wed Jul 07, 2010 4:12 pm

IDIPOL=3 and slab calculations

#1 Post by CStar » Fri Jul 23, 2010 2:51 pm

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?
Last edited by CStar on Fri Jul 23, 2010 2:51 pm, edited 1 time in total.

CStar
Newbie
Newbie
Posts: 12
Joined: Wed Jul 07, 2010 4:12 pm

IDIPOL=3 and slab calculations

#2 Post by CStar » Wed Aug 11, 2010 11:36 pm

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.

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

IDIPOL=3 and slab calculations

#3 Post by admin » Thu Aug 26, 2010 1:44 pm

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.

Post Reply