Dear all,
the purpose of my calculations is to get dE/dk of charge carriers for specific 
band number in point (x,y,z).
For this after common wien calculations I use two methods (I do calculations 
for Bi2Te3 in symmetric point (x,x,x)):

1. I create klist file with points with coordinates:
        (x+dx,y,z)
        (x-dx,y,z)
        (x,y+dx,z)
        (x,y-dx,z)
        (x,y,z+dx)
        (x,y,z-dx)
   Do lapw1 (and lapwso), and get E(...)'s
   Calculate dE/dk: dE/dk_x = (E(x+dx,y,z)-E(x-dx,y,z))/2dx and e.t.c.

2. I clear klist file and put there point with coordinate (x,y,z)
   Do optic (with pmat flag)

I have few problems:
1. In first method results dramatically depend on dx, for dx=1e-6 and dx=1e-7 
results differ by two times... no convergence.
2. I get completely different results by this two methods.
I understand that I could lose some factor, but in first method I have 
symmetric vector, like (x,x,x) as expected, in second I have not symmetric 
result, I get vector like (x,y,z).

I have few questions:
1. What method is right (optimal)?
2. Why results are not equalt? 
3. Why I haven't convergence in first method?
4. Why in second method I receive non symmetric vector? (Even if result is in 
cartesian coordinates I expect vector like (0,0,x) for this point in Bi2Te3).

Hope somebody can help me or give advice.
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html

Reply via email to