[Pw_forum] PWscf can not find out symmetry operations!!!
英利 牛
niuyingli at yahoo.com.cn
Tue Aug 29 18:16:29 CEST 2006
Dear all,
I am calculate properties of diamond
surface(001),the lattice parameters of the slab is :
2.5221A,5.0442A,16.0253A,90,90,90
after atoms relax ,the atom positions(crystal) are:
C 0.24978 -0.11238 0.21129
C 0.24980 0.60988 0.21129
C 0.74978 -0.01422 0.17583
C 0.74980 0.51175 0.17582
C 0.74981 0.24879 0.13029
C 0.74981 0.74881 0.11624
C 0.24982 0.24883 0.07780
C 0.24982 0.74885 0.07061
C 0.24984 -0.00014 0.02461
C 0.24984 0.49787 0.02461
C 0.74985 -0.00211 -0.02461
C 0.74986 0.49988 -0.02461
C 0.74987 0.24889 -0.07061
C 0.74987 0.74888 -0.07780
C 0.24987 0.24889 -0.11624
C 0.24988 0.74888 -0.13029
C 0.24990 0.01186 -0.17582
C 0.24987 0.48589 -0.17583
C 0.74990 0.11000 -0.21129
C 0.74987 0.38773 -0.21129
but PWscf can not find any symmetry operations of the
system.
I setup the struct and use Material Studio find 8
space group operations include inversion! The
tolerance is 0.01A.
Then I read PW manual:
the number of significant figures in the atomic
positions is not large enough. In file PW/eqvect.f90,
the variable accep is used to decide whether a
rotation is a symmetry operation. Its current value
(10-5 ) is quite strict: a rotated atom must coincide
with another atom to 5 significant digits. You may
change the value of accep and recompile.
I change accep from 1.0d-5 to 1.0d-2 then to 1.0d-1
,but PWscf still find no symmetry.
What should I do?
Why pwscf choose 1.0d-5,so strict?
___________________________________________________________
Mp3疯狂搜-新歌热歌高速下
http://music.yahoo.com.cn/?source=mail_mailbox_footer
More information about the Pw_forum
mailing list