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:25
|
|
|
|
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.6412 a.u.
|
|
unit-cell volume = 236.0493 (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 = 30.0000 Ry
|
|
charge density cutoff = 180.0000 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.641170 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.11987523 Ry
|
|
-----------------------------------------------
|
|
|
|
Starting wfcs are 16 randomized atomic wfcs
|
|
|
|
total cpu time spent up to now is 25.4 secs
|
|
|
|
Self-consistent Calculation
|
|
|
|
iteration # 1 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.00E-02, avg # of iterations = 2.0
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15290138 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.5 secs
|
|
|
|
total energy = -45.81464016 Ry
|
|
Harris-Foulkes estimate = -46.06086379 Ry
|
|
estimated scf accuracy < 0.43954209 Ry
|
|
|
|
iteration # 2 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 2.75E-03, avg # of iterations = 2.0
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15250949 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.6 secs
|
|
|
|
total energy = -45.88020136 Ry
|
|
Harris-Foulkes estimate = -45.87895069 Ry
|
|
estimated scf accuracy < 0.00557486 Ry
|
|
|
|
iteration # 3 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 3.48E-05, avg # of iterations = 2.2
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15272467 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.7 secs
|
|
|
|
total energy = -45.88103678 Ry
|
|
Harris-Foulkes estimate = -45.88083774 Ry
|
|
estimated scf accuracy < 0.00041630 Ry
|
|
|
|
iteration # 4 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 2.60E-06, avg # of iterations = 1.9
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15279605 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.7 secs
|
|
|
|
total energy = -45.88108659 Ry
|
|
Harris-Foulkes estimate = -45.88108423 Ry
|
|
estimated scf accuracy < 0.00000290 Ry
|
|
|
|
iteration # 5 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 1.81E-08, avg # of iterations = 3.1
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15279534 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.8 secs
|
|
|
|
total energy = -45.88108911 Ry
|
|
Harris-Foulkes estimate = -45.88108910 Ry
|
|
estimated scf accuracy < 0.00000012 Ry
|
|
|
|
iteration # 6 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 7.42E-10, avg # of iterations = 3.1
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15279564 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 25.9 secs
|
|
|
|
total energy = -45.88108913 Ry
|
|
Harris-Foulkes estimate = -45.88108915 Ry
|
|
estimated scf accuracy < 0.00000005 Ry
|
|
|
|
iteration # 7 ecut= 30.00 Ry beta= 0.70
|
|
Davidson diagonalization with overlap
|
|
ethr = 3.23E-10, avg # of iterations = 2.5
|
|
|
|
|
|
-----------------------------------------------
|
|
Non-local corr. energy = 0.15279273 Ry
|
|
-----------------------------------------------
|
|
|
|
|
|
total cpu time spent up to now is 26.0 secs
|
|
|
|
End of self-consistent calculation
|
|
|
|
k = 0.1250 0.2165 0.0458 ( 646 PWs) bands (ev):
|
|
|
|
-11.5264 -11.2692 -0.0598 0.6944 0.7336 1.6561 1.7746 1.8020
|
|
|
|
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.4540 -11.3476 0.3857 0.7057 0.7231 1.0882 1.7820 1.7934
|
|
|
|
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.0012 -7.8097 -5.0824 -4.9390 -0.5105 -0.4386 3.9429 5.0425
|
|
|
|
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.9465 -7.8673 -5.0417 -4.9823 -0.4898 -0.4600 4.2595 4.7178
|
|
|
|
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.0763 -9.8438 -2.0915 -1.9926 0.2592 0.3191 1.6680 3.2327
|
|
|
|
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.0106 -9.9143 -2.0629 -2.0219 0.2764 0.3012 2.0797 2.7210
|
|
|
|
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.2634 -11.9936 -0.9538 0.8227 2.4753 2.5180 3.1443 3.1754
|
|
|
|
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.1876 -12.0760 -0.4938 0.2319 2.4878 2.5055 3.1534 3.1663
|
|
|
|
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.3966 -6.3031 -5.4650 -5.4278 -2.7781 -2.6882 5.6602 6.2663
|
|
|
|
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.3661 -6.3271 -5.4593 -5.4436 -2.7532 -2.7159 5.8919 6.1558
|
|
|
|
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.3653 -9.1451 -3.7809 -3.6460 0.8438 0.8919 2.4638 3.8973
|
|
|
|
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.3029 -9.2118 -3.7422 -3.6864 0.8579 0.8777 2.8479 3.4375
|
|
|
|
occupation numbers
|
|
1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
|
|
|
|
highest occupied level (ev): 6.2663
|
|
|
|
! total energy = -45.88108914 Ry
|
|
Harris-Foulkes estimate = -45.88108914 Ry
|
|
estimated scf accuracy < 4.5E-09 Ry
|
|
|
|
The total energy is the sum of the following terms:
|
|
|
|
one-electron contribution = -11.17965974 Ry
|
|
hartree contribution = 13.63785457 Ry
|
|
xc contribution = -14.43015002 Ry
|
|
ewald contribution = -33.90913395 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= 34.81
|
|
0.00028286 -0.00000000 0.00000000 41.61 -0.00 0.00
|
|
-0.00000000 0.00028286 0.00000000 -0.00 41.61 0.00
|
|
0.00000000 0.00000000 0.00014410 0.00 0.00 21.20
|
|
|
|
kinetic stress (kbar) 10748.59 -0.00 0.00
|
|
-0.00 10748.59 0.00
|
|
0.00 0.00 12437.07
|
|
|
|
local stress (kbar) -18533.87 -0.00 -0.00
|
|
-0.00 -18533.87 -0.00
|
|
-0.00 -0.00 25693.82
|
|
|
|
nonloc. stress (kbar) -473.26 -0.00 0.00
|
|
-0.00 -473.26 0.00
|
|
0.00 0.00 -469.23
|
|
|
|
hartree stress (kbar) 7438.50 0.00 -0.00
|
|
0.00 7438.50 -0.00
|
|
-0.00 -0.00 -6377.94
|
|
|
|
exc-cor stress (kbar) -2751.57 0.00 0.00
|
|
0.00 -2751.57 0.00
|
|
0.00 0.00 -2996.00
|
|
|
|
corecor stress (kbar) 0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
0.00 0.00 0.00
|
|
|
|
ewald stress (kbar) 3606.00 0.00 -0.00
|
|
0.00 3606.00 -0.00
|
|
-0.00 -0.00 -28344.05
|
|
|
|
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.53
|
|
|
|
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.26s CPU 25.29s WALL ( 1 calls)
|
|
electrons : 0.49s CPU 0.55s 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.03s CPU 0.03s WALL ( 1 calls)
|
|
wfcinit:atom : 0.00s CPU 0.00s WALL ( 12 calls)
|
|
wfcinit:wfcr : 0.03s CPU 0.03s WALL ( 12 calls)
|
|
potinit : 25.21s CPU 25.24s WALL ( 1 calls)
|
|
hinit0 : 0.01s CPU 0.01s WALL ( 1 calls)
|
|
|
|
Called by electrons:
|
|
c_bands : 0.31s CPU 0.33s WALL ( 7 calls)
|
|
sum_band : 0.06s CPU 0.06s WALL ( 7 calls)
|
|
v_of_rho : 25.31s CPU 25.36s WALL ( 8 calls)
|
|
v_h : 0.00s CPU 0.00s WALL ( 8 calls)
|
|
v_xc : 25.31s CPU 25.35s 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.18s CPU 25.21s WALL ( 1 calls)
|
|
|
|
Called by c_bands:
|
|
init_us_2 : 0.01s CPU 0.01s WALL ( 204 calls)
|
|
cegterg : 0.29s CPU 0.31s 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.24s CPU 0.25s 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.04s 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.02s CPU 0.02s WALL ( 298 calls)
|
|
vloc_psi : 0.21s CPU 0.22s WALL ( 298 calls)
|
|
add_vuspsi : 0.01s CPU 0.01s WALL ( 298 calls)
|
|
|
|
General routines
|
|
calbec : 0.02s CPU 0.03s WALL ( 442 calls)
|
|
fft : 0.07s CPU 0.08s WALL ( 554 calls)
|
|
ffts : 0.00s CPU 0.00s WALL ( 15 calls)
|
|
fftw : 0.22s CPU 0.23s WALL ( 4904 calls)
|
|
interpolate : 0.00s CPU 0.00s WALL ( 8 calls)
|
|
|
|
Parallel routines
|
|
fft_scatt_xy : 0.04s CPU 0.04s WALL ( 5473 calls)
|
|
fft_scatt_yz : 0.12s CPU 0.13s WALL ( 5473 calls)
|
|
|
|
PWSCF : 25.93s CPU 26.14s WALL
|
|
|
|
|
|
This run was terminated on: 12:51:51 23Sep2019
|
|
|
|
=------------------------------------------------------------------------------=
|
|
JOB DONE.
|
|
=------------------------------------------------------------------------------=
|