[Wien] question on how to rotate the axis in qtl propram
王静
wangjingjing at ciac.jl.cn
Thu Mar 1 08:19:04 CET 2012
Hi, there,
I have a system in space group P-1 (No. 2) and I want to calculate the partial DOS for d orbitals by qtl program. I would like a local coordinate with axis pointing towards the oxygens in the octahedron instead of the defaulted. What I know is the value of “locrot” in the .inq file should be set to 2. What I want is the xy plane to change 45° from the basal plane. Then I set the z axis along [0, 0, 1] which I though the default value since I would not like it to be changed. Meanwhile x axis is set to along [1,1,0]. My .inq file is as the following:
___________________________________
-9.0 3.0 Emin Emax
12 number of atoms
1 -2 0 2 iatom,qsplit,symmetrize,locrot
4 0 1 2 3 nL, l-values
2 -2 0 2 iatom,qsplit,symmetrize,locrot
4 0 1 2 3 nL, l-values
3 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
4 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
5 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
6 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
7 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
8 -2 0 2 iatom,qsplit,symmetrize,locrot
2 0 1 nL, l-values
9 -2 0 2 iatom,qsplit,symmetrize,locrot
3 0 1 2 nL, l-values
10 -2 0 2 iatom,qsplit,symmetrize,locrot
3 0 1 2 nL, l-values
11 -2 0 2 iatom,qsplit,symmetrize,locrot
3 0 1 2 nL, l-values
12 -2 0 2 iatom,qsplit,symmetrize,locrot
3 0 1 2 nL, l-values
0. 0. 1. new axis z
1. 1. 0. new axis x
___________________________________
When “x qtl -up -p” is executed, it generated the following message:
___________________________________
calculating QTL's from parallel vectors
new x and z axes are not orthogonal for atom 1
angle: 48.4255235866551 degrees
x,z vectors not orthogonal
0.025u 0.043s 0:00.10 60.0% 0+0k 0+0io 0pf+0w
___________________________________
My question is
(1) I doubt the default lattice vector is [1,0,0] for x axis, [0,1,0] for y axis and [0,0,1] for z axis. Am I right? If not, which file or where could I get the related information?
(2) Any suggestion about the solutions for the case mentioned above would be appreciated. The structure parameters are a=10.535601, b=10.631599, c=14.925057, α=90.02, β=90.00, γ= 90.00.
Thanks in advance!
Wang J.
wangjingjing at ciac.jl.cn
2012-03-01
More information about the Wien
mailing list