quantum-espresso/test-suite/epw_mob_ibte_sym/benchmark.out.git.inp=epw2....

331 lines
17 KiB
Plaintext

``:oss/
`.+s+. .+ys--yh+ `./ss+.
-sh//yy+` +yy +yy -+h+-oyy
-yh- .oyy/.-sh. .syo-.:sy- /yh
`.-.` `yh+ -oyyyo. `/syys: oys `.`
`/+ssys+-` `sh+ ` oys` .:osyo`
-yh- ./syyooyo` .sys+/oyo--yh/
`yy+ .-:-. `-/+/:` -sh-
/yh. oys
``..---hho---------` .---------..` `.-----.` -hd+---.
`./osmNMMMMMMMMMMMMMMMs. +NNMMMMMMMMNNmh+. yNMMMMMNm- oNMMMMMNmo++:`
+sy--/sdMMMhyyyyyyyNMMh- .oyNMMmyyyyyhNMMm+` -yMMMdyyo:` .oyyNMMNhs+syy`
-yy/ /MMM+.`-+/``mMMy- `mMMh:`````.dMMN:` `MMMy-`-dhhy```mMMy:``+hs
-yy+` /MMMo:-mMM+`-oo/. mMMh: `dMMN/` dMMm:`dMMMMy..MMMo-.+yo`
.sys`/MMMMNNMMMs- mMMmyooooymMMNo: oMMM/sMMMMMM++MMN//oh:
`sh+/MMMhyyMMMs- `-` mMMMMMMMMMNmy+-` -MMMhMMMsmMMmdMMd/yy+
`-/+++oyy-/MMM+.`/hh/.`mNm:` mMMd+/////:-.` NMMMMMd/:NMMMMMy:/yyo/:.`
+os+//:-..-oMMMo:--:::-/MMMo. .-mMMd+---` hMMMMN+. oMMMMMo. `-+osyso:`
syo `mNMMMMMNNNNNNNNMMMo.oNNMMMMMNNNN:` +MMMMs:` dMMMN/` ``:syo
/yh` :syyyyyyyyyyyyyyyy+.`+syyyyyyyyo:` .oyys:` .oyys:` +yh
-yh- ```````````````` ````````` `` `` oys
-+h/------------------------::::::::://////++++++++++++++++++++++///////::::/yd:
shdddddddddddddddddddddddddddddhhhhhhhhyyyyyssssssssssssssssyyyyyyyhhhhhhhddddh`
Lee, H., Poncé, S., Bushick, K., Hajinazar, S., Lafuente-Bartolome, J.,Leveillee, J.,
Lian, C., Lihm, J., Macheda, F., Mori, H., Paudyal, H., Sio, W., Tiwari, S.,
Zacharias, M., Zhang, X., Bonini, N., Kioupakis, E., Margine, E.R., and Giustino F.,
npj Comput Mater 9, 156 (2023)
Program EPW v.5.8 starts on 9Jan2024 at 13:40:44
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);
"P. Giannozzi et al., J. Chem. Phys. 152 154105 (2020);
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 4 processors
MPI processes distributed on 1 nodes
K-points division: npool = 4
35222 MiB available memory on the printing compute node when the environment starts
Reading input from epw2.in
Reading supplied temperature list.
------------------------------------------------------------------------
RESTART - RESTART - RESTART - RESTART
Restart is done without reading PWSCF save file.
Be aware that some consistency checks are therefore not done.
------------------------------------------------------------------------
--
bravais-lattice index = 0
lattice parameter (a_0) = 0.0000 a.u.
unit-cell volume = 0.0000 (a.u.)^3
number of atoms/cell = 0
number of atomic types = 0
kinetic-energy cut-off = 0.0000 Ry
charge density cut-off = 0.0000 Ry
Exchange-correlation= not set
( -1 -1 -1 -1 -1 -1 -1)
celldm(1)= 0.00000 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.0000 0.0000 0.0000 )
a(2) = ( 0.0000 0.0000 0.0000 )
a(3) = ( 0.0000 0.0000 0.0000 )
reciprocal axes: (cart. coord. in units 2 pi/a_0)
b(1) = ( 0.0000 0.0000 0.0000 )
b(2) = ( 0.0000 0.0000 0.0000 )
b(3) = ( 0.0000 0.0000 0.0000 )
Atoms inside the unit cell:
Cartesian axes
site n. atom mass positions (a_0 units)
No symmetry!
G cutoff = 0.0000 ( 0 G-vectors) FFT grid: ( 0, 0, 0)
number of k points= 0
cart. coord. in units 2pi/a_0
EPW : 0.00s CPU 0.00s WALL
EPW : 0.00s CPU 0.00s WALL
-------------------------------------------------------------------
Using si.ukk from disk
-------------------------------------------------------------------
Symmetries of Bravais lattice: 48
Symmetries of crystal: 48
Do not need to read .epb files; read .fmt files
Band disentanglement is used: nbndsub = 8
Use zone-centred Wigner-Seitz cells
Number of WS vectors for electrons 279
Number of WS vectors for phonons 19
Number of WS vectors for electron-phonon 19
Maximum number of cores for efficient parallelization 114
Results may improve by using use_ws == .TRUE.
Reading Hamiltonian, Dynamical matrix and EP vertex in Wann rep from file
Reading interatomic force constants
Read Z* and epsilon
IFC last -0.0026125
Norm of the difference between old and new effective charges: 0.0000000
Norm of the difference between old and new force-constants: 0.0000051
Imposed crystal ASR
Finished reading ifcs
Finished reading Wann rep data from file
===================================================================
Memory usage: VmHWM = 52Mb
VmPeak = 3740Mb
===================================================================
Using uniform q-mesh: 7 7 7
Size of q point mesh for interpolation: 343
Using uniform MP k-mesh: 7 7 7
Size of k point mesh for interpolation: 40
Max number of k points per pool: 10
Fermi energy coarse grid = 6.255484 eV
Fermi energy is calculated from the fine k-mesh: Ef = 6.602987 eV
Warning: check if difference with Fermi level fine grid makes sense
===================================================================
Applying a scissor shift of 0.70000 eV to the CB 5
ibndmin = 2 ebndmin = 2.620 eV
ibndmax = 7 ebndmax = 10.503 eV
Number of ep-matrix elements per pool : 1080 ~= 8.44 Kb (@ 8 bytes/ DP)
A selecq.fmt file was found but re-created because selecqread == .FALSE.
Number selected, total 50 50
Number selected, total 100 100
Number selected, total 150 150
Number selected, total 200 200
Number selected, total 250 250
Number selected, total 300 300
We only need to compute 343 q-points
Valence band maximum = 6.255484 eV
Conduction band minimum = 7.660797 eV
Temperature 400.000 K
Mobility VB Fermi level = 6.861898 eV
Mobility CB Fermi level = 7.030495 eV
Temperature 500.000 K
Mobility VB Fermi level = 7.013557 eV
Mobility CB Fermi level = 6.872901 eV
===================================================================
Scattering rate for IBTE
===================================================================
Restart and restart_step inputs deactivated (restart point at every q-points).
No intermediate mobility will be shown.
Fermi Surface thickness = 4.000000 eV
This is computed with respect to the fine Fermi level 6.602987 eV
Only states between 2.602987 eV and 10.602987 eV will be included
Save matrix elements larger than threshold: 0.236107215342E-22
Progression iq (fine) = 50/ 343
Progression iq (fine) = 100/ 343
Progression iq (fine) = 150/ 343
Progression iq (fine) = 200/ 343
Progression iq (fine) = 250/ 343
Progression iq (fine) = 300/ 343
400.000 7.0305 0.100001E+14
500.000 6.8729 0.100000E+14
epmatkqread automatically changed to .TRUE. as all scattering have been computed.
===================================================================
Memory usage: VmHWM = 62Mb
VmPeak = 3861Mb
===================================================================
Number of elements per core 265
Symmetry mapping finished
=============================================================================================
BTE in the self-energy relaxation time approximation (SERTA)
=============================================================================================
=============================================================================================
Temp Fermi Elec density Population SR Drift Elec mobility
[K] [eV] [cm^-3] [e per cell] [cm^2/Vs]
=============================================================================================
400.000 7.0305 0.10000E+14 0.00000E+00 0.283441E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.283441E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.283441E+02
500.000 6.8729 0.10000E+14 0.00000E+00 0.192631E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.192631E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.192631E+02
=============================================================================================
Start solving iterative Boltzmann Transport Equation
=============================================================================================
Iteration number: 1
=============================================================================================
Temp Fermi Elec density Population SR Drift Elec mobility
[K] [eV] [cm^-3] [e per cell] [cm^2/Vs]
=============================================================================================
400.000 7.0305 0.10000E+14 0.00000E+00 0.241569E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.241569E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.241569E+02
500.000 6.8729 0.10000E+14 0.00000E+00 0.164613E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.164613E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.164613E+02
0.241569E+02 Max error
Iteration number: 2
=============================================================================================
Temp Fermi Elec density Population SR Drift Elec mobility
[K] [eV] [cm^-3] [e per cell] [cm^2/Vs]
=============================================================================================
400.000 7.0305 0.10000E+14 0.00000E+00 0.247755E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.247755E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.247755E+02
500.000 6.8729 0.10000E+14 0.00000E+00 0.168690E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.168690E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.168690E+02
0.618605E+00 Max error
Iteration number: 3
=============================================================================================
Temp Fermi Elec density Population SR Drift Elec mobility
[K] [eV] [cm^-3] [e per cell] [cm^2/Vs]
=============================================================================================
400.000 7.0305 0.10000E+14 0.00000E+00 0.246841E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.246841E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.246841E+02
500.000 6.8729 0.10000E+14 0.00000E+00 0.168097E+02 0.000000E+00 0.000000E+00
0.00000E+00 0.000000E+00 0.168097E+02 0.000000E+00
0.00000E+00 0.000000E+00 0.000000E+00 0.168097E+02
0.913899E-01 Max error
Iteration number: 4
=============================================================================================
The iteration reached the maximum but did not converge.
=============================================================================================
Unfolding on the coarse grid
elphon_wrap : 0.01s CPU 0.01s WALL ( 1 calls)
INITIALIZATION:
Electron-Phonon interpolation
ephwann : 1.71s CPU 2.48s WALL ( 1 calls)
ep-interp : 1.44s CPU 2.20s WALL ( 343 calls)
DynW2B : 0.01s CPU 0.02s WALL ( 343 calls)
HamW2B : 0.14s CPU 0.15s WALL ( 3445 calls)
ephW2Bp : 0.73s CPU 1.26s WALL ( 343 calls)
ephW2B : 0.14s CPU 0.15s WALL ( 1715 calls)
print_ibte : 0.20s CPU 0.38s WALL ( 343 calls)
Total program execution
EPW : 1.72s CPU 2.49s WALL
% Copyright (C) 2016-2023 EPW-Collaboration
===============================================================================
Please consider citing the following papers.
% Paper describing the method on which EPW relies
F. Giustino and M. L. Cohen and S. G. Louie, Phys. Rev. B 76, 165108 (2007)
% Papers describing the EPW software
H. Lee et al., npj Comput. Mater. 9, 156 (2023)
S. Ponc\'e, E.R. Margine, C. Verdi and F. Giustino, Comput. Phys. Commun. 209, 116 (2016)
J. Noffsinger et al., Comput. Phys. Commun. 181, 2140 (2010)
% Since you used the [scattering/iterative_bte] input, please consider also citing
S. Ponc\'e, E. R. Margine and F. Giustino, Phys. Rev. B 97, 121201 (2018)
F. Macheda and N. Bonini, Phys. Rev. B 98, 201201 (2018)
For your convenience, this information is also reported in the
functionality-dependent EPW.bib file.
===============================================================================