Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 20:55:30 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 57 47 13 1953 1453 215 Max 58 48 14 1956 1469 218 Sum 2061 1701 481 70319 52635 7809 bravais-lattice index = 14 lattice parameter (alat) = 8.9290 a.u. unit-cell volume = 711.8721 (a.u.)^3 number of atoms/cell = 4 number of atomic types = 2 number of electrons = 32.00 number of Kohn-Sham states= 40 kinetic-energy cutoff = 67.0000 Ry charge density cutoff = 324.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)= 8.928955 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 Na read from file: /users/gautes/Pseudo/Na.rel-pbe-spnl-rrkjus_psl.1.0.0.UPF MD5 check sum: 13bcaa2858b974004ed12db349b14e76 Pseudo is Ultrasoft + core correction, Zval = 9.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1121 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 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 Na 9.00 22.98980 Na( 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= 10 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 Dense grid: 70319 G-vectors FFT dimensions: ( 54, 54, 54) Smooth grid: 52635 G-vectors FFT dimensions: ( 48, 48, 48) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.23 Mb ( 380, 40) NL pseudopotentials 0.16 Mb ( 190, 56) Each V/rho on FFT grid 0.09 Mb ( 5832) Each G-vector array 0.01 Mb ( 1953) G-vector shells 0.00 Mb ( 427) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.93 Mb ( 380, 160) Each subspace H/S matrix 0.02 Mb ( 40, 40) Each matrix 0.07 Mb ( 56, 2, 40) Arrays for rho mixing 0.71 Mb ( 5832, 8) Check: negative/imaginary core charge= -0.000003 0.000000 Initial potential from superposition of free atoms starting charge 31.98112, renormalised to 32.00000 Starting wfc are 38 randomized atomic wfcs + 2 random wfc total cpu time spent up to now is 1.1 secs per-process dynamical memory: 12.5 Mb Self-consistent Calculation iteration # 1 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 4.6 total cpu time spent up to now is 2.0 secs total energy = -311.13336266 Ry Harris-Foulkes estimate = -312.87305250 Ry estimated scf accuracy < 2.19668977 Ry iteration # 2 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.86E-03, avg # of iterations = 3.3 total cpu time spent up to now is 2.9 secs total energy = -311.12764058 Ry Harris-Foulkes estimate = -314.77283295 Ry estimated scf accuracy < 11.12632353 Ry iteration # 3 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.86E-03, avg # of iterations = 3.0 total cpu time spent up to now is 3.8 secs total energy = -312.38034408 Ry Harris-Foulkes estimate = -312.38757549 Ry estimated scf accuracy < 0.08642373 Ry iteration # 4 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.70E-04, avg # of iterations = 1.7 total cpu time spent up to now is 4.5 secs total energy = -312.35709794 Ry Harris-Foulkes estimate = -312.38225529 Ry estimated scf accuracy < 0.06686826 Ry iteration # 5 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.09E-04, avg # of iterations = 2.0 total cpu time spent up to now is 5.2 secs total energy = -312.36278672 Ry Harris-Foulkes estimate = -312.36277258 Ry estimated scf accuracy < 0.00070084 Ry iteration # 6 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.19E-06, avg # of iterations = 3.9 total cpu time spent up to now is 6.1 secs total energy = -312.36293468 Ry Harris-Foulkes estimate = -312.36296457 Ry estimated scf accuracy < 0.00091021 Ry iteration # 7 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.19E-06, avg # of iterations = 1.0 total cpu time spent up to now is 6.8 secs total energy = -312.36281823 Ry Harris-Foulkes estimate = -312.36294335 Ry estimated scf accuracy < 0.00042261 Ry iteration # 8 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.32E-06, avg # of iterations = 2.0 total cpu time spent up to now is 7.5 secs total energy = -312.36284634 Ry Harris-Foulkes estimate = -312.36284651 Ry estimated scf accuracy < 0.00000295 Ry iteration # 9 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.21E-09, avg # of iterations = 2.4 total cpu time spent up to now is 8.3 secs total energy = -312.36284624 Ry Harris-Foulkes estimate = -312.36284668 Ry estimated scf accuracy < 0.00000150 Ry iteration # 10 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.69E-09, avg # of iterations = 2.0 total cpu time spent up to now is 9.0 secs total energy = -312.36284638 Ry Harris-Foulkes estimate = -312.36284639 Ry estimated scf accuracy < 0.00000003 Ry iteration # 11 ecut= 67.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.96E-11, avg # of iterations = 2.8 total cpu time spent up to now is 9.8 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 6619 PWs) bands (ev): -49.8004 -49.8004 -49.7911 -49.7911 -49.7911 -49.7911 -21.9022 -21.9022 -21.9022 -21.9022 -21.8818 -21.8818 -21.7556 -21.7556 -21.7556 -21.7556 -21.7426 -21.7426 -21.5989 -21.5989 -21.5643 -21.5643 -21.5643 -21.5643 -9.7088 -9.7088 0.6135 0.6135 1.1924 1.1924 1.2019 1.2019 1.2019 1.2019 6.3948 6.3948 6.3948 6.3948 7.5901 7.5919 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 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2000 ( 6615 PWs) bands (ev): -49.7997 -49.7997 -49.7919 -49.7919 -49.7911 -49.7911 -21.9019 -21.9019 -21.8997 -21.8997 -21.8828 -21.8828 -21.7553 -21.7553 -21.7541 -21.7541 -21.7436 -21.7436 -21.5979 -21.5979 -21.5701 -21.5701 -21.5654 -21.5654 -9.6744 -9.6744 0.2972 0.2972 1.2213 1.2213 1.2277 1.2277 1.6105 1.6105 5.2627 5.2627 6.6265 6.6265 7.5046 7.5081 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 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.9999 0.9999 0.9999 0.9999 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4000 ( 6572 PWs) bands (ev): -49.7978 -49.7978 -49.7937 -49.7937 -49.7911 -49.7911 -21.9013 -21.9013 -21.8941 -21.8941 -21.8859 -21.8859 -21.7548 -21.7548 -21.7504 -21.7504 -21.7462 -21.7462 -21.5932 -21.5932 -21.5827 -21.5827 -21.5671 -21.5671 -9.6177 -9.6177 -0.0476 -0.0476 1.2646 1.2646 1.2709 1.2709 2.4550 2.4550 3.7690 3.7690 7.1584 7.1584 7.3385 7.3419 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 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.9984 0.9984 0.9974 0.9974 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.2000 ( 6629 PWs) bands (ev): -49.7990 -49.7990 -49.7922 -49.7922 -49.7915 -49.7915 -21.9019 -21.9019 -21.8987 -21.8987 -21.8833 -21.8833 -21.7575 -21.7575 -21.7525 -21.7525 -21.7421 -21.7421 -21.6024 -21.6024 -21.5755 -21.5755 -21.5609 -21.5609 -9.6405 -9.6405 0.2839 0.2839 0.9003 0.9003 1.2475 1.2475 2.0695 2.0695 5.3667 5.3667 5.4223 5.4223 7.5450 7.5451 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 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.9995 0.9995 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.4000 ( 6597 PWs) bands (ev): -49.7972 -49.7972 -49.7937 -49.7937 -49.7917 -49.7917 -21.9024 -21.9024 -21.8954 -21.8954 -21.8846 -21.8846 -21.7591 -21.7591 -21.7514 -21.7514 -21.7397 -21.7397 -21.6094 -21.6094 -21.5795 -21.5795 -21.5579 -21.5579 -9.5848 -9.5848 0.0887 0.0887 0.8707 0.8707 1.2858 1.2858 2.7999 2.7999 3.9123 3.9123 5.6660 5.6660 7.5616 7.5616 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 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.9923 0.9923 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4000 0.4000 ( 6586 PWs) bands (ev): -49.7957 -49.7957 -49.7939 -49.7939 -49.7930 -49.7930 -21.9055 -21.9055 -21.8941 -21.8941 -21.8844 -21.8844 -21.7630 -21.7630 -21.7508 -21.7508 -21.7338 -21.7338 -21.6244 -21.6244 -21.5806 -21.5806 -21.5479 -21.5479 -9.5306 -9.5306 0.3334 0.3334 0.4630 0.4630 1.3157 1.3157 2.9334 2.9334 4.4033 4.4033 4.8603 4.8603 6.6403 6.6403 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 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.9347 0.9347 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.2000 ( 6601 PWs) bands (ev): -49.7983 -49.7983 -49.7922 -49.7922 -49.7922 -49.7922 -21.9022 -21.9022 -21.8991 -21.8991 -21.8827 -21.8827 -21.7623 -21.7623 -21.7478 -21.7478 -21.7406 -21.7406 -21.6128 -21.6128 -21.5709 -21.5709 -21.5592 -21.5592 -9.6072 -9.6072 0.3219 0.3219 0.9489 0.9489 0.9551 0.9551 2.3956 2.3956 5.3995 5.3995 5.4013 5.4013 6.0383 6.0383 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 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 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.4000 ( 6586 PWs) bands (ev): -49.7967 -49.7967 -49.7937 -49.7937 -49.7922 -49.7922 -21.9033 -21.9033 -21.8988 -21.8988 -21.8822 -21.8822 -21.7662 -21.7662 -21.7446 -21.7446 -21.7377 -21.7377 -21.6251 -21.6251 -21.5688 -21.5688 -21.5556 -21.5556 -9.5523 -9.5523 0.2133 0.2133 0.8466 0.8466 1.0403 1.0403 2.9707 2.9707 4.1538 4.1538 5.3428 5.3428 6.6231 6.6231 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 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 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.4000 0.4000 ( 6575 PWs) bands (ev): -49.7954 -49.7954 -49.7939 -49.7939 -49.7932 -49.7932 -21.9064 -21.9064 -21.9006 -21.9006 -21.8803 -21.8803 -21.7720 -21.7720 -21.7411 -21.7411 -21.7330 -21.7330 -21.6412 -21.6412 -21.5647 -21.5647 -21.5479 -21.5479 -9.4989 -9.4989 0.4309 0.4309 0.5457 0.5457 1.0500 1.0500 3.1137 3.1137 4.5092 4.5092 4.6979 4.6979 6.5284 6.5284 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 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 0.0000 0.0000 0.0000 0.0000 k = 0.4000 0.4000 0.4000 ( 6614 PWs) bands (ev): -49.7948 -49.7948 -49.7939 -49.7939 -49.7939 -49.7939 -21.9081 -21.9081 -21.9073 -21.9073 -21.8779 -21.8779 -21.7797 -21.7797 -21.7327 -21.7327 -21.7317 -21.7317 -21.6587 -21.6587 -21.5493 -21.5493 -21.5473 -21.5473 -9.4467 -9.4467 0.5999 0.5999 0.6905 0.6905 0.6975 0.6975 3.5454 3.5454 4.4046 4.4046 4.4073 4.4073 5.6280 5.6280 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 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 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 1.3519 ev ! total energy = -312.36284639 Ry Harris-Foulkes estimate = -312.36284639 Ry estimated scf accuracy < 3.3E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -297.35654929 Ry hartree contribution = 154.80870361 Ry xc contribution = -50.04695639 Ry ewald contribution = -119.76797682 Ry smearing contrib. (-TS) = -0.00006750 Ry convergence has been achieved in 11 iterations Writing output data file Na3N.save init_run : 0.33s CPU 0.36s WALL ( 1 calls) electrons : 8.62s CPU 8.77s WALL ( 1 calls) Called by init_run: wfcinit : 0.23s CPU 0.25s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 7.11s CPU 7.23s WALL ( 11 calls) sum_band : 1.29s CPU 1.30s WALL ( 11 calls) v_of_rho : 0.05s CPU 0.05s WALL ( 12 calls) v_h : 0.00s CPU 0.00s WALL ( 12 calls) v_xc : 0.05s CPU 0.05s WALL ( 12 calls) newd : 0.10s CPU 0.10s WALL ( 12 calls) mix_rho : 0.04s CPU 0.04s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.02s WALL ( 230 calls) cegterg : 6.92s CPU 7.02s WALL ( 110 calls) Called by sum_band: sum_band:bec : 0.04s CPU 0.03s WALL ( 110 calls) addusdens : 0.09s CPU 0.09s WALL ( 11 calls) Called by *egterg: h_psi : 5.01s CPU 5.10s WALL ( 407 calls) s_psi : 0.09s CPU 0.08s WALL ( 407 calls) g_psi : 0.01s CPU 0.01s WALL ( 287 calls) cdiaghg : 1.34s CPU 1.33s WALL ( 397 calls) cegterg:over : 0.17s CPU 0.18s WALL ( 287 calls) cegterg:upda : 0.18s CPU 0.17s WALL ( 287 calls) cegterg:last : 0.08s CPU 0.08s WALL ( 110 calls) cdiaghg:chol : 0.08s CPU 0.08s WALL ( 397 calls) cdiaghg:inve : 0.05s CPU 0.03s WALL ( 397 calls) cdiaghg:para : 0.07s CPU 0.09s WALL ( 794 calls) Called by h_psi: h_psi:vloc : 4.73s CPU 4.81s WALL ( 407 calls) h_psi:vnl : 0.24s CPU 0.27s WALL ( 407 calls) add_vuspsi : 0.11s CPU 0.10s WALL ( 407 calls) General routines calbec : 0.18s CPU 0.22s WALL ( 517 calls) fft : 0.12s CPU 0.11s WALL ( 356 calls) ffts : 0.02s CPU 0.02s WALL ( 92 calls) fftw : 5.39s CPU 5.47s WALL ( 50776 calls) interpolate : 0.05s CPU 0.05s WALL ( 92 calls) Parallel routines fft_scatter : 2.00s CPU 2.06s WALL ( 51224 calls) PWSCF : 10.10s CPU 10.99s WALL This run was terminated on: 20:55:41 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=