EFIELD and problem with convergence

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
spakinfo
Newbie
Newbie
Posts: 19
Joined: Mon Sep 29, 2008 3:10 am

EFIELD and problem with convergence

#1 Post by spakinfo » Wed Dec 15, 2010 8:38 pm

Dear all,

I'm trying to do electronic structure calculations for a slab under a finite electric field. My system has hexagonal symmetry, and I use the following tags in INCAR to switch on a saw-like electric field of 0.3 eV/A along c axis:

ISMEAR =0
NSW = 1
IBRION = -1
ISPIN = 1
EDIFF = 0.00001
EDIFFG = -0.005
ICHARG = 2
ISYM=0
LMAXMIX=6
ENCUT=500
EFIELD =0.3
LDIPOL=.TRUE.
IDIPOL= 3
NELM=200
PREC=A

With the above settings, however, electronic convergence seems to be out of reach. Without them, on the other hand, there seems to be no such a problem. Now I'm wondering, If this issue is due to
1) The hexagonal symmetry of my structure
2) The amplitude of electric field
3) mixing parameters

Any help is greatly appreciated.
Spak
Last edited by spakinfo on Wed Dec 15, 2010 8:38 pm, edited 1 time in total.

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

EFIELD and problem with convergence

#2 Post by admin » Thu Dec 30, 2010 1:44 pm

please first of all try to fix the origin of the dipol. Usually the bad electronic convergence is caused by a permanent shift of the origin of the dipolemoment
Last edited by admin on Thu Dec 30, 2010 1:44 pm, edited 1 time in total.

spakinfo
Newbie
Newbie
Posts: 19
Joined: Mon Sep 29, 2008 3:10 am

EFIELD and problem with convergence

#3 Post by spakinfo » Thu Dec 30, 2010 6:18 pm

Thanks admin.
Does it matter where to fix the dipole? Should I fix it somewhere inside the slab (e.g. at the center) or elsewhere in the vacuum region?
Last edited by spakinfo on Thu Dec 30, 2010 6:18 pm, edited 1 time in total.

Post Reply