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

Hongsheng Zhao zhaohscas at yahoo.com.cn
Tue Jun 28 15:22:23 CEST 2011


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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: graphene.inp
Type: chemical/x-gamess-input
Size: 731 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110628/b14cdcfd/attachment.bin>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: graphene.pwscf
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110628/b14cdcfd/attachment.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: supportInfo.kpath
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110628/b14cdcfd/attachment-0001.ksh>


More information about the users mailing list