Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 17:27:55 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 47 47 13 1489 1477 227 Max 48 48 14 1492 1486 234 Sum 1725 1725 497 53643 53355 8289 bravais-lattice index = 14 lattice parameter (alat) = 8.2014 a.u. unit-cell volume = 551.6526 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 14.00 number of Kohn-Sham states= 22 kinetic-energy cutoff = 80.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)= 8.201411 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 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 PseudoPot. # 2 for Cs read from file: /users/gautes/Pseudo/Cs.rel-pbe-oncvpsp.UPF MD5 check sum: b26d75608c6aa7e713949ad9f7a49707 Pseudo is Norm-conserving, Zval = 9.0 Generated using ONCVPSP code by D. R. Hamann Using radial grid of 2414 points, 14 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 l(11) = 3 l(12) = 3 l(13) = 3 l(14) = 3 atomic species valence mass pseudopotential N 5.00 14.00670 N( 1.00) Cs 9.00 132.90550 Cs( 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: 53643 G-vectors FFT dimensions: ( 48, 48, 48) Smooth grid: 53355 G-vectors FFT dimensions: ( 48, 48, 48) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.13 Mb ( 380, 22) NL pseudopotentials 0.22 Mb ( 190, 76) Each V/rho on FFT grid 0.07 Mb ( 4608) Each G-vector array 0.01 Mb ( 1492) G-vector shells 0.00 Mb ( 350) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.51 Mb ( 380, 88) Each subspace H/S matrix 0.01 Mb ( 22, 22) Each matrix 0.05 Mb ( 76, 2, 22) Arrays for rho mixing 0.56 Mb ( 4608, 8) Check: negative/imaginary core charge= -0.000003 0.000000 Initial potential from superposition of free atoms starting charge 13.96180, renormalised to 14.00000 Starting wfc are 18 randomized atomic wfcs + 4 random wfc total cpu time spent up to now is 1.2 secs per-process dynamical memory: 38.1 Mb Self-consistent Calculation iteration # 1 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 5.5 total cpu time spent up to now is 2.6 secs total energy = -59.49388661 Ry Harris-Foulkes estimate = -59.67138422 Ry estimated scf accuracy < 0.22514863 Ry iteration # 2 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.61E-03, avg # of iterations = 3.2 total cpu time spent up to now is 3.9 secs total energy = -59.50583595 Ry Harris-Foulkes estimate = -59.67429826 Ry estimated scf accuracy < 0.38821842 Ry iteration # 3 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.61E-03, avg # of iterations = 2.0 total cpu time spent up to now is 4.9 secs total energy = -59.59262651 Ry Harris-Foulkes estimate = -59.59560414 Ry estimated scf accuracy < 0.00931508 Ry iteration # 4 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.65E-05, avg # of iterations = 2.8 total cpu time spent up to now is 6.0 secs total energy = -59.59396786 Ry Harris-Foulkes estimate = -59.59405398 Ry estimated scf accuracy < 0.00038266 Ry iteration # 5 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap c_bands: 1 eigenvalues not converged c_bands: 1 eigenvalues not converged ethr = 2.73E-06, avg # of iterations = 6.3 total cpu time spent up to now is 7.6 secs total energy = -59.59406756 Ry Harris-Foulkes estimate = -59.59410289 Ry estimated scf accuracy < 0.00009252 Ry iteration # 6 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.61E-07, avg # of iterations = 2.8 total cpu time spent up to now is 8.6 secs total energy = -59.59407955 Ry Harris-Foulkes estimate = -59.59408052 Ry estimated scf accuracy < 0.00000221 Ry iteration # 7 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.58E-08, avg # of iterations = 3.1 total cpu time spent up to now is 9.9 secs total energy = -59.59408066 Ry Harris-Foulkes estimate = -59.59408066 Ry estimated scf accuracy < 0.00000001 Ry iteration # 8 ecut= 80.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.32E-11, avg # of iterations = 4.2 total cpu time spent up to now is 11.4 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 6667 PWs) bands (ev): -17.7017 -17.7017 -8.8788 -8.8788 -6.1853 -6.1853 -4.6702 -4.6702 -4.6702 -4.6702 2.5193 2.5193 2.5640 2.5640 2.5640 2.5640 5.3403 5.3403 6.2536 6.2536 6.2536 6.2536 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 0.0000 0.0000 k = 0.0000 0.0000 0.1667 ( 6679 PWs) bands (ev): -17.6890 -17.6890 -8.8866 -8.8866 -6.2108 -6.2108 -4.7541 -4.7541 -4.6510 -4.6510 2.2899 2.2899 2.5108 2.5108 2.5354 2.5354 5.8471 5.8471 6.3762 6.3762 6.7445 6.7445 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9988 0.9988 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.3333 ( 6676 PWs) bands (ev): -17.6634 -17.6634 -8.9050 -8.9050 -6.2712 -6.2712 -4.8907 -4.8907 -4.6119 -4.6119 1.8873 1.8873 2.4585 2.4585 2.4751 2.4751 6.6663 6.6663 6.7675 6.7675 7.9975 7.9975 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.0036 0.0036 0.0011 0.0011 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000-0.5000 ( 6696 PWs) bands (ev): -17.6504 -17.6504 -8.9157 -8.9157 -6.3057 -6.3057 -4.9450 -4.9450 -4.5920 -4.5920 1.7181 1.7181 2.4315 2.4315 2.4434 2.4434 6.8452 6.8452 7.2413 7.2413 8.7781 8.7781 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.0255 0.0255 0.0108 0.0108 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.1667 ( 6675 PWs) bands (ev): -17.6765 -17.6765 -8.8871 -8.8871 -6.2400 -6.2400 -4.8010 -4.8010 -4.6839 -4.6839 2.2140 2.2140 2.4361 2.4361 2.5131 2.5131 6.1979 6.1979 6.6931 6.6931 7.0246 7.0246 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.0183 0.0183 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667 0.3333 ( 6670 PWs) bands (ev): -17.6513 -17.6513 -8.8903 -8.8903 -6.3079 -6.3079 -4.9414 -4.9414 -4.6713 -4.6713 1.9715 1.9715 2.3886 2.3886 2.4840 2.4840 6.7334 6.7334 7.1627 7.1627 7.9483 7.9483 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.3796 0.3796 0.0006 0.0006 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1667-0.5000 ( 6678 PWs) bands (ev): -17.6386 -17.6386 -8.8932 -8.8932 -6.3462 -6.3462 -5.0034 -5.0034 -4.6590 -4.6590 1.8524 1.8524 2.3725 2.3725 2.4682 2.4682 6.9416 6.9416 7.6576 7.6576 7.7823 7.7823 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.6664 0.6664 0.0017 0.0017 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333 0.3333 ( 6670 PWs) bands (ev): -17.6269 -17.6269 -8.8611 -8.8611 -6.3920 -6.3920 -5.0877 -5.0877 -4.7118 -4.7118 2.0847 2.0847 2.2776 2.2776 2.5154 2.5154 6.7372 6.7372 7.1307 7.1307 8.8990 8.8992 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.9995 0.9995 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.3333-0.5000 ( 6658 PWs) bands (ev): -17.6146 -17.6146 -8.8467 -8.8467 -6.4393 -6.4393 -5.1634 -5.1634 -4.7146 -4.7146 2.0966 2.0966 2.2622 2.2622 2.5310 2.5310 6.8065 6.8065 6.8426 6.8426 8.8749 8.8749 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.9998 0.9998 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000-0.5000-0.5000 ( 6632 PWs) bands (ev): -17.6026 -17.6026 -8.8228 -8.8228 -6.4920 -6.4920 -5.2474 -5.2474 -4.7227 -4.7227 2.1935 2.1935 2.2330 2.2330 2.5635 2.5635 6.4958 6.4958 6.7146 6.7146 8.7396 8.7396 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.1667 ( 6700 PWs) bands (ev): -17.6642 -17.6642 -8.8828 -8.8828 -6.2692 -6.2692 -4.8037 -4.8037 -4.7707 -4.7707 2.2331 2.2331 2.3841 2.3841 2.4157 2.4157 6.5805 6.5805 7.1016 7.1016 7.1337 7.1337 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.4592 0.4592 0.0769 0.0769 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667 0.3333 ( 6654 PWs) bands (ev): -17.6394 -17.6394 -8.8759 -8.8759 -6.3383 -6.3383 -4.9580 -4.9580 -4.7621 -4.7621 2.0452 2.0452 2.3181 2.3181 2.4271 2.4271 7.1706 7.1706 7.3368 7.3368 7.9691 7.9691 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.9909 0.9909 0.0348 0.0348 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.1667-0.5000 ( 6642 PWs) bands (ev): -17.6268 -17.6268 -8.8735 -8.8735 -6.3779 -6.3779 -5.0276 -5.0276 -4.7503 -4.7503 1.9330 1.9330 2.2762 2.2762 2.4663 2.4663 7.4665 7.4665 7.6239 7.6239 7.6800 7.6800 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.9996 0.9996 0.0020 0.0020 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333 0.3333 ( 6664 PWs) bands (ev): -17.6153 -17.6153 -8.8473 -8.8473 -6.4111 -6.4111 -5.0699 -5.0699 -4.8302 -4.8302 2.0952 2.0952 2.2026 2.2026 2.3787 2.3787 7.2248 7.2248 7.5067 7.5067 8.6348 8.6348 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 0.5580 0.5580 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667 0.3333-0.5000 ( 6652 PWs) bands (ev): -17.6031 -17.6031 -8.8331 -8.8331 -6.4537 -6.4537 -5.1430 -5.1430 -4.8315 -4.8315 2.0622 2.0622 2.1541 2.1541 2.4014 2.4014 7.0962 7.0962 7.4514 7.4514 8.4556 8.4557 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 0.1922 0.1922 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1667-0.5000-0.5000 ( 6656 PWs) bands (ev): -17.5912 -17.5912 -8.8124 -8.8124 -6.4981 -6.4981 -5.2179 -5.2179 -4.8392 -4.8392 2.1022 2.1022 2.1305 2.1305 2.3465 2.3465 6.8654 6.8654 7.3895 7.3895 8.4397 8.4413 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 0.9308 0.9308 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.3333 0.3333 0.3333 ( 6665 PWs) bands (ev): -17.5918 -17.5918 -8.8195 -8.8195 -6.4621 -6.4621 -5.0367 -5.0367 -5.0303 -5.0303 2.0383 2.0383 2.0579 2.0579 2.1858 2.1858 7.5588 7.5588 8.2092 8.2092 8.2441 8.2441 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 0.0000 0.0000 k = 0.3333 0.3333-0.5000 ( 6672 PWs) bands (ev): -17.5799 -17.5799 -8.8055 -8.8055 -6.4947 -6.4947 -5.1196 -5.1196 -5.0145 -5.0145 1.9676 1.9676 1.9761 1.9761 2.1432 2.1432 7.5087 7.5087 8.0739 8.0739 8.2447 8.2447 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 0.0000 0.0000 k = 0.3333-0.5000-0.5000 ( 6692 PWs) bands (ev): -17.5682 -17.5682 -8.7916 -8.7916 -6.5223 -6.5223 -5.1571 -5.1571 -5.0406 -5.0406 1.9194 1.9194 1.9240 1.9240 1.9830 1.9830 7.5846 7.5846 8.0706 8.0706 8.1006 8.1006 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 0.0000 0.0000 k =-0.5000-0.5000-0.5000 ( 6704 PWs) bands (ev): -17.5565 -17.5565 -8.7810 -8.7810 -6.5413 -6.5413 -5.1259 -5.1259 -5.1259 -5.1259 1.8176 1.8176 1.8279 1.8279 1.8279 1.8279 7.9453 7.9453 7.9453 7.9453 7.9841 7.9841 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 0.0000 0.0000 the Fermi energy is 2.3819 ev ! total energy = -59.59408068 Ry Harris-Foulkes estimate = -59.59408068 Ry estimated scf accuracy < 6.7E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -13.70444990 Ry hartree contribution = 11.15376884 Ry xc contribution = -11.57159395 Ry ewald contribution = -45.47117314 Ry smearing contrib. (-TS) = -0.00063253 Ry convergence has been achieved in 8 iterations Writing output data file CsN.save init_run : 0.48s CPU 0.55s WALL ( 1 calls) electrons : 9.98s CPU 10.25s WALL ( 1 calls) Called by init_run: wfcinit : 0.27s CPU 0.29s WALL ( 1 calls) potinit : 0.02s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 8.19s CPU 8.39s WALL ( 8 calls) sum_band : 0.97s CPU 0.99s WALL ( 8 calls) v_of_rho : 0.03s CPU 0.03s WALL ( 9 calls) v_h : 0.00s CPU 0.00s WALL ( 9 calls) v_xc : 0.03s CPU 0.03s WALL ( 9 calls) newd : 0.83s CPU 0.83s WALL ( 9 calls) mix_rho : 0.03s CPU 0.02s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.06s CPU 0.04s WALL ( 340 calls) cegterg : 7.14s CPU 7.28s WALL ( 160 calls) Called by sum_band: sum_band:bec : 0.01s CPU 0.01s WALL ( 160 calls) addusdens : 0.03s CPU 0.03s WALL ( 8 calls) Called by *egterg: h_psi : 5.07s CPU 5.22s WALL ( 780 calls) s_psi : 0.12s CPU 0.10s WALL ( 780 calls) g_psi : 0.02s CPU 0.02s WALL ( 600 calls) cdiaghg : 1.59s CPU 1.62s WALL ( 760 calls) cegterg:over : 0.12s CPU 0.17s WALL ( 600 calls) cegterg:upda : 0.19s CPU 0.19s WALL ( 600 calls) cegterg:last : 0.06s CPU 0.06s WALL ( 164 calls) cdiaghg:chol : 0.11s CPU 0.10s WALL ( 760 calls) cdiaghg:inve : 0.01s CPU 0.02s WALL ( 760 calls) cdiaghg:para : 0.12s CPU 0.12s WALL ( 1520 calls) Called by h_psi: h_psi:vloc : 4.56s CPU 4.73s WALL ( 780 calls) h_psi:vnl : 0.50s CPU 0.47s WALL ( 780 calls) add_vuspsi : 0.28s CPU 0.28s WALL ( 780 calls) General routines calbec : 0.28s CPU 0.24s WALL ( 940 calls) fft : 0.02s CPU 0.06s WALL ( 263 calls) ffts : 0.02s CPU 0.01s WALL ( 68 calls) fftw : 4.99s CPU 5.24s WALL ( 48016 calls) interpolate : 0.03s CPU 0.03s WALL ( 68 calls) Parallel routines fft_scatter : 1.86s CPU 1.94s WALL ( 48347 calls) PWSCF : 11.58s CPU 12.85s WALL This run was terminated on: 17:28: 8 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=