Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 17: 4:51 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 Parallel version (MPI), running on 36 processors R & G space division: proc/nbgrp/npool/nimage = 36 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: one sub-group per k-point group (pool) will be used scalapack distributed-memory algorithm (size of sub-group: 4* 4 procs) Parallelization info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Min 38 22 6 1066 450 78 Max 39 23 7 1069 476 83 Sum 1369 793 241 38401 16663 2897 bravais-lattice index = 14 lattice parameter (alat) = 7.3363 a.u. unit-cell volume = 394.8483 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 12.00 number of Kohn-Sham states= 20 kinetic-energy cutoff = 46.0000 Ry charge density cutoff = 321.0000 Ry convergence threshold = 1.0E-08 mixing beta = 0.7000 number of iterations used = 8 plain mixing Exchange-correlation = PBE ( 1 4 3 4 0 0) Non magnetic calculation with spin-orbit celldm(1)= 7.336294 celldm(2)= 1.000000 celldm(3)= 1.000000 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.000000 1.000000 0.000000 ) a(3) = ( 0.000000 0.000000 1.000000 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 1.000000 -0.000000 -0.000000 ) b(2) = ( 0.000000 1.000000 -0.000000 ) b(3) = ( 0.000000 0.000000 1.000000 ) PseudoPot. # 1 for Cl read from file: /users/gautes/Pseudo/Cl.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 689bc0dc0caaad1d61892f37b0005eca Pseudo is Ultrasoft + core correction, Zval = 7.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1157 points, 10 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 l(7) = 2 l(8) = 2 l(9) = 2 l(10) = 2 Q(r) pseudized with 0 coefficients PseudoPot. # 2 for N read from file: /users/gautes/Pseudo/N.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: e87e56825df8daeb07642eb4d268bcfd Pseudo is Ultrasoft + core correction, Zval = 5.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1085 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 1 l(6) = 1 Q(r) pseudized with 0 coefficients atomic species valence mass pseudopotential Cl 7.00 35.45300 Cl( 1.00) N 5.00 14.00670 N( 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) = ( -1 0 0 ) ( 0 -1 0 ) ( 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 ) ( 0 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 ) ( 0 -1 0 ) ( 0 0 -1 ) 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) = ( 0 -1 0 ) ( -1 0 0 ) ( 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 ) ( 1 0 0 ) ( 0 0 1 ) 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) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 1 ) 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) = ( 0 0 1 ) ( 0 -1 0 ) ( 1 0 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) = ( 0 0 -1 ) ( 0 -1 0 ) ( -1 0 0 ) 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 0 -1 ) ( 0 1 0 ) ( 1 0 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 0 1 ) ( 0 1 0 ) ( -1 0 0 ) 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 0 0 ) ( 0 0 1 ) ( 0 1 0 ) 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) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 -1 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) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) 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 0 0 ) ( 0 0 1 ) ( 0 -1 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) = ( 0 0 1 ) ( 1 0 0 ) ( 0 1 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 0 -1 ) ( -1 0 0 ) ( 0 1 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) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 -1 0 ) 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 0 1 ) ( -1 0 0 ) ( 0 -1 0 ) 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 1 0 ) ( 0 0 1 ) ( 1 0 0 ) 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) = ( 0 -1 0 ) ( 0 0 -1 ) ( 1 0 0 ) 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 -1 0 ) ( 0 0 1 ) ( -1 0 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) = ( 0 1 0 ) ( 0 0 -1 ) ( -1 0 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) = ( 1 0 0 ) ( 0 1 0 ) ( 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 ) ( 0 -1 0 ) ( 0 0 1 ) 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) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) 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) = ( 0 1 0 ) ( 1 0 0 ) ( 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 ) ( -1 0 0 ) ( 0 0 -1 ) 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) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 -1 ) 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) = ( 0 0 -1 ) ( 0 1 0 ) ( -1 0 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) = ( 0 0 1 ) ( 0 1 0 ) ( 1 0 0 ) 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 0 1 ) ( 0 -1 0 ) ( -1 0 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 0 -1 ) ( 0 -1 0 ) ( 1 0 0 ) 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 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) 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) = ( 1 0 0 ) ( 0 0 1 ) ( 0 1 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) = ( -1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) 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 0 0 ) ( 0 0 -1 ) ( 0 1 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) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 -1 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 0 1 ) ( 1 0 0 ) ( 0 -1 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) = ( 0 0 1 ) ( -1 0 0 ) ( 0 1 0 ) 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 0 -1 ) ( 1 0 0 ) ( 0 1 0 ) 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 -1 0 ) ( 0 0 -1 ) ( -1 0 0 ) 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) = ( 0 1 0 ) ( 0 0 1 ) ( -1 0 0 ) 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 1 0 ) ( 0 0 -1 ) ( 1 0 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) = ( 0 -1 0 ) ( 0 0 1 ) ( 1 0 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 number of k points= 20 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.1666667), wk = 0.0277778 k( 3) = ( 0.0000000 0.0000000 0.3333333), wk = 0.0277778 k( 4) = ( 0.0000000 0.0000000 -0.5000000), wk = 0.0138889 k( 5) = ( 0.0000000 0.1666667 0.1666667), wk = 0.0555556 k( 6) = ( 0.0000000 0.1666667 0.3333333), wk = 0.1111111 k( 7) = ( 0.0000000 0.1666667 -0.5000000), wk = 0.0555556 k( 8) = ( 0.0000000 0.3333333 0.3333333), wk = 0.0555556 k( 9) = ( 0.0000000 0.3333333 -0.5000000), wk = 0.0555556 k( 10) = ( 0.0000000 -0.5000000 -0.5000000), wk = 0.0138889 k( 11) = ( 0.1666667 0.1666667 0.1666667), wk = 0.0370370 k( 12) = ( 0.1666667 0.1666667 0.3333333), wk = 0.1111111 k( 13) = ( 0.1666667 0.1666667 -0.5000000), wk = 0.0555556 k( 14) = ( 0.1666667 0.3333333 0.3333333), wk = 0.1111111 k( 15) = ( 0.1666667 0.3333333 -0.5000000), wk = 0.1111111 k( 16) = ( 0.1666667 -0.5000000 -0.5000000), wk = 0.0277778 k( 17) = ( 0.3333333 0.3333333 0.3333333), wk = 0.0370370 k( 18) = ( 0.3333333 0.3333333 -0.5000000), wk = 0.0555556 k( 19) = ( 0.3333333 -0.5000000 -0.5000000), wk = 0.0277778 k( 20) = ( -0.5000000 -0.5000000 -0.5000000), wk = 0.0046296 Dense grid: 38401 G-vectors FFT dimensions: ( 45, 45, 45) Smooth grid: 16663 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.04 Mb ( 132, 20) NL pseudopotentials 0.05 Mb ( 66, 48) Each V/rho on FFT grid 0.06 Mb ( 4050) Each G-vector array 0.01 Mb ( 1068) G-vector shells 0.00 Mb ( 253) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.16 Mb ( 132, 80) Each subspace H/S matrix 0.01 Mb ( 20, 20) Each matrix 0.03 Mb ( 48, 2, 20) Arrays for rho mixing 0.49 Mb ( 4050, 8) Check: negative/imaginary core charge= -0.000004 0.000000 Initial potential from superposition of free atoms starting charge 11.99999, renormalised to 12.00000 Starting wfc are 16 randomized atomic wfcs + 4 random wfc total cpu time spent up to now is 1.6 secs per-process dynamical memory: 18.9 Mb Self-consistent Calculation iteration # 1 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 6.1 total cpu time spent up to now is 2.3 secs total energy = -54.31976627 Ry Harris-Foulkes estimate = -54.49423601 Ry estimated scf accuracy < 0.21188180 Ry iteration # 2 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.77E-03, avg # of iterations = 2.3 total cpu time spent up to now is 2.7 secs total energy = -54.15761168 Ry Harris-Foulkes estimate = -55.57122543 Ry estimated scf accuracy < 8.60288611 Ry iteration # 3 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.77E-03, avg # of iterations = 1.9 total cpu time spent up to now is 3.1 secs total energy = -54.47466394 Ry Harris-Foulkes estimate = -54.47474189 Ry estimated scf accuracy < 0.00418140 Ry iteration # 4 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.48E-05, avg # of iterations = 3.6 total cpu time spent up to now is 3.6 secs total energy = -54.47479794 Ry Harris-Foulkes estimate = -54.47488436 Ry estimated scf accuracy < 0.00316848 Ry iteration # 5 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.64E-05, avg # of iterations = 1.0 total cpu time spent up to now is 3.9 secs total energy = -54.47406813 Ry Harris-Foulkes estimate = -54.47486311 Ry estimated scf accuracy < 0.00561070 Ry iteration # 6 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.64E-05, avg # of iterations = 1.0 total cpu time spent up to now is 4.2 secs total energy = -54.47438007 Ry Harris-Foulkes estimate = -54.47447878 Ry estimated scf accuracy < 0.00032942 Ry iteration # 7 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.75E-06, avg # of iterations = 1.9 total cpu time spent up to now is 4.6 secs total energy = -54.47435210 Ry Harris-Foulkes estimate = -54.47451560 Ry estimated scf accuracy < 0.00164058 Ry iteration # 8 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.75E-06, avg # of iterations = 1.0 total cpu time spent up to now is 5.0 secs total energy = -54.47443693 Ry Harris-Foulkes estimate = -54.47444570 Ry estimated scf accuracy < 0.00003721 Ry iteration # 9 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.10E-07, avg # of iterations = 1.0 total cpu time spent up to now is 5.3 secs total energy = -54.47444078 Ry Harris-Foulkes estimate = -54.47444262 Ry estimated scf accuracy < 0.00003895 Ry iteration # 10 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.10E-07, avg # of iterations = 1.0 total cpu time spent up to now is 5.6 secs total energy = -54.47444015 Ry Harris-Foulkes estimate = -54.47444110 Ry estimated scf accuracy < 0.00001564 Ry iteration # 11 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.30E-07, avg # of iterations = 1.0 total cpu time spent up to now is 6.0 secs total energy = -54.47444020 Ry Harris-Foulkes estimate = -54.47444042 Ry estimated scf accuracy < 0.00000229 Ry iteration # 12 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.91E-08, avg # of iterations = 1.0 total cpu time spent up to now is 6.3 secs total energy = -54.47444029 Ry Harris-Foulkes estimate = -54.47444029 Ry estimated scf accuracy < 0.00000001 Ry iteration # 13 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.96E-11, avg # of iterations = 2.1 total cpu time spent up to now is 6.8 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2103 PWs) bands (ev): -15.3266 -15.3266 -13.6706 -13.6706 -3.0699 -3.0699 -2.9759 -2.9759 -2.9759 -2.9759 -1.6093 -1.6093 -1.5486 -1.5486 -1.5486 -1.5486 5.8113 5.8114 8.1111 8.1111 occupation numbers 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.0000 0.0000 0.1667 ( 2068 PWs) bands (ev): -15.2761 -15.2761 -13.6982 -13.6982 -3.1356 -3.1356 -2.9429 -2.9429 -2.8896 -2.8896 -1.8235 -1.8235 -1.6615 -1.6615 -1.6272 -1.6272 6.2960 6.2960 8.5443 8.5459 occupation numbers 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.0000 0.0000 0.3333 ( 2056 PWs) bands (ev): -15.1659 -15.1659 -13.7628 -13.7628 -3.3359 -3.3359 -2.7366 -2.7366 -2.6671 -2.6671 -2.2759 -2.2759 -1.8639 -1.8639 -1.8351 -1.8351 7.3433 7.3433 8.9286 8.9286 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.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000-0.5000 ( 2088 PWs) bands (ev): -15.1050 -15.1050 -13.8012 -13.8012 -3.4247 -3.4247 -2.5751 -2.5751 -2.4901 -2.4901 -2.4895 -2.4895 -2.0195 -2.0195 -2.0059 -2.0059 8.0360 8.0361 9.3777 9.3779 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.9144 0.9144 0.7970 0.7970 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 2063 PWs) bands (ev): -15.2308 -15.2308 -13.7159 -13.7159 -3.2817 -3.2817 -2.8894 -2.8894 -2.8185 -2.8185 -2.0948 -2.0948 -1.7097 -1.7097 -1.5805 -1.5805 6.7121 6.7121 8.6153 8.6153 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9996 0.9996 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 2064 PWs) bands (ev): -15.1336 -15.1336 -13.7582 -13.7582 -3.4970 -3.4970 -3.0046 -3.0046 -2.6484 -2.6484 -2.1960 -2.1960 -1.8709 -1.8709 -1.6972 -1.6972 7.6121 7.6121 9.0366 9.0367 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.0002 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667-0.5000 ( 2080 PWs) bands (ev): -15.0809 -15.0809 -13.7836 -13.7836 -3.5792 -3.5792 -3.1648 -3.1648 -2.5093 -2.5093 -2.0129 -2.0129 -1.9955 -1.9955 -1.9020 -1.9020 8.2755 8.2755 8.9914 8.9914 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.8683 0.8683 0.6465 0.6465 0.0019 0.0019 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 2073 PWs) bands (ev): -15.0666 -15.0666 -13.7506 -13.7506 -3.7561 -3.7561 -3.5004 -3.5004 -2.5333 -2.5333 -1.9529 -1.9529 -1.9246 -1.9246 -1.6142 -1.6142 8.0167 8.0167 8.6916 8.6916 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0740 0.0740 0.0099 0.0099 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 2096 PWs) bands (ev): -15.0317 -15.0317 -13.7478 -13.7478 -3.8511 -3.8511 -3.7338 -3.7338 -2.4618 -2.4618 -1.9644 -1.9644 -1.7298 -1.7298 -1.7106 -1.7106 8.1700 8.1700 8.4612 8.4612 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.1571 0.1571 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 2096 PWs) bands (ev): -15.0067 -15.0067 -13.7296 -13.7296 -3.9908 -3.9908 -3.9430 -3.9430 -2.4373 -2.4373 -1.9479 -1.9479 -1.6405 -1.6405 -1.6217 -1.6217 7.8117 7.8117 8.4804 8.4804 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0525 0.0525 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 2064 PWs) bands (ev): -15.1901 -15.1901 -13.7261 -13.7261 -3.2152 -3.2152 -3.1480 -3.1480 -2.7996 -2.7996 -2.1285 -2.1285 -1.6970 -1.6970 -1.6594 -1.6594 7.1117 7.1117 9.0397 9.0428 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 k = 0.1667 0.1667 0.3333 ( 2078 PWs) bands (ev): -15.1036 -15.1036 -13.7513 -13.7513 -3.5052 -3.5052 -3.0996 -3.0996 -2.9205 -2.9205 -1.9998 -1.9998 -1.9266 -1.9266 -1.6973 -1.6973 7.9543 7.9543 9.4750 9.4776 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.7156 0.7156 0.0114 0.0114 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 2084 PWs) bands (ev): -15.0574 -15.0574 -13.7668 -13.7668 -3.6339 -3.6339 -3.2624 -3.2624 -2.7263 -2.7263 -2.1229 -2.1229 -1.8559 -1.8559 -1.7797 -1.7797 8.7025 8.7025 8.9533 8.9534 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.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333 0.3333 ( 2090 PWs) bands (ev): -15.0416 -15.0416 -13.7385 -13.7385 -3.6424 -3.6424 -3.4896 -3.4896 -2.8616 -2.8616 -2.0275 -2.0275 -1.8201 -1.8201 -1.7429 -1.7429 8.4077 8.4077 9.1871 9.1871 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9505 0.9505 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 2088 PWs) bands (ev): -15.0096 -15.0096 -13.7327 -13.7327 -3.7721 -3.7721 -3.6400 -3.6400 -2.7219 -2.7219 -2.1345 -2.1345 -1.7759 -1.7759 -1.7428 -1.7428 8.4366 8.4366 9.1400 9.1400 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 k = 0.1667-0.5000-0.5000 ( 2080 PWs) bands (ev): -14.9853 -14.9853 -13.7155 -13.7155 -3.8684 -3.8684 -3.8178 -3.8178 -2.6580 -2.6580 -2.0781 -2.0781 -1.7863 -1.7863 -1.7696 -1.7696 8.1386 8.1386 9.2184 9.2185 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9987 0.9987 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 2100 PWs) bands (ev): -14.9907 -14.9907 -13.7143 -13.7143 -3.5270 -3.5270 -3.3540 -3.3540 -3.2984 -3.2984 -2.0913 -2.0913 -2.0626 -2.0626 -1.8195 -1.8195 8.9026 8.9026 9.6627 9.6652 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9995 0.9995 0.9961 0.9961 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 2078 PWs) bands (ev): -14.9646 -14.9646 -13.7023 -13.7023 -3.6076 -3.6076 -3.4713 -3.4713 -3.0811 -3.0811 -2.3374 -2.3374 -1.9827 -1.9827 -1.9415 -1.9415 8.8827 8.8827 9.4213 9.4213 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.4163 0.4163 0.0335 0.0335 0.0000 0.0000 0.0000 0.0000 k = 0.3333-0.5000-0.5000 ( 2096 PWs) bands (ev): -14.9418 -14.9418 -13.6869 -13.6869 -3.5506 -3.5506 -3.4865 -3.4865 -3.0591 -3.0591 -2.3249 -2.3249 -2.1548 -2.1548 -2.1479 -2.1479 8.8341 8.8342 9.4784 9.4799 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 k =-0.5000-0.5000-0.5000 ( 2008 PWs) bands (ev): -14.9197 -14.9197 -13.6724 -13.6724 -3.3292 -3.3292 -3.2112 -3.2112 -3.2112 -3.2112 -2.4541 -2.4541 -2.4353 -2.4353 -2.4353 -2.4353 9.2112 9.2156 9.2165 9.2165 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is -1.9873 ev ! total energy = -54.47444029 Ry Harris-Foulkes estimate = -54.47444029 Ry estimated scf accuracy < 2.1E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -12.57404673 Ry hartree contribution = 10.02706721 Ry xc contribution = -15.65583670 Ry ewald contribution = -36.27110948 Ry smearing contrib. (-TS) = -0.00051459 Ry convergence has been achieved in 13 iterations Writing output data file NCl.save init_run : 0.22s CPU 0.30s WALL ( 1 calls) electrons : 4.78s CPU 5.30s WALL ( 1 calls) Called by init_run: wfcinit : 0.10s CPU 0.12s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 3.46s CPU 3.82s WALL ( 13 calls) sum_band : 0.93s CPU 0.96s WALL ( 13 calls) v_of_rho : 0.05s CPU 0.04s WALL ( 14 calls) v_h : 0.01s CPU 0.00s WALL ( 14 calls) v_xc : 0.04s CPU 0.04s WALL ( 14 calls) newd : 0.30s CPU 0.32s WALL ( 14 calls) mix_rho : 0.02s CPU 0.03s WALL ( 13 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.02s WALL ( 540 calls) cegterg : 3.26s CPU 3.34s WALL ( 260 calls) Called by sum_band: sum_band:bec : 0.26s CPU 0.27s WALL ( 260 calls) addusdens : 0.19s CPU 0.19s WALL ( 13 calls) Called by *egterg: h_psi : 1.88s CPU 1.89s WALL ( 779 calls) s_psi : 0.08s CPU 0.07s WALL ( 779 calls) g_psi : 0.00s CPU 0.01s WALL ( 499 calls) cdiaghg : 1.16s CPU 1.19s WALL ( 759 calls) cegterg:over : 0.08s CPU 0.08s WALL ( 499 calls) cegterg:upda : 0.05s CPU 0.07s WALL ( 499 calls) cegterg:last : 0.01s CPU 0.03s WALL ( 260 calls) cdiaghg:chol : 0.05s CPU 0.08s WALL ( 759 calls) cdiaghg:inve : 0.02s CPU 0.02s WALL ( 759 calls) cdiaghg:para : 0.12s CPU 0.09s WALL ( 1518 calls) Called by h_psi: h_psi:vloc : 1.69s CPU 1.71s WALL ( 779 calls) h_psi:vnl : 0.18s CPU 0.18s WALL ( 779 calls) add_vuspsi : 0.09s CPU 0.09s WALL ( 779 calls) General routines calbec : 0.13s CPU 0.12s WALL ( 1039 calls) fft : 0.08s CPU 0.08s WALL ( 418 calls) ffts : 0.00s CPU 0.01s WALL ( 108 calls) fftw : 1.88s CPU 1.89s WALL ( 58432 calls) interpolate : 0.02s CPU 0.03s WALL ( 108 calls) Parallel routines fft_scatter : 1.25s CPU 1.24s WALL ( 58958 calls) PWSCF : 6.47s CPU 8.83s WALL This run was terminated on: 17: 5: 0 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=