Program PWSCF v.5.3.0 (svn rev. 11974) starts on 5Jan2017 at 16:47:44 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 61 30 9 2171 758 119 Max 62 31 10 2174 777 123 Sum 2217 1101 325 78205 27633 4337 bravais-lattice index = 14 lattice parameter (alat) = 9.2181 a.u. unit-cell volume = 783.2888 (a.u.)^3 number of atoms/cell = 5 number of atomic types = 3 number of electrons = 40.00 number of Kohn-Sham states= 48 kinetic-energy cutoff = 41.0000 Ry charge density cutoff = 327.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.218083 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 C read from file: /users/gautes/Pseudo/C.rel-pbe-n-rrkjus_psl.1.0.0.UPF MD5 check sum: de79cff2a1a990998107dfc0a7a38bf3 Pseudo is Ultrasoft + core correction, Zval = 4.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9656 Using radial grid of 1073 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. # 3 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) C 4.00 12.01070 C( 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: 78205 G-vectors FFT dimensions: ( 54, 54, 54) Smooth grid: 27633 G-vectors FFT dimensions: ( 40, 40, 40) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 0.14 Mb ( 196, 48) NL pseudopotentials 0.22 Mb ( 98, 150) Each V/rho on FFT grid 0.09 Mb ( 5832) Each G-vector array 0.02 Mb ( 2174) G-vector shells 0.00 Mb ( 429) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 0.57 Mb ( 196, 192) Each subspace H/S matrix 0.04 Mb ( 48, 48) Each matrix 0.22 Mb ( 150, 2, 48) Arrays for rho mixing 0.71 Mb ( 5832, 8) Check: negative/imaginary core charge= -0.000003 0.000000 Initial potential from superposition of free atoms starting charge 39.98549, renormalised to 40.00000 Starting wfc are 94 randomized atomic wfcs total cpu time spent up to now is 3.6 secs per-process dynamical memory: 50.6 Mb Self-consistent Calculation iteration # 1 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 2.0 total cpu time spent up to now is 4.7 secs total energy = -267.32907941 Ry Harris-Foulkes estimate = -268.47851986 Ry estimated scf accuracy < 1.44898196 Ry iteration # 2 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.62E-03, avg # of iterations = 4.1 total cpu time spent up to now is 6.0 secs total energy = -267.32298670 Ry Harris-Foulkes estimate = -269.30840757 Ry estimated scf accuracy < 5.22267397 Ry iteration # 3 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.62E-03, avg # of iterations = 3.5 total cpu time spent up to now is 7.2 secs total energy = -268.14510247 Ry Harris-Foulkes estimate = -268.15597150 Ry estimated scf accuracy < 0.03989269 Ry iteration # 4 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.97E-05, avg # of iterations = 2.9 total cpu time spent up to now is 8.2 secs total energy = -268.14458414 Ry Harris-Foulkes estimate = -268.15044395 Ry estimated scf accuracy < 0.01437087 Ry iteration # 5 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.59E-05, avg # of iterations = 2.5 total cpu time spent up to now is 9.3 secs total energy = -268.14632200 Ry Harris-Foulkes estimate = -268.14637605 Ry estimated scf accuracy < 0.00018319 Ry iteration # 6 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.58E-07, avg # of iterations = 6.4 total cpu time spent up to now is 11.0 secs total energy = -268.14637961 Ry Harris-Foulkes estimate = -268.14644539 Ry estimated scf accuracy < 0.00013219 Ry iteration # 7 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.30E-07, avg # of iterations = 2.8 total cpu time spent up to now is 12.1 secs total energy = -268.14640105 Ry Harris-Foulkes estimate = -268.14640223 Ry estimated scf accuracy < 0.00000269 Ry iteration # 8 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.72E-09, avg # of iterations = 3.8 total cpu time spent up to now is 13.3 secs total energy = -268.14640187 Ry Harris-Foulkes estimate = -268.14640210 Ry estimated scf accuracy < 0.00000050 Ry iteration # 9 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.25E-09, avg # of iterations = 2.8 total cpu time spent up to now is 14.4 secs total energy = -268.14640193 Ry Harris-Foulkes estimate = -268.14640195 Ry estimated scf accuracy < 0.00000002 Ry iteration # 10 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.93E-11, avg # of iterations = 4.0 total cpu time spent up to now is 15.6 secs total energy = -268.14640194 Ry Harris-Foulkes estimate = -268.14640195 Ry estimated scf accuracy < 0.00000001 Ry iteration # 11 ecut= 41.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.59E-11, avg # of iterations = 2.7 total cpu time spent up to now is 16.7 secs End of self-consistent calculation k = 0.0000 0.0000 0.0000 ( 3431 PWs) bands (ev): -32.7667 -32.7667 -32.6556 -32.6556 -32.6556 -32.6556 -13.8562 -13.8562 -13.8562 -13.8562 -13.5510 -13.5510 -12.5464 -12.5464 -12.4482 -12.4482 -12.4482 -12.4482 -12.2853 -12.2853 -12.0213 -12.0213 -12.0213 -12.0213 -0.4678 -0.4678 3.9220 3.9220 8.7470 8.7470 8.7470 8.7470 8.9400 8.9400 8.9400 8.9400 9.0598 9.0598 9.2392 9.2392 9.2442 9.2442 9.2442 9.2442 11.1009 11.1009 11.1739 11.1739 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.0305 0.0305 0.0214 0.0214 0.0214 0.0214 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.2000 ( 3499 PWs) bands (ev): -32.7574 -32.7574 -32.6650 -32.6650 -32.6557 -32.6557 -13.8530 -13.8530 -13.8293 -13.8293 -13.5744 -13.5744 -12.5250 -12.5250 -12.4406 -12.4406 -12.4299 -12.4299 -12.2946 -12.2946 -12.0796 -12.0796 -12.0400 -12.0400 -0.2453 -0.2453 4.1510 4.1510 7.8250 7.8250 8.6840 8.6840 8.7267 8.7267 8.7712 8.7712 9.0870 9.0870 9.1838 9.1838 9.1919 9.1919 9.3830 9.3830 10.6823 10.6823 11.5177 11.5178 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.9996 0.9996 0.6501 0.6501 0.5051 0.5051 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.0000 0.4000 ( 3482 PWs) bands (ev): -32.7358 -32.7358 -32.6866 -32.6866 -32.6556 -32.6556 -13.8479 -13.8479 -13.7651 -13.7651 -13.6334 -13.6334 -12.5076 -12.5076 -12.3832 -12.3832 -12.3656 -12.3656 -12.3253 -12.3253 -12.2181 -12.2181 -12.0718 -12.0718 0.1670 0.1670 4.6283 4.6283 6.9028 6.9028 7.4921 7.4921 8.4417 8.4417 8.5211 8.5211 8.9788 8.9788 9.0080 9.0080 9.6843 9.6843 9.8724 9.8724 11.2221 11.2221 11.6569 11.6569 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.2000 ( 3485 PWs) bands (ev): -32.7489 -32.7489 -32.6684 -32.6684 -32.6607 -32.6607 -13.8408 -13.8408 -13.8195 -13.8195 -13.5985 -13.5985 -12.5135 -12.5135 -12.4442 -12.4442 -12.4090 -12.4090 -12.2879 -12.2879 -12.1345 -12.1345 -12.0440 -12.0440 -0.0257 -0.0257 4.3635 4.3635 7.6159 7.6159 8.0859 8.0859 8.3450 8.3450 8.4752 8.4752 8.7864 8.7864 9.0722 9.0722 9.5986 9.5986 9.7096 9.7096 11.2062 11.2062 11.3632 11.3632 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.9999 0.9999 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.2000 0.4000 ( 3474 PWs) bands (ev): -32.7285 -32.7285 -32.6867 -32.6867 -32.6625 -32.6625 -13.8342 -13.8342 -13.7719 -13.7719 -13.6568 -13.6568 -12.5050 -12.5050 -12.4280 -12.4280 -12.3611 -12.3611 -12.2782 -12.2782 -12.2241 -12.2241 -12.0740 -12.0740 0.3850 0.3850 4.7963 4.7963 6.8526 6.8526 7.2222 7.2222 7.9955 7.9955 8.2372 8.2372 8.6782 8.6782 8.9041 8.9041 10.1178 10.1178 10.1966 10.1966 11.2971 11.2971 11.4075 11.4075 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.0000 0.4000 0.4000 ( 3468 PWs) bands (ev): -32.7105 -32.7105 -32.6889 -32.6889 -32.6779 -32.6779 -13.8043 -13.8043 -13.7756 -13.7756 -13.7099 -13.7099 -12.4945 -12.4945 -12.4635 -12.4635 -12.3550 -12.3550 -12.2745 -12.2745 -12.2054 -12.2054 -12.0751 -12.0751 0.7998 0.7998 5.1617 5.1617 6.3647 6.3647 6.5202 6.5202 7.9788 7.9788 8.0665 8.0665 8.2725 8.2725 8.7621 8.7621 10.7510 10.7510 10.8440 10.8440 11.0828 11.0828 11.3873 11.3874 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.2000 ( 3476 PWs) bands (ev): -32.7411 -32.7411 -32.6683 -32.6683 -32.6683 -32.6683 -13.8362 -13.8362 -13.8105 -13.8105 -13.6235 -13.6235 -12.5171 -12.5171 -12.4472 -12.4472 -12.3904 -12.3904 -12.2723 -12.2723 -12.1723 -12.1723 -12.0418 -12.0418 0.1913 0.1913 4.5517 4.5517 7.5589 7.5589 7.6077 7.6077 7.7543 7.7543 8.4737 8.4737 8.7545 8.7545 8.7589 8.7589 9.9376 9.9376 9.9390 9.9390 11.4028 11.4028 11.6844 11.6844 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.2000 0.4000 ( 3481 PWs) bands (ev): -32.7225 -32.7225 -32.6868 -32.6868 -32.6682 -32.6682 -13.8320 -13.8320 -13.7766 -13.7766 -13.6806 -13.6806 -12.5209 -12.5209 -12.4570 -12.4570 -12.3450 -12.3450 -12.2698 -12.2698 -12.1956 -12.1956 -12.0650 -12.0650 0.6008 0.6008 4.9164 4.9164 6.7028 6.7028 7.0950 7.0950 7.4667 7.4667 8.1314 8.1314 8.5677 8.5677 8.8837 8.8837 10.4046 10.4046 10.4339 10.4339 11.3542 11.3542 11.5234 11.5234 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.2000 0.4000 0.4000 ( 3457 PWs) bands (ev): -32.7073 -32.7073 -32.6887 -32.6887 -32.6808 -32.6808 -13.8260 -13.8260 -13.7763 -13.7763 -13.7295 -13.7295 -12.5329 -12.5329 -12.5020 -12.5020 -12.3162 -12.3162 -12.2711 -12.2711 -12.1444 -12.1444 -12.0588 -12.0588 1.0159 1.0159 5.1478 5.1478 6.2956 6.2956 6.3769 6.3769 7.2121 7.2121 8.2445 8.2445 8.4268 8.4268 8.8355 8.8355 10.8922 10.8922 10.9124 10.9124 11.2316 11.2316 11.3979 11.3979 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 k = 0.4000 0.4000 0.4000 ( 3464 PWs) bands (ev): -32.6992 -32.6992 -32.6885 -32.6885 -32.6885 -32.6885 -13.8591 -13.8591 -13.7778 -13.7778 -13.7613 -13.7613 -12.5751 -12.5751 -12.5426 -12.5426 -12.2771 -12.2771 -12.2616 -12.2616 -12.0670 -12.0670 -12.0350 -12.0350 1.4411 1.4411 5.1092 5.1092 6.1645 6.1645 6.1726 6.1726 6.1954 6.1954 8.5685 8.5685 8.8387 8.8387 8.8821 8.8821 10.8644 10.8644 10.8890 10.8890 11.0033 11.0033 11.3497 11.3498 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 1.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 the Fermi energy is 9.1922 ev ! total energy = -268.14640195 Ry Harris-Foulkes estimate = -268.14640195 Ry estimated scf accuracy < 2.5E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -74.92783197 Ry hartree contribution = 50.28098839 Ry xc contribution = -60.84823231 Ry ewald contribution = -182.65119121 Ry smearing contrib. (-TS) = -0.00013485 Ry convergence has been achieved in 11 iterations Writing output data file Y3AlC.save init_run : 0.87s CPU 1.26s WALL ( 1 calls) electrons : 12.28s CPU 13.29s WALL ( 1 calls) Called by init_run: wfcinit : 0.47s CPU 0.65s WALL ( 1 calls) potinit : 0.02s CPU 0.02s WALL ( 1 calls) Called by electrons: c_bands : 9.06s CPU 9.66s WALL ( 11 calls) sum_band : 2.22s CPU 2.25s WALL ( 11 calls) v_of_rho : 0.05s CPU 0.05s WALL ( 12 calls) v_h : 0.00s CPU 0.00s WALL ( 12 calls) v_xc : 0.05s CPU 0.05s WALL ( 12 calls) newd : 0.96s CPU 0.98s WALL ( 12 calls) mix_rho : 0.03s CPU 0.03s WALL ( 11 calls) Called by c_bands: init_us_2 : 0.01s CPU 0.03s WALL ( 230 calls) cegterg : 8.63s CPU 8.75s WALL ( 110 calls) Called by sum_band: sum_band:bec : 0.46s CPU 0.47s WALL ( 110 calls) addusdens : 0.64s CPU 0.64s WALL ( 11 calls) Called by *egterg: h_psi : 5.61s CPU 5.70s WALL ( 495 calls) s_psi : 0.33s CPU 0.37s WALL ( 495 calls) g_psi : 0.00s CPU 0.01s WALL ( 375 calls) cdiaghg : 2.42s CPU 2.41s WALL ( 485 calls) cegterg:over : 0.20s CPU 0.20s WALL ( 375 calls) cegterg:upda : 0.20s CPU 0.18s WALL ( 375 calls) cegterg:last : 0.08s CPU 0.06s WALL ( 110 calls) cdiaghg:chol : 0.12s CPU 0.14s WALL ( 485 calls) cdiaghg:inve : 0.11s CPU 0.09s WALL ( 485 calls) cdiaghg:para : 0.16s CPU 0.14s WALL ( 970 calls) Called by h_psi: h_psi:vloc : 4.73s CPU 4.84s WALL ( 495 calls) h_psi:vnl : 0.86s CPU 0.86s WALL ( 495 calls) add_vuspsi : 0.46s CPU 0.44s WALL ( 495 calls) General routines calbec : 0.52s CPU 0.54s WALL ( 605 calls) fft : 0.14s CPU 0.12s WALL ( 356 calls) ffts : 0.00s CPU 0.02s WALL ( 92 calls) fftw : 5.26s CPU 5.35s WALL ( 67896 calls) interpolate : 0.04s CPU 0.05s WALL ( 92 calls) Parallel routines fft_scatter : 1.98s CPU 2.11s WALL ( 68344 calls) PWSCF : 15.58s CPU 18.89s WALL This run was terminated on: 16:48: 3 5Jan2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=