Program PWSCF v.5.3.0 (svn rev. 11974) starts on 24Dec2016 at 19:41: 2 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 Y.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s) 4S 4P 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 47 28 8 1461 684 106 Max 48 29 9 1464 697 111 Sum 1701 1033 293 52635 24837 3887 bravais-lattice index = 14 lattice parameter (alat) = 9.1047 a.u. unit-cell volume = 754.7392 (a.u.)^3 number of atoms/cell = 4 number of atomic types = 2 number of electrons = 36.00 number of Kohn-Sham states= 44 kinetic-energy cutoff = 39.0000 Ry charge density cutoff = 258.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)= 9.104700 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 Al read from file: /users/gautes/Pseudo/Al.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: d344aa71f13a435d235890da883ec65e Pseudo is Ultrasoft + core correction, Zval = 3.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1135 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 Y read from file: /users/gautes/Pseudo/Y.rel-pbe-spn-rrkjus_psl.1.0.0.UPF MD5 check sum: 25a826a599bbcba7ebf619f9598d90b1 Pseudo is Ultrasoft + core correction, Zval = 11.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1223 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 Al 3.00 26.98150 Al( 1.00) Y 11.00 88.90590 Y( 1.00) 48 Sym. Ops., with inversion, found s frac. trans. isym = 1 identity cryst. s( 1) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s( 1) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 2 180 deg rotation - cart. axis [0,0,1] cryst. s( 2) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 3 180 deg rotation - cart. axis [0,1,0] cryst. s( 3) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s( 3) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 4 180 deg rotation - cart. axis [1,0,0] cryst. s( 4) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s( 4) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 5 180 deg rotation - cart. axis [1,1,0] cryst. s( 5) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s( 5) = ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 6 180 deg rotation - cart. axis [1,-1,0] cryst. s( 6) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s( 6) = ( -0.0000000 -1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 7 90 deg rotation - cart. axis [0,0,-1] cryst. s( 7) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s( 7) = ( -0.0000000 1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 8 90 deg rotation - cart. axis [0,0,1] cryst. s( 8) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s( 8) = ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 9 180 deg rotation - cart. axis [1,0,1] cryst. s( 9) = ( 0 0 1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s( 9) = ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 10 180 deg rotation - cart. axis [-1,0,1] cryst. s(10) = ( 0 0 -1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(10) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 11 90 deg rotation - cart. axis [0,1,0] cryst. s(11) = ( 0 0 -1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(11) = ( -0.0000000 0.0000000 1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 12 90 deg rotation - cart. axis [0,-1,0] cryst. s(12) = ( 0 0 1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(12) = ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 13 180 deg rotation - cart. axis [0,1,1] cryst. s(13) = ( -1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(13) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 14 180 deg rotation - cart. axis [0,1,-1] cryst. s(14) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(14) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 15 90 deg rotation - cart. axis [-1,0,0] cryst. s(15) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(15) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 16 90 deg rotation - cart. axis [1,0,0] cryst. s(16) = ( 1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(16) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 17 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(17) = ( 0 0 1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(17) = ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 18 120 deg rotation - cart. axis [-1,1,1] cryst. s(18) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(18) = ( -0.0000000 -1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 19 120 deg rotation - cart. axis [1,1,-1] cryst. s(19) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(19) = ( -0.0000000 1.0000000 -0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 20 120 deg rotation - cart. axis [1,-1,1] cryst. s(20) = ( 0 0 1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(20) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 21 120 deg rotation - cart. axis [1,1,1] cryst. s(21) = ( 0 1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(21) = ( 0.0000000 0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 22 120 deg rotation - cart. axis [-1,1,-1] cryst. s(22) = ( 0 -1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(22) = ( -0.0000000 -0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 23 120 deg rotation - cart. axis [1,-1,-1] cryst. s(23) = ( 0 -1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(23) = ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 24 120 deg rotation - cart. axis [-1,-1,1] cryst. s(24) = ( 0 1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(24) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 25 inversion cryst. s(25) = ( -1 0 0 ) ( 0 -1 0 ) ( 0 0 -1 ) cart. s(25) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 26 inv. 180 deg rotation - cart. axis [0,0,1] cryst. s(26) = ( 1 0 0 ) ( 0 1 0 ) ( 0 0 -1 ) cart. s(26) = ( 1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 27 inv. 180 deg rotation - cart. axis [0,1,0] cryst. s(27) = ( 1 0 0 ) ( 0 -1 0 ) ( 0 0 1 ) cart. s(27) = ( 1.0000000 -0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 28 inv. 180 deg rotation - cart. axis [1,0,0] cryst. s(28) = ( -1 0 0 ) ( 0 1 0 ) ( 0 0 1 ) cart. s(28) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 29 inv. 180 deg rotation - cart. axis [1,1,0] cryst. s(29) = ( 0 -1 0 ) ( -1 0 0 ) ( 0 0 1 ) cart. s(29) = ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 30 inv. 180 deg rotation - cart. axis [1,-1,0] cryst. s(30) = ( 0 1 0 ) ( 1 0 0 ) ( 0 0 1 ) cart. s(30) = ( 0.0000000 1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) isym = 31 inv. 90 deg rotation - cart. axis [0,0,-1] cryst. s(31) = ( 0 1 0 ) ( -1 0 0 ) ( 0 0 -1 ) cart. s(31) = ( 0.0000000 -1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 32 inv. 90 deg rotation - cart. axis [0,0,1] cryst. s(32) = ( 0 -1 0 ) ( 1 0 0 ) ( 0 0 -1 ) cart. s(32) = ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) isym = 33 inv. 180 deg rotation - cart. axis [1,0,1] cryst. s(33) = ( 0 0 -1 ) ( 0 1 0 ) ( -1 0 0 ) cart. s(33) = ( -0.0000000 0.0000000 -1.0000000 ) ( -0.0000000 1.0000000 -0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 34 inv. 180 deg rotation - cart. axis [-1,0,1] cryst. s(34) = ( 0 0 1 ) ( 0 1 0 ) ( 1 0 0 ) cart. s(34) = ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 35 inv. 90 deg rotation - cart. axis [0,1,0] cryst. s(35) = ( 0 0 1 ) ( 0 -1 0 ) ( -1 0 0 ) cart. s(35) = ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 36 inv. 90 deg rotation - cart. axis [0,-1,0] cryst. s(36) = ( 0 0 -1 ) ( 0 -1 0 ) ( 1 0 0 ) cart. s(36) = ( -0.0000000 -0.0000000 1.0000000 ) ( -0.0000000 -1.0000000 0.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 37 inv. 180 deg rotation - cart. axis [0,1,1] cryst. s(37) = ( 1 0 0 ) ( 0 0 -1 ) ( 0 -1 0 ) cart. s(37) = ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 38 inv. 180 deg rotation - cart. axis [0,1,-1] cryst. s(38) = ( 1 0 0 ) ( 0 0 1 ) ( 0 1 0 ) cart. s(38) = ( 1.0000000 0.0000000 0.0000000 ) ( 0.0000000 0.0000000 1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 39 inv. 90 deg rotation - cart. axis [-1,0,0] cryst. s(39) = ( -1 0 0 ) ( 0 0 1 ) ( 0 -1 0 ) cart. s(39) = ( -1.0000000 0.0000000 -0.0000000 ) ( 0.0000000 0.0000000 -1.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 40 inv. 90 deg rotation - cart. axis [1,0,0] cryst. s(40) = ( -1 0 0 ) ( 0 0 -1 ) ( 0 1 0 ) cart. s(40) = ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 41 inv. 120 deg rotation - cart. axis [-1,-1,-1] cryst. s(41) = ( 0 0 -1 ) ( -1 0 0 ) ( 0 -1 0 ) cart. s(41) = ( -0.0000000 -1.0000000 -0.0000000 ) ( -0.0000000 0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 42 inv. 120 deg rotation - cart. axis [-1,1,1] cryst. s(42) = ( 0 0 1 ) ( 1 0 0 ) ( 0 -1 0 ) cart. s(42) = ( 0.0000000 1.0000000 -0.0000000 ) ( 0.0000000 -0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 43 inv. 120 deg rotation - cart. axis [1,1,-1] cryst. s(43) = ( 0 0 1 ) ( -1 0 0 ) ( 0 1 0 ) cart. s(43) = ( 0.0000000 -1.0000000 0.0000000 ) ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) isym = 44 inv. 120 deg rotation - cart. axis [1,-1,1] cryst. s(44) = ( 0 0 -1 ) ( 1 0 0 ) ( 0 1 0 ) cart. s(44) = ( -0.0000000 1.0000000 0.0000000 ) ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) isym = 45 inv. 120 deg rotation - cart. axis [1,1,1] cryst. s(45) = ( 0 -1 0 ) ( 0 0 -1 ) ( -1 0 0 ) cart. s(45) = ( -0.0000000 -0.0000000 -1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 46 inv. 120 deg rotation - cart. axis [-1,1,-1] cryst. s(46) = ( 0 1 0 ) ( 0 0 1 ) ( -1 0 0 ) cart. s(46) = ( 0.0000000 0.0000000 -1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) isym = 47 inv. 120 deg rotation - cart. axis [1,-1,-1] cryst. s(47) = ( 0 1 0 ) ( 0 0 -1 ) ( 1 0 0 ) cart. s(47) = ( 0.0000000 -0.0000000 1.0000000 ) ( 1.0000000 -0.0000000 -0.0000000 ) ( 0.0000000 -1.0000000 0.0000000 ) isym = 48 inv. 120 deg rotation - cart. axis [-1,-1,1] cryst. s(48) = ( 0 -1 0 ) ( 0 0 1 ) ( 1 0 0 ) cart. s(48) = ( -0.0000000 0.0000000 1.0000000 ) ( -1.0000000 0.0000000 0.0000000 ) ( 0.0000000 1.0000000 -0.0000000 ) double point group O_h (m-3m) there are 16 classes and 6 irreducible representations the character table: E -E 8C3 -8C3 3C2 6C4 -6C4 6C2' i -i 8S6 -8S6 -3C2 -6C2' G_6+ 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 2.00 -2.00 1.00 -1.00 G_7+ 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 2.00 -2.00 1.00 -1.00 G_8+ 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 4.00 -4.00 -1.00 1.00 G_6- 2.00 -2.00 1.00 -1.00 0.00 1.41 -1.41 0.00 -2.00 2.00 -1.00 1.00 G_7- 2.00 -2.00 1.00 -1.00 0.00 -1.41 1.41 0.00 -2.00 2.00 -1.00 1.00 G_8- 4.00 -4.00 -1.00 1.00 0.00 0.00 0.00 0.00 -4.00 4.00 1.00 -1.00 3s_h 6S4 -6S4 6s_d -3s_h -6s_d G_6+ 0.00 1.41 -1.41 0.00 G_7+ 0.00 -1.41 1.41 0.00 G_8+ 0.00 0.00 0.00 0.00 G_6- 0.00 -1.41 1.41 0.00 G_7- 0.00 1.41 -1.41 0.00 G_8- 0.00 0.00 0.00 0.00 the symmetry operations in each class and the name of the first element: E 1 identity -E -1 identity E 8C3 17 19 20 18 24 21 22 23 120 deg rotation - cart. axis [-1,-1,-1] -8C3 -17 -19 -20 -18 -24 -21 -22 -23 120 deg rotation - cart. axis [-1,-1,-1] E 3C2 -3C2 2 -2 4 -4 3 -3 180 deg rotation - cart. axis [0,0,1] 6C4 7 8 15 16 12 11 90 deg rotation - cart. axis [0,0,-1] -6C4 -7 -8 -15 -16 -12 -11 90 deg rotation - cart. axis [0,0,-1] E 6C2'-6C2' 5 -5 6 -6 14 -13 -14 13 -10 -9 10 9 180 deg rotation - cart. axis [1,1,0] i 25 inversion -i -25 inversion E 8S6 41 43 44 42 48 45 46 47 inv. 120 deg rotation - cart. axis [-1,-1,-1] -8S6 -41 -43 -44 -42 -48 -45 -46 -47 inv. 120 deg rotation - cart. axis [-1,-1,-1] E 3s_h-3s_h 26 -26 28 -28 27 -27 inv. 180 deg rotation - cart. axis [0,0,1] 6S4 31 32 39 40 36 35 inv. 90 deg rotation - cart. axis [0,0,-1] -6S4 -31 -32 -39 -40 -36 -35 inv. 90 deg rotation - cart. axis [0,0,-1] E 6s_d-6s_d 29 -29 30 -30 38 -37 -38 37 -34 -33 34 33 inv. 180 deg rotation - cart. axis [1,1,0] Cartesian axes number of k points= 10 Fermi-Dirac smearing, width (Ry)= 0.0010 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 cryst. coord. k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0080000 k( 2) = ( 0.0000000 0.0000000 0.2000000), wk = 0.0480000 k( 3) = ( 0.0000000 0.0000000 0.4000000), wk = 0.0480000 k( 4) = ( 0.0000000 0.2000000 0.2000000), wk = 0.0960000 k( 5) = ( 0.0000000 0.2000000 0.4000000), wk = 0.1920000 k( 6) = ( 0.0000000 0.4000000 0.4000000), wk = 0.0960000 k( 7) = ( 0.2000000 0.2000000 0.2000000), wk = 0.0640000 k( 8) = ( 0.2000000 0.2000000 0.4000000), wk = 0.1920000 k( 9) = ( 0.2000000 0.4000000 0.4000000), wk = 0.1920000 k( 10) = ( 0.4000000 0.4000000 0.4000000), wk = 0.0640000 Dense grid: 52635 G-vectors FFT dimensions: ( 48, 48, 48) Smooth grid: 24837 G-vectors FFT dimensions: ( 40, 40, 40) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.12 Mb ( 178, 44) NL pseudopotentials 0.18 Mb ( 89, 136) Each V/rho on FFT grid 0.07 Mb ( 4608) Each G-vector array 0.01 Mb ( 1463) G-vector shells 0.00 Mb ( 345) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.48 Mb ( 178, 176) Each subspace H/S matrix 0.03 Mb ( 44, 44) Each matrix 0.18 Mb ( 136, 2, 44) 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 35.98553, renormalised to 36.00000 Starting wfc are 86 randomized atomic wfcs total cpu time spent up to now is 2.8 secs per-process dynamical memory: 34.9 Mb Self-consistent Calculation iteration # 1 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 2.0 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 8.07E-05, avg # of iterations = 3.0 total cpu time spent up to now is 4.4 secs total energy = -255.60401117 Ry Harris-Foulkes estimate = -255.61879241 Ry estimated scf accuracy < 0.03550515 Ry iteration # 2 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.86E-05, avg # of iterations = 2.0 total cpu time spent up to now is 5.3 secs total energy = -255.60870727 Ry Harris-Foulkes estimate = -255.61078869 Ry estimated scf accuracy < 0.00391124 Ry iteration # 3 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.09E-05, avg # of iterations = 3.1 total cpu time spent up to now is 6.2 secs total energy = -255.60960599 Ry Harris-Foulkes estimate = -255.61040246 Ry estimated scf accuracy < 0.00173661 Ry iteration # 4 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.82E-06, avg # of iterations = 2.4 total cpu time spent up to now is 7.0 secs total energy = -255.60998689 Ry Harris-Foulkes estimate = -255.60999789 Ry estimated scf accuracy < 0.00005939 Ry iteration # 5 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.65E-07, avg # of iterations = 3.4 total cpu time spent up to now is 8.0 secs total energy = -255.61000136 Ry Harris-Foulkes estimate = -255.61000175 Ry estimated scf accuracy < 0.00000093 Ry iteration # 6 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.59E-09, avg # of iterations = 4.2 total cpu time spent up to now is 9.0 secs total energy = -255.61000189 Ry Harris-Foulkes estimate = -255.61000190 Ry estimated scf accuracy < 0.00000003 Ry iteration # 7 ecut= 39.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.95E-11, avg # of iterations = 4.0 total cpu time spent up to now is 10.2 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 3071 PWs) bands (ev): -33.0961 -33.0961 -32.9745 -32.9745 -32.9745 -32.9745 -14.1645 -14.1645 -14.1645 -14.1645 -13.8241 -13.8241 -12.9025 -12.9025 -12.7183 -12.7183 -12.7183 -12.7183 -12.5322 -12.5322 -12.3193 -12.3193 -12.3193 -12.3193 3.3343 3.3343 6.7648 6.7648 8.4768 8.4768 8.4768 8.4768 8.9742 8.9742 8.9792 8.9792 8.9792 8.9792 10.2243 10.2243 10.2243 10.2243 10.2458 10.2458 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2000 ( 3086 PWs) bands (ev): -33.0862 -33.0862 -32.9846 -32.9846 -32.9743 -32.9743 -14.1606 -14.1606 -14.1335 -14.1335 -13.8505 -13.8505 -12.8734 -12.8734 -12.7140 -12.7140 -12.6992 -12.6992 -12.5429 -12.5429 -12.3759 -12.3759 -12.3393 -12.3393 3.6320 3.6320 6.9144 6.9144 7.9452 7.9452 8.8291 8.8291 8.9009 8.9009 8.9150 8.9150 9.1235 9.1235 10.2493 10.2493 10.2801 10.2801 10.4095 10.4095 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4000 ( 3086 PWs) bands (ev): -33.0625 -33.0625 -33.0084 -33.0084 -32.9740 -32.9740 -14.1544 -14.1544 -14.0598 -14.0598 -13.9171 -13.9171 -12.8290 -12.8290 -12.6607 -12.6607 -12.6494 -12.6494 -12.5784 -12.5784 -12.5107 -12.5107 -12.3724 -12.3724 4.3626 4.3626 6.2717 6.2717 8.0010 8.0010 8.7430 8.7430 8.7738 8.7738 9.4386 9.4386 9.6469 9.6469 10.5387 10.5387 10.5963 10.5963 11.0882 11.0883 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9999 0.9999 0.0016 0.0016 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.2000 ( 3097 PWs) bands (ev): -33.0768 -33.0768 -32.9884 -32.9884 -32.9798 -32.9798 -14.1450 -14.1450 -14.1229 -14.1229 -13.8771 -13.8771 -12.8484 -12.8484 -12.7134 -12.7134 -12.6889 -12.6889 -12.5368 -12.5368 -12.4278 -12.4278 -12.3499 -12.3499 3.9052 3.9052 7.0659 7.0659 7.9300 7.9300 8.3418 8.3418 8.7583 8.7583 9.0640 9.0640 9.4652 9.4652 10.0514 10.0514 10.4038 10.4038 10.6514 10.6514 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9990 0.9990 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.4000 ( 3105 PWs) bands (ev): -33.0545 -33.0545 -33.0085 -33.0085 -32.9817 -32.9817 -14.1368 -14.1368 -14.0672 -14.0672 -13.9420 -13.9420 -12.8109 -12.8109 -12.6983 -12.6983 -12.6478 -12.6478 -12.5537 -12.5537 -12.4980 -12.4980 -12.3892 -12.3892 4.5566 4.5566 6.5263 6.5263 7.9809 7.9809 8.3788 8.3788 8.6476 8.6476 9.6399 9.6399 9.6530 9.6530 10.3716 10.3716 10.4026 10.4027 11.1171 11.1172 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0026 0.0026 0.0010 0.0010 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4000 0.4000 ( 3116 PWs) bands (ev): -33.0348 -33.0348 -33.0109 -33.0109 -32.9988 -32.9988 -14.0962 -14.0962 -14.0736 -14.0736 -13.9998 -13.9998 -12.7733 -12.7733 -12.7371 -12.7371 -12.6284 -12.6284 -12.5836 -12.5836 -12.4539 -12.4539 -12.4170 -12.4170 4.9724 4.9724 7.0119 7.0119 7.2607 7.2607 8.4956 8.4956 8.5778 8.5778 9.6903 9.6903 9.8336 9.8336 10.0876 10.0877 10.5632 10.5632 11.1069 11.1070 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.2000 ( 3101 PWs) bands (ev): -33.0683 -33.0683 -32.9882 -32.9882 -32.9882 -32.9882 -14.1353 -14.1353 -14.1150 -14.1150 -13.9042 -13.9042 -12.8287 -12.8287 -12.7298 -12.7298 -12.6681 -12.6681 -12.5227 -12.5227 -12.4736 -12.4736 -12.3483 -12.3483 4.1466 4.1466 7.0970 7.0970 8.1703 8.1703 8.1936 8.1936 8.3463 8.3463 9.4226 9.4226 9.4284 9.4284 9.9741 9.9741 10.6139 10.6139 10.6307 10.6307 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9999 0.9999 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.4000 ( 3111 PWs) bands (ev): -33.0478 -33.0478 -33.0086 -33.0086 -32.9881 -32.9881 -14.1294 -14.1294 -14.0735 -14.0735 -13.9673 -13.9673 -12.8027 -12.8027 -12.7376 -12.7376 -12.6271 -12.6271 -12.5573 -12.5573 -12.4770 -12.4770 -12.3784 -12.3784 4.6917 4.6917 6.7051 6.7051 7.7754 7.7754 8.5376 8.5376 8.6423 8.6423 9.4947 9.4947 9.6866 9.6866 10.1537 10.1537 10.4564 10.4564 10.6930 10.6930 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9912 0.9912 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.4000 0.4000 ( 3100 PWs) bands (ev): -33.0311 -33.0311 -33.0107 -33.0107 -33.0020 -33.0020 -14.1102 -14.1102 -14.0768 -14.0768 -14.0220 -14.0220 -12.8006 -12.8006 -12.7675 -12.7675 -12.6041 -12.6041 -12.5756 -12.5756 -12.4251 -12.4251 -12.3789 -12.3789 4.9374 4.9374 7.0661 7.0661 7.4876 7.4876 8.4698 8.4698 8.7372 8.7372 9.3153 9.3153 9.5999 9.5999 9.9876 9.9876 10.8145 10.8145 10.8835 10.8835 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0470 0.0470 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4000 0.4000 0.4000 ( 3112 PWs) bands (ev): -33.0222 -33.0222 -33.0106 -33.0106 -33.0106 -33.0106 -14.1304 -14.1304 -14.0802 -14.0802 -14.0604 -14.0604 -12.8545 -12.8545 -12.7789 -12.7789 -12.5775 -12.5775 -12.5688 -12.5688 -12.3659 -12.3659 -12.3402 -12.3402 4.8612 4.8612 7.4051 7.4051 8.0096 8.0096 8.0273 8.0273 8.9604 8.9604 8.9844 8.9844 9.0701 9.0701 10.2925 10.2925 11.0561 11.0561 11.0790 11.0791 occupation numbers 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 9.5589 ev ! total energy = -255.61000190 Ry Harris-Foulkes estimate = -255.61000191 Ry estimated scf accuracy < 5.5E-09 Ry The total energy is the sum of the following terms: one-electron contribution = -68.10004306 Ry hartree contribution = 43.72187269 Ry xc contribution = -56.18756308 Ry ewald contribution = -175.04416290 Ry smearing contrib. (-TS) = -0.00010555 Ry convergence has been achieved in 7 iterations Writing output data file Y3Al.save init_run : 0.66s CPU 0.73s WALL ( 1 calls) electrons : 7.21s CPU 7.51s WALL ( 1 calls) Called by init_run: wfcinit : 0.40s CPU 0.43s WALL ( 1 calls) potinit : 0.01s CPU 0.01s WALL ( 1 calls) Called by electrons: c_bands : 5.44s CPU 5.59s WALL ( 8 calls) sum_band : 1.28s CPU 1.30s WALL ( 8 calls) v_of_rho : 0.03s CPU 0.03s WALL ( 8 calls) v_h : 0.01s CPU 0.00s WALL ( 8 calls) v_xc : 0.02s CPU 0.02s WALL ( 8 calls) newd : 0.52s CPU 0.53s WALL ( 8 calls) mix_rho : 0.01s CPU 0.01s WALL ( 8 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.01s WALL ( 170 calls) cegterg : 5.18s CPU 5.27s WALL ( 80 calls) Called by sum_band: sum_band:bec : 0.32s CPU 0.33s WALL ( 80 calls) addusdens : 0.27s CPU 0.27s WALL ( 8 calls) Called by *egterg: h_psi : 3.54s CPU 3.66s WALL ( 331 calls) s_psi : 0.22s CPU 0.24s WALL ( 331 calls) g_psi : 0.00s CPU 0.01s WALL ( 241 calls) cdiaghg : 1.36s CPU 1.36s WALL ( 311 calls) cegterg:over : 0.09s CPU 0.12s WALL ( 241 calls) cegterg:upda : 0.11s CPU 0.10s WALL ( 241 calls) cegterg:last : 0.05s CPU 0.04s WALL ( 80 calls) cdiaghg:chol : 0.10s CPU 0.08s WALL ( 311 calls) cdiaghg:inve : 0.06s CPU 0.04s WALL ( 311 calls) cdiaghg:para : 0.11s CPU 0.08s WALL ( 622 calls) Called by h_psi: h_psi:vloc : 3.00s CPU 3.14s WALL ( 331 calls) h_psi:vnl : 0.52s CPU 0.51s WALL ( 331 calls) add_vuspsi : 0.28s CPU 0.27s WALL ( 331 calls) General routines calbec : 0.32s CPU 0.30s WALL ( 411 calls) fft : 0.06s CPU 0.05s WALL ( 242 calls) ffts : 0.00s CPU 0.01s WALL ( 64 calls) fftw : 3.27s CPU 3.47s WALL ( 46728 calls) interpolate : 0.02s CPU 0.03s WALL ( 64 calls) Parallel routines fft_scatter : 1.31s CPU 1.33s WALL ( 47034 calls) PWSCF : 9.91s CPU 10.93s WALL This run was terminated on: 19:41:13 24Dec2016 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=