Problem with interpolation of band with ypp
Posted: Wed Oct 23, 2019 12:03 pm
Dear developers,
I am a new user of Yambo code. I want to obtain the band structure of Si within the GW approximation. The calculation is done successfully. But the following problem occurs when the post processing for the band interpolation is done with the help of ypp.
<---> [01] CPU structure, Files & I/O Directories
<---> [02] Y(ambo) P(ost)/(re) P(rocessor)
<---> [03] Core DB
<---> :: Electrons : 8.000000
<---> :: Temperature [ev]: 0.000000
<---> :: Lattice factors [a.u.]: 5.100000 5.100000 5.100000
<---> :: K points : 120
<---> :: Bands : 80
<---> :: Symmetries : 24
<---> :: RL vectors : 2733
<---> [04] K-point grid
<---> :: Q-points (IBZ): 125
<---> :: X K-points (IBZ): 120
<---> [05] CORE Variables Setup
<---> [05.01] Unit cells
<---> [05.02] Symmetries
<---> [WARNING]Spatial Inversion not found among the given symmetry list
<---> [05.03] RL shells
<---> [05.04] K-grid lattice
<---> [05.05] Energies [ev] & Occupations
<---> [06] Interpolation tool
<---> [06.01] Loading special Points for the FCC lattice
<---> [06.02] External/Internal QP corrections
<---> E<SAVE/ndb.QP[ E-E PPA E-E PPA: 27.21138 XG: 1 Xb: 1 80 Scb: 1 80]
<---> [dE_from_DB-Nearest K] Exact matches [o/o]: 100.0000
<---> [06.03] Interpolation@work: Circuit
<---> [06.03.01] Special K-points list [rlu]
0.0000000 0.0000000 0.0000000 G
0.9999999 0.9999999 0.0000000 X
0.0000000 0.9999999 0.0000000 L
0.4999999 0.9999999 -0.5000000 W
-0.7499999 0.7499999 0.0000000 K
-0.7499999 0.7499999 0.0000000 U
[ERROR] STOP signal received while in :[06.03.01] Special K-points list [rlu]
[ERROR]No points given in input
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0.
The input file(ypp.in) looks like as follows
electrons # [R] Electrons (and holes)
bnds # [R] Bands
INTERP_mode= "NN" # Interpolation mode (NN=nearest point, BOLTZ=boltztrap aproach)
cooIn= "rlu" # Points coordinates (in) cc/rlu/iku/alat
cooOut= "rlu" # Points coordinates (out) cc/rlu/iku/alat
% BANDS_bands
3 | 8 | # Number of bands
%
% INTERP_Grid
-1 |-1 |-1 | # Interpolation BZ Grid
%
INTERP_Shell_Fac= 20.00000 # The bigger it is a higher number of shells is used
CIRCUIT_E_DB_path= "none" # SAVE obtained from the QE `bands` run (alternative to %BANDS_kpts)
BANDS_path= "" # BANDS path points labels (G,M,K,L...)
BANDS_steps= 10 # Number of divisions
#BANDS_built_in # Print the bands of the generating points of the circuit using the nearest internal point
GfnQPdb= "E < SAVE/ndb.QP" # [EXTQP G] Database
GfnQP_N= 1 # [EXTQP G] Interpolation neighbours
% GfnQP_E
0.000000 | 1.000000 | 1.000000 | # [EXTQP G] E parameters (c/v) eV|adim|adim
%
GfnQP_Z= ( 1.000000 , 0.000000 ) # [EXTQP G] Z factor (c/v)
GfnQP_Wv_E= 0.000000 eV # [EXTQP G] W Energy reference (valence)
% GfnQP_Wv
0.00 | 0.00 | 0.00 | # [EXTQP G] W parameters (valence) eV| 1|eV^-1
%
GfnQP_Wv_dos= 0.000000 eV # [EXTQP G] W dos pre-factor (valence)
GfnQP_Wc_E= 0.000000 eV # [EXTQP G] W Energy reference (conduction)
% GfnQP_Wc
0.00 | 0.00 | 0.00 | # [EXTQP G] W parameters (conduction) eV| 1 |eV^-1
%
GfnQP_Wc_dos= 0.000000 eV # [EXTQP G] W dos pre-factor (conduction)
%BANDS_kpts # K points of the bands circuit
%
| %BANDS_kpts # K points of the bands circuit
0.00000 |0.00000 |0.00000|
0.00000 |0.50000 |0.50000|
0.50000 |0.50000 |0.50000|
| %
Regards,
Bramha
I am a new user of Yambo code. I want to obtain the band structure of Si within the GW approximation. The calculation is done successfully. But the following problem occurs when the post processing for the band interpolation is done with the help of ypp.
<---> [01] CPU structure, Files & I/O Directories
<---> [02] Y(ambo) P(ost)/(re) P(rocessor)
<---> [03] Core DB
<---> :: Electrons : 8.000000
<---> :: Temperature [ev]: 0.000000
<---> :: Lattice factors [a.u.]: 5.100000 5.100000 5.100000
<---> :: K points : 120
<---> :: Bands : 80
<---> :: Symmetries : 24
<---> :: RL vectors : 2733
<---> [04] K-point grid
<---> :: Q-points (IBZ): 125
<---> :: X K-points (IBZ): 120
<---> [05] CORE Variables Setup
<---> [05.01] Unit cells
<---> [05.02] Symmetries
<---> [WARNING]Spatial Inversion not found among the given symmetry list
<---> [05.03] RL shells
<---> [05.04] K-grid lattice
<---> [05.05] Energies [ev] & Occupations
<---> [06] Interpolation tool
<---> [06.01] Loading special Points for the FCC lattice
<---> [06.02] External/Internal QP corrections
<---> E<SAVE/ndb.QP[ E-E PPA E-E PPA: 27.21138 XG: 1 Xb: 1 80 Scb: 1 80]
<---> [dE_from_DB-Nearest K] Exact matches [o/o]: 100.0000
<---> [06.03] Interpolation@work: Circuit
<---> [06.03.01] Special K-points list [rlu]
0.0000000 0.0000000 0.0000000 G
0.9999999 0.9999999 0.0000000 X
0.0000000 0.9999999 0.0000000 L
0.4999999 0.9999999 -0.5000000 W
-0.7499999 0.7499999 0.0000000 K
-0.7499999 0.7499999 0.0000000 U
[ERROR] STOP signal received while in :[06.03.01] Special K-points list [rlu]
[ERROR]No points given in input
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0.
The input file(ypp.in) looks like as follows
electrons # [R] Electrons (and holes)
bnds # [R] Bands
INTERP_mode= "NN" # Interpolation mode (NN=nearest point, BOLTZ=boltztrap aproach)
cooIn= "rlu" # Points coordinates (in) cc/rlu/iku/alat
cooOut= "rlu" # Points coordinates (out) cc/rlu/iku/alat
% BANDS_bands
3 | 8 | # Number of bands
%
% INTERP_Grid
-1 |-1 |-1 | # Interpolation BZ Grid
%
INTERP_Shell_Fac= 20.00000 # The bigger it is a higher number of shells is used
CIRCUIT_E_DB_path= "none" # SAVE obtained from the QE `bands` run (alternative to %BANDS_kpts)
BANDS_path= "" # BANDS path points labels (G,M,K,L...)
BANDS_steps= 10 # Number of divisions
#BANDS_built_in # Print the bands of the generating points of the circuit using the nearest internal point
GfnQPdb= "E < SAVE/ndb.QP" # [EXTQP G] Database
GfnQP_N= 1 # [EXTQP G] Interpolation neighbours
% GfnQP_E
0.000000 | 1.000000 | 1.000000 | # [EXTQP G] E parameters (c/v) eV|adim|adim
%
GfnQP_Z= ( 1.000000 , 0.000000 ) # [EXTQP G] Z factor (c/v)
GfnQP_Wv_E= 0.000000 eV # [EXTQP G] W Energy reference (valence)
% GfnQP_Wv
0.00 | 0.00 | 0.00 | # [EXTQP G] W parameters (valence) eV| 1|eV^-1
%
GfnQP_Wv_dos= 0.000000 eV # [EXTQP G] W dos pre-factor (valence)
GfnQP_Wc_E= 0.000000 eV # [EXTQP G] W Energy reference (conduction)
% GfnQP_Wc
0.00 | 0.00 | 0.00 | # [EXTQP G] W parameters (conduction) eV| 1 |eV^-1
%
GfnQP_Wc_dos= 0.000000 eV # [EXTQP G] W dos pre-factor (conduction)
%BANDS_kpts # K points of the bands circuit
%
| %BANDS_kpts # K points of the bands circuit
0.00000 |0.00000 |0.00000|
0.00000 |0.50000 |0.50000|
0.50000 |0.50000 |0.50000|
| %
Regards,
Bramha