Dear developers and users,
I'm trying to calculate the epw1.in of Pb according to the EPW manual, but an error was reported in the epw1.out file. And I couldn't find any mistake in the calculation, the files I used in the following attachments. I would appreciate it if someone could give me some suggestions on dealing with this problem.
Best wishes,
Xiangyue
Error in routine elphel2_shuffle (1) when calculated epw1.in of Pb
Moderator: stiwari
Error in routine elphel2_shuffle (1) when calculated epw1.in of Pb
- Attachments
-
- scf.in文件
&control
calculation='scf'
restart_mode='from_scratch',
prefix='pb',
pseudo_dir = '../',
outdir='./'
/
&system
ibrav= 2,
celldm(1) = 9.2225583816,
nat= 1,
ntyp= 1,
ecutwfc = 30.0
occupations='smearing',
smearing='marzari-vanderbilt',
degauss=0.05
nbnd=10
/
&electrons
conv_thr = 1.0d-10
mixing_beta = 0.7
/
ATOMIC_SPECIES
Pb 207.2 pb_s.UPF
ATOMIC_POSITIONS
Pb 0.00 0.00 0.00
K_POINTS {automatic}
6 6 6 0 0 0
ph.in文件
--
&inputph
prefix = 'pb',
fildyn = 'pb.dyn',
fildvscf = 'dvscf',
ldisp = .true.,
nq1 = 3,
nq2 = 3,
nq3 = 3,
tr2_ph = 1.0d-12
/
nscf.in 文件
&control
calculation='nscf'
restart_mode='from_scratch',
prefix='pb',
pseudo_dir = '../',
outdir='./'
/
&system
ibrav= 2,
celldm(1) = 9.2225583816,
nat= 1,
ntyp= 1,
ecutwfc = 30.0
occupations='smearing',
smearing='marzari-vanderbilt',
degauss=0.05
nbnd=10
/
&electrons
conv_thr = 1.0d-10
mixing_beta = 0.7
/
ATOMIC_SPECIES
Pb 207.2 pb_s.UPF
ATOMIC_POSITIONS
Pb 0.00 0.00 0.00
K_POINTS crystal
216
0.00000000 0.00000000 0.00000000 4.629630e-03
0.00000000 0.00000000 0.16666667 4.629630e-03
0.00000000 0.00000000 0.33333333 4.629630e-03
0.00000000 0.00000000 0.50000000 4.629630e-03
0.00000000 0.00000000 0.66666667 4.629630e-03
0.00000000 0.00000000 0.83333333 4.629630e-03
0.00000000 0.16666667 0.00000000 4.629630e-03
0.00000000 0.16666667 0.16666667 4.629630e-03
0.00000000 0.16666667 0.33333333 4.629630e-03
0.00000000 0.16666667 0.50000000 4.629630e-03
0.00000000 0.16666667 0.66666667 4.629630e-03
0.00000000 0.16666667 0.83333333 4.629630e-03
0.00000000 0.33333333 0.00000000 4.629630e-03
0.00000000 0.33333333 0.16666667 4.629630e-03
0.00000000 0.33333333 0.33333333 4.629630e-03
0.00000000 0.33333333 0.50000000 4.629630e-03
0.00000000 0.33333333 0.66666667 4.629630e-03
0.00000000 0.33333333 0.83333333 4.629630e-03
0.00000000 0.50000000 0.00000000 4.629630e-03
0.00000000 0.50000000 0.16666667 4.629630e-03
0.00000000 0.50000000 0.33333333 4.629630e-03
0.00000000 0.50000000 0.50000000 4.629630e-03
0.00000000 0.50000000 0.66666667 4.629630e-03
0.00000000 0.50000000 0.83333333 4.629630e-03
0.00000000 0.66666667 0.00000000 4.629630e-03
0.00000000 0.66666667 0.16666667 4.629630e-03
0.00000000 0.66666667 0.33333333 4.629630e-03
0.00000000 0.66666667 0.50000000 4.629630e-03
0.00000000 0.66666667 0.66666667 4.629630e-03
0.00000000 0.66666667 0.83333333 4.629630e-03
0.00000000 0.83333333 0.00000000 4.629630e-03
0.00000000 0.83333333 0.16666667 4.629630e-03
0.00000000 0.83333333 0.33333333 4.629630e-03
0.00000000 0.83333333 0.50000000 4.629630e-03
0.00000000 0.83333333 0.66666667 4.629630e-03
0.00000000 0.83333333 0.83333333 4.629630e-03
0.16666667 0.00000000 0.00000000 4.629630e-03
0.16666667 0.00000000 0.16666667 4.629630e-03
0.16666667 0.00000000 0.33333333 4.629630e-03
0.16666667 0.00000000 0.50000000 4.629630e-03
0.16666667 0.00000000 0.66666667 4.629630e-03
0.16666667 0.00000000 0.83333333 4.629630e-03
0.16666667 0.16666667 0.00000000 4.629630e-03
0.16666667 0.16666667 0.16666667 4.629630e-03
0.16666667 0.16666667 0.33333333 4.629630e-03
0.16666667 0.16666667 0.50000000 4.629630e-03
0.16666667 0.16666667 0.66666667 4.629630e-03
0.16666667 0.16666667 0.83333333 4.629630e-03
0.16666667 0.33333333 0.00000000 4.629630e-03
0.16666667 0.33333333 0.16666667 4.629630e-03
...
epw1.in 文件
--
&inputepw
prefix = 'pb',
amass(1) = 207.2
outdir = './'
dvscf_dir = '../phonon/save'
elph = .true.
epwwrite = .true.
epwread = .false.
wannierize = .true.
nbndsub = 4
bands_skipped = 'exclude_bands = 1-5'
num_iter = 300
dis_win_max = 21
dis_froz_max= 13.5
proj(1) = 'Pb:sp3'
wannier_plot= .true.
band_plot = .true.
filkf = 'path2.dat'
filqf = 'path2.dat'
nk1 = 6
nk2 = 6
nk3 = 6
nq1 = 3
nq2 = 3
nq3 = 3
/
epw1.out文件
Program EPW v.5.3.0 starts on 5Mar2022 at 5:14:58
This program is part of the open-source Quantum ESPRESSO suite
for quantum simulation of materials; please cite
"P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
"P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);
URL http://www.quantum-espresso.org",
in publications or presentations arising from this work. More details at
http://www.quantum-espresso.org/quote
Parallel version (MPI), running on 40 processors
MPI processes distributed on 1 nodes
R & G space division: proc/nbgrp/npool/nimage = 40
Fft bands division: nmany = 1
No temperature supplied. Setting temps(:) to 300 K.
Reading xml data from directory:
./pb.save/
IMPORTANT: XC functional enforced from input :
Exchange-correlation= PZ
( 1 1 0 0 0 0 0)
Any further DFT definition will be discarded
Please, verify this is what you really want
Parallelization info
--------------------
sticks: dense smooth PW G-vecs: dense smooth PW
Min 8 8 3 105 105 23
Max 9 9 4 109 109 27
Sum 349 349 121 4285 4285 941
Reading collected, re-writing distributed wavefunctions
--
bravais-lattice index = 2
lattice parameter (a_0) = 9.2226 a.u.
unit-cell volume = 196.1075 (a.u.)^3
number of atoms/cell = 1
number of atomic types = 1
kinetic-energy cut-off = 30.0000 Ry
charge density cut-off = 120.0000 Ry
Exchange-correlation= PZ
( 1 1 0 0 0 0 0)
celldm(1)= 9.22256 celldm(2)= 0.00000 celldm(3)= 0.00000
celldm(4)= 0.00000 celldm(5)= 0.00000 celldm(6)= 0.00000
crystal axes: (cart. coord. in units of a_0)
a(1) = ( -0.5000 0.0000 0.5000 )
a(2) = ( 0.0000 0.5000 0.5000 )
a(3) = ( -0.5000 0.5000 0.0000 )
reciprocal axes: (cart. coord. in units 2 pi/a_0)
b(1) = ( -1.0000 -1.0000 1.0000 )
b(2) = ( 1.0000 1.0000 1.0000 )
b(3) = ( -1.0000 1.0000 -1.0000 )
Atoms inside the unit cell:
Cartesian axes
site n. atom mass positions (a_0 units)
1 Pb 207.2000 tau( 1) = ( 0.00000 0.00000 0.00000 )
49 Sym.Ops. (with q -> -q+G )
Writing out Wannier function cube files
nr1s = 24, nr2s = 24, nr3s = 24
write_plot: wannier_plot_supercell = 5 5 5
Wannier Function Num: 1 Maximum Im/Re Ratio = 0.000000
Wannier Function Num: 2 Maximum Im/Re Ratio = 0.000000
Wannier Function Num: 3 Maximum Im/Re Ratio = 0.000000
Wannier Function Num: 4 Maximum Im/Re Ratio = 0.000000
cube files written
-------------------------------------------------------------------
WANNIER : 16.88s CPU 75.35s WALL ( 1 calls)
-------------------------------------------------------------------
Calculating kgmap
Progress kgmap: ########################################
kmaps : 0.00s CPU 0.00s WALL ( 1 calls)
Symmetries of Bravais lattice: 48
Symmetries of crystal: 48
===================================================================
irreducible q point # 1
===================================================================
Symmetries of small group of q: 48
in addition sym. q -> -q+G:
Number of q in the star = 1
List of q in the star:
1 0.000000000 0.000000000 0.000000000
Imposing acoustic sum rule on the dynamical matrix
q( 1 ) = ( 0.0000000 0.0000000 0.0000000 )
Error in routine elphel2_shuffle (1):
only one proc per pool in shuffle mode
stopping ... - epw1.out.png (43.93 KiB) Viewed 2682 times
Re: Error in routine elphel2_shuffle (1) when calculated epw1.in of Pb
Dear Xiangyue:
Sincerely,
H. Lee
Currently, you need to use the same number of pools as the number of cores when running EPW, that is, if you use 48 cores, you need to specify the same number of pools by adding to the run script "-nk 48" .Error in routine elphel2_shuffle (1):
only one proc per pool in shuffle mode
stopping ...
Sincerely,
H. Lee
Re: Error in routine elphel2_shuffle (1) when calculated epw1.in of Pb
I can run successfully according to the method you said. Thanks!