[Wien] error while doing LDA+U calculations

M G madhura1802 at gmail.com
Sat Apr 18 09:32:06 CEST 2009


Hello all,
I am running wien2k08.2 (release21/4/2008) on a intel xeons with operating
system linuxif9, fortran compiler  10 and math libraries 10.
I am facing some problems while doing calculations for GGA+U with Spin Orbit
coupling for pure ZnO and then Cu doped ZnO.
with rkmax 6 and tetra method.
The procedure I am following for wortzite ZnO unit cell (2 atoms with mult
2) is:
1. run normal GGA scf.
2. Then apply U = 0.5 Ry for Zn 3d and no RLO but mag field along Z
direction in SO.
   I am attaching the inorb inso and indmc along with.
3. Then run scf with "runsp -fc 1 -orb -so".
4. Then calculate DOS or SO (lapw1 -orb up/dn ; lapwso -orb up; lapw2 -so
   up/dn; and then tetra / spaghetti -so up/dn)

The problem is Zn 3d (observed around 5-6 eV) instead of shifting downwards
in E shift upwards and thus pushing O-p up, making ZnO metallic.
Now if I apply a U = -0.5 Ry (-ve ?) Zn 3d shift downwards but Ef is
shifting down as well making conduction bands cross the Ef.
Also the shift observed (hardly 1eV)is not as much as expected in published
papers. Also no band gap opening is seen even though I am using SIC method.

I had faced same problem for GGA+U without SO as well.

I  have searched wien-mailing list and userguide  as well but could not
figure out the possible reason for this strange behavior. Can you tell
probably where I am going wrong.

   With Best Regards,
   Madhura

***************************************************************************
Ms. Madhura Gokhale.
Research Student
Department of Physics,
University of Pune                          Tel.(O)+91-20-25692678 ext 429
Pune-411007                                  Fax.(O)+91-20-25691684
MAHARASHTRA, INDIA.             Email : madhura1802 at yahoo.com
***************************************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20090418/191fdb06/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ZnO-2atom.indmc
Type: application/octet-stream
Size: 211 bytes
Desc: not available
Url : http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20090418/191fdb06/ZnO-2atom.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ZnO-2atom.inorb
Type: application/octet-stream
Size: 286 bytes
Desc: not available
Url : http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20090418/191fdb06/ZnO-2atom-0001.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ZnO-2atom.inso
Type: application/octet-stream
Size: 416 bytes
Desc: not available
Url : http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20090418/191fdb06/ZnO-2atom-0002.obj


More information about the Wien mailing list