Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 16:33:49 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 Cd.rel-pbe-n-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4D 4D renormalized file Ca.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 3S 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 74 21 6 2892 441 70 Max 75 22 7 2895 458 75 Sum 2693 777 233 104199 16135 2601 bravais-lattice index = 14 lattice parameter (alat) = 7.2547 a.u. unit-cell volume = 381.8131 (a.u.)^3 number of atoms/cell = 2 number of atomic types = 2 number of electrons = 22.00 number of Kohn-Sham states= 30 kinetic-energy cutoff = 46.0000 Ry charge density cutoff = 639.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.254658 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 Cd read from file: /users/gautes/Pseudo/Cd.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: 1d41bce1441572f114805845d9d822b0 Pseudo is Ultrasoft + core correction, Zval = 12.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1239 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 Ca read from file: /users/gautes/Pseudo/Ca.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 60bd286d514994a6a955285ae841e6ef Pseudo is Ultrasoft + core correction, Zval = 10.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1169 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 atomic species valence mass pseudopotential Cd 12.00 112.41100 Cd( 1.00) Ca 10.00 40.07800 Ca( 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.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0029155 k( 2) = ( 0.0000000 0.0000000 0.1428571), wk = 0.0174927 k( 3) = ( 0.0000000 0.0000000 0.2857143), wk = 0.0174927 k( 4) = ( 0.0000000 0.0000000 0.4285714), wk = 0.0174927 k( 5) = ( 0.0000000 0.1428571 0.1428571), wk = 0.0349854 k( 6) = ( 0.0000000 0.1428571 0.2857143), wk = 0.0699708 k( 7) = ( 0.0000000 0.1428571 0.4285714), wk = 0.0699708 k( 8) = ( 0.0000000 0.2857143 0.2857143), wk = 0.0349854 k( 9) = ( 0.0000000 0.2857143 0.4285714), wk = 0.0699708 k( 10) = ( 0.0000000 0.4285714 0.4285714), wk = 0.0349854 k( 11) = ( 0.1428571 0.1428571 0.1428571), wk = 0.0233236 k( 12) = ( 0.1428571 0.1428571 0.2857143), wk = 0.0699708 k( 13) = ( 0.1428571 0.1428571 0.4285714), wk = 0.0699708 k( 14) = ( 0.1428571 0.2857143 0.2857143), wk = 0.0699708 k( 15) = ( 0.1428571 0.2857143 0.4285714), wk = 0.1399417 k( 16) = ( 0.1428571 0.4285714 0.4285714), wk = 0.0699708 k( 17) = ( 0.2857143 0.2857143 0.2857143), wk = 0.0233236 k( 18) = ( 0.2857143 0.2857143 0.4285714), wk = 0.0699708 k( 19) = ( 0.2857143 0.4285714 0.4285714), wk = 0.0699708 k( 20) = ( 0.4285714 0.4285714 0.4285714), wk = 0.0233236 Dense grid: 104199 G-vectors FFT dimensions: ( 60, 60, 60) Smooth grid: 16135 G-vectors FFT dimensions: ( 32, 32, 32) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.05 Mb ( 110, 30) NL pseudopotentials 0.06 Mb ( 55, 68) Each V/rho on FFT grid 0.11 Mb ( 7200) Each G-vector array 0.02 Mb ( 2894) G-vector shells 0.00 Mb ( 578) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.20 Mb ( 110, 120) Each subspace H/S matrix 0.01 Mb ( 30, 30) Each matrix 0.06 Mb ( 68, 2, 30) Arrays for rho mixing 0.88 Mb ( 7200, 8) Initial potential from superposition of free atoms starting charge 21.99128, renormalised to 22.00000 Starting wfc are 34 randomized atomic wfcs total cpu time spent up to now is 2.8 secs per-process dynamical memory: 51.7 Mb Self-consistent Calculation iteration # 1 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 8.7 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.08E-04, avg # of iterations = 8.4 total cpu time spent up to now is 5.7 secs total energy = -175.48449803 Ry Harris-Foulkes estimate = -175.49604258 Ry estimated scf accuracy < 0.01758128 Ry iteration # 2 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.99E-05, avg # of iterations = 5.1 total cpu time spent up to now is 6.9 secs total energy = -175.48863692 Ry Harris-Foulkes estimate = -175.49469983 Ry estimated scf accuracy < 0.01051648 Ry iteration # 3 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.78E-05, avg # of iterations = 2.7 total cpu time spent up to now is 7.7 secs total energy = -175.49123257 Ry Harris-Foulkes estimate = -175.49120253 Ry estimated scf accuracy < 0.00010998 Ry iteration # 4 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.00E-07, avg # of iterations = 5.5 total cpu time spent up to now is 9.0 secs total energy = -175.49144802 Ry Harris-Foulkes estimate = -175.49147868 Ry estimated scf accuracy < 0.00008066 Ry iteration # 5 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.67E-07, avg # of iterations = 2.1 total cpu time spent up to now is 9.7 secs total energy = -175.49145021 Ry Harris-Foulkes estimate = -175.49145485 Ry estimated scf accuracy < 0.00000973 Ry iteration # 6 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.42E-08, avg # of iterations = 5.5 total cpu time spent up to now is 10.8 secs total energy = -175.49145707 Ry Harris-Foulkes estimate = -175.49145793 Ry estimated scf accuracy < 0.00000186 Ry iteration # 7 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.44E-09, avg # of iterations = 2.2 total cpu time spent up to now is 11.5 secs total energy = -175.49145685 Ry Harris-Foulkes estimate = -175.49145721 Ry estimated scf accuracy < 0.00000047 Ry iteration # 8 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.14E-09, avg # of iterations = 4.8 total cpu time spent up to now is 12.5 secs total energy = -175.49145716 Ry Harris-Foulkes estimate = -175.49145724 Ry estimated scf accuracy < 0.00000016 Ry iteration # 9 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.16E-10, avg # of iterations = 2.4 total cpu time spent up to now is 13.3 secs total energy = -175.49145716 Ry Harris-Foulkes estimate = -175.49145718 Ry estimated scf accuracy < 0.00000003 Ry iteration # 10 ecut= 46.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.35E-10, avg # of iterations = 5.0 total cpu time spent up to now is 14.3 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 2007 PWs) bands (ev): -33.3336 -33.3336 -14.5561 -14.5561 -14.1768 -14.1768 -14.1768 -14.1768 -0.8410 -0.8410 -0.8410 -0.8410 -0.1667 -0.1667 -0.1667 -0.1667 -0.0450 -0.0450 1.9345 1.9345 9.6074 9.6074 9.6074 9.6074 12.0879 12.0879 12.4639 12.4641 12.4642 12.4642 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 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.1429 ( 2007 PWs) bands (ev): -33.3324 -33.3324 -14.5616 -14.5616 -14.1891 -14.1891 -14.1754 -14.1754 -0.8385 -0.8385 -0.8258 -0.8258 -0.1648 -0.1648 -0.1542 -0.1542 -0.0521 -0.0521 2.1681 2.1681 9.1161 9.1161 9.6531 9.6531 11.3477 11.3478 11.4681 11.4682 11.8776 11.8776 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 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.2857 ( 2004 PWs) bands (ev): -33.3298 -33.3298 -14.5754 -14.5754 -14.2145 -14.2145 -14.1721 -14.1721 -0.8334 -0.8334 -0.7990 -0.7990 -0.1626 -0.1626 -0.1203 -0.1203 -0.0654 -0.0654 2.8564 2.8564 7.5632 7.5632 9.7638 9.7638 10.7284 10.7284 10.8223 10.8223 12.1004 12.1004 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4286 ( 2008 PWs) bands (ev): -33.3277 -33.3277 -14.5885 -14.5885 -14.2327 -14.2327 -14.1694 -14.1694 -0.8297 -0.8297 -0.7848 -0.7848 -0.1636 -0.1636 -0.0882 -0.0882 -0.0727 -0.0727 3.8925 3.8925 5.8854 5.8854 9.8618 9.8618 10.3964 10.3964 10.4826 10.4826 12.2809 12.7884 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.1429 ( 1997 PWs) bands (ev): -33.3313 -33.3313 -14.5665 -14.5665 -14.1945 -14.1945 -14.1804 -14.1804 -0.8314 -0.8314 -0.8181 -0.8181 -0.1589 -0.1589 -0.1470 -0.1470 -0.0616 -0.0616 2.3971 2.3971 9.1347 9.1347 9.3221 9.3221 10.7049 10.7049 11.2501 11.2502 12.4894 12.4896 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.2857 ( 2012 PWs) bands (ev): -33.3287 -33.3287 -14.5794 -14.5794 -14.2177 -14.2177 -14.1793 -14.1793 -0.8271 -0.8271 -0.7986 -0.7986 -0.1576 -0.1576 -0.1162 -0.1162 -0.0784 -0.0784 3.0670 3.0670 7.7266 7.7266 9.4721 9.4721 10.4341 10.4341 10.8694 10.8694 12.0684 12.0684 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.1429 0.4286 ( 2001 PWs) bands (ev): -33.3265 -33.3265 -14.5914 -14.5914 -14.2356 -14.2356 -14.1766 -14.1766 -0.8252 -0.8252 -0.7887 -0.7887 -0.1613 -0.1613 -0.0959 -0.0959 -0.0746 -0.0746 4.0526 4.0526 6.1195 6.1195 9.7878 9.7878 10.1320 10.1320 10.6059 10.6059 12.0994 12.0995 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.2857 ( 2007 PWs) bands (ev): -33.3261 -33.3261 -14.5897 -14.5897 -14.2321 -14.2321 -14.1861 -14.1861 -0.8224 -0.8224 -0.7978 -0.7978 -0.1502 -0.1502 -0.1258 -0.1258 -0.0720 -0.0720 3.6715 3.6715 8.0428 8.0428 8.3199 8.3199 10.3473 10.3473 11.0204 11.0204 11.7770 11.7771 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 0.2501 0.2501 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2857 0.4286 ( 2020 PWs) bands (ev): -33.3240 -33.3240 -14.5998 -14.5998 -14.2472 -14.2472 -14.1862 -14.1862 -0.8238 -0.8238 -0.8014 -0.8014 -0.1592 -0.1592 -0.1361 -0.1361 -0.0400 -0.0400 4.5037 4.5037 6.7773 6.7773 8.7545 8.7545 10.0884 10.0884 11.0163 11.0163 11.7562 11.7562 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4286 0.4286 ( 2013 PWs) bands (ev): -33.3220 -33.3220 -14.6079 -14.6079 -14.2599 -14.2599 -14.1884 -14.1884 -0.8275 -0.8275 -0.8158 -0.8158 -0.1757 -0.1757 -0.1496 -0.1496 -0.0003 -0.0003 5.0768 5.0768 7.5325 7.5325 7.7698 7.7698 9.6250 9.6250 11.4006 11.4006 11.5419 11.5419 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.1429 ( 1986 PWs) bands (ev): -33.3301 -33.3301 -14.5710 -14.5710 -14.1938 -14.1938 -14.1916 -14.1916 -0.8208 -0.8208 -0.8159 -0.8159 -0.1480 -0.1480 -0.1440 -0.1440 -0.0733 -0.0733 2.6198 2.6198 9.1814 9.1814 9.2415 9.2415 10.1604 10.1604 11.3599 11.3599 11.4069 11.4070 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.2857 ( 2009 PWs) bands (ev): -33.3275 -33.3275 -14.5831 -14.5831 -14.2185 -14.2185 -14.1888 -14.1888 -0.8192 -0.8192 -0.7992 -0.7992 -0.1487 -0.1487 -0.1098 -0.1098 -0.0942 -0.0942 3.2691 3.2691 7.8906 7.8906 9.4091 9.4091 10.2282 10.2282 10.6793 10.6793 11.6694 11.6694 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.1429 0.4286 ( 2013 PWs) bands (ev): -33.3254 -33.3254 -14.5945 -14.5945 -14.2370 -14.2370 -14.1856 -14.1856 -0.8200 -0.8200 -0.7916 -0.7916 -0.1570 -0.1570 -0.1032 -0.1032 -0.0741 -0.0741 4.2090 4.2090 6.3421 6.3421 9.6555 9.6555 10.1991 10.1991 10.5350 10.5350 11.8219 11.8219 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.2857 0.2857 ( 1998 PWs) bands (ev): -33.3249 -33.3249 -14.5929 -14.5929 -14.2300 -14.2300 -14.1986 -14.1986 -0.8135 -0.8135 -0.7987 -0.7987 -0.1365 -0.1365 -0.1281 -0.1281 -0.0779 -0.0779 3.8529 3.8529 8.2491 8.2491 8.3748 8.3748 10.3195 10.3195 10.6034 10.6034 11.4235 11.4235 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.2857 0.4286 ( 2016 PWs) bands (ev): -33.3229 -33.3229 -14.6026 -14.6026 -14.2457 -14.2457 -14.1983 -14.1983 -0.8169 -0.8169 -0.8000 -0.8000 -0.1604 -0.1604 -0.1239 -0.1239 -0.0416 -0.0416 4.6584 4.6584 6.9519 6.9519 8.7468 8.7468 10.2784 10.2784 10.8268 10.8268 11.4836 11.4836 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.1429 0.4286 0.4286 ( 2030 PWs) bands (ev): -33.3209 -33.3209 -14.6107 -14.6107 -14.2580 -14.2580 -14.2015 -14.2015 -0.8194 -0.8194 -0.8115 -0.8115 -0.1824 -0.1824 -0.1252 -0.1252 -0.0014 -0.0014 5.2467 5.2467 7.6073 7.6073 7.7709 7.7709 9.8566 9.8566 11.3279 11.3279 11.5642 11.5642 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.2857 0.2857 ( 2028 PWs) bands (ev): -33.3224 -33.3224 -14.6018 -14.6018 -14.2264 -14.2264 -14.2244 -14.2244 -0.8063 -0.8063 -0.7996 -0.7996 -0.1568 -0.1568 -0.0878 -0.0878 -0.0765 -0.0765 4.3832 4.3832 8.3748 8.3748 8.4960 8.4960 9.5303 9.5303 11.0047 11.0047 11.0078 11.0078 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.2857 0.4286 ( 2007 PWs) bands (ev): -33.3203 -33.3203 -14.6102 -14.6102 -14.2448 -14.2448 -14.2215 -14.2215 -0.8109 -0.8109 -0.7981 -0.7981 -0.1790 -0.1790 -0.0805 -0.0805 -0.0339 -0.0339 5.1181 5.1181 7.3711 7.3711 8.4897 8.4897 10.3436 10.3436 10.8062 10.8062 11.2303 11.2303 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2857 0.4286 0.4286 ( 2011 PWs) bands (ev): -33.3183 -33.3183 -14.6179 -14.6179 -14.2535 -14.2535 -14.2292 -14.2292 -0.8099 -0.8099 -0.8021 -0.8021 -0.1976 -0.1976 -0.0598 -0.0598 -0.0009 -0.0009 5.7622 5.7622 7.6726 7.6726 7.8060 7.8060 10.3616 10.3616 10.7518 10.7518 11.6711 11.6711 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4286 0.4286 0.4286 ( 1995 PWs) bands (ev): -33.3162 -33.3162 -14.6249 -14.6249 -14.2500 -14.2500 -14.2497 -14.2497 -0.8050 -0.8050 -0.8007 -0.8007 -0.2104 -0.2104 -0.0002 -0.0002 0.0022 0.0022 6.5452 6.5452 7.6149 7.6149 7.7416 7.7416 9.1727 9.1727 11.6619 11.6619 11.8090 11.8092 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 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 8.0278 ev ! total energy = -175.49145717 Ry Harris-Foulkes estimate = -175.49145718 Ry estimated scf accuracy < 5.3E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -60.46171484 Ry hartree contribution = 41.88447685 Ry xc contribution = -34.95595526 Ry ewald contribution = -121.95822450 Ry smearing contrib. (-TS) = -0.00003942 Ry convergence has been achieved in 10 iterations Writing output data file CaCd.save init_run : 0.57s CPU 0.78s WALL ( 1 calls) electrons : 10.83s CPU 11.65s WALL ( 1 calls) Called by init_run: wfcinit : 0.20s CPU 0.26s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 8.04s CPU 8.68s WALL ( 11 calls) sum_band : 1.95s CPU 1.96s WALL ( 11 calls) v_of_rho : 0.06s CPU 0.05s WALL ( 11 calls) v_h : 0.00s CPU 0.00s WALL ( 11 calls) v_xc : 0.05s CPU 0.05s WALL ( 11 calls) newd : 0.78s CPU 0.83s WALL ( 11 calls) mix_rho : 0.02s CPU 0.03s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.03s CPU 0.02s WALL ( 460 calls) cegterg : 7.77s CPU 7.97s WALL ( 220 calls) Called by sum_band: sum_band:bec : 0.43s CPU 0.44s WALL ( 220 calls) addusdens : 0.80s CPU 0.81s WALL ( 11 calls) Called by *egterg: h_psi : 3.26s CPU 3.39s WALL ( 1290 calls) s_psi : 0.20s CPU 0.19s WALL ( 1290 calls) g_psi : 0.02s CPU 0.01s WALL ( 1050 calls) cdiaghg : 3.84s CPU 3.95s WALL ( 1250 calls) cegterg:over : 0.20s CPU 0.22s WALL ( 1050 calls) cegterg:upda : 0.15s CPU 0.17s WALL ( 1050 calls) cegterg:last : 0.07s CPU 0.06s WALL ( 234 calls) cdiaghg:chol : 0.18s CPU 0.20s WALL ( 1250 calls) cdiaghg:inve : 0.04s CPU 0.08s WALL ( 1250 calls) cdiaghg:para : 0.26s CPU 0.25s WALL ( 2500 calls) Called by h_psi: h_psi:vloc : 2.85s CPU 2.95s WALL ( 1290 calls) h_psi:vnl : 0.41s CPU 0.43s WALL ( 1290 calls) add_vuspsi : 0.23s CPU 0.25s WALL ( 1290 calls) General routines calbec : 0.21s CPU 0.23s WALL ( 1510 calls) fft : 0.11s CPU 0.11s WALL ( 335 calls) ffts : 0.01s CPU 0.01s WALL ( 88 calls) fftw : 2.98s CPU 3.01s WALL ( 91692 calls) interpolate : 0.04s CPU 0.04s WALL ( 88 calls) Parallel routines fft_scatter : 2.03s CPU 1.96s WALL ( 92115 calls) PWSCF : 13.58s CPU 15.96s WALL This run was terminated on: 16:34: 5 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=