manoj at phys.ufl.edu
Wed Apr 22 05:44:30 CEST 2009
Dear PWSCF users,
I have a question regarding real band structure calculated from
PWCOND. Complex band code requires kx,ky and E as input and gives kz
as output. I have noticed that the answer for kz depends whether you have
unit cell vectors as (a1,a2,a3) or (a1,a2,-a3). I did a calculation for
Cu(001) with 2 atoms tetragonal unit cell, with the following input file -
pseudo_dir = '/home/manoj/espresso-4.0.1/pseudo',
ibrav = 0,
conv_thr = 1.0e-8
Cu 63.55 Cu.pz-d-rrkjus.UPF
Cu 0.0 0.0 0.0
Cu 0.5 0.5 0.707106781
8 8 8 1 1 1
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.414213562
0.0 0.0 1.0
This gives me reasonable answer and it matches up with existing results.
Now when I changed a3 to -a3, so that the unit cell now has become
(a1,a2,-a3), the answer does not match up with the previous one.
Just to give you one concrete example-
at Fermi energy for kx,ky (0,0) for unit cell (a1,a2,a3) we get kz=-0.1571
while for (a1,a2,-a3) we get kz=0.07323.
I understand that for two different set of unit cell vectors we should not
in general expect k_z to match, but this is a very special example, as it
is just mirror symmetry about xy plane. so shouldn't kz in one set have
value -kz in other?
Thanks for your attention.
Department of Physics
University of Florida
More information about the users