Error on k+q does not fall on k-grid

Post here questions linked with issue while running the EPW code

Moderator: stiwari

Post Reply
jameszt
Posts: 1
Joined: Wed Aug 17, 2022 11:41 am
Affiliation: The University of Adelaide

Error on k+q does not fall on k-grid

Post by jameszt »

Dear all,
I wanted to use EPW to calculate the carrier mobility of graphene. However, I meet some problems on using EPW v5.3.1, since this is the version on our school HPC. The most recent problem is called:
Error in routine ktokpmq (1):
k+q does not fall on k-grid
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Can anyone help me resolve this problem, appreciate!
This is my scf.in
&control
calculation='scf',
prefix='graphene',
restart_mode='from_scratch'
pseudo_dir = './',
outdir='./',
wf_collect =.true.
!verbosity=.high.
tprnfor = .true.
tstress = .true.
/
&system
ibrav= 0,
a = 4.64872629,
nat= 2,
ntyp= 1,
ecutwfc = 50.0,
occupations= 'smearing'
smearing= 'mp'
degauss= 0.01
nbnd = 4
/
&electrons
diagonalization='david'
conv_thr = 1.0d-10
mixing_beta = 0.7
/
ATOMIC_SPECIES
C 12.0107 C.pbesol-n-kjpaw_psl.1.0.0.UPF
CELL_PARAMETERS (alat= 4.64872629)
1.002304867 0.000000000 0.000000000
-0.501152434 0.868021478 0.000000000
0.000000000 0.000000000 8.130081301
ATOMIC_POSITIONS (crystal)
C 0.000000 0.000000 0.25000000
C 0.333333 0.666667 0.25000000
K_POINTS {automatic}
6 6 6 0 0 0

This is my ph.in
--
&inputph
prefix = 'graphene',
epsil = .false.,
fildyn = 'graphene.dyn.xml',
ldisp = .true.
fildvscf = 'dvscf'
nq1=3,
nq2=3,
nq3=3,
tr2_ph = 1.0d-12,

This is my nscf.in
/
&control
calculation='nscf'
prefix='graphene',
pseudo_dir = './',
outdir='./'
tstress = .true.
tprnfor = .true.
verbosity = 'high'
wf_collect = .false.
/
&system
ibrav= 0,
a = 4.6595,
nat= 2,
ntyp= 1,
ecutwfc = 50.0,
occupations = 'smearing'
smearing = 'mp'
degauss = 0.01
nbnd = 10
/
&electrons
conv_thr = 1.0d-8
mixing_beta = 0.7
/
ATOMIC_SPECIES
C 12.0107 C.pbesol-n-kjpaw_psl.1.0.0.UPF
CELL_PARAMETERS (alat= 4.64872629)
1.000000000 -0.000000000 0.000000000
-0.500000000 0.8670000000 0.000000000
0.000000000 0.000000000 8.130081301
ATOMIC_POSITIONS (crystal)
C 0.000000 0.000000 0.25000000
C 0.333333 0.666667 0.25000000
K_POINTS crystal
216
0.0000000000000 0.0000000000000 0.0000000000000 4.629630e-03
0.0000000000000 0.0000000000000 0.1666666666667 4.629630e-03
0.0000000000000 0.0000000000000 0.3333333333333 4.629630e-03
0.0000000000000 0.0000000000000 0.5000000000000 4.629630e-03
0.0000000000000 0.0000000000000 0.6666666666667 4.629630e-03
0.0000000000000 0.0000000000000 0.8333333333333 4.629630e-03
0.0000000000000 0.1666666666667 0.0000000000000 4.629630e-03
0.0000000000000 0.1666666666667 0.1666666666667 4.629630e-03
0.0000000000000 0.1666666666667 0.3333333333333 4.629630e-03
0.0000000000000 0.1666666666667 0.5000000000000 4.629630e-03
0.0000000000000 0.1666666666667 0.6666666666667 4.629630e-03
0.0000000000000 0.1666666666667 0.8333333333333 4.629630e-03
0.0000000000000 0.3333333333333 0.0000000000000 4.629630e-03
0.0000000000000 0.3333333333333 0.1666666666667 4.629630e-03
0.0000000000000 0.3333333333333 0.3333333333333 4.629630e-03
0.0000000000000 0.3333333333333 0.5000000000000 4.629630e-03
0.0000000000000 0.3333333333333 0.6666666666667 4.629630e-03
0.0000000000000 0.3333333333333 0.8333333333333 4.629630e-03
0.0000000000000 0.5000000000000 0.0000000000000 4.629630e-03
0.0000000000000 0.5000000000000 0.1666666666667 4.629630e-03
0.0000000000000 0.5000000000000 0.3333333333333 4.629630e-03
0.0000000000000 0.5000000000000 0.5000000000000 4.629630e-03
0.0000000000000 0.5000000000000 0.6666666666667 4.629630e-03
0.0000000000000 0.5000000000000 0.8333333333333 4.629630e-03
0.0000000000000 0.6666666666667 0.0000000000000 4.629630e-03
0.0000000000000 0.6666666666667 0.1666666666667 4.629630e-03
0.0000000000000 0.6666666666667 0.3333333333333 4.629630e-03
0.0000000000000 0.6666666666667 0.5000000000000 4.629630e-03
0.0000000000000 0.6666666666667 0.6666666666667 4.629630e-03
0.0000000000000 0.6666666666667 0.8333333333333 4.629630e-03
0.0000000000000 0.8333333333333 0.0000000000000 4.629630e-03
0.0000000000000 0.8333333333333 0.1666666666667 4.629630e-03
0.0000000000000 0.8333333333333 0.3333333333333 4.629630e-03
0.0000000000000 0.8333333333333 0.5000000000000 4.629630e-03
0.0000000000000 0.8333333333333 0.6666666666667 4.629630e-03
0.0000000000000 0.8333333333333 0.8333333333333 4.629630e-03
0.1666666666667 0.0000000000000 0.0000000000000 4.629630e-03
0.1666666666667 0.0000000000000 0.1666666666667 4.629630e-03
0.1666666666667 0.0000000000000 0.3333333333333 4.629630e-03
0.1666666666667 0.0000000000000 0.5000000000000 4.629630e-03
0.1666666666667 0.0000000000000 0.6666666666667 4.629630e-03
0.1666666666667 0.0000000000000 0.8333333333333 4.629630e-03
0.1666666666667 0.1666666666667 0.0000000000000 4.629630e-03
0.1666666666667 0.1666666666667 0.1666666666667 4.629630e-03
0.1666666666667 0.1666666666667 0.3333333333333 4.629630e-03
0.1666666666667 0.1666666666667 0.5000000000000 4.629630e-03
0.1666666666667 0.1666666666667 0.6666666666667 4.629630e-03
0.1666666666667 0.1666666666667 0.8333333333333 4.629630e-03
0.1666666666667 0.3333333333333 0.0000000000000 4.629630e-03
0.1666666666667 0.3333333333333 0.1666666666667 4.629630e-03
0.1666666666667 0.3333333333333 0.3333333333333 4.629630e-03
0.1666666666667 0.3333333333333 0.5000000000000 4.629630e-03
0.1666666666667 0.3333333333333 0.6666666666667 4.629630e-03
0.1666666666667 0.3333333333333 0.8333333333333 4.629630e-03
0.1666666666667 0.5000000000000 0.0000000000000 4.629630e-03
0.1666666666667 0.5000000000000 0.1666666666667 4.629630e-03
0.1666666666667 0.5000000000000 0.3333333333333 4.629630e-03
0.1666666666667 0.5000000000000 0.5000000000000 4.629630e-03
0.1666666666667 0.5000000000000 0.6666666666667 4.629630e-03
0.1666666666667 0.5000000000000 0.8333333333333 4.629630e-03
0.1666666666667 0.6666666666667 0.0000000000000 4.629630e-03
0.1666666666667 0.6666666666667 0.1666666666667 4.629630e-03
0.1666666666667 0.6666666666667 0.3333333333333 4.629630e-03
0.1666666666667 0.6666666666667 0.5000000000000 4.629630e-03
0.1666666666667 0.6666666666667 0.6666666666667 4.629630e-03
0.1666666666667 0.6666666666667 0.8333333333333 4.629630e-03
0.1666666666667 0.8333333333333 0.0000000000000 4.629630e-03
0.1666666666667 0.8333333333333 0.1666666666667 4.629630e-03
0.1666666666667 0.8333333333333 0.3333333333333 4.629630e-03
0.1666666666667 0.8333333333333 0.5000000000000 4.629630e-03
0.1666666666667 0.8333333333333 0.6666666666667 4.629630e-03
0.1666666666667 0.8333333333333 0.8333333333333 4.629630e-03
0.3333333333333 0.0000000000000 0.0000000000000 4.629630e-03
0.3333333333333 0.0000000000000 0.1666666666667 4.629630e-03
0.3333333333333 0.0000000000000 0.3333333333333 4.629630e-03
0.3333333333333 0.0000000000000 0.5000000000000 4.629630e-03
0.3333333333333 0.0000000000000 0.6666666666667 4.629630e-03
0.3333333333333 0.0000000000000 0.8333333333333 4.629630e-03
0.3333333333333 0.1666666666667 0.0000000000000 4.629630e-03
0.3333333333333 0.1666666666667 0.1666666666667 4.629630e-03
0.3333333333333 0.1666666666667 0.3333333333333 4.629630e-03
0.3333333333333 0.1666666666667 0.5000000000000 4.629630e-03
0.3333333333333 0.1666666666667 0.6666666666667 4.629630e-03
0.3333333333333 0.1666666666667 0.8333333333333 4.629630e-03
0.3333333333333 0.3333333333333 0.0000000000000 4.629630e-03
0.3333333333333 0.3333333333333 0.1666666666667 4.629630e-03
0.3333333333333 0.3333333333333 0.3333333333333 4.629630e-03
0.3333333333333 0.3333333333333 0.5000000000000 4.629630e-03
0.3333333333333 0.3333333333333 0.6666666666667 4.629630e-03
0.3333333333333 0.3333333333333 0.8333333333333 4.629630e-03
0.3333333333333 0.5000000000000 0.0000000000000 4.629630e-03
0.3333333333333 0.5000000000000 0.1666666666667 4.629630e-03
0.3333333333333 0.5000000000000 0.3333333333333 4.629630e-03
0.3333333333333 0.5000000000000 0.5000000000000 4.629630e-03
0.3333333333333 0.5000000000000 0.6666666666667 4.629630e-03
0.3333333333333 0.5000000000000 0.8333333333333 4.629630e-03
0.3333333333333 0.6666666666667 0.0000000000000 4.629630e-03
0.3333333333333 0.6666666666667 0.1666666666667 4.629630e-03
0.3333333333333 0.6666666666667 0.3333333333333 4.629630e-03
0.3333333333333 0.6666666666667 0.5000000000000 4.629630e-03
0.3333333333333 0.6666666666667 0.6666666666667 4.629630e-03
0.3333333333333 0.6666666666667 0.8333333333333 4.629630e-03
0.3333333333333 0.8333333333333 0.0000000000000 4.629630e-03
0.3333333333333 0.8333333333333 0.1666666666667 4.629630e-03
0.3333333333333 0.8333333333333 0.3333333333333 4.629630e-03
0.3333333333333 0.8333333333333 0.5000000000000 4.629630e-03
0.3333333333333 0.8333333333333 0.6666666666667 4.629630e-03
0.3333333333333 0.8333333333333 0.8333333333333 4.629630e-03
0.5000000000000 0.0000000000000 0.0000000000000 4.629630e-03
0.5000000000000 0.0000000000000 0.1666666666667 4.629630e-03
0.5000000000000 0.0000000000000 0.3333333333333 4.629630e-03
0.5000000000000 0.0000000000000 0.5000000000000 4.629630e-03
0.5000000000000 0.0000000000000 0.6666666666667 4.629630e-03
0.5000000000000 0.0000000000000 0.8333333333333 4.629630e-03
0.5000000000000 0.1666666666667 0.0000000000000 4.629630e-03
0.5000000000000 0.1666666666667 0.1666666666667 4.629630e-03
0.5000000000000 0.1666666666667 0.3333333333333 4.629630e-03
0.5000000000000 0.1666666666667 0.5000000000000 4.629630e-03
0.5000000000000 0.1666666666667 0.6666666666667 4.629630e-03
0.5000000000000 0.1666666666667 0.8333333333333 4.629630e-03
0.5000000000000 0.3333333333333 0.0000000000000 4.629630e-03
0.5000000000000 0.3333333333333 0.1666666666667 4.629630e-03
0.5000000000000 0.3333333333333 0.3333333333333 4.629630e-03
0.5000000000000 0.3333333333333 0.5000000000000 4.629630e-03
0.5000000000000 0.3333333333333 0.6666666666667 4.629630e-03
0.5000000000000 0.3333333333333 0.8333333333333 4.629630e-03
0.5000000000000 0.5000000000000 0.0000000000000 4.629630e-03
0.5000000000000 0.5000000000000 0.1666666666667 4.629630e-03
0.5000000000000 0.5000000000000 0.3333333333333 4.629630e-03
0.5000000000000 0.5000000000000 0.5000000000000 4.629630e-03
0.5000000000000 0.5000000000000 0.6666666666667 4.629630e-03
0.5000000000000 0.5000000000000 0.8333333333333 4.629630e-03
0.5000000000000 0.6666666666667 0.0000000000000 4.629630e-03
0.5000000000000 0.6666666666667 0.1666666666667 4.629630e-03
0.5000000000000 0.6666666666667 0.3333333333333 4.629630e-03
0.5000000000000 0.6666666666667 0.5000000000000 4.629630e-03
0.5000000000000 0.6666666666667 0.6666666666667 4.629630e-03
0.5000000000000 0.6666666666667 0.8333333333333 4.629630e-03
0.5000000000000 0.8333333333333 0.0000000000000 4.629630e-03
0.5000000000000 0.8333333333333 0.1666666666667 4.629630e-03
0.5000000000000 0.8333333333333 0.3333333333333 4.629630e-03
0.5000000000000 0.8333333333333 0.5000000000000 4.629630e-03
0.5000000000000 0.8333333333333 0.6666666666667 4.629630e-03
0.5000000000000 0.8333333333333 0.8333333333333 4.629630e-03
0.6666666666667 0.0000000000000 0.0000000000000 4.629630e-03
0.6666666666667 0.0000000000000 0.1666666666667 4.629630e-03
0.6666666666667 0.0000000000000 0.3333333333333 4.629630e-03
0.6666666666667 0.0000000000000 0.5000000000000 4.629630e-03
0.6666666666667 0.0000000000000 0.6666666666667 4.629630e-03
0.6666666666667 0.0000000000000 0.8333333333333 4.629630e-03
0.6666666666667 0.1666666666667 0.0000000000000 4.629630e-03
0.6666666666667 0.1666666666667 0.1666666666667 4.629630e-03
0.6666666666667 0.1666666666667 0.3333333333333 4.629630e-03
0.6666666666667 0.1666666666667 0.5000000000000 4.629630e-03
0.6666666666667 0.1666666666667 0.6666666666667 4.629630e-03
0.6666666666667 0.1666666666667 0.8333333333333 4.629630e-03
0.6666666666667 0.3333333333333 0.0000000000000 4.629630e-03
0.6666666666667 0.3333333333333 0.1666666666667 4.629630e-03
0.6666666666667 0.3333333333333 0.3333333333333 4.629630e-03
0.6666666666667 0.3333333333333 0.5000000000000 4.629630e-03
0.6666666666667 0.3333333333333 0.6666666666667 4.629630e-03
0.6666666666667 0.3333333333333 0.8333333333333 4.629630e-03
0.6666666666667 0.5000000000000 0.0000000000000 4.629630e-03
0.6666666666667 0.5000000000000 0.1666666666667 4.629630e-03
0.6666666666667 0.5000000000000 0.3333333333333 4.629630e-03
0.6666666666667 0.5000000000000 0.5000000000000 4.629630e-03
0.6666666666667 0.5000000000000 0.6666666666667 4.629630e-03
0.6666666666667 0.5000000000000 0.8333333333333 4.629630e-03
0.6666666666667 0.6666666666667 0.0000000000000 4.629630e-03
0.6666666666667 0.6666666666667 0.1666666666667 4.629630e-03
0.6666666666667 0.6666666666667 0.3333333333333 4.629630e-03
0.6666666666667 0.6666666666667 0.5000000000000 4.629630e-03
0.6666666666667 0.6666666666667 0.6666666666667 4.629630e-03
0.6666666666667 0.6666666666667 0.8333333333333 4.629630e-03
0.6666666666667 0.8333333333333 0.0000000000000 4.629630e-03
0.6666666666667 0.8333333333333 0.1666666666667 4.629630e-03
0.6666666666667 0.8333333333333 0.3333333333333 4.629630e-03
0.6666666666667 0.8333333333333 0.5000000000000 4.629630e-03
0.6666666666667 0.8333333333333 0.6666666666667 4.629630e-03
0.6666666666667 0.8333333333333 0.8333333333333 4.629630e-03
0.8333333333333 0.0000000000000 0.0000000000000 4.629630e-03
0.8333333333333 0.0000000000000 0.1666666666667 4.629630e-03
0.8333333333333 0.0000000000000 0.3333333333333 4.629630e-03
0.8333333333333 0.0000000000000 0.5000000000000 4.629630e-03
0.8333333333333 0.0000000000000 0.6666666666667 4.629630e-03
0.8333333333333 0.0000000000000 0.8333333333333 4.629630e-03
0.8333333333333 0.1666666666667 0.0000000000000 4.629630e-03
0.8333333333333 0.1666666666667 0.1666666666667 4.629630e-03
0.8333333333333 0.1666666666667 0.3333333333333 4.629630e-03
0.8333333333333 0.1666666666667 0.5000000000000 4.629630e-03
0.8333333333333 0.1666666666667 0.6666666666667 4.629630e-03
0.8333333333333 0.1666666666667 0.8333333333333 4.629630e-03
0.8333333333333 0.3333333333333 0.0000000000000 4.629630e-03
0.8333333333333 0.3333333333333 0.1666666666667 4.629630e-03
0.8333333333333 0.3333333333333 0.3333333333333 4.629630e-03
0.8333333333333 0.3333333333333 0.5000000000000 4.629630e-03
0.8333333333333 0.3333333333333 0.6666666666667 4.629630e-03
0.8333333333333 0.3333333333333 0.8333333333333 4.629630e-03
0.8333333333333 0.5000000000000 0.0000000000000 4.629630e-03
0.8333333333333 0.5000000000000 0.1666666666667 4.629630e-03
0.8333333333333 0.5000000000000 0.3333333333333 4.629630e-03
0.8333333333333 0.5000000000000 0.5000000000000 4.629630e-03
0.8333333333333 0.5000000000000 0.6666666666667 4.629630e-03
0.8333333333333 0.5000000000000 0.8333333333333 4.629630e-03
0.8333333333333 0.6666666666667 0.0000000000000 4.629630e-03
0.8333333333333 0.6666666666667 0.1666666666667 4.629630e-03
0.8333333333333 0.6666666666667 0.3333333333333 4.629630e-03
0.8333333333333 0.6666666666667 0.5000000000000 4.629630e-03
0.8333333333333 0.6666666666667 0.6666666666667 4.629630e-03
0.8333333333333 0.6666666666667 0.8333333333333 4.629630e-03
0.8333333333333 0.8333333333333 0.0000000000000 4.629630e-03
0.8333333333333 0.8333333333333 0.1666666666667 4.629630e-03
0.8333333333333 0.8333333333333 0.3333333333333 4.629630e-03
0.8333333333333 0.8333333333333 0.5000000000000 4.629630e-03
0.8333333333333 0.8333333333333 0.6666666666667 4.629630e-03
0.8333333333333 0.8333333333333 0.8333333333333 4.629630e-03

This is my epw.in
--
&inputepw
prefix = 'graphene'
amass(1) = 12.01078
outdir = './'

iverbosity = 0

elph = .true.
epbwrite = .true.
epbread = .false.

epwwrite = .true.
epwread = .false.

nbndsub = 4

wannierize = .true.
num_iter = 300
iprint = 2
!dis_win_max = 2
!dis_froz_max= 7
dis_win_min= -11.73
dis_win_max= -4
!dis_froz_min = -9
!dis_froz_max = -4
proj(1) = 'f=0,0,0:l=-3'
!proj(1) = 'C:sp2'
elecselfen = .true.
phonselfen = .true.
a2f = .true.


fsthick = 1.36056981 ! eV
temps = 300 ! K (same as PRB 76, 165108)
degaussw = 0.1 ! eV

dvscf_dir = './save'

nkf1 = 18
nkf2 = 18
nkf3 = 18

nqf1 = 3
nqf2 = 3
nqf3 = 3

nk1 = 6
nk2 = 6
nk3 = 6

nq1 = 3
nq2 = 3
nq3 = 3
/

hlee
Posts: 415
Joined: Thu Aug 03, 2017 12:24 pm
Affiliation: The University of Texas at Austin

Re: Error on k+q does not fall on k-grid

Post by hlee »

Dear jameszt:

This is a well-known issue with EPW which usually arises when (1) we use relaxed structures, (2) there are differences in the number of significant digits of floating-point numbers, etc. and basically it is due to the floating-point representation.

I am working on addressing it in the next major release of EPW, but for the current and previous versions of EPW, we have little choice; for example, in order to work around it, you can modify the code as described in the following post:
viewtopic.php?p=4368#p4368

PS) You are using PAW dataset for Carbon and rigorously speaking, you have to use norm-conserving pseudopotentials instead of ultrasoft or PAW since there is an issue with the definition of electron-phonon vertex in these cases.

Sincerelly,

H. Lee

Post Reply