Porblem on GaN Tutorial

[yeoh]

Dear moderator,

I am new to EPW. I am trying to run the example for GaN.
I successfully run the scf and nscf calculation as per the tutorial.
However when I run the epw.x, I encounter the following error:

Error: dis_spheres_first_wann is larger than num_bands-num_wann+1
Error: examine the output/error file for details

Can you please help?

Best regards
Yeoh

[sponce]

Dear yeoh,

If you ran "make" inside EPW/src, you need to run it a second time.

Can you try to do "make" inside EPW/src twice and then try again?

Best,

Samuel.

PS: If you do it from the root in /espresso, this should be done automatically by "make epw".

[yeoh]

Dear Dr Samuel,

Thank you very much for your response. I try ran "make" inside EPW/src twice but it did not work.
However, when I ran "make" inside EPW folder twice, it works.

So I proceed with the GaN tutorial using:

epw.x <epw.in>epw.out

At the end of the epw.out file, I get:

===================================================================
irreducible q point # 28
===================================================================

Symmetries of small group of q: 6

Number of q in the star = 2
List of q in the star:
1 0.333333333 0.577350269 -0.306727120
2 -0.333333333 0.577350269 -0.306727120

q( 215 ) = ( 0.3333333 0.5773503 -0.3067271 )
BMN calculated

q( 216 ) = ( -0.3333333 0.5773503 -0.3067271 )
BMN calculated


Writing epmatq on .epb files


===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= PID 9021 RUNNING AT localhost.localdomain
= EXIT CODE: 9
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)

Any idea, why I get the above message?

Best regards
Yeoh

[sponce]

This error has been reported previously on this forum.

Its because of PGI

Please see viewtopic.php?f=6&t=83

Best,

Samuel

[yeoh]

Dear Dr Sponce,

As recommended, I change all the occurrences in EPW/src from

myfmt = *(3x,E15.5)) to myfmt = "(1000(3x,E15.5))"

and then i ran "make" twice in EPW/src.

I still received the same error:
===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= PID 23983 RUNNING AT localhost.localdomain
= EXIT CODE: 9
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)

Did I miss out anything?

Thanks

Regards
Yeoh

[sponce]

Dear Yeoh,

Did you modify all occurrence of 'myfmt' in the code ? There should be 3 or 4.

If so I will need more info (all input files and the place where the code crash).

Best,

Samuel

[yeoh]

Dear Dr Ponce,

I modified all occurrence of 'myfmt' in the code, which located in
1. ephwann_shuffle.f90
2. selfen_phon.f90
Both are inside the src folder.

The scf input are as follow:

&control
calculation='scf'
prefix='gan'
restart_mode='from_scratch'
pseudo_dir='/home/common/gan/pp'
outdir='./'
tprnfor=.true.
tstress=.true.
wf_collect=.true.
/
&system
ibrav = 4
celldm(1)=5.95484286816
celldm(3)=1.63011343669
nat = 4
ntyp = 2
ecutwfc = 160
/
&electrons
diagonalization='david'
mixing_beta=0.7
conv_thr=1.0d-10
/
ATOMIC_SPECIES
Ga 69.723 Ga_ONCV_LDA-1.0.upf
N 14.007 N_ONCV_LDA-1.0.upf
ATOMIC_POSITIONS crystal
Ga 0.666666666667 0.333333333333 0.0
N 0.666666666667 0.333333333333 0.376429222
Ga 0.333333333333 0.666666666667 0.5
N 0.333333333333 0.666666666667 0.876429222
K_POINTS automatic
6 6 6 0 0 0

The nscf input:

&control
calculation='nscf'
prefix='gan'
restart_mode='from_scratch'
pseudo_dir='/home/common/ganm/pp'
outdir='./'
tprnfor=.true.
tstress=.true.
verbosity='high'
/
&system
ibrav = 4
celldm(1)=5.95484286816
celldm(3)=1.63011343669
nat = 4
ntyp = 2
ecutwfc = 160
/
&electrons
diagonalization='david'
mixing_beta=0.7
conv_thr=1.0d-10
/
ATOMIC_SPECIES
Ga 69.723 Ga_ONCV_LDA-1.0.upf
N 14.007 N_ONCV_LDA-1.0.upf
ATOMIC_POSITIONS crystal
Ga 0.666666666667 0.333333333333 0.0
N 0.666666666667 0.333333333333 0.376429222
Ga 0.333333333333 0.666666666667 0.5
N 0.333333333333 0.666666666667 0.876429222
K_POINTS crystal
216
0.0 0.0 0.0 0.00462962962963
0.0 0.0 0.166666666667 0.00462962962963
0.0 0.0 0.333333333333 0.00462962962963
0.0 0.0 0.5 0.00462962962963
0.0 0.0 0.666666666667 0.00462962962963
0.0 0.0 0.833333333333 0.00462962962963
0.0 0.166666666667 0.0 0.00462962962963
0.0 0.166666666667 0.166666666667 0.00462962962963
0.0 0.166666666667 0.333333333333 0.00462962962963
0.0 0.166666666667 0.5 0.00462962962963
0.0 0.166666666667 0.666666666667 0.00462962962963
0.0 0.166666666667 0.833333333333 0.00462962962963
0.0 0.333333333333 0.0 0.00462962962963
0.0 0.333333333333 0.166666666667 0.00462962962963
0.0 0.333333333333 0.333333333333 0.00462962962963
0.0 0.333333333333 0.5 0.00462962962963
0.0 0.333333333333 0.666666666667 0.00462962962963
0.0 0.333333333333 0.833333333333 0.00462962962963
0.0 0.5 0.0 0.00462962962963
0.0 0.5 0.166666666667 0.00462962962963
0.0 0.5 0.333333333333 0.00462962962963
0.0 0.5 0.5 0.00462962962963
0.0 0.5 0.666666666667 0.00462962962963
0.0 0.5 0.833333333333 0.00462962962963
0.0 0.666666666667 0.0 0.00462962962963
0.0 0.666666666667 0.166666666667 0.00462962962963
0.0 0.666666666667 0.333333333333 0.00462962962963
0.0 0.666666666667 0.5 0.00462962962963
0.0 0.666666666667 0.666666666667 0.00462962962963
0.0 0.666666666667 0.833333333333 0.00462962962963
0.0 0.833333333333 0.0 0.00462962962963
0.0 0.833333333333 0.166666666667 0.00462962962963
0.0 0.833333333333 0.333333333333 0.00462962962963
0.0 0.833333333333 0.5 0.00462962962963
0.0 0.833333333333 0.666666666667 0.00462962962963
0.0 0.833333333333 0.833333333333 0.00462962962963
0.166666666667 0.0 0.0 0.00462962962963
0.166666666667 0.0 0.166666666667 0.00462962962963
0.166666666667 0.0 0.333333333333 0.00462962962963
0.166666666667 0.0 0.5 0.00462962962963
0.166666666667 0.0 0.666666666667 0.00462962962963
0.166666666667 0.0 0.833333333333 0.00462962962963
0.166666666667 0.166666666667 0.0 0.00462962962963
0.166666666667 0.166666666667 0.166666666667 0.00462962962963
0.166666666667 0.166666666667 0.333333333333 0.00462962962963
0.166666666667 0.166666666667 0.5 0.00462962962963
0.166666666667 0.166666666667 0.666666666667 0.00462962962963
0.166666666667 0.166666666667 0.833333333333 0.00462962962963
0.166666666667 0.333333333333 0.0 0.00462962962963
0.166666666667 0.333333333333 0.166666666667 0.00462962962963
0.166666666667 0.333333333333 0.333333333333 0.00462962962963
0.166666666667 0.333333333333 0.5 0.00462962962963
0.166666666667 0.333333333333 0.666666666667 0.00462962962963
0.166666666667 0.333333333333 0.833333333333 0.00462962962963
0.166666666667 0.5 0.0 0.00462962962963
0.166666666667 0.5 0.166666666667 0.00462962962963
0.166666666667 0.5 0.333333333333 0.00462962962963
0.166666666667 0.5 0.5 0.00462962962963
0.166666666667 0.5 0.666666666667 0.00462962962963
0.166666666667 0.5 0.833333333333 0.00462962962963
0.166666666667 0.666666666667 0.0 0.00462962962963
0.166666666667 0.666666666667 0.166666666667 0.00462962962963
0.166666666667 0.666666666667 0.333333333333 0.00462962962963
0.166666666667 0.666666666667 0.5 0.00462962962963
0.166666666667 0.666666666667 0.666666666667 0.00462962962963
0.166666666667 0.666666666667 0.833333333333 0.00462962962963
0.166666666667 0.833333333333 0.0 0.00462962962963
0.166666666667 0.833333333333 0.166666666667 0.00462962962963
0.166666666667 0.833333333333 0.333333333333 0.00462962962963
0.166666666667 0.833333333333 0.5 0.00462962962963
0.166666666667 0.833333333333 0.666666666667 0.00462962962963
0.166666666667 0.833333333333 0.833333333333 0.00462962962963
0.333333333333 0.0 0.0 0.00462962962963
0.333333333333 0.0 0.166666666667 0.00462962962963
0.333333333333 0.0 0.333333333333 0.00462962962963
0.333333333333 0.0 0.5 0.00462962962963
0.333333333333 0.0 0.666666666667 0.00462962962963
0.333333333333 0.0 0.833333333333 0.00462962962963
0.333333333333 0.166666666667 0.0 0.00462962962963
0.333333333333 0.166666666667 0.166666666667 0.00462962962963
0.333333333333 0.166666666667 0.333333333333 0.00462962962963
0.333333333333 0.166666666667 0.5 0.00462962962963
0.333333333333 0.166666666667 0.666666666667 0.00462962962963
0.333333333333 0.166666666667 0.833333333333 0.00462962962963
0.333333333333 0.333333333333 0.0 0.00462962962963
0.333333333333 0.333333333333 0.166666666667 0.00462962962963
0.333333333333 0.333333333333 0.333333333333 0.00462962962963
0.333333333333 0.333333333333 0.5 0.00462962962963
0.333333333333 0.333333333333 0.666666666667 0.00462962962963
0.333333333333 0.333333333333 0.833333333333 0.00462962962963
0.333333333333 0.5 0.0 0.00462962962963
0.333333333333 0.5 0.166666666667 0.00462962962963
0.333333333333 0.5 0.333333333333 0.00462962962963
0.333333333333 0.5 0.5 0.00462962962963
0.333333333333 0.5 0.666666666667 0.00462962962963
0.333333333333 0.5 0.833333333333 0.00462962962963
0.333333333333 0.666666666667 0.0 0.00462962962963
0.333333333333 0.666666666667 0.166666666667 0.00462962962963
0.333333333333 0.666666666667 0.333333333333 0.00462962962963
0.333333333333 0.666666666667 0.5 0.00462962962963
0.333333333333 0.666666666667 0.666666666667 0.00462962962963
0.333333333333 0.666666666667 0.833333333333 0.00462962962963
0.333333333333 0.833333333333 0.0 0.00462962962963
0.333333333333 0.833333333333 0.166666666667 0.00462962962963
0.333333333333 0.833333333333 0.333333333333 0.00462962962963
0.333333333333 0.833333333333 0.5 0.00462962962963
0.333333333333 0.833333333333 0.666666666667 0.00462962962963
0.333333333333 0.833333333333 0.833333333333 0.00462962962963
0.5 0.0 0.0 0.00462962962963
0.5 0.0 0.166666666667 0.00462962962963
0.5 0.0 0.333333333333 0.00462962962963
0.5 0.0 0.5 0.00462962962963
0.5 0.0 0.666666666667 0.00462962962963
0.5 0.0 0.833333333333 0.00462962962963
0.5 0.166666666667 0.0 0.00462962962963
0.5 0.166666666667 0.166666666667 0.00462962962963
0.5 0.166666666667 0.333333333333 0.00462962962963
0.5 0.166666666667 0.5 0.00462962962963
0.5 0.166666666667 0.666666666667 0.00462962962963
0.5 0.166666666667 0.833333333333 0.00462962962963
0.5 0.333333333333 0.0 0.00462962962963
0.5 0.333333333333 0.166666666667 0.00462962962963
0.5 0.333333333333 0.333333333333 0.00462962962963
0.5 0.333333333333 0.5 0.00462962962963
0.5 0.333333333333 0.666666666667 0.00462962962963
0.5 0.333333333333 0.833333333333 0.00462962962963
0.5 0.5 0.0 0.00462962962963
0.5 0.5 0.166666666667 0.00462962962963
0.5 0.5 0.333333333333 0.00462962962963
0.5 0.5 0.5 0.00462962962963
0.5 0.5 0.666666666667 0.00462962962963
0.5 0.5 0.833333333333 0.00462962962963
0.5 0.666666666667 0.0 0.00462962962963
0.5 0.666666666667 0.166666666667 0.00462962962963
0.5 0.666666666667 0.333333333333 0.00462962962963
0.5 0.666666666667 0.5 0.00462962962963
0.5 0.666666666667 0.666666666667 0.00462962962963
0.5 0.666666666667 0.833333333333 0.00462962962963
0.5 0.833333333333 0.0 0.00462962962963
0.5 0.833333333333 0.166666666667 0.00462962962963
0.5 0.833333333333 0.333333333333 0.00462962962963
0.5 0.833333333333 0.5 0.00462962962963
0.5 0.833333333333 0.666666666667 0.00462962962963
0.5 0.833333333333 0.833333333333 0.00462962962963
0.666666666667 0.0 0.0 0.00462962962963
0.666666666667 0.0 0.166666666667 0.00462962962963
0.666666666667 0.0 0.333333333333 0.00462962962963
0.666666666667 0.0 0.5 0.00462962962963
0.666666666667 0.0 0.666666666667 0.00462962962963
0.666666666667 0.0 0.833333333333 0.00462962962963
0.666666666667 0.166666666667 0.0 0.00462962962963
0.666666666667 0.166666666667 0.166666666667 0.00462962962963
0.666666666667 0.166666666667 0.333333333333 0.00462962962963
0.666666666667 0.166666666667 0.5 0.00462962962963
0.666666666667 0.166666666667 0.666666666667 0.00462962962963
0.666666666667 0.166666666667 0.833333333333 0.00462962962963
0.666666666667 0.333333333333 0.0 0.00462962962963
0.666666666667 0.333333333333 0.166666666667 0.00462962962963
0.666666666667 0.333333333333 0.333333333333 0.00462962962963
0.666666666667 0.333333333333 0.5 0.00462962962963
0.666666666667 0.333333333333 0.666666666667 0.00462962962963
0.666666666667 0.333333333333 0.833333333333 0.00462962962963
0.666666666667 0.5 0.0 0.00462962962963
0.666666666667 0.5 0.166666666667 0.00462962962963
0.666666666667 0.5 0.333333333333 0.00462962962963
0.666666666667 0.5 0.5 0.00462962962963
0.666666666667 0.5 0.666666666667 0.00462962962963
0.666666666667 0.5 0.833333333333 0.00462962962963
0.666666666667 0.666666666667 0.0 0.00462962962963
0.666666666667 0.666666666667 0.166666666667 0.00462962962963
0.666666666667 0.666666666667 0.333333333333 0.00462962962963
0.666666666667 0.666666666667 0.5 0.00462962962963
0.666666666667 0.666666666667 0.666666666667 0.00462962962963
0.666666666667 0.666666666667 0.833333333333 0.00462962962963
0.666666666667 0.833333333333 0.0 0.00462962962963
0.666666666667 0.833333333333 0.166666666667 0.00462962962963
0.666666666667 0.833333333333 0.333333333333 0.00462962962963
0.666666666667 0.833333333333 0.5 0.00462962962963
0.666666666667 0.833333333333 0.666666666667 0.00462962962963
0.666666666667 0.833333333333 0.833333333333 0.00462962962963
0.833333333333 0.0 0.0 0.00462962962963
0.833333333333 0.0 0.166666666667 0.00462962962963
0.833333333333 0.0 0.333333333333 0.00462962962963
0.833333333333 0.0 0.5 0.00462962962963
0.833333333333 0.0 0.666666666667 0.00462962962963
0.833333333333 0.0 0.833333333333 0.00462962962963
0.833333333333 0.166666666667 0.0 0.00462962962963
0.833333333333 0.166666666667 0.166666666667 0.00462962962963
0.833333333333 0.166666666667 0.333333333333 0.00462962962963
0.833333333333 0.166666666667 0.5 0.00462962962963
0.833333333333 0.166666666667 0.666666666667 0.00462962962963
0.833333333333 0.166666666667 0.833333333333 0.00462962962963
0.833333333333 0.333333333333 0.0 0.00462962962963
0.833333333333 0.333333333333 0.166666666667 0.00462962962963
0.833333333333 0.333333333333 0.333333333333 0.00462962962963
0.833333333333 0.333333333333 0.5 0.00462962962963
0.833333333333 0.333333333333 0.666666666667 0.00462962962963
0.833333333333 0.333333333333 0.833333333333 0.00462962962963
0.833333333333 0.5 0.0 0.00462962962963
0.833333333333 0.5 0.166666666667 0.00462962962963
0.833333333333 0.5 0.333333333333 0.00462962962963
0.833333333333 0.5 0.5 0.00462962962963
0.833333333333 0.5 0.666666666667 0.00462962962963
0.833333333333 0.5 0.833333333333 0.00462962962963
0.833333333333 0.666666666667 0.0 0.00462962962963
0.833333333333 0.666666666667 0.166666666667 0.00462962962963
0.833333333333 0.666666666667 0.333333333333 0.00462962962963
0.833333333333 0.666666666667 0.5 0.00462962962963
0.833333333333 0.666666666667 0.666666666667 0.00462962962963
0.833333333333 0.666666666667 0.833333333333 0.00462962962963
0.833333333333 0.833333333333 0.0 0.00462962962963
0.833333333333 0.833333333333 0.166666666667 0.00462962962963
0.833333333333 0.833333333333 0.333333333333 0.00462962962963
0.833333333333 0.833333333333 0.5 0.00462962962963
0.833333333333 0.833333333333 0.666666666667 0.00462962962963
0.833333333333 0.833333333333 0.833333333333 0.00462962962963

The epw input are as follow:
--
&inputepw
prefix = 'gan'
amass(1) = 69.723
amass(2) = 14.007
outdir = './'

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

lpolar = .true.

epwwrite = .true.
epwread = .false.

nbndsub = 18
nbndskip = 0

wannierize = .true.
num_iter = 500
!dis_win_max = 11
!dis_froz_max= 11
proj(1) = 'N:sp3'
proj(2) = 'Ga:d'

elinterp = .true.
phinterp = .true.

tshuffle2 = .true.
tphases = .false.

elecselfen = .false.
phonselfen = .true.
a2f = .false.

parallel_k = .true.
parallel_q = .false.

fsthick = 3.0
eptemp = 300.0
degaussw = 0.1
efermi_read = .true.
fermi_energy= 9.9775

dvscf_dir = '/home/common/gan/save'
filqf = 'gan_band.qpt'

nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 6
nq2 = 6
nq3 = 6

nkf1 = 10
nkf2 = 10
nkf3 = 10
/
28 cartesian
0.000000000 0.000000000 0.000000000 1
0.000000000 0.000000000 0.102242373 1
0.000000000 0.000000000 0.204484747 1
0.000000000 0.000000000 -0.306727120 1
0.000000000 0.192450090 0.000000000 1
0.000000000 0.192450090 0.102242373 1
0.000000000 0.192450090 0.204484747 1
0.000000000 0.192450090 -0.306727120 1
0.000000000 0.384900179 0.000000000 1
0.000000000 0.384900179 0.102242373 1
0.000000000 0.384900179 0.204484747 1
0.000000000 0.384900179 -0.306727120 1
0.000000000 -0.577350269 0.000000000 1
0.000000000 -0.577350269 0.102242373 1
0.000000000 -0.577350269 0.204484747 1
0.000000000 -0.577350269 -0.306727120 1
0.166666667 0.288675135 0.000000000 1
0.166666667 0.288675135 0.102242373 1
0.166666667 0.288675135 0.204484747 1
0.166666667 0.288675135 -0.306727120 1
0.166666667 0.481125224 0.000000000 1
0.166666667 0.481125224 0.102242373 1
0.166666667 0.481125224 0.204484747 1
0.166666667 0.481125224 -0.306727120 1
0.333333333 0.577350269 0.000000000 1
0.333333333 0.577350269 0.102242373 1
0.333333333 0.577350269 0.204484747 1
0.333333333 0.577350269 -0.306727120 1

The error occur is captured in the epw output file:

===================================================================
irreducible q point # 28
===================================================================

Symmetries of small group of q: 6

Number of q in the star = 2
List of q in the star:
1 0.333333333 0.577350269 -0.306727120
2 -0.333333333 0.577350269 -0.306727120

q( 215 ) = ( 0.3333333 0.5773503 -0.3067271 )
BMN calculated

q( 216 ) = ( -0.3333333 0.5773503 -0.3067271 )
BMN calculated


Writing epmatq on .epb files


===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= PID 23983 RUNNING AT localhost.localdomain
= EXIT CODE: 9
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)

P/S: The scf, nscf and epw.in file is taken same as in the tutorial.

Thanks!

Best regards
Yeoh

[sponce]

Dear Yeoh,

Nothing obvious come to my mind as the source of error.

However, I did a lot of debugging since the version 5.4 and I'm quite convince that it should now work.

There will be the release of QE 6.0 (that will include EPW 4.1) in the next couple of days (end of Sept likely) that should fix it.

If it still do not work then, I will investigate more (currently this example works perfectly for me on intel compiler).

Best,

Samuel

[yeoh]

Dear Dr Sponce,

Thank you very much!

Best regards
Yeoh