[Pw_forum] The issue of k-points in CARTESIAN coordinates obtained by XCrySDen-1.5.24.

Hongsheng Zhao zhaohscas at yahoo.com.cn
Wed Jun 29 02:08:21 CEST 2011


On 06/28/2011 09:45 PM, jia chen wrote:
> Hi Hongsheng,
>
> The difference is 2*pi, I guess it is just a matter of convention.

Thanks a lot, I've got it.

>
> Regards
> jia
>
> On Tue, Jun 28, 2011 at 9:22 AM, Hongsheng Zhao<zhaohscas at yahoo.com.cn>  wrote:
>> Dear Prof. Tone and others,
>>
>> I use pre-compiled xc-1.5.24-linux_x86-semishared.tar.gz on my Debian box.
>>
>> Now I use the k-path selection function to prepare k-point file used for
>> pwscf.  My system is graphene.  The graphene.inp file is as follows:
>>
>> -------------
>>      &CONTROL
>> calculation = 'scf'
>> prefix='graphene'
>> pseudo_dir = './'
>> outdir='./'
>> tstress=.true.
>> tprnfor = .true.
>> forc_conv_thr=1.0d-4
>> nstep=200
>> /
>>      &SYSTEM
>> ibrav=0
>> celldm(1)=1.889726
>> nat=2
>> ntyp=1
>> ecutwfc=36
>> /
>>      &ELECTRONS
>> mixing_beta=0.7
>> conv_thr =  1.0d-8
>> electron_maxstep=200
>> /
>>      &IONS
>> trust_radius_max=0.2
>> /
>>      &CELL
>> cell_dynamics='bfgs'
>> /
>> ATOMIC_SPECIES
>>     C  12.0110   C.UPF
>>
>> CELL_PARAMETERS (alat)
>>     2.1314617179   -1.2306000020    0.0000000000
>>     0.0000000000    2.4612000040    0.0000000000
>>     0.0000000000    0.0000000000   36.9180000600
>>
>> ATOMIC_POSITIONS (crystal)
>>   C       0.000000000     0.000000000     0.000000000
>>   C       0.666666671     0.333333335     0.000000000
>> K_POINTS {automatic}
>>   4 4 4  1 1 1
>> -------------
>>
>> I select the k-path as M -->  G -->K -->  M , and then generate the
>> corresponding files, i.e.,  graphene.pwscf and supportInfo.kpath.
>>
>> In the file supportInfo.kpath, you can find the following information:
>>
>> -------------
>>     Selected k-points in CARTESIAN coordinates:
>>        0.23458    0.00000    0.00000    M
>>        0.00000    0.00000    0.00000    G
>>        0.23458    0.13543    0.00000    K
>>        0.23458    0.00000    0.00000    M
>>
>>     Selected k-points in crystal coordinates:
>>     (with respect to RECIPROCAL-PRIMITIVE vectors)
>>        0.50000    0.00000    0.00000    M
>>        0.00000    0.00000    0.00000    G
>>        0.33333    0.33333    0.00000    K
>>        0.50000    0.00000    0.00000    M
>>
>>     Selected k-points in crystal coordinates:
>>     (with respect to RECIPROCAL-CONVENTIONAL vectors)
>>        0.50000    0.00000    0.00000    M
>>        0.00000    0.00000    0.00000    G
>>        0.33333    0.33333    0.00000    K
>>        0.50000    0.00000    0.00000    M
>> ---------------
>>
>> Here, I have some puzzles on the k-points in CARTESIAN coordinates obtained
>> by XCrySDen.  For example, the K point in BZ has the fractional coordinates
>> as this: ( 0.33333    0.33333    0.00000 ).
>>
>> And in my case, I've the lattice vectors as follows:
>>
>> ----------
>> CELL_PARAMETERS (alat)
>>     2.1314617179   -1.2306000020    0.0000000000
>>     0.0000000000    2.4612000040    0.0000000000
>>     0.0000000000    0.0000000000   36.9180000600
>> -----------
>>
>> So, I can obtain the reciprocal lattice like this:
>>
>> B ^ T = 2*Pi* A ^ (-1),  where, A and B are 3*3 matrix whose column vectors
>>   are corresponding to the real space lattice vectors and reciprocal lattice
>> vectors, respectively.
>>
>> Thus, I obtain B as follows from within matlab:
>>
>>   >>    B= 2*pi*inv([2.1314617179   -1.2306000020    0.0000000000;
>>     0.0000000000    2.4612000040    0.0000000000;
>>     0.0000000000    0.0000000000   36.9180000600])'
>>
>> B =
>>
>>      2.9478         0         0
>>      1.4739    2.5529         0
>>           0         0    0.1702
>>
>> So the CARTESIAN coordinates for K point ( 0.33333    0.33333    0.00000 )
>> should be:
>>
>>   >>  [0.33333    0.33333    0.00000]*B
>>
>> ans =
>>
>>      1.4739    0.8510         0
>>
>> As you can see, your XCrySDen gives the CARTESIAN coordinates for K point (
>> 0.33333    0.33333    0.00000 ) as follows:
>>
>>        0.23458    0.13543    0.00000    K
>>
>> I cann't figure out what's the reason.  Could you please give me some hints
>> on this issue?
>>
>> Attached please find all the files mentioned above.
>>
>> Regards.
>> --
>> Hongsheng Zhao<zhaohscas at yahoo.com.cn>
>> School of Physics and Electrical Information Science,
>> Ningxia University, Yinchuan 750021, China
>>
>> _______________________________________________
>> Pw_forum mailing list
>> Pw_forum at pwscf.org
>> http://www.democritos.it/mailman/listinfo/pw_forum
>>
>>
>
>
>


-- 
Hongsheng Zhao <zhaohscas at yahoo.com.cn>
School of Physics and Electrical Information Science,
Ningxia University, Yinchuan 750021, China



More information about the users mailing list