forum@abinit.org
Subject:
The ABINIT Users Mailing List ( CLOSED )
List archive
- From: Xenophon Krokidis <Xenophon.Krokidis@scienomics.com>
- To: forum@abinit.org
- Subject: Re: [abinit-forum] problem with elastic properties calculation
- Date: Thu, 05 Nov 2009 14:44:52 +0100
- Organization: Scienomics
Title: Signature text
Hi Don,
thanks a lot, I completely missed this one !!!!
Xenophon
Xenophon Krokidis, PhD
Scienomics
T: + 33 (0)1 40 07 56 60
F: + 33 (0)1 71 19 75 83
M: + 33 (0)6 76 68 06 47 |
D. R. Hamann wrote:
4AEF1E87.3070609@mat-simresearch.com" type="cite">Dear
Xenophon,
This is a feature -- not an error. Your shiftk set is wrong for fcc.
Look at the input file for exercise 6 of the tutorial on elastic
properties. Alternatively, you can use nshiftk 1, shiftk 0 0 0, and a
bigger ngkpt set for comparable convergence.
Don Hamann
xenophon.krokidis@scienomics.com wrote:
Dear All,
I am using the input below to compute elastic properties but I get the
following error message when the 3rd dataset is executed (note that if
I use
nshiftk = 1 and shiftk 0.0 0.0 0.0 I don't get this error).
Thak you for the help.
Xenophon
ERROR<<<<<
symkchk : found identity, with number 1
symkchk : ERROR -
k-point set must have full space-group symmetry
there is no match for kpt 1 transformed by symmetry 3
Action : change kptopt to 2 or 3 and/or change or use shiftk
shiftk = 0 0 0 is always a safe choice.
leave_new : decision taken to exit ...
INPUT<<<<<
# Diamond
acell 3.5668 3.5668 3.5668 Angstr
rprim 0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
chkprim 0
ecut 150.0 eV
ecutsm 13.6057 eV
iscf 5
nkpt 0
kptopt 1
ngkpt 4 4 4
nshiftk 4
shiftk 0.500000 0.500000 0.500000
0.500000 0.000000 0.000000
0.000000 0.500000 0.500000
0.000000 0.000000 0.500000
nstep 50
ixc 1
nband 6
nsppol 1
natom 2
natrd 2
ntypat 1
occopt 1
prtden 1
typat 1 1
xred 0 0 0
0.25 0.25 0.25
znucl 6
# Calculation of elastic constants
iscf1 5
kptopt1 1
tolvrs1 1.0d-18
getwfk2 -1
iscf2 -3
kptopt2 2
nqpt2 1
qpt2 0.0 0.0 0.0
rfelfd2 2
rfdir2 1 1 1
tolwfr2 1.0d-20
getddk3 -1
getwfk3 -2
iscf3 5
kptopt3 2
nqpt3 1
qpt3 0.0 0.0 0.0
rfphon3 1
rfatpol3 1 2
rfstrs3 3
rfdir3 1 1 1
tolvrs3 1.0d-10
ndtset 3
|
begin:vcard
fn:Xenophon Krokidis
n:Krokidis;Xenophon
org:Scienomics
adr:;;17 square Edouard VII;Paris;;75009;France
email;internet:Xenophon.Krokidis@Scienomics.com
title:Managing Director
tel;work:+33 (0)1 53 43 51 05
tel;fax:+ 33 (0)1 53 43 92 92
tel;cell:+33 (0)6 76 68 06 47
version:2.1
end:vcard
Archive powered by MHonArc 2.6.16.