Re: THG
Posted: Thu Dec 19, 2024 9:26 am
Dear Sanoj
sorry I do not agree with you, I use QE 7.3 and when I set a ibrav/=0 in the output I get something like
celldm(1)= 6.015395 celldm(2)= 0.000000 celldm(3)= 7.435501
celldm(4)= 7.435501 celldm(5)= 7.435501 celldm(6)= 7.435501
crystal axes: (cart. coord. in units of alat)
a(1) = ( 1.000000 0.000000 0.000000 )
a(2) = ( -0.500000 0.866025 0.000000 )
a(3) = ( 0.000000 0.000000 7.435501 )
reciprocal axes: (cart. coord. in units 2 pi/alat)
b(1) = ( 1.000000 0.577350 0.000000 )
b(2) = ( 0.000000 1.154701 0.000000 )
b(3) = ( 0.000000 0.000000 0.134490 )
this is exactly the cell. If you multiply the crystal axes for alat [ celldm(1)]
they are your new unit cell.
let me know what you get
Claudio
sorry I do not agree with you, I use QE 7.3 and when I set a ibrav/=0 in the output I get something like
celldm(1)= 6.015395 celldm(2)= 0.000000 celldm(3)= 7.435501
celldm(4)= 7.435501 celldm(5)= 7.435501 celldm(6)= 7.435501
crystal axes: (cart. coord. in units of alat)
a(1) = ( 1.000000 0.000000 0.000000 )
a(2) = ( -0.500000 0.866025 0.000000 )
a(3) = ( 0.000000 0.000000 7.435501 )
reciprocal axes: (cart. coord. in units 2 pi/alat)
b(1) = ( 1.000000 0.577350 0.000000 )
b(2) = ( 0.000000 1.154701 0.000000 )
b(3) = ( 0.000000 0.000000 0.134490 )
this is exactly the cell. If you multiply the crystal axes for alat [ celldm(1)]
they are your new unit cell.
let me know what you get
Claudio