Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] problem with elastic properties calculation

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] problem with elastic properties calculation


Chronological Thread 
  • 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.

Top of Page