mirror of https://gitlab.com/QEF/q-e.git
915 lines
38 KiB
Plaintext
915 lines
38 KiB
Plaintext
|
|
Program PWSCF v.6.4.1 starts on 23Sep2019 at 12:51:52
|
|
|
|
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 4 processors
|
|
|
|
MPI processes distributed on 1 nodes
|
|
R & G space division: proc/nbgrp/npool/nimage = 4
|
|
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
|
|
|
|
IMPORTANT: XC functional enforced from input :
|
|
Exchange-correlation= VDW-DF
|
|
( 1 4 4 0 1 0 0)
|
|
Any further DFT definition will be discarded
|
|
Please, verify this is what you really want
|
|
|
|
file C.pbe-rrkjus.UPF: wavefunction(s) 2S 2P renormalized
|
|
|
|
Subspace diagonalization in iterative solution of the eigenvalue problem:
|
|
a serial algorithm will be used
|
|
|
|
|
|
Parallelization info
|
|
--------------------
|
|
sticks: dense smooth PW G-vecs: dense smooth PW
|
|
Min 65 46 15 2392 1302 275
|
|
Max 67 48 16 2397 1304 278
|
|
Sum 265 187 61 9583 5211 1107
|
|
|
|
|
|
|
|
bravais-lattice index = 4
|
|
lattice parameter (alat) = 4.6512 a.u.
|
|
unit-cell volume = 237.5783 (a.u.)^3
|
|
number of atoms/cell = 4
|
|
number of atomic types = 1
|
|
number of electrons = 16.00
|
|
number of Kohn-Sham states= 8
|
|
kinetic-energy cutoff = 29.8711 Ry
|
|
charge density cutoff = 179.2268 Ry
|
|
convergence threshold = 1.0E-08
|
|
mixing beta = 0.7000
|
|
number of iterations used = 8 plain mixing
|
|
Exchange-correlation= VDW-DF
|
|
( 1 4 4 0 1 0 0)
|
|
|
|
celldm(1)= 4.651170 celldm(2)= 0.000000 celldm(3)= 2.726400
|
|
celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000
|
|
|
|
crystal axes: (cart. coord. in units of alat)
|
|
a(1) = ( 1.000000 0.000000 0.000000 )
|
|
a(2) = ( -0.500000 0.866025 0.000000 )
|
|
a(3) = ( 0.000000 0.000000 2.726400 )
|
|
|
|
reciprocal axes: (cart. coord. in units 2 pi/alat)
|
|
b(1) = ( 1.000000 0.577350 0.000000 )
|
|
b(2) = ( 0.000000 1.154701 0.000000 )
|
|
b(3) = ( 0.000000 0.000000 0.366784 )
|
|
|
|
|
|
PseudoPot. # 1 for C read from file:
|
|
/benchmarks/pseudo/C.pbe-rrkjus.UPF
|
|
MD5 check sum: c9ac5a99bc85b198593446162950cd17
|
|
Pseudo is Ultrasoft, Zval = 4.0
|
|
Generated by new atomic code, or converted to UPF format
|
|
Using radial grid of 627 points, 4 beta functions with:
|
|
l(1) = 0
|
|
l(2) = 0
|
|
l(3) = 1
|
|
l(4) = 1
|
|
Q(r) pseudized with 0 coefficients
|
|
|
|
|
|
atomic species valence mass pseudopotential
|
|
C 4.00 12.00000 C ( 1.00)
|
|
|
|
24 Sym. Ops., with inversion, found (12 have fractional translation)
|
|
|
|
|
|
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) = ( -1 0 0 ) f =( 0.0000000 )
|
|
( 0 -1 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 3 180 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s( 3) = ( -1 0 0 )
|
|
( 1 1 0 )
|
|
( 0 0 -1 )
|
|
|
|
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) = ( 1 0 0 ) f =( 0.0000000 )
|
|
( -1 -1 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s( 4) = ( 1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 5 60 deg rotation - cryst. axis [0,0,1]
|
|
|
|
cryst. s( 5) = ( 1 1 0 ) f =( 0.0000000 )
|
|
( -1 0 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s( 5) = ( 0.5000000 -0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( 0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 6 60 deg rotation - cryst. axis [0,0,-1]
|
|
|
|
cryst. s( 6) = ( 0 -1 0 ) f =( 0.0000000 )
|
|
( 1 1 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s( 6) = ( 0.5000000 0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( -0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 7 120 deg rotation - cryst. axis [0,0,1]
|
|
|
|
cryst. s( 7) = ( 0 1 0 )
|
|
( -1 -1 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s( 7) = ( -0.5000000 -0.8660254 0.0000000 )
|
|
( 0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 8 120 deg rotation - cryst. axis [0,0,-1]
|
|
|
|
cryst. s( 8) = ( -1 -1 0 )
|
|
( 1 0 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s( 8) = ( -0.5000000 0.8660254 0.0000000 )
|
|
( -0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 9 180 deg rotation - cryst. axis [1,-1,0]
|
|
|
|
cryst. s( 9) = ( 0 -1 0 )
|
|
( -1 0 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s( 9) = ( 0.5000000 -0.8660254 0.0000000 )
|
|
( -0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 10 180 deg rotation - cryst. axis [2,1,0]
|
|
|
|
cryst. s(10) = ( 1 1 0 )
|
|
( 0 -1 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(10) = ( 0.5000000 0.8660254 0.0000000 )
|
|
( 0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 11 180 deg rotation - cryst. axis [0,1,0]
|
|
|
|
cryst. s(11) = ( -1 -1 0 ) f =( 0.0000000 )
|
|
( 0 1 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s(11) = ( -0.5000000 -0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( -0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 12 180 deg rotation - cryst. axis [1,1,0]
|
|
|
|
cryst. s(12) = ( 0 1 0 ) f =( 0.0000000 )
|
|
( 1 0 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s(12) = ( -0.5000000 0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( 0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 13 inversion
|
|
|
|
cryst. s(13) = ( -1 0 0 ) f =( 0.0000000 )
|
|
( 0 -1 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s(13) = ( -1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 14 inv. 180 deg rotation - cart. axis [0,0,1]
|
|
|
|
cryst. s(14) = ( 1 0 0 )
|
|
( 0 1 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(14) = ( 1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 15 inv. 180 deg rotation - cart. axis [0,1,0]
|
|
|
|
cryst. s(15) = ( 1 0 0 ) f =( 0.0000000 )
|
|
( -1 -1 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s(15) = ( 1.0000000 0.0000000 0.0000000 ) f =( 0.0000000 )
|
|
( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 16 inv. 180 deg rotation - cart. axis [1,0,0]
|
|
|
|
cryst. s(16) = ( -1 0 0 )
|
|
( 1 1 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s(16) = ( -1.0000000 0.0000000 0.0000000 )
|
|
( 0.0000000 1.0000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 17 inv. 60 deg rotation - cryst. axis [0,0,1]
|
|
|
|
cryst. s(17) = ( -1 -1 0 )
|
|
( 1 0 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(17) = ( -0.5000000 0.8660254 0.0000000 )
|
|
( -0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 18 inv. 60 deg rotation - cryst. axis [0,0,-1]
|
|
|
|
cryst. s(18) = ( 0 1 0 )
|
|
( -1 -1 0 )
|
|
( 0 0 -1 )
|
|
|
|
cart. s(18) = ( -0.5000000 -0.8660254 0.0000000 )
|
|
( 0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 )
|
|
|
|
|
|
isym = 19 inv. 120 deg rotation - cryst. axis [0,0,1]
|
|
|
|
cryst. s(19) = ( 0 -1 0 ) f =( 0.0000000 )
|
|
( 1 1 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s(19) = ( 0.5000000 0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( -0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 20 inv. 120 deg rotation - cryst. axis [0,0,-1]
|
|
|
|
cryst. s(20) = ( 1 1 0 ) f =( 0.0000000 )
|
|
( -1 0 0 ) ( 0.0000000 )
|
|
( 0 0 -1 ) ( 0.5000000 )
|
|
|
|
cart. s(20) = ( 0.5000000 -0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( 0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 -1.0000000 ) ( 1.3632000 )
|
|
|
|
|
|
isym = 21 inv. 180 deg rotation - cryst. axis [1,-1,0]
|
|
|
|
cryst. s(21) = ( 0 1 0 ) f =( 0.0000000 )
|
|
( 1 0 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s(21) = ( -0.5000000 0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( 0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 22 inv. 180 deg rotation - cryst. axis [2,1,0]
|
|
|
|
cryst. s(22) = ( -1 -1 0 ) f =( 0.0000000 )
|
|
( 0 1 0 ) ( 0.0000000 )
|
|
( 0 0 1 ) ( -0.5000000 )
|
|
|
|
cart. s(22) = ( -0.5000000 -0.8660254 0.0000000 ) f =( 0.0000000 )
|
|
( -0.8660254 0.5000000 0.0000000 ) ( 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 ) ( -1.3632000 )
|
|
|
|
|
|
isym = 23 inv. 180 deg rotation - cryst. axis [0,1,0]
|
|
|
|
cryst. s(23) = ( 1 1 0 )
|
|
( 0 -1 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s(23) = ( 0.5000000 0.8660254 0.0000000 )
|
|
( 0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
isym = 24 inv. 180 deg rotation - cryst. axis [1,1,0]
|
|
|
|
cryst. s(24) = ( 0 -1 0 )
|
|
( -1 0 0 )
|
|
( 0 0 1 )
|
|
|
|
cart. s(24) = ( 0.5000000 -0.8660254 0.0000000 )
|
|
( -0.8660254 -0.5000000 0.0000000 )
|
|
( 0.0000000 0.0000000 1.0000000 )
|
|
|
|
|
|
point group D_6h(6/mmm)
|
|
there are 12 classes
|
|
the character table:
|
|
|
|
E 2C6 2C3 C2 3C2' 3C2'' i 2S3 2S6 s_h 3s_d 3s_v
|
|
A_1g 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00
|
|
A_2g 1.00 1.00 1.00 1.00 -1.00 -1.00 1.00 1.00 1.00 1.00 -1.00 -1.00
|
|
B_1g 1.00 -1.00 1.00 -1.00 1.00 -1.00 1.00 -1.00 1.00 -1.00 1.00 -1.00
|
|
B_2g 1.00 -1.00 1.00 -1.00 -1.00 1.00 1.00 -1.00 1.00 -1.00 -1.00 1.00
|
|
E_1g 2.00 1.00 -1.00 -2.00 0.00 0.00 2.00 1.00 -1.00 -2.00 0.00 0.00
|
|
E_2g 2.00 -1.00 -1.00 2.00 0.00 0.00 2.00 -1.00 -1.00 2.00 0.00 0.00
|
|
A_1u 1.00 1.00 1.00 1.00 1.00 1.00 -1.00 -1.00 -1.00 -1.00 -1.00 -1.00
|
|
A_2u 1.00 1.00 1.00 1.00 -1.00 -1.00 -1.00 -1.00 -1.00 -1.00 1.00 1.00
|
|
B_1u 1.00 -1.00 1.00 -1.00 1.00 -1.00 -1.00 1.00 -1.00 1.00 -1.00 1.00
|
|
B_2u 1.00 -1.00 1.00 -1.00 -1.00 1.00 -1.00 1.00 -1.00 1.00 1.00 -1.00
|
|
E_1u 2.00 1.00 -1.00 -2.00 0.00 0.00 -2.00 -1.00 1.00 2.00 0.00 0.00
|
|
E_2u 2.00 -1.00 -1.00 2.00 0.00 0.00 -2.00 1.00 1.00 -2.00 0.00 0.00
|
|
|
|
the symmetry operations in each class and the name of the first element:
|
|
|
|
E 1
|
|
identity
|
|
2C6 5 6
|
|
60 deg rotation - cryst. axis [0,0,1]
|
|
2C3 7 8
|
|
120 deg rotation - cryst. axis [0,0,1]
|
|
C2 2
|
|
180 deg rotation - cart. axis [0,0,1]
|
|
3C2' 4 12 11
|
|
180 deg rotation - cart. axis [1,0,0]
|
|
3C2'' 3 9 10
|
|
180 deg rotation - cart. axis [0,1,0]
|
|
i 13
|
|
inversion
|
|
2S3 17 18
|
|
inv. 60 deg rotation - cryst. axis [0,0,1]
|
|
2S6 19 20
|
|
inv. 120 deg rotation - cryst. axis [0,0,1]
|
|
s_h 14
|
|
inv. 180 deg rotation - cart. axis [0,0,1]
|
|
3s_d 16 24 23
|
|
inv. 180 deg rotation - cart. axis [1,0,0]
|
|
3s_v 15 21 22
|
|
inv. 180 deg rotation - cart. axis [0,1,0]
|
|
|
|
Cartesian axes
|
|
|
|
site n. atom positions (alat units)
|
|
1 C tau( 1) = ( 0.0000000 0.0000000 0.0000000 )
|
|
2 C tau( 2) = ( 0.0000000 0.5773503 0.0000000 )
|
|
3 C tau( 3) = ( 0.0000000 0.0000000 1.3632000 )
|
|
4 C tau( 4) = ( 0.5000000 0.2886751 1.3632000 )
|
|
|
|
Crystallographic axes
|
|
|
|
site n. atom positions (cryst. coord.)
|
|
1 C tau( 1) = ( 0.0000000 0.0000000 0.0000000 )
|
|
2 C tau( 2) = ( 0.3333333 0.6666667 0.0000000 )
|
|
3 C tau( 3) = ( 0.0000000 0.0000000 0.5000000 )
|
|
4 C tau( 4) = ( 0.6666667 0.3333333 0.5000000 )
|
|
|
|
number of k points= 12
|
|
cart. coord. in units 2pi/alat
|
|
k( 1) = ( 0.1250000 0.2165064 0.0458480), wk = 0.1250000
|
|
k( 2) = ( 0.1250000 0.2165064 0.1375440), wk = 0.1250000
|
|
k( 3) = ( 0.1250000 0.5051815 0.0458480), wk = 0.2500000
|
|
k( 4) = ( 0.1250000 0.5051815 0.1375440), wk = 0.2500000
|
|
k( 5) = ( 0.1250000 -0.3608439 0.0458480), wk = 0.2500000
|
|
k( 6) = ( 0.1250000 -0.3608439 0.1375440), wk = 0.2500000
|
|
k( 7) = ( 0.1250000 -0.0721688 0.0458480), wk = 0.1250000
|
|
k( 8) = ( 0.1250000 -0.0721688 0.1375440), wk = 0.1250000
|
|
k( 9) = ( 0.3750000 0.6495191 0.0458480), wk = 0.1250000
|
|
k( 10) = ( 0.3750000 0.6495191 0.1375440), wk = 0.1250000
|
|
k( 11) = ( 0.3750000 -0.2165064 0.0458480), wk = 0.1250000
|
|
k( 12) = ( 0.3750000 -0.2165064 0.1375440), wk = 0.1250000
|
|
|
|
cryst. coord.
|
|
k( 1) = ( 0.1250000 0.1250000 0.1250000), wk = 0.1250000
|
|
k( 2) = ( 0.1250000 0.1250000 0.3750000), wk = 0.1250000
|
|
k( 3) = ( 0.1250000 0.3750000 0.1250000), wk = 0.2500000
|
|
k( 4) = ( 0.1250000 0.3750000 0.3750000), wk = 0.2500000
|
|
k( 5) = ( 0.1250000 -0.3750000 0.1250000), wk = 0.2500000
|
|
k( 6) = ( 0.1250000 -0.3750000 0.3750000), wk = 0.2500000
|
|
k( 7) = ( 0.1250000 -0.1250000 0.1250000), wk = 0.1250000
|
|
k( 8) = ( 0.1250000 -0.1250000 0.3750000), wk = 0.1250000
|
|
k( 9) = ( 0.3750000 0.3750000 0.1250000), wk = 0.1250000
|
|
k( 10) = ( 0.3750000 0.3750000 0.3750000), wk = 0.1250000
|
|
k( 11) = ( 0.3750000 -0.3750000 0.1250000), wk = 0.1250000
|
|
k( 12) = ( 0.3750000 -0.3750000 0.3750000), wk = 0.1250000
|
|
|
|
Dense grid: 9583 G-vectors FFT dimensions: ( 20, 20, 60)
|
|
|
|
Smooth grid: 5211 G-vectors FFT dimensions: ( 18, 18, 48)
|
|
|
|
Dynamical RAM for wfc: 0.02 MB
|
|
|
|
Dynamical RAM for wfc (w. buffer): 0.26 MB
|
|
|
|
Dynamical RAM for str. fact: 0.04 MB
|
|
|
|
Dynamical RAM for local pot: 0.00 MB
|
|
|
|
Dynamical RAM for nlocal pot: 0.08 MB
|
|
|
|
Dynamical RAM for qrad: 0.23 MB
|
|
|
|
Dynamical RAM for rho,v,vnew: 0.25 MB
|
|
|
|
Dynamical RAM for rhoin: 0.08 MB
|
|
|
|
Dynamical RAM for rho*nmix: 0.59 MB
|
|
|
|
Dynamical RAM for G-vectors: 0.15 MB
|
|
|
|
Dynamical RAM for h,s,v(r/c): 0.05 MB
|
|
|
|
Dynamical RAM for <psi|beta>: 0.00 MB
|
|
|
|
Dynamical RAM for psi: 0.08 MB
|
|
|
|
Dynamical RAM for hpsi: 0.08 MB
|
|
|
|
Dynamical RAM for spsi: 0.08 MB
|
|
|
|
Dynamical RAM for wfcinit/wfcrot: 0.09 MB
|
|
|
|
Dynamical RAM for addusdens: 1.72 MB
|
|
|
|
Dynamical RAM for addusforce: 1.96 MB
|
|
|
|
Dynamical RAM for addusstress: 1.83 MB
|
|
|
|
Estimated static dynamical RAM per process > 1.97 MB
|
|
|
|
Estimated max dynamical RAM per process > 3.93 MB
|
|
|
|
Estimated total dynamical RAM > 15.71 MB
|
|
|
|
Initial potential from superposition of free atoms
|
|
|
|
starting charge 15.99979, renormalised to 16.00000
|
|
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% %
|
|
% You are using vdW-DF, which was implemented by the Thonhauser group. %
|
|
% Please cite the following two papers that made this development %
|
|
% possible and the two reviews that describe the various versions: %
|
|
% %
|
|
% T. Thonhauser et al., PRL 115, 136402 (2015). %
|
|
% T. Thonhauser et al., PRB 76, 125112 (2007). %
|
|
% K. Berland et al., Rep. Prog. Phys. 78, 066501 (2015). %
|
|
% D.C. Langreth et al., J. Phys.: Condens. Matter 21, 084203 (2009). %
|
|
% %
|
|
% %
|
|
% If you are calculating the stress with vdW-DF, please also cite: %
|
|
% %
|
|
% R. Sabatini et al., J. Phys.: Condens. Matter 24, 424209 (2012). %
|
|
% %
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
|
|
Carrying out vdW-DF run using the following parameters:
|
|
Nqs = 20 Npoints = 1024 r_max = 100.000
|
|
q_mesh = 0.00001000 0.04494208 0.09755937 0.15916263
|
|
0.23128650 0.31572767 0.41458969 0.53033537
|
|
0.66584808 0.82450364 1.01025438 1.22772762
|
|
1.48234092 1.78043706 2.12944203 2.53805004
|
|
3.01644009 3.57652955 4.23227104 5.00000000
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.12042970 Ry
|
|
-----------------------------------------------
|
|
|
|
Starting wfcs are 16 randomized atomic wfcs
|
|
|
|
total cpu time spent up to now is 25.2 secs
|
|
|
|
Self-consistent Calculation
|
|
|
|
iteration # 1 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.00E-02, avg # of iterations = 2.0
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15323944 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.3 secs
|
|
|
|
total energy = -45.81600632 Ry
|
|
Harris-Foulkes estimate = -46.05767994 Ry
|
|
estimated scf accuracy < 0.43244218 Ry
|
|
|
|
iteration # 2 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 2.70E-03, avg # of iterations = 2.0
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15284658 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.4 secs
|
|
|
|
total energy = -45.88048395 Ry
|
|
Harris-Foulkes estimate = -45.87924325 Ry
|
|
estimated scf accuracy < 0.00555529 Ry
|
|
|
|
iteration # 3 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 3.47E-05, avg # of iterations = 2.2
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15306398 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.5 secs
|
|
|
|
total energy = -45.88130574 Ry
|
|
Harris-Foulkes estimate = -45.88110842 Ry
|
|
estimated scf accuracy < 0.00041559 Ry
|
|
|
|
iteration # 4 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 2.60E-06, avg # of iterations = 1.9
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15313476 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.5 secs
|
|
|
|
total energy = -45.88135474 Ry
|
|
Harris-Foulkes estimate = -45.88135242 Ry
|
|
estimated scf accuracy < 0.00000282 Ry
|
|
|
|
iteration # 5 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.77E-08, avg # of iterations = 3.1
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15313404 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.6 secs
|
|
|
|
total energy = -45.88135724 Ry
|
|
Harris-Foulkes estimate = -45.88135723 Ry
|
|
estimated scf accuracy < 0.00000012 Ry
|
|
|
|
iteration # 6 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 7.25E-10, avg # of iterations = 3.1
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15313435 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.7 secs
|
|
|
|
total energy = -45.88135727 Ry
|
|
Harris-Foulkes estimate = -45.88135728 Ry
|
|
estimated scf accuracy < 0.00000005 Ry
|
|
|
|
iteration # 7 ecut= 29.87 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 3.16E-10, avg # of iterations = 2.5
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15313147 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.7 secs
|
|
|
|
End of self-consistent calculation
|
|
|
|
k = 0.1250 0.2165 0.0458 ( 646 PWs) bands (ev):
|
|
|
|
-11.5504 -11.2965 -0.0904 0.6272 0.6661 1.6106 1.7135 1.7406
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250 0.2165 0.1375 ( 654 PWs) bands (ev):
|
|
|
|
-11.4789 -11.3739 0.3513 0.6385 0.6557 1.0478 1.7208 1.7320
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250 0.5052 0.0458 ( 662 PWs) bands (ev):
|
|
|
|
-8.0411 -7.8523 -5.1256 -4.9844 -0.5612 -0.4901 3.8891 4.9781
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250 0.5052 0.1375 ( 662 PWs) bands (ev):
|
|
|
|
-7.9872 -7.9091 -5.0855 -5.0270 -0.5407 -0.5112 4.2027 4.6566
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250-0.3608 0.0458 ( 661 PWs) bands (ev):
|
|
|
|
-10.1062 -9.8766 -2.1486 -2.0510 0.2055 0.2646 1.6280 3.1784
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250-0.3608 0.1375 ( 657 PWs) bands (ev):
|
|
|
|
-10.0413 -9.9463 -2.1204 -2.0800 0.2224 0.2469 2.0361 2.6716
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250-0.0722 0.0458 ( 639 PWs) bands (ev):
|
|
|
|
-12.2844 -12.0183 -0.9797 0.7816 2.4044 2.4467 3.0746 3.1054
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.1250-0.0722 0.1375 ( 635 PWs) bands (ev):
|
|
|
|
-12.2096 -12.0995 -0.5235 0.1960 2.4168 2.4343 3.0836 3.0964
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.3750 0.6495 0.0458 ( 647 PWs) bands (ev):
|
|
|
|
-6.4481 -6.3558 -5.5085 -5.4725 -2.8146 -2.7265 5.5958 6.1966
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.3750 0.6495 0.1375 ( 662 PWs) bands (ev):
|
|
|
|
-6.4181 -6.3796 -5.5031 -5.4878 -2.7903 -2.7537 5.8251 6.0866
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.3750-0.2165 0.0458 ( 658 PWs) bands (ev):
|
|
|
|
-9.3983 -9.1811 -3.8297 -3.6966 0.7878 0.8353 2.4191 3.8391
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
k = 0.3750-0.2165 0.1375 ( 656 PWs) bands (ev):
|
|
|
|
-9.3367 -9.2468 -3.7916 -3.7365 0.8016 0.8213 2.7998 3.3839
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
highest occupied level (ev): 6.1966
|
|
|
|
! total energy = -45.88135727 Ry
|
|
Harris-Foulkes estimate = -45.88135728 Ry
|
|
estimated scf accuracy < 4.4E-09 Ry
|
|
|
|
The total energy is the sum of the following terms:
|
|
|
|
one-electron contribution = -11.28163530 Ry
|
|
hartree contribution = 13.64793017 Ry
|
|
xc contribution = -14.41142272 Ry
|
|
ewald contribution = -33.83622943 Ry
|
|
|
|
convergence has been achieved in 7 iterations
|
|
|
|
Forces acting on atoms (cartesian axes, Ry/au):
|
|
|
|
atom 1 type 1 force = 0.00000000 -0.00000000 0.00000000
|
|
atom 2 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 3 type 1 force = 0.00000000 -0.00000000 0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 3 type 1 force = -0.00000000 -0.00000000 0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = 0.00000000 -0.00000000 0.00000000
|
|
atom 3 type 1 force = 0.00000000 0.00000000 -0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = 0.00000000 -0.00000000 0.00000000
|
|
atom 3 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 3 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 3 type 1 force = 0.00000000 0.00000000 0.00000000
|
|
atom 4 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
|
|
atom 2 type 1 force = -0.00000000 -0.00000000 -0.00000000
|
|
atom 3 type 1 force = -0.00000000 0.00000000 0.00000000
|
|
atom 4 type 1 force = 0.00000000 -0.00000000 0.00000000
|
|
|
|
Total force = 0.000000 Total SCF correction = 0.000000
|
|
|
|
|
|
Computing stress (Cartesian axis) and pressure
|
|
|
|
total stress (Ry/bohr**3) (kbar) P= 16.79
|
|
0.00010287 0.00000000 -0.00000000 15.13 0.00 -0.00
|
|
0.00000000 0.00010287 -0.00000000 0.00 15.13 -0.00
|
|
-0.00000000 -0.00000000 0.00013673 -0.00 -0.00 20.11
|
|
|
|
kinetic stress (kbar) 10653.05 -0.00 -0.00
|
|
-0.00 10653.05 -0.00
|
|
-0.00 -0.00 12340.52
|
|
|
|
local stress (kbar) -18417.57 -0.00 -0.00
|
|
-0.00 -18417.57 -0.00
|
|
-0.00 -0.00 25473.58
|
|
|
|
nonloc. stress (kbar) -464.98 0.00 0.00
|
|
0.00 -464.98 0.00
|
|
0.00 0.00 -462.29
|
|
|
|
hartree stress (kbar) 7392.66 -0.00 0.00
|
|
-0.00 7392.66 -0.00
|
|
0.00 -0.00 -6334.73
|
|
|
|
exc-cor stress (kbar) -2730.35 0.00 -0.00
|
|
0.00 -2730.35 0.00
|
|
-0.00 0.00 -2972.99
|
|
|
|
corecor stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
ewald stress (kbar) 3575.09 0.00 -0.00
|
|
0.00 3575.09 0.00
|
|
-0.00 0.00 -28101.07
|
|
|
|
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
|
|
|
|
DFT-D3 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) 7.23 -0.00 0.00
|
|
-0.00 7.23 0.00
|
|
0.00 0.00 77.09
|
|
|
|
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 /benchmarks/tempdir/graphite.save/
|
|
|
|
init_run : 25.06s CPU 25.09s WALL ( 1 calls)
|
|
electrons : 0.43s CPU 0.50s WALL ( 1 calls)
|
|
forces : 0.01s CPU 0.01s WALL ( 1 calls)
|
|
stress : 0.05s CPU 0.05s WALL ( 1 calls)
|
|
|
|
Called by init_run:
|
|
wfcinit : 0.02s CPU 0.03s WALL ( 1 calls)
|
|
wfcinit:atom : 0.00s CPU 0.00s WALL ( 12 calls)
|
|
wfcinit:wfcr : 0.02s CPU 0.03s WALL ( 12 calls)
|
|
potinit : 25.02s CPU 25.04s WALL ( 1 calls)
|
|
hinit0 : 0.01s CPU 0.01s WALL ( 1 calls)
|
|
|
|
Called by electrons:
|
|
c_bands : 0.27s CPU 0.31s WALL ( 7 calls)
|
|
sum_band : 0.05s CPU 0.06s WALL ( 7 calls)
|
|
v_of_rho : 25.11s CPU 25.15s WALL ( 8 calls)
|
|
v_h : 0.00s CPU 0.00s WALL ( 8 calls)
|
|
v_xc : 25.11s CPU 25.15s WALL ( 8 calls)
|
|
newd : 0.01s CPU 0.02s WALL ( 8 calls)
|
|
mix_rho : 0.00s CPU 0.00s WALL ( 7 calls)
|
|
vdW_kernel : 25.00s CPU 25.02s WALL ( 1 calls)
|
|
|
|
Called by c_bands:
|
|
init_us_2 : 0.01s CPU 0.01s WALL ( 204 calls)
|
|
cegterg : 0.25s CPU 0.29s WALL ( 84 calls)
|
|
|
|
Called by sum_band:
|
|
sum_band:bec : 0.00s CPU 0.00s WALL ( 84 calls)
|
|
addusdens : 0.01s CPU 0.01s WALL ( 7 calls)
|
|
|
|
Called by *egterg:
|
|
h_psi : 0.20s CPU 0.23s WALL ( 298 calls)
|
|
s_psi : 0.01s CPU 0.01s WALL ( 298 calls)
|
|
g_psi : 0.00s CPU 0.00s WALL ( 202 calls)
|
|
cdiaghg : 0.03s CPU 0.04s WALL ( 286 calls)
|
|
cegterg:over : 0.01s CPU 0.01s WALL ( 202 calls)
|
|
cegterg:upda : 0.00s CPU 0.01s WALL ( 202 calls)
|
|
cegterg:last : 0.00s CPU 0.00s WALL ( 84 calls)
|
|
|
|
Called by h_psi:
|
|
h_psi:calbec : 0.01s CPU 0.02s WALL ( 298 calls)
|
|
vloc_psi : 0.18s CPU 0.20s WALL ( 298 calls)
|
|
add_vuspsi : 0.01s CPU 0.01s WALL ( 298 calls)
|
|
|
|
General routines
|
|
calbec : 0.02s CPU 0.02s WALL ( 442 calls)
|
|
fft : 0.07s CPU 0.07s WALL ( 554 calls)
|
|
ffts : 0.00s CPU 0.00s WALL ( 15 calls)
|
|
fftw : 0.18s CPU 0.21s WALL ( 4908 calls)
|
|
interpolate : 0.00s CPU 0.00s WALL ( 8 calls)
|
|
|
|
Parallel routines
|
|
fft_scatt_xy : 0.04s CPU 0.04s WALL ( 5477 calls)
|
|
fft_scatt_yz : 0.09s CPU 0.10s WALL ( 5477 calls)
|
|
|
|
PWSCF : 25.69s CPU 25.90s WALL
|
|
|
|
|
|
This run was terminated on: 12:52:17 23Sep2019
|
|
|
|
=------------------------------------------------------------------------------=
|
|
JOB DONE.
|
|
=------------------------------------------------------------------------------=
|