error in gw

Moderators: gmatteo, bruneval, rigna

Post Reply
yaser
Posts: 19
Joined: Wed Mar 10, 2010 6:03 am

error in gw

Post by yaser » Wed Mar 10, 2010 6:30 am

Dear Abiniters,
I want to calculate the band gap with GW approximation (first part of lesson
one ) for a nonsymmorphic material whit space group 206 .
The question is : when i get symmorphic=0 i will face with error in log file:
( Simple Lattice Grid
symkpt : found identity, with number 1

inread : ERROR -
Attempted to read ndig= 10 integer digits,
from string(1:ndig)= OPTDRIVER2, to initialize an integer variable

inread : WARNING -
Note that this string contains the letter O.
It is likely that this letter should be replaced by the number 0.

inarray : ERROR -
It occurred reading data for keyword " ISTWFK1

",
looking for 12 array elements.
There is a problem with the input file : maybe
a disagreement between the declared dimension of the array,
and the number of data actually provided.
Action : correct your input file !

leave_new : decision taken to exit ...1).


-and when symmorphic=1 : ( chkinp: Checking input parameters for consistency,
jdtset= 1.

chkinp : ERROR -
When kssform==1 or 2, one cannot use pseudopotentials with
more than one projector per angular momentum channel.
Note that the ABINIT Troullier-Martins psps (pspcod=6),
are suitable for use with kssform==1 or 2.
Action : change kssform, nbandkss, or the problematic pseudopotential.

chkinp: ERROR -
When nbandkss/=0, all the components of tnons must be zero.
However, for the symmetry operation number 3, tnons = 0.5000 0.0000
0.5000.
Action : use the symmetry finder (nsym=0) with symmorphi==0.

leave_new : decision taken to exit ...1).

My Dear Abiniters ,please help me .
I am wating to hear you as soon as possible.
tank you
asadi.

gonze
Posts: 391
Joined: Fri Aug 14, 2009 8:29 pm

Re: error in gw

Post by gonze » Wed Mar 10, 2010 7:40 pm

Hi,
There must be a typographical error in your input file.
Could you post it ?
X.

yaser
Posts: 19
Joined: Wed Mar 10, 2010 6:03 am

Re: error in gw

Post by yaser » Thu Mar 11, 2010 7:35 am

Dear gonze,
Hi, Thank you for your reply. the input file is:
--------------------------------------------------------------------------------------------------------------------------

ndtset 3


# Definition of parameters for the calculation of the KSS file

nbandkss1 -1 # Number of bands in KSS file (-1 means the maximum possible)
nband1 100 # Number of (occ and empty) bands to be computed
istwfk1 10*1

# Calculation of the screening (epsilon^-1 matrix)
optdriver2 3 # Screening calculation
getkss2 -1 # Obtain KSS file from previous dataset
nband2 150 # Bands to be used in the screening calculation
ecutwfn2 5 # Cut-off energy of the planewave set to represent the wavefunctions
ecuteps2 6.6 # Cut-off energy of the planewave set to represent the dielectric matrix
ppmfrq2 16.8 # Imaginary frequency where to calculate the screening

# Calculation of the Self-Energy matrix elements (GW corrections)
optdriver3 4 # Self-Energy calculation
getkss3 -2 # Obtain KSS file from dataset 1
getscr3 -1 # Obtain SCR file from previous dataset
nband3 200 # Bands to be used in the Self-Energy calculation
ecutwfn3 8 # Planewaves to be used to represent the wavefunctions
ecutsigx3 9 # Dimension of the G sum in Sigma_x
# (the dimension in Sigma_c is controlled by npweps)
nkptgw3 1 # number of k-point where to calculate the GW correction
kptgw3 # k-points
0.000 0.000 0.000
bdgw3 96 97 # calculate GW corrections for bands from 4 to 5


# Data common to the three different datasets

acell 3*18.852
ntypat 2
znucl 49 8
natom 40
natrd 40
spgroup 206

brvltt -1
typat 16*1 24*2
xred
0.25000000 0.25000000 0.25000000
0.25000000 0.75000000 0.75000000
0.75000000 0.25000000 0.75000000
0.25000000 0.25000000 0.75000000
...
...
...
...

kptopt 1

nshiftk 1
shiftk 0.5 0.5 0.5

ngkpt 4 4 4
# Use only symmorphic operations
symmorphi 0




# Definition of the planewave basis set (at convergence 16 Rydberg 8 Hartree)
ecut 16 # Maximal kinetic energy cut-off, in Hartree

# Definition of the SCF procedure
nstep 1500 # Maximal number of SCF cycles
toldfe 1.0d-6 # Will stop when this tolerance is achieved on total energy
diemac 12.0 # Although this is not mandatory, it is worth to
# precondition the SCF cycle. The model dielectric
# function used as the standard preconditioner
# is described in the "dielng" input variable section.
# Here, we follow the prescription for bulk silicon.

mverstra
Posts: 646
Joined: Wed Aug 19, 2009 12:01 pm

Re: error in gw

Post by mverstra » Thu Mar 11, 2010 11:11 am

> istwfk1 10*1

is the problem, and trying to read into the next fields looking for data.
You have more than 10 k-points, so adjust the 10.
On most platforms you can use
istwfk *1
to tell the parser to fill all needed values with 1.

Matthieu
Matthieu Verstraete
University of Liege, Belgium

yaser
Posts: 19
Joined: Wed Mar 10, 2010 6:03 am

Re: error in gw

Post by yaser » Thu Mar 11, 2010 4:01 pm

Dear mverstra,
Hi,thank you for your reply.
When i let istwfk*1 in input fille then i have the following error.
----------------------------------------------------------------------------------------------------------
chkinp: machine precision is 2.2204460492503131E-16

chkinp: Checking input parameters for consistency, jdtset= 1.

chkinp : ERROR -
When kssform==1 or 2, one cannot use pseudopotentials with
more than one projector per angular momentum channel.
Note that the ABINIT Troullier-Martins psps (pspcod=6),
are suitable for use with kssform==1 or 2.
Action : change kssform, nbandkss, or the problematic pseudopotential.

chkinp: Checking input parameters for consistency, jdtset= 2.

chkinp: Checking input parameters for consistency, jdtset= 3.

leave_new : decision taken to exit ...
1
----------------------------------------------------------
thank you.
asadi.

User avatar
gmatteo
Posts: 269
Joined: Sun Aug 16, 2009 5:40 pm

Re: error in gw

Post by gmatteo » Thu Mar 11, 2010 4:26 pm

Dear Asadi,

you are probably using an old version of the code.
Could you try the latest stable version available on the web site?
Moreover I suggest you use kssform=3 for generating the KSS file
instead of the direct diagonalization.
This topic has been already discussed in this forum.

An additional comment: the treatment of the optical limit (q-->0) of the
heads and of the wings of the dielectric matrix is not coded
when at least one of the NC pseudos has more than one projector per angular channel.
GW results in bulk systems should tend to the correct value when the
q-sampling of the BZ is increased as the contribution given by q=\Gamma will decrease.
On the other hand the optical spectra produced in the screening part might be affected
by important errors when multi-projector pseudopotentials are used.

Hope it helps,
Matteo

yaser
Posts: 19
Joined: Wed Mar 10, 2010 6:03 am

Re: error in gw

Post by yaser » Thu Mar 11, 2010 6:47 pm

Dear gonze,Matthieu and Matteo,

Thanks for your reply. That's helpful!

Best regards,

Asadi.

Post Reply