[0] DAPL startup(): trying to open DAPL provider from I_MPI_DAPL_PROVIDER: ofa-v2-mlx4_0-1 [0] MPI startup(): DAPL provider ofa-v2-mlx4_0-1 [0] MPI startup(): shm and dapl data transfer modes Message from routine get_command_line: unexpected argument # 2 :-inp Program PWSCF v.5.1.1 starts on 14Jul2015 at 15: 0:15 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 1 processors Waiting for input... Reading input from standard input Warning: card &CELL ignored Warning: card CELL_DYNAMICS ='BFGS' ignored Warning: card WMASS = 0.004, ignored Warning: card / ignored Current dimensions of program PWSCF are: Max number of different atomic species (ntypx) = 10 Max number of k-points (npk) = 40000 Max angular momentum in pseudopotentials (lmaxx) = 3 Subspace diagonalization in iterative solution of the eigenvalue problem: a serial algorithm will be used G-vector sticks info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Sum 5871 5871 1591 346487 346487 48877 bravais-lattice index = 0 lattice parameter (alat) = 1.8897 a.u. unit-cell volume = 5517.4909 (a.u.)^3 number of atoms/cell = 72 number of atomic types = 3 number of electrons = 404.00 number of Kohn-Sham states= 242 kinetic-energy cutoff = 60.0000 Ry charge density cutoff = 240.0000 Ry convergence threshold = 1.0E-02 mixing beta = 0.7000 number of iterations used = 4 plain mixing Exchange-correlation = LDA ( 1 1 0 0 0 0) nstep = 1000 celldm(1)= 1.889726 celldm(2)= 0.000000 celldm(3)= 0.000000 celldm(4)= 0.000000 celldm(5)= 0.000000 celldm(6)= 0.000000 crystal axes: (cart. coord. in units of alat) a(1) = ( 7.569000 0.000000 0.000000 ) a(2) = ( 0.000000 11.353500 0.000000 ) a(3) = ( 0.000000 0.000000 9.514300 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 0.132118 0.000000 0.000000 ) b(2) = ( 0.000000 0.088079 0.000000 ) b(3) = ( 0.000000 0.000000 0.105105 ) PseudoPot. # 1 for Ti read from file: /export/home/wmulwa/scratch5/RESEARCH/PAW/Ti.pz-n-nc.UPF MD5 check sum: 7d96eb281fc152189740e93005b9c377 Pseudo is Norm-conserving + core correction, Zval = 4.0 Generated using "atomic" code by A. Dal Corso v.5.0.99 svn rev. 10869 Using radial grid of 1177 points, 2 beta functions with: l(1) = 1 l(2) = 2 PseudoPot. # 2 for O read from file: /export/home/wmulwa/scratch5/RESEARCH/PAW/O.pz-n-nc.UPF MD5 check sum: 1f9736618f826250aad27c6ecbad5d4d Pseudo is Norm-conserving, Zval = 6.0 Generated using "atomic" code by A. Dal Corso v.5.0.99 svn rev. 10869 Using radial grid of 1095 points, 1 beta functions with: l(1) = 1 PseudoPot. # 3 for Er read from file: /export/home/wmulwa/scratch5/RESEARCH/PAW/Er.pz-n-nc.UPF MD5 check sum: 222954db37a7f3c2a88a493f31090e26 Pseudo is Norm-conserving + core correction, Zval = 14.0 Generated using "atomic" code by A. Dal Corso v.5.0.99 svn rev. 10869 Using radial grid of 1267 points, 3 beta functions with: l(1) = 1 l(2) = 2 l(3) = 3 atomic species valence mass pseudopotential Ti 4.00 47.86700 Ti( 1.00) O 6.00 15.99940 O( 1.00) Er 14.00 168.93420 Er( 1.00) Simplified LDA+U calculation (l_max = 2) with parameters (eV): atomic species L U alpha J0 beta Ti 2 10.0000 0.0000 0.0000 0.0000 O 1 5.0000 0.0000 0.0000 0.0000 No symmetry found Cartesian axes site n. atom positions (alat units) 1 Ti tau( 1) = ( -0.0079167 -0.0166166 -0.0162763 ) 2 O tau( 2) = ( 0.0004923 -0.0072591 1.9387376 ) 3 O tau( 3) = ( 1.8937752 1.8737803 2.6538090 ) 4 O tau( 4) = ( 1.8871220 -0.0309827 9.1068904 ) 5 O tau( 5) = ( 1.8899683 -0.0966719 5.2208269 ) 6 Ti tau( 6) = ( 1.8842798 -0.0265774 7.1684321 ) 7 O tau( 7) = ( 1.8829039 1.8866006 6.8343755 ) 8 Er tau( 8) = ( 1.8953666 1.8928920 4.7467126 ) 9 O tau( 9) = ( -0.0038087 -0.0121536 7.5345397 ) 10 O tau( 10) = ( -0.0807362 1.9022568 4.2437149 ) 11 O tau( 11) = ( -0.0078257 1.8913343 0.3560987 ) 12 Ti tau( 12) = ( -0.0133771 1.8863570 2.2959407 ) 13 Ti tau( 13) = ( 3.7925486 7.5938188 -0.0321625 ) 14 O tau( 14) = ( 3.7864584 7.6645130 1.9150520 ) 15 O tau( 15) = ( 5.6776212 9.4827504 2.7881184 ) 16 O tau( 16) = ( 5.6805295 7.5812657 9.1143907 ) 17 O tau( 17) = ( 5.6761296 7.5763724 5.1957337 ) 18 Ti tau( 18) = ( 5.6845210 7.5858798 7.1511627 ) 19 O tau( 19) = ( 5.6853757 9.4684379 6.7186277 ) 20 Ti tau( 20) = ( 5.6753289 9.4579383 4.7507415 ) 21 O tau( 21) = ( 3.7896398 7.5998743 7.5423000 ) 22 O tau( 22) = ( 3.7850354 9.4468616 4.3327561 ) 23 O tau( 23) = ( 3.7928310 9.4765597 0.3683840 ) 24 Ti tau( 24) = ( 3.7834085 9.4999624 2.3647485 ) 25 Ti tau( 25) = ( -0.0065117 7.5672160 0.0263758 ) 26 O tau( 26) = ( -0.0000557 7.5729192 1.9917384 ) 27 O tau( 27) = ( 1.8915939 9.4759153 2.8023974 ) 28 O tau( 28) = ( 1.8896256 7.5811417 9.1372889 ) 29 O tau( 29) = ( 1.8916002 7.5785201 5.2217898 ) 30 Ti tau( 30) = ( 1.8843037 7.5871549 7.1922409 ) 31 O tau( 31) = ( 1.8835730 9.4455952 6.7664482 ) 32 Ti tau( 32) = ( 1.8933346 9.4221781 4.7696519 ) 33 O tau( 33) = ( -0.0049752 7.5776334 7.5777243 ) 34 O tau( 34) = ( -0.0009056 9.4399907 4.3469454 ) 35 O tau( 35) = ( -0.0081713 9.4541209 0.4165785 ) 36 Ti tau( 36) = ( 0.0012386 9.4649489 2.3844439 ) 37 Ti tau( 37) = ( 3.7942263 3.7730641 -0.0537437 ) 38 O tau( 38) = ( 3.7929350 3.6954290 1.8772653 ) 39 O tau( 39) = ( 5.7585079 5.6669355 2.8917114 ) 40 O tau( 40) = ( 5.6804532 3.7789646 9.0922920 ) 41 O tau( 41) = ( 5.6789851 3.7836184 5.1742885 ) 42 Ti tau( 42) = ( 5.6842224 3.7670695 7.1278855 ) 43 O tau( 43) = ( 5.6846775 5.6776234 6.7787207 ) 44 Ti tau( 44) = ( 5.6917174 5.6825442 4.8403488 ) 45 O tau( 45) = ( 3.7871553 3.7670737 7.5189957 ) 46 O tau( 46) = ( 3.7829439 5.6949292 4.4870739 ) 47 O tau( 47) = ( 3.7937639 5.6822769 0.3066756 ) 48 Er tau( 48) = ( 3.7814686 5.6758918 2.3977336 ) 49 Ti tau( 49) = ( -0.0077347 3.8014987 0.0070369 ) 50 O tau( 50) = ( -0.0020658 3.7853974 1.9604559 ) 51 O tau( 51) = ( 1.8081492 5.6702382 2.9100388 ) 52 O tau( 52) = ( 1.8896421 3.8018219 9.1301374 ) 53 O tau( 53) = ( 1.8834252 3.8745174 5.2588698 ) 54 Ti tau( 54) = ( 1.8828143 3.7972731 7.1903330 ) 55 O tau( 55) = ( 1.8843503 5.7029946 6.8247333 ) 56 Ti tau( 56) = ( 1.8743315 5.7192807 4.8634133 ) 57 O tau( 57) = ( -0.0037676 3.7904436 7.5569067 ) 58 O tau( 58) = ( -0.0023180 5.7004436 4.4460036 ) 59 O tau( 59) = ( -0.0082843 5.6832726 0.4507572 ) 60 Ti tau( 60) = ( 0.0007606 5.6782883 2.4300225 ) 61 Ti tau( 61) = ( 3.7926879 -0.0179595 -0.0569045 ) 62 O tau( 62) = ( 3.7852601 -0.0094843 1.9127335 ) 63 O tau( 63) = ( 5.6791212 1.8683838 2.6886321 ) 64 O tau( 64) = ( 5.6817874 -0.0090977 9.0718423 ) 65 O tau( 65) = ( 5.6768929 -0.0038714 5.1436263 ) 66 Ti tau( 66) = ( 5.6832091 0.0020616 7.1087364 ) 67 O tau( 67) = ( 5.6851485 1.8857825 6.6836630 ) 68 Ti tau( 68) = ( 5.6759931 1.8906217 4.7036619 ) 69 O tau( 69) = ( 3.7871926 -0.0122676 7.5117632 ) 70 O tau( 70) = ( 3.8676789 1.8990267 4.2259094 ) 71 O tau( 71) = ( 3.7924550 1.8663707 0.3104066 ) 72 Ti tau( 72) = ( 3.8008923 1.8500145 2.2725961 ) number of k points= 8 Marzari-Vanderbilt smearing, width (Ry)= 0.0000 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.2500000 k( 2) = ( 0.0000000 0.0000000 -0.0525525), wk = 0.2500000 k( 3) = ( 0.0000000 -0.0440393 0.0000000), wk = 0.2500000 k( 4) = ( 0.0000000 -0.0440393 -0.0525525), wk = 0.2500000 k( 5) = ( -0.0660589 0.0000000 0.0000000), wk = 0.2500000 k( 6) = ( -0.0660589 0.0000000 -0.0525525), wk = 0.2500000 k( 7) = ( -0.0660589 -0.0440393 0.0000000), wk = 0.2500000 k( 8) = ( -0.0660589 -0.0440393 -0.0525525), wk = 0.2500000 Dense grid: 346487 G-vectors FFT dimensions: ( 72, 108, 90) Largest allocated arrays est. size (Mb) dimensions Kohn-Sham Wavefunctions 160.22 Mb ( 43388, 242) Atomic Hubbard wavefuncts 168.16 Mb ( 43388, 254) NL pseudopotentials 231.72 Mb ( 43388, 350) Each V/rho on FFT grid 10.68 Mb ( 699840) Each G-vector array 2.64 Mb ( 346487) G-vector shells 0.25 Mb ( 33370) Largest temporary arrays est. size (Mb) dimensions Auxiliary wavefunctions 640.86 Mb ( 43388, 968) Each subspace H/S matrix 14.30 Mb ( 968, 968) Each matrix 1.29 Mb ( 350, 242) Arrays for rho mixing 42.71 Mb ( 699840, 4) Check: negative/imaginary core charge= -0.000069 0.000000 Initial potential from superposition of free atoms starting charge 403.96532, renormalised to 404.00000 Number of +U iterations with fixed ns = 0 Starting occupations: --- enter write_ns --- LDA+U parameters: U( 1) = 10.00000000 alpha( 1) = 0.00000000 U( 2) = 5.00000000 alpha( 2) = 0.00000000 atom 1 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 2 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 3 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 4 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 5 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 6 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 7 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 9 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 10 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 11 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 12 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 13 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 14 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 15 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 16 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 17 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 18 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 19 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 20 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 21 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 22 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 23 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 24 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 25 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 26 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 27 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 28 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 29 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 30 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 31 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 32 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 33 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 34 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 35 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 36 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 37 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 38 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 39 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 40 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 41 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 42 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 43 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 44 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 45 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 46 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 47 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 49 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 50 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 51 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 52 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 53 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 54 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 55 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 56 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 57 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 58 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 59 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 60 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 61 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 62 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 63 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 64 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 65 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 66 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 67 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 68 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 atom 69 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 70 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 71 Tr[ns(na)] = 4.00000 eigenvalues: 0.667 0.667 0.667 eigenvectors: 1.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 1.000 occupations: 0.667 0.000 0.000 0.000 0.667 0.000 0.000 0.000 0.667 atom 72 Tr[ns(na)] = 2.00000 eigenvalues: 0.200 0.200 0.200 0.200 0.200 eigenvectors: 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 occupations: 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 0.000 0.000 0.000 0.000 0.000 0.200 N of occupied +U levels = 236.000000 --- exit write_ns --- Atomic wfc used for LDA+U Projector are NOT orthogonalized Starting wfc are 422 atomic wfcs total cpu time spent up to now is 639.2 secs Self-consistent Calculation iteration # 1 ecut= 60.00 Ry beta=0.70 Davidson diagonalization with overlap APPLICATION TERMINATED WITH THE EXIT STRING: Killed (signal 9)