Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 16: 2: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); 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 file Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S renormalized 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 40 24 7 1158 560 93 Max 41 25 8 1163 573 96 Sum 1457 885 277 41755 20341 3407 bravais-lattice index = 14 lattice parameter (alat) = 7.5362 a.u. unit-cell volume = 428.0179 (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 = 50.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.536227 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 Br read from file: /users/gautes/Pseudo/Br.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 5d639767687de88634f0d235e6526498 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 1213 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 Br 7.00 79.90400 Br( 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: 41755 G-vectors FFT dimensions: ( 45, 45, 45) Smooth grid: 20341 G-vectors FFT dimensions: ( 36, 36, 36) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.05 Mb ( 156, 20) NL pseudopotentials 0.03 Mb ( 78, 28) Each V/rho on FFT grid 0.06 Mb ( 4050) Each G-vector array 0.01 Mb ( 1160) G-vector shells 0.00 Mb ( 258) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.19 Mb ( 156, 80) Each subspace H/S matrix 0.01 Mb ( 20, 20) Each matrix 0.02 Mb ( 28, 2, 20) Arrays for rho mixing 0.49 Mb ( 4050, 8) Check: negative/imaginary core charge= -0.000003 0.000000 Initial potential from superposition of free atoms starting charge 11.99997, renormalised to 12.00000 Starting wfc are 16 randomized atomic wfcs + 4 random wfc total cpu time spent up to now is 1.1 secs per-process dynamical memory: 10.0 Mb Self-consistent Calculation iteration # 1 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 4.8 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.22E-04, avg # of iterations = 4.0 total cpu time spent up to now is 2.1 secs total energy = -49.71319418 Ry Harris-Foulkes estimate = -49.72361457 Ry estimated scf accuracy < 0.01480939 Ry iteration # 2 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.23E-04, avg # of iterations = 2.0 total cpu time spent up to now is 2.6 secs total energy = -49.69526514 Ry Harris-Foulkes estimate = -49.76238621 Ry estimated scf accuracy < 0.38389297 Ry iteration # 3 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.23E-04, avg # of iterations = 1.9 total cpu time spent up to now is 3.0 secs total energy = -49.72114512 Ry Harris-Foulkes estimate = -49.72148786 Ry estimated scf accuracy < 0.00187068 Ry iteration # 4 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.56E-05, avg # of iterations = 1.6 total cpu time spent up to now is 3.4 secs total energy = -49.72128012 Ry Harris-Foulkes estimate = -49.72143141 Ry estimated scf accuracy < 0.00080887 Ry iteration # 5 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.74E-06, avg # of iterations = 1.5 total cpu time spent up to now is 3.9 secs total energy = -49.72134232 Ry Harris-Foulkes estimate = -49.72144099 Ry estimated scf accuracy < 0.00081784 Ry iteration # 6 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.74E-06, avg # of iterations = 1.1 total cpu time spent up to now is 4.2 secs total energy = -49.72139920 Ry Harris-Foulkes estimate = -49.72140090 Ry estimated scf accuracy < 0.00003422 Ry iteration # 7 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.85E-07, avg # of iterations = 1.1 total cpu time spent up to now is 4.6 secs total energy = -49.72139893 Ry Harris-Foulkes estimate = -49.72139961 Ry estimated scf accuracy < 0.00001002 Ry iteration # 8 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.35E-08, avg # of iterations = 1.0 total cpu time spent up to now is 5.0 secs total energy = -49.72139901 Ry Harris-Foulkes estimate = -49.72139915 Ry estimated scf accuracy < 0.00000120 Ry iteration # 9 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.96E-09, avg # of iterations = 1.0 total cpu time spent up to now is 5.4 secs total energy = -49.72139909 Ry Harris-Foulkes estimate = -49.72139909 Ry estimated scf accuracy < 0.00000004 Ry iteration # 10 ecut= 50.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.07E-10, avg # of iterations = 1.0 total cpu time spent up to now is 5.7 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2517 PWs) bands (ev): -14.6940 -14.6940 -12.7043 -12.7043 -2.3588 -2.3588 -1.9932 -1.9932 -1.9932 -1.9932 -0.6925 -0.6925 -0.4809 -0.4809 -0.4809 -0.4809 6.0013 6.0013 7.9744 7.9744 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 ( 2557 PWs) bands (ev): -14.6467 -14.6467 -12.7283 -12.7283 -2.3995 -2.3995 -2.0162 -2.0162 -1.8933 -1.8933 -0.8831 -0.8831 -0.6700 -0.6700 -0.5731 -0.5731 6.5267 6.5267 8.1895 8.5473 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 ( 2561 PWs) bands (ev): -14.5454 -14.5454 -12.7834 -12.7834 -2.6214 -2.6214 -1.8503 -1.8503 -1.6234 -1.6234 -1.3386 -1.3386 -0.9270 -0.9270 -0.8286 -0.8286 7.5305 7.5305 8.7658 8.7658 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.0000-0.5000 ( 2560 PWs) bands (ev): -14.4907 -14.4907 -12.8151 -12.8151 -2.7481 -2.7481 -1.6708 -1.6708 -1.5518 -1.5518 -1.3499 -1.3499 -1.1087 -1.1087 -1.0955 -1.0955 8.1011 8.1011 9.2225 9.2225 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.9921 0.9921 0.9795 0.9795 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 2558 PWs) bands (ev): -14.6036 -14.6036 -12.7433 -12.7433 -2.5133 -2.5133 -2.0262 -2.0262 -1.8096 -1.8096 -1.1686 -1.1686 -0.7294 -0.7294 -0.5450 -0.5450 6.9740 6.9740 8.5578 8.5578 occupation numbers 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 2564 PWs) bands (ev): -14.5124 -14.5124 -12.7784 -12.7784 -2.7681 -2.7681 -2.1421 -2.1421 -1.6553 -1.6553 -1.2534 -1.2534 -0.9179 -0.9179 -0.6867 -0.6867 7.8148 7.8148 9.0649 9.0649 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.1667-0.5000 ( 2544 PWs) bands (ev): -14.4637 -14.4637 -12.7990 -12.7990 -2.8958 -2.8958 -2.2796 -2.2796 -1.4994 -1.4994 -1.0859 -1.0859 -0.9970 -0.9970 -0.9508 -0.9508 8.3366 8.3366 9.2916 9.2917 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9595 0.9595 0.0331 0.0331 0.0011 0.0011 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 2562 PWs) bands (ev): -14.4440 -14.4440 -12.7699 -12.7699 -3.0414 -3.0414 -2.6849 -2.6849 -1.5381 -1.5381 -1.0179 -1.0179 -0.9558 -0.9558 -0.5993 -0.5993 8.1306 8.1306 8.6646 8.6646 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.1371 0.1371 0.0017 0.0017 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 2560 PWs) bands (ev): -14.4087 -14.4087 -12.7664 -12.7664 -3.1796 -3.1796 -2.9034 -2.9034 -1.4453 -1.4453 -1.0519 -1.0519 -0.7584 -0.7584 -0.6836 -0.6836 8.1500 8.1500 8.4545 8.4545 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.6596 0.6596 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 2544 PWs) bands (ev): -14.3805 -14.3805 -12.7499 -12.7499 -3.3407 -3.3407 -3.1193 -3.1193 -1.4152 -1.4152 -1.0348 -1.0348 -0.6623 -0.6623 -0.5831 -0.5831 7.7085 7.7085 8.4566 8.4566 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.3549 0.3549 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 2557 PWs) bands (ev): -14.5641 -14.5641 -12.7515 -12.7515 -2.5105 -2.5105 -2.2042 -2.2042 -1.8499 -1.8499 -1.2026 -1.2026 -0.7354 -0.7354 -0.6186 -0.6186 7.4279 7.4279 9.0564 9.0741 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 ( 2572 PWs) bands (ev): -14.4811 -14.4811 -12.7716 -12.7716 -2.7860 -2.7860 -2.2451 -2.2451 -1.9738 -1.9738 -1.0529 -1.0529 -0.9573 -0.9573 -0.6892 -0.6892 8.2401 8.2401 9.4999 9.5013 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.6766 0.6766 0.0019 0.0019 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 2556 PWs) bands (ev): -14.4374 -14.4374 -12.7838 -12.7838 -2.9486 -2.9486 -2.3864 -2.3864 -1.7734 -1.7734 -1.2114 -1.2114 -0.8603 -0.8603 -0.7869 -0.7869 8.8654 8.8654 9.1930 9.1930 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.3333 0.3333 ( 2560 PWs) bands (ev): -14.4168 -14.4168 -12.7589 -12.7589 -2.9489 -2.9489 -2.6490 -2.6490 -1.9371 -1.9371 -1.0897 -1.0897 -0.8472 -0.8472 -0.7354 -0.7354 8.5713 8.5713 9.2409 9.2409 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9689 0.9689 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 2552 PWs) bands (ev): -14.3836 -14.3836 -12.7530 -12.7530 -3.1000 -3.1000 -2.8028 -2.8028 -1.7773 -1.7773 -1.2207 -1.2207 -0.8084 -0.8084 -0.7262 -0.7262 8.4109 8.4109 9.2851 9.2851 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 ( 2556 PWs) bands (ev): -14.3561 -14.3561 -12.7374 -12.7374 -3.2202 -3.2202 -2.9837 -2.9837 -1.6996 -1.6996 -1.1656 -1.1656 -0.8292 -0.8292 -0.7520 -0.7520 8.0414 8.0414 9.3437 9.3437 occupation numbers 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.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 2533 PWs) bands (ev): -14.3610 -14.3610 -12.7369 -12.7369 -2.8693 -2.8693 -2.5187 -2.5187 -2.3970 -2.3970 -1.1772 -1.1772 -1.0921 -1.0921 -0.8303 -0.8303 8.9713 8.9713 9.6387 9.6413 occupation numbers 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.9739 0.9739 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 2536 PWs) bands (ev): -14.3325 -14.3325 -12.7260 -12.7260 -2.9523 -2.9523 -2.6189 -2.6189 -2.2081 -2.2081 -1.4241 -1.4241 -1.0342 -1.0342 -0.9619 -0.9619 8.8409 8.8409 9.3643 9.3645 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.3461 0.3461 0.0026 0.0026 0.0000 0.0000 0.0000 0.0000 k = 0.3333-0.5000-0.5000 ( 2516 PWs) bands (ev): -14.3064 -14.3064 -12.7123 -12.7123 -2.9307 -2.9307 -2.6205 -2.6205 -2.1781 -2.1781 -1.4117 -1.4117 -1.2327 -1.2327 -1.1842 -1.1842 8.7568 8.7568 9.4369 9.4379 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 ( 2512 PWs) bands (ev): -14.2811 -14.2811 -12.6996 -12.6996 -2.7833 -2.7833 -2.3071 -2.3071 -2.3071 -2.3071 -1.5397 -1.5397 -1.5217 -1.5217 -1.5217 -1.5217 9.1515 9.1525 9.1577 9.1577 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.0429 ev ! total energy = -49.72139908 Ry Harris-Foulkes estimate = -49.72139909 Ry estimated scf accuracy < 4.5E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -9.36845381 Ry hartree contribution = 8.51690301 Ry xc contribution = -13.56056441 Ry ewald contribution = -35.30885186 Ry smearing contrib. (-TS) = -0.00043201 Ry convergence has been achieved in 10 iterations Writing output data file BrN.save init_run : 0.24s CPU 0.31s WALL ( 1 calls) electrons : 3.77s CPU 4.84s WALL ( 1 calls) Called by init_run: wfcinit : 0.13s CPU 0.17s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 3.10s CPU 3.91s WALL ( 11 calls) sum_band : 0.56s CPU 0.57s WALL ( 11 calls) v_of_rho : 0.02s CPU 0.03s WALL ( 11 calls) v_h : 0.00s CPU 0.00s WALL ( 11 calls) v_xc : 0.02s CPU 0.03s WALL ( 11 calls) newd : 0.04s CPU 0.04s WALL ( 11 calls) mix_rho : 0.02s CPU 0.02s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.02s WALL ( 460 calls) cegterg : 3.02s CPU 3.08s WALL ( 220 calls) Called by sum_band: sum_band:bec : 0.04s CPU 0.03s WALL ( 220 calls) addusdens : 0.04s CPU 0.04s WALL ( 11 calls) Called by *egterg: h_psi : 1.90s CPU 1.96s WALL ( 661 calls) s_psi : 0.02s CPU 0.03s WALL ( 661 calls) g_psi : 0.00s CPU 0.01s WALL ( 421 calls) cdiaghg : 1.05s CPU 0.96s WALL ( 621 calls) cegterg:over : 0.04s CPU 0.06s WALL ( 421 calls) cegterg:upda : 0.02s CPU 0.06s WALL ( 421 calls) cegterg:last : 0.01s CPU 0.03s WALL ( 220 calls) cdiaghg:chol : 0.05s CPU 0.06s WALL ( 621 calls) cdiaghg:inve : 0.02s CPU 0.01s WALL ( 621 calls) cdiaghg:para : 0.09s CPU 0.08s WALL ( 1242 calls) Called by h_psi: h_psi:vloc : 1.80s CPU 1.87s WALL ( 661 calls) h_psi:vnl : 0.10s CPU 0.08s WALL ( 661 calls) add_vuspsi : 0.04s CPU 0.03s WALL ( 661 calls) General routines calbec : 0.09s CPU 0.07s WALL ( 881 calls) fft : 0.09s CPU 0.07s WALL ( 335 calls) ffts : 0.00s CPU 0.01s WALL ( 88 calls) fftw : 2.04s CPU 2.12s WALL ( 48560 calls) interpolate : 0.04s CPU 0.03s WALL ( 88 calls) Parallel routines fft_scatter : 1.07s CPU 1.09s WALL ( 48983 calls) PWSCF : 4.96s CPU 7.88s WALL This run was terminated on: 16: 2:59 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=