Personal tools
You are here: Home Support Help Center FAQs kptgen : accuracy of atomic positions

kptgen : accuracy of atomic positions

Up to Table of Contents

When I used KPTGEN today (May 27, 1999), I had a curious experience.
(In fact, no curious... Only I might have not realized the following.)

If I use the following input for "kptgen", it fails.

 


Example unit cell

7.1431642 7.1431642 12.402272 acell
0.8660254037844386 0.5 0 -0.8660254037844386 0.5 0 0 0 1 rprim
5 3 natom, ntype
1 2 2 3 3
0.00000 0.00000 0.00000
0.33333 0.66667 0.29000
0.66667 0.33333 0.71000
0.33333 0.66667 0.64000
0.66667 0.33333 0.36000
0 0 ibrav, inpgd
0 nkpt
4 4 4 0 0 0 Monkhorst-Pack indices and origin for grid

 


The error message in standard output is as follows:
cpointg error: atom 2 doesn't fit symmetry 2
rotated coord -2.062033266329 -3.571617815821 3.596658880000
other atom 2.062095127946 3.571582100000 8.805613120000
difference -4.124128394275 -7.143199915821 -5.208954240000
phases -.500036275549 -.499981861897 -.876306680044
rotation matrix: -.50 -.87 .00 .87 -.50 .00 .00 .00 1.00
nonsymmorphic transl: .000000000000 .000000000000 .000000000000

 


But If I use the following input (the digits of "xred" are increased),

it did work !

 


Example unit cell

7.1431642 7.1431642 12.402272 acell
0.8660254037844386 0.5 0 -0.8660254037844386 0.5 0 0 0 1 rprim
5 3 natom, ntype
1 2 2 3 3
0.000000000 0.000000000 0.000000000
0.333333333 0.666666667 0.290000000
0.666666667 0.333333333 0.710000000
0.333333333 0.666666667 0.640000000
0.666666667 0.333333333 0.360000000
0 0 ibrav, inpgd
0 nkpt
4 4 4 0 0 0 Monkhorst-Pack indices and origin for grid

 


Thus I really realized that the digits of "xred" is very important and that enough digits of "xred" are prepared in all the input examples

(in.*). (The criterion condition (the parameter of "tolerance") in cpointg.f seems to be set severely.)

+-------------------------------+
Masayoshi Mikami
Mitsubishi Chemical Corporation (Research and Development Division, Yokohama Research Center)
Phone: +81-45-963-3265 Fax:+81-45-963-3947
E-mail: mmikami@rc.m-kagaku.co.jp (office), mmikami@yk.rim.or.jp (home)
+-------------------------------+

Former website

http://wwwold.abinit.org

« July 2010 »
July
MoTuWeThFrSaSu
1234
567891011
12131415161718
19202122232425
262728293031
Site status
Stable