mirror of https://gitlab.com/QEF/q-e.git
979 lines
37 KiB
Plaintext
979 lines
37 KiB
Plaintext
|
|
Program PWSCF v.5.4.0 (svn rev. 12493) starts on 21Jun2016 at 16:49: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);
|
|
URL http://www.quantum-espresso.org",
|
|
in publications or presentations arising from this work. More details at
|
|
http://www.quantum-espresso.org/quote
|
|
|
|
HOST : @host@
|
|
ARCH : x86_64
|
|
CC : cc
|
|
CPP : cpp
|
|
F90 : mpif90
|
|
F77 : gfortran
|
|
DFLAGS : -D__GFORTRAN -D__STD_F95 -D__FFTW -D__MPI -D__PARA
|
|
BLAS LIBS : /home/sponce/program/espresso/BLAS/blas.a
|
|
LAPACK LIBS : /home/sponce/program/espresso/lapack-3.2/lapack.a
|
|
FFT LIBS :
|
|
MASS LIBS :
|
|
|
|
Parallel version (MPI), running on 1 processors
|
|
Waiting for input...
|
|
Reading input from standard input
|
|
|
|
Current dimensions of program PWSCF are:
|
|
Max number of different atomic species (ntypx) = 10
|
|
Max number of k-points (npk) = 40000
|
|
Max angular momentum in pseudopotentials (lmaxx) = 3
|
|
|
|
Subspace diagonalization in iterative solution of the eigenvalue problem:
|
|
a serial algorithm will be used
|
|
|
|
|
|
G-vector sticks info
|
|
--------------------
|
|
sticks: dense smooth PW G-vecs: dense smooth PW
|
|
Sum 703 703 211 12387 12387 1989
|
|
|
|
|
|
|
|
bravais-lattice index = 2
|
|
lattice parameter (alat) = 9.2700 a.u.
|
|
unit-cell volume = 199.1495 (a.u.)^3
|
|
number of atoms/cell = 1
|
|
number of atomic types = 1
|
|
number of electrons = 14.00
|
|
number of Kohn-Sham states= 22
|
|
kinetic-energy cutoff = 60.0000 Ry
|
|
charge density cutoff = 240.0000 Ry
|
|
convergence threshold = 1.0E-10
|
|
mixing beta = 0.7000
|
|
number of iterations used = 8 plain mixing
|
|
Exchange-correlation = PZ ( 1 1 0 0 0 0)
|
|
Non magnetic calculation with spin-orbit
|
|
|
|
|
|
celldm(1)= 9.270000 celldm(2)= 0.000000 celldm(3)= 0.000000
|
|
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
|
|
|
|
crystal axes: (cart. coord. in units of alat)
|
|
a(1) = ( -0.500000 0.000000 0.500000 )
|
|
a(2) = ( 0.000000 0.500000 0.500000 )
|
|
a(3) = ( -0.500000 0.500000 0.000000 )
|
|
|
|
reciprocal axes: (cart. coord. in units 2 pi/alat)
|
|
b(1) = ( -1.000000 -1.000000 1.000000 )
|
|
b(2) = ( 1.000000 1.000000 1.000000 )
|
|
b(3) = ( -1.000000 1.000000 -1.000000 )
|
|
|
|
|
|
PseudoPot. # 1 for Pb read from file:
|
|
./pb_s.UPF
|
|
MD5 check sum: 7774f596943a2055f0d8d19eea6f2ae8
|
|
Pseudo is Norm-conserving, Zval = 14.0
|
|
Generated using "atomic" code by A. Dal Corso (Quantum ESPRESSO distribution)
|
|
Using radial grid of 1281 points, 4 beta functions with:
|
|
l(1) = 2
|
|
l(2) = 2
|
|
l(3) = 1
|
|
l(4) = 1
|
|
|
|
atomic species valence mass pseudopotential
|
|
Pb 14.00 207.20000 Pb( 1.00)
|
|
|
|
48 Sym. Ops., with inversion, found
|
|
|
|
|
|
s frac. trans.
|
|
|
|
isym = 1 identity
|
|
|
|
cryst. s( 1) = ( 1 0 0 )
|
|
( 0 1 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 2 180 deg rotation - cart. axis [0,0,1]
|
|
|
|
cryst. s( 2) = ( 0 1 -1 )
|
|
( 1 0 -1 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 3 180 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s( 3) = ( -1 0 0 )
|
|
( -1 0 1 )
|
|
( -1 1 0 )
|
|
|
|
cart. s( 3) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 4 180 deg rotation - cart. axis [1,0,0]
|
|
|
|
cryst. s( 4) = ( 0 -1 1 )
|
|
( 0 -1 0 )
|
|
( 1 -1 0 )
|
|
|
|
cart. s( 4) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 5 180 deg rotation - cart. axis [1,1,0]
|
|
|
|
cryst. s( 5) = ( 0 -1 0 )
|
|
( -1 0 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s( 5) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 6 180 deg rotation - cart. axis [1,-1,0]
|
|
|
|
cryst. s( 6) = ( -1 0 1 )
|
|
( 0 -1 1 )
|
|
( 0 0 1 )
|
|
|
|
cart. s( 6) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 7 90 deg rotation - cart. axis [0,0,-1]
|
|
|
|
cryst. s( 7) = ( 0 1 0 )
|
|
( 0 1 -1 )
|
|
( -1 1 0 )
|
|
|
|
cart. s( 7) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 8 90 deg rotation - cart. axis [0,0,1]
|
|
|
|
cryst. s( 8) = ( 1 0 -1 )
|
|
( 1 0 0 )
|
|
( 1 -1 0 )
|
|
|
|
cart. s( 8) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 9 180 deg rotation - cart. axis [1,0,1]
|
|
|
|
cryst. s( 9) = ( -1 0 0 )
|
|
( 0 0 -1 )
|
|
( 0 -1 0 )
|
|
|
|
cart. s( 9) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 10 180 deg rotation - cart. axis [-1,0,1]
|
|
|
|
cryst. s(10) = ( 1 0 0 )
|
|
( 1 -1 0 )
|
|
( 1 0 -1 )
|
|
|
|
cart. s(10) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 11 90 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s(11) = ( 0 1 -1 )
|
|
( -1 1 0 )
|
|
( 0 1 0 )
|
|
|
|
cart. s(11) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 12 90 deg rotation - cart. axis [0,-1,0]
|
|
|
|
cryst. s(12) = ( 0 -1 1 )
|
|
( 0 0 1 )
|
|
( -1 0 1 )
|
|
|
|
cart. s(12) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 13 180 deg rotation - cart. axis [0,1,1]
|
|
|
|
cryst. s(13) = ( -1 1 0 )
|
|
( 0 1 0 )
|
|
( 0 1 -1 )
|
|
|
|
cart. s(13) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 14 180 deg rotation - cart. axis [0,1,-1]
|
|
|
|
cryst. s(14) = ( 0 0 -1 )
|
|
( 0 -1 0 )
|
|
( -1 0 0 )
|
|
|
|
cart. s(14) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 15 90 deg rotation - cart. axis [-1,0,0]
|
|
|
|
cryst. s(15) = ( 0 0 1 )
|
|
( -1 0 1 )
|
|
( 0 -1 1 )
|
|
|
|
cart. s(15) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 16 90 deg rotation - cart. axis [1,0,0]
|
|
|
|
cryst. s(16) = ( 1 -1 0 )
|
|
( 1 0 -1 )
|
|
( 1 0 0 )
|
|
|
|
cart. s(16) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 17 120 deg rotation - cart. axis [-1,-1,-1]
|
|
|
|
cryst. s(17) = ( -1 0 1 )
|
|
( -1 1 0 )
|
|
( -1 0 0 )
|
|
|
|
cart. s(17) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 18 120 deg rotation - cart. axis [-1,1,1]
|
|
|
|
cryst. s(18) = ( 0 1 0 )
|
|
( 0 0 1 )
|
|
( 1 0 0 )
|
|
|
|
cart. s(18) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 19 120 deg rotation - cart. axis [1,1,-1]
|
|
|
|
cryst. s(19) = ( 1 0 -1 )
|
|
( 0 0 -1 )
|
|
( 0 1 -1 )
|
|
|
|
cart. s(19) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 20 120 deg rotation - cart. axis [1,-1,1]
|
|
|
|
cryst. s(20) = ( 0 -1 0 )
|
|
( 1 -1 0 )
|
|
( 0 -1 1 )
|
|
|
|
cart. s(20) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 21 120 deg rotation - cart. axis [1,1,1]
|
|
|
|
cryst. s(21) = ( 0 0 -1 )
|
|
( 0 1 -1 )
|
|
( 1 0 -1 )
|
|
|
|
cart. s(21) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 22 120 deg rotation - cart. axis [-1,1,-1]
|
|
|
|
cryst. s(22) = ( -1 1 0 )
|
|
( -1 0 0 )
|
|
( -1 0 1 )
|
|
|
|
cart. s(22) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 23 120 deg rotation - cart. axis [1,-1,-1]
|
|
|
|
cryst. s(23) = ( 0 0 1 )
|
|
( 1 0 0 )
|
|
( 0 1 0 )
|
|
|
|
cart. s(23) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 24 120 deg rotation - cart. axis [-1,-1,1]
|
|
|
|
cryst. s(24) = ( 1 -1 0 )
|
|
( 0 -1 1 )
|
|
( 0 -1 0 )
|
|
|
|
cart. s(24) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 25 inversion
|
|
|
|
cryst. s(25) = ( -1 0 0 )
|
|
( 0 -1 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(25) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 26 inv. 180 deg rotation - cart. axis [0,0,1]
|
|
|
|
cryst. s(26) = ( 0 -1 1 )
|
|
( -1 0 1 )
|
|
( 0 0 1 )
|
|
|
|
cart. s(26) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 27 inv. 180 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s(27) = ( 1 0 0 )
|
|
( 1 0 -1 )
|
|
( 1 -1 0 )
|
|
|
|
cart. s(27) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 28 inv. 180 deg rotation - cart. axis [1,0,0]
|
|
|
|
cryst. s(28) = ( 0 1 -1 )
|
|
( 0 1 0 )
|
|
( -1 1 0 )
|
|
|
|
cart. s(28) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 29 inv. 180 deg rotation - cart. axis [1,1,0]
|
|
|
|
cryst. s(29) = ( 0 1 0 )
|
|
( 1 0 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s(29) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 30 inv. 180 deg rotation - cart. axis [1,-1,0]
|
|
|
|
cryst. s(30) = ( 1 0 -1 )
|
|
( 0 1 -1 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(30) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 31 inv. 90 deg rotation - cart. axis [0,0,-1]
|
|
|
|
cryst. s(31) = ( 0 -1 0 )
|
|
( 0 -1 1 )
|
|
( 1 -1 0 )
|
|
|
|
cart. s(31) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 32 inv. 90 deg rotation - cart. axis [0,0,1]
|
|
|
|
cryst. s(32) = ( -1 0 1 )
|
|
( -1 0 0 )
|
|
( -1 1 0 )
|
|
|
|
cart. s(32) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 33 inv. 180 deg rotation - cart. axis [1,0,1]
|
|
|
|
cryst. s(33) = ( 1 0 0 )
|
|
( 0 0 1 )
|
|
( 0 1 0 )
|
|
|
|
cart. s(33) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 34 inv. 180 deg rotation - cart. axis [-1,0,1]
|
|
|
|
cryst. s(34) = ( -1 0 0 )
|
|
( -1 1 0 )
|
|
( -1 0 1 )
|
|
|
|
cart. s(34) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 35 inv. 90 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s(35) = ( 0 -1 1 )
|
|
( 1 -1 0 )
|
|
( 0 -1 0 )
|
|
|
|
cart. s(35) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 36 inv. 90 deg rotation - cart. axis [0,-1,0]
|
|
|
|
cryst. s(36) = ( 0 1 -1 )
|
|
( 0 0 -1 )
|
|
( 1 0 -1 )
|
|
|
|
cart. s(36) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 37 inv. 180 deg rotation - cart. axis [0,1,1]
|
|
|
|
cryst. s(37) = ( 1 -1 0 )
|
|
( 0 -1 0 )
|
|
( 0 -1 1 )
|
|
|
|
cart. s(37) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 38 inv. 180 deg rotation - cart. axis [0,1,-1]
|
|
|
|
cryst. s(38) = ( 0 0 1 )
|
|
( 0 1 0 )
|
|
( 1 0 0 )
|
|
|
|
cart. s(38) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 39 inv. 90 deg rotation - cart. axis [-1,0,0]
|
|
|
|
cryst. s(39) = ( 0 0 -1 )
|
|
( 1 0 -1 )
|
|
( 0 1 -1 )
|
|
|
|
cart. s(39) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 40 inv. 90 deg rotation - cart. axis [1,0,0]
|
|
|
|
cryst. s(40) = ( -1 1 0 )
|
|
( -1 0 1 )
|
|
( -1 0 0 )
|
|
|
|
cart. s(40) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 41 inv. 120 deg rotation - cart. axis [-1,-1,-1]
|
|
|
|
cryst. s(41) = ( 1 0 -1 )
|
|
( 1 -1 0 )
|
|
( 1 0 0 )
|
|
|
|
cart. s(41) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 42 inv. 120 deg rotation - cart. axis [-1,1,1]
|
|
|
|
cryst. s(42) = ( 0 -1 0 )
|
|
( 0 0 -1 )
|
|
( -1 0 0 )
|
|
|
|
cart. s(42) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 43 inv. 120 deg rotation - cart. axis [1,1,-1]
|
|
|
|
cryst. s(43) = ( -1 0 1 )
|
|
( 0 0 1 )
|
|
( 0 -1 1 )
|
|
|
|
cart. s(43) = ( 0.0000000 -1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 44 inv. 120 deg rotation - cart. axis [1,-1,1]
|
|
|
|
cryst. s(44) = ( 0 1 0 )
|
|
( -1 1 0 )
|
|
( 0 1 -1 )
|
|
|
|
cart. s(44) = ( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
|
|
|
|
isym = 45 inv. 120 deg rotation - cart. axis [1,1,1]
|
|
|
|
cryst. s(45) = ( 0 0 1 )
|
|
( 0 -1 1 )
|
|
( -1 0 1 )
|
|
|
|
cart. s(45) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 46 inv. 120 deg rotation - cart. axis [-1,1,-1]
|
|
|
|
cryst. s(46) = ( 1 -1 0 )
|
|
( 1 0 0 )
|
|
( 1 0 -1 )
|
|
|
|
cart. s(46) = ( 0.0000000 0.0000000 -1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 47 inv. 120 deg rotation - cart. axis [1,-1,-1]
|
|
|
|
cryst. s(47) = ( 0 0 -1 )
|
|
( -1 0 0 )
|
|
( 0 -1 0 )
|
|
|
|
cart. s(47) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 )
|
|
|
|
|
|
isym = 48 inv. 120 deg rotation - cart. axis [-1,-1,1]
|
|
|
|
cryst. s(48) = ( -1 1 0 )
|
|
( 0 1 -1 )
|
|
( 0 1 0 )
|
|
|
|
cart. s(48) = ( 0.0000000 0.0000000 1.0000000 )
|
|
( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
|
|
|
|
double point group O_h (m-3m)
|
|
there are 16 classes and 6 irreducible representations
|
|
the character table:
|
|
|
|
E -E 8C3 -8C3 3C2 6C4 -6C4 6C2' i -i 8S6 -8S6
|
|
-3C2 -6C2'
|
|
G_6+ 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 2.00 -2.00 1.00 -1.00
|
|
G_7+ 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 2.00 -2.00 1.00 -1.00
|
|
G_8+ 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 4.00 -4.00 -1.00 1.00
|
|
G_6- 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 -2.00 2.00 -1.00 1.00
|
|
G_7- 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 -2.00 2.00 -1.00 1.00
|
|
G_8- 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 -4.00 4.00 1.00 -1.00
|
|
|
|
3s_h 6S4 -6S4 6s_d
|
|
-3s_h -6s_d
|
|
G_6+ 0.00 1.41 -1.41 0.00
|
|
G_7+ 0.00 -1.41 1.41 0.00
|
|
G_8+ 0.00 0.00 0.00 0.00
|
|
G_6- 0.00 -1.41 1.41 0.00
|
|
G_7- 0.00 1.41 -1.41 0.00
|
|
G_8- 0.00 0.00 0.00 0.00
|
|
|
|
the symmetry operations in each class and the name of the first element:
|
|
|
|
E 1
|
|
identity
|
|
-E -1
|
|
identity E
|
|
8C3 17 19 20 18 24 21 22 23
|
|
120 deg rotation - cart. axis [-1,-1,-1]
|
|
-8C3 -17 -19 -20 -18 -24 -21 -22 -23
|
|
120 deg rotation - cart. axis [-1,-1,-1] E
|
|
3C2 -3C2 2 -2 4 -4 3 -3
|
|
180 deg rotation - cart. axis [0,0,1]
|
|
6C4 7 8 15 16 12 11
|
|
90 deg rotation - cart. axis [0,0,-1]
|
|
-6C4 -7 -8 -15 -16 -12 -11
|
|
90 deg rotation - cart. axis [0,0,-1] E
|
|
6C2'-6C2' 5 -5 6 -6 14 -13 -14 13 -10 -9 10 9
|
|
180 deg rotation - cart. axis [1,1,0]
|
|
i 25
|
|
inversion
|
|
-i -25
|
|
inversion E
|
|
8S6 41 43 44 42 48 45 46 47
|
|
inv. 120 deg rotation - cart. axis [-1,-1,-1]
|
|
-8S6 -41 -43 -44 -42 -48 -45 -46 -47
|
|
inv. 120 deg rotation - cart. axis [-1,-1,-1] E
|
|
3s_h-3s_h 26 -26 28 -28 27 -27
|
|
inv. 180 deg rotation - cart. axis [0,0,1]
|
|
6S4 31 32 39 40 36 35
|
|
inv. 90 deg rotation - cart. axis [0,0,-1]
|
|
-6S4 -31 -32 -39 -40 -36 -35
|
|
inv. 90 deg rotation - cart. axis [0,0,-1] E
|
|
6s_d-6s_d 29 -29 30 -30 38 -37 -38 37 -34 -33 34 33
|
|
inv. 180 deg rotation - cart. axis [1,1,0]
|
|
|
|
Cartesian axes
|
|
|
|
site n. atom positions (alat units)
|
|
1 Pb tau( 1) = ( 0.0000000 0.0000000 0.0000000 )
|
|
|
|
Crystallographic axes
|
|
|
|
site n. atom positions (cryst. coord.)
|
|
1 Pb tau( 1) = ( 0.0000000 0.0000000 0.0000000 )
|
|
|
|
number of k points= 4 Methfessel-Paxton smearing, width (Ry)= 0.0250
|
|
cart. coord. in units 2pi/alat
|
|
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0370370
|
|
k( 2) = ( -0.3333333 0.3333333 -0.3333333), wk = 0.2962963
|
|
k( 3) = ( 0.0000000 0.6666667 0.0000000), wk = 0.2222222
|
|
k( 4) = ( 0.6666667 -0.0000000 0.6666667), wk = 0.4444444
|
|
|
|
cryst. coord.
|
|
k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0370370
|
|
k( 2) = ( 0.0000000 0.0000000 0.3333333), wk = 0.2962963
|
|
k( 3) = ( 0.0000000 0.3333333 0.3333333), wk = 0.2222222
|
|
k( 4) = ( 0.0000000 0.3333333 -0.3333333), wk = 0.4444444
|
|
|
|
Dense grid: 12387 G-vectors FFT dimensions: ( 36, 36, 36)
|
|
|
|
Largest allocated arrays est. size (Mb) dimensions
|
|
Kohn-Sham Wavefunctions 1.06 Mb ( 3160, 22)
|
|
NL pseudopotentials 0.39 Mb ( 1580, 16)
|
|
Each V/rho on FFT grid 0.71 Mb ( 46656)
|
|
Each G-vector array 0.09 Mb ( 12387)
|
|
G-vector shells 0.00 Mb ( 175)
|
|
Largest temporary arrays est. size (Mb) dimensions
|
|
Auxiliary wavefunctions 4.24 Mb ( 3160, 88)
|
|
Each subspace H/S matrix 0.12 Mb ( 88, 88)
|
|
Each <psi_i|beta_j> matrix 0.01 Mb ( 16, 2, 22)
|
|
Arrays for rho mixing 5.70 Mb ( 46656, 8)
|
|
|
|
Initial potential from superposition of free atoms
|
|
|
|
starting charge 13.99923, renormalised to 14.00000
|
|
Starting wfc are 18 randomized atomic wfcs + 4 random wfc
|
|
|
|
total cpu time spent up to now is 1.1 secs
|
|
|
|
per-process dynamical memory: 24.7 Mb
|
|
|
|
Self-consistent Calculation
|
|
|
|
iteration # 1 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.00E-02, avg # of iterations = 4.2
|
|
|
|
Threshold (ethr) on eigenvalues was too large:
|
|
Diagonalizing with lowered threshold
|
|
|
|
Davidson diagonalization with overlap
|
|
ethr = 4.59E-05, avg # of iterations = 5.5
|
|
|
|
total cpu time spent up to now is 4.9 secs
|
|
|
|
total energy = -119.59560943 Ry
|
|
Harris-Foulkes estimate = -119.59933786 Ry
|
|
estimated scf accuracy < 0.00640183 Ry
|
|
|
|
iteration # 2 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 4.57E-05, avg # of iterations = 2.0
|
|
|
|
total cpu time spent up to now is 6.7 secs
|
|
|
|
total energy = -119.59667894 Ry
|
|
Harris-Foulkes estimate = -119.59683963 Ry
|
|
estimated scf accuracy < 0.00032296 Ry
|
|
|
|
iteration # 3 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 2.31E-06, avg # of iterations = 2.8
|
|
|
|
total cpu time spent up to now is 8.6 secs
|
|
|
|
total energy = -119.59672281 Ry
|
|
Harris-Foulkes estimate = -119.59672661 Ry
|
|
estimated scf accuracy < 0.00001460 Ry
|
|
|
|
iteration # 4 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.04E-07, avg # of iterations = 1.0
|
|
|
|
total cpu time spent up to now is 10.1 secs
|
|
|
|
total energy = -119.59672360 Ry
|
|
Harris-Foulkes estimate = -119.59672359 Ry
|
|
estimated scf accuracy < 0.00000007 Ry
|
|
|
|
iteration # 5 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 4.91E-10, avg # of iterations = 2.2
|
|
|
|
total cpu time spent up to now is 12.0 secs
|
|
|
|
total energy = -119.59672363 Ry
|
|
Harris-Foulkes estimate = -119.59672364 Ry
|
|
estimated scf accuracy < 0.00000001 Ry
|
|
|
|
iteration # 6 ecut= 60.00 Ry beta=0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 8.53E-11, avg # of iterations = 2.0
|
|
|
|
total cpu time spent up to now is 13.8 secs
|
|
|
|
End of self-consistent calculation
|
|
|
|
k = 0.0000 0.0000 0.0000 ( 1471 PWs) bands (ev):
|
|
|
|
-6.9829 -6.9829 -6.9829 -6.9829 -4.4809 -4.4809 -4.3693 -4.3693
|
|
-4.3693 -4.3693 -0.5407 -0.5407 16.4807 16.4807 19.9086 19.9086
|
|
19.9086 19.9086 21.9538 21.9538 23.0227 23.0227
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
1.0000 1.0000 1.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000
|
|
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
|
|
|
|
k =-0.3333 0.3333-0.3333 ( 1540 PWs) bands (ev):
|
|
|
|
-6.9717 -6.9717 -6.9212 -6.9212 -4.4494 -4.4494 -4.3804 -4.3804
|
|
-4.2501 -4.2501 1.6151 1.6151 9.6769 9.6769 16.3684 16.3684
|
|
18.3349 18.3349 20.9084 20.9084 22.5924 22.5924
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
1.0000 1.0000 1.0000 1.0000 1.0343 1.0343 -0.0000 -0.0000
|
|
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
|
|
|
|
k = 0.0000 0.6667 0.0000 ( 1580 PWs) bands (ev):
|
|
|
|
-6.9619 -6.9619 -6.9142 -6.9142 -4.4458 -4.4458 -4.3762 -4.3762
|
|
-4.1926 -4.1926 2.3307 2.3307 11.0145 11.0145 13.8283 13.8283
|
|
15.0491 15.0491 18.1329 18.1329 19.0118 19.0118
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
1.0000 1.0000 1.0000 1.0000 -0.0007 -0.0007 -0.0000 -0.0000
|
|
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
|
|
|
|
k = 0.6667-0.0000 0.6667 ( 1568 PWs) bands (ev):
|
|
|
|
-6.9428 -6.9428 -6.8996 -6.8996 -4.4177 -4.4177 -4.3345 -4.3345
|
|
-4.2406 -4.2406 4.2194 4.2194 9.0297 9.0297 10.0997 10.0997
|
|
15.1164 15.1164 18.7741 18.7741 25.8614 25.8614
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.5607 0.5607
|
|
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
|
|
|
|
the Fermi energy is 10.1242 ev
|
|
|
|
! total energy = -119.59672363 Ry
|
|
Harris-Foulkes estimate = -119.59672363 Ry
|
|
estimated scf accuracy < 1.6E-11 Ry
|
|
|
|
The total energy is the sum of the following terms:
|
|
|
|
one-electron contribution = -33.13354618 Ry
|
|
hartree contribution = 25.02819815 Ry
|
|
xc contribution = -14.54931403 Ry
|
|
ewald contribution = -96.93991011 Ry
|
|
smearing contrib. (-TS) = -0.00215146 Ry
|
|
|
|
convergence has been achieved in 6 iterations
|
|
|
|
Forces acting on atoms (Ry/au):
|
|
|
|
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
The non-local contrib. to forces
|
|
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
The ionic contribution to forces
|
|
atom 1 type 1 force = -0.00000000 -0.00000000 -0.00000000
|
|
The local contribution to forces
|
|
atom 1 type 1 force = 0.00000000 -0.00000000 -0.00000000
|
|
The core correction contribution to forces
|
|
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
The Hubbard contrib. to forces
|
|
atom 1 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
The SCF correction term to forces
|
|
atom 1 type 1 force = -0.00000000 0.00000000 -0.00000000
|
|
|
|
Total force = 0.000000 Total SCF correction = 0.000000
|
|
|
|
|
|
entering subroutine stress ...
|
|
|
|
total stress (Ry/bohr**3) (kbar) P= -3.68
|
|
-0.00002504 -0.00000000 0.00000000 -3.68 -0.00 0.00
|
|
-0.00000000 -0.00002504 -0.00000000 -0.00 -3.68 -0.00
|
|
0.00000000 -0.00000000 -0.00002504 0.00 -0.00 -3.68
|
|
|
|
kinetic stress (kbar) 34590.62 -0.00 0.00
|
|
-0.00 34590.62 0.00
|
|
0.00 0.00 34590.62
|
|
|
|
local stress (kbar) 32802.63 0.00 0.00
|
|
0.00 32802.63 -0.00
|
|
0.00 -0.00 32802.63
|
|
|
|
nonloc. stress (kbar) -46356.98 0.00 0.00
|
|
0.00 -46356.98 -0.00
|
|
0.00 -0.00 -46356.98
|
|
|
|
hartree stress (kbar) 6162.50 -0.00 -0.00
|
|
-0.00 6162.50 0.00
|
|
-0.00 0.00 6162.50
|
|
|
|
exc-cor stress (kbar) -3333.70 0.00 0.00
|
|
0.00 -3333.70 0.00
|
|
0.00 0.00 -3333.70
|
|
|
|
corecor stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
ewald stress (kbar) -23868.75 -0.00 0.00
|
|
-0.00 -23868.75 0.00
|
|
0.00 0.00 -23868.75
|
|
|
|
hubbard stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
london stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
XDM stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
dft-nl stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
TS-vdW stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
|
|
|
|
Writing output data file pb.save
|
|
|
|
init_run : 1.94s CPU 0.97s WALL ( 1 calls)
|
|
electrons : 25.56s CPU 12.80s WALL ( 1 calls)
|
|
forces : 0.09s CPU 0.05s WALL ( 1 calls)
|
|
stress : 0.47s CPU 0.24s WALL ( 1 calls)
|
|
|
|
Called by init_run:
|
|
wfcinit : 1.12s CPU 0.56s WALL ( 1 calls)
|
|
potinit : 0.07s CPU 0.03s WALL ( 1 calls)
|
|
|
|
Called by electrons:
|
|
c_bands : 21.21s CPU 10.62s WALL ( 7 calls)
|
|
sum_band : 4.08s CPU 2.04s WALL ( 7 calls)
|
|
v_of_rho : 0.10s CPU 0.05s WALL ( 7 calls)
|
|
v_h : 0.03s CPU 0.02s WALL ( 7 calls)
|
|
v_xc : 0.07s CPU 0.03s WALL ( 7 calls)
|
|
mix_rho : 0.15s CPU 0.07s WALL ( 7 calls)
|
|
|
|
Called by c_bands:
|
|
init_us_2 : 0.13s CPU 0.06s WALL ( 68 calls)
|
|
cegterg : 20.70s CPU 10.36s WALL ( 28 calls)
|
|
|
|
Called by sum_band:
|
|
|
|
Called by *egterg:
|
|
h_psi : 19.60s CPU 9.80s WALL ( 111 calls)
|
|
g_psi : 0.12s CPU 0.06s WALL ( 79 calls)
|
|
cdiaghg : 0.20s CPU 0.10s WALL ( 103 calls)
|
|
cegterg:over : 0.76s CPU 0.38s WALL ( 79 calls)
|
|
cegterg:upda : 0.42s CPU 0.21s WALL ( 79 calls)
|
|
cegterg:last : 0.27s CPU 0.13s WALL ( 28 calls)
|
|
|
|
Called by h_psi:
|
|
h_psi:vloc : 19.13s CPU 9.57s WALL ( 111 calls)
|
|
h_psi:vnl : 0.44s CPU 0.22s WALL ( 111 calls)
|
|
add_vuspsi : 0.22s CPU 0.11s WALL ( 111 calls)
|
|
|
|
General routines
|
|
calbec : 0.28s CPU 0.14s WALL ( 131 calls)
|
|
fft : 0.33s CPU 0.17s WALL ( 79 calls)
|
|
fftw : 20.87s CPU 10.44s WALL ( 7784 calls)
|
|
|
|
Parallel routines
|
|
fft_scatter : 0.85s CPU 0.42s WALL ( 7863 calls)
|
|
|
|
PWSCF : 28.30s CPU 14.20s WALL
|
|
|
|
|
|
This run was terminated on: 16:49:59 21Jun2016
|
|
|
|
=------------------------------------------------------------------------------=
|
|
JOB DONE.
|
|
=------------------------------------------------------------------------------=
|