Program PWSCF v.6.0 (svn rev. 13079) starts on 4Jun2020 at 4:28: 0 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 / 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 file Fe.pz-nd-rrkjus.UPF: wavefunction(s) 4S renormalized file Fe.pz-nd-rrkjus.UPF: wavefunction(s) 4S renormalized 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 3115 3115 847 279295 279295 39555 bravais-lattice index = 0 lattice parameter (alat) = 11.2415 a.u. unit-cell volume = 2424.0806 (a.u.)^3 number of atoms/cell = 24 number of atomic types = 3 number of electrons = 168.00 number of Kohn-Sham states= 101 kinetic-energy cutoff = 90.0000 Ry charge density cutoff = 360.0000 Ry convergence threshold = 1.0E-06 mixing beta = 0.7000 number of iterations used = 8 plain mixing Exchange-correlation = LDA ( 1 1 0 0 0 0) nstep = 50 celldm(1)= 11.241507 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) = ( 1.000000 0.000000 0.000000 ) a(2) = ( -0.500000 0.866025 0.000000 ) a(3) = ( 0.000000 0.000000 1.970348 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 1.000000 0.577350 0.000000 ) b(2) = ( 0.000000 1.154701 0.000000 ) b(3) = ( 0.000000 0.000000 0.507525 ) PseudoPot. # 1 for Fe read from file: /home/poonam/fes2/pseudo/Fe.pz-nd-rrkjus.UPF MD5 check sum: 2e083728ad07023434bc1cc596eb954d Pseudo is Ultrasoft + core correction, Zval = 8.0 Generated by new atomic code, or converted to UPF format Using radial grid of 957 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 2 l(6) = 2 Q(r) pseudized with 0 coefficients PseudoPot. # 2 for Fe read from file: /home/poonam/fes2/pseudo/Fe.pz-nd-rrkjus.UPF MD5 check sum: 2e083728ad07023434bc1cc596eb954d Pseudo is Ultrasoft + core correction, Zval = 8.0 Generated by new atomic code, or converted to UPF format Using radial grid of 957 points, 6 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 l(5) = 2 l(6) = 2 Q(r) pseudized with 0 coefficients PseudoPot. # 3 for S read from file: /home/poonam/fes2/pseudo/S.pz-n-rrkjus_psl.0.1.UPF MD5 check sum: fce4376d1d3ac8b10d990fc1f69cce3a Pseudo is Ultrasoft + core correction, Zval = 6.0 Generated using "atomic" code by A. Dal Corso v.5.0.2 svn rev. 9415 Using radial grid of 1151 points, 4 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 Q(r) pseudized with 0 coefficients atomic species valence mass pseudopotential Fe1 8.00 1.00000 Fe( 1.00) Fe2 8.00 1.00000 Fe( 1.00) S 6.00 1.00000 S( 1.00) Starting magnetic structure atomic species magnetization Fe1 0.500 Fe2 -0.500 S 0.000 Simplified LDA+U calculation (l_max = 2) with parameters (eV): atomic species L U alpha J0 beta Fe1 2 7.1000 0.0000 0.0000 0.0000 Fe2 2 7.1000 0.0000 0.0000 0.0000 3 Sym. Ops. (no inversion) found (note: 3 additional sym.ops. were found but ignored their fractional translations are incommensurate with FFT grid) Cartesian axes site n. atom positions (alat units) 1 Fe1 tau( 1) = ( 0.3594264 0.0624582 0.2399340 ) 2 Fe1 tau( 2) = ( 0.7661964 0.2800433 0.2399340 ) 3 Fe1 tau( 3) = ( 0.3743772 0.5235239 0.2399340 ) 4 Fe2 tau( 4) = ( 0.3594264 0.0624582 0.7452398 ) 5 Fe2 tau( 5) = ( 0.7661964 0.2800433 0.7452398 ) 6 Fe2 tau( 6) = ( 0.3743772 0.5235239 0.7452398 ) 7 Fe1 tau( 7) = ( -0.1256228 0.3425015 1.2251078 ) 8 Fe1 tau( 8) = ( -0.1405736 0.8035672 1.2251078 ) 9 Fe1 tau( 9) = ( 0.2661964 0.5859821 1.2251078 ) 10 Fe2 tau( 10) = ( -0.1256228 0.3425015 1.7304136 ) 11 Fe2 tau( 11) = ( -0.1405736 0.8035672 1.7304136 ) 12 Fe2 tau( 12) = ( 0.2661964 0.5859821 1.7304136 ) 13 S tau( 13) = ( 0.0000000 0.0000000 0.0000000 ) 14 S tau( 14) = ( 0.0000000 0.0000000 0.9851738 ) 15 S tau( 15) = ( -0.0000000 0.5773503 0.0495379 ) 16 S tau( 16) = ( -0.0000000 0.5773503 0.9356359 ) 17 S tau( 17) = ( 0.5000000 0.2886751 1.9208097 ) 18 S tau( 18) = ( 0.5000000 0.2886751 1.0347117 ) 19 S tau( 19) = ( 0.1643325 0.8659802 0.4925869 ) 20 S tau( 20) = ( -0.3321271 0.5753515 0.4925869 ) 21 S tau( 21) = ( 0.1677945 0.2907192 0.4925869 ) 22 S tau( 22) = ( 0.6677945 0.5753062 1.4777607 ) 23 S tau( 23) = ( 0.6643325 0.0000453 1.4777607 ) 24 S tau( 24) = ( 0.1678729 0.2906739 1.4777607 ) number of k points= 44 Marzari-Vanderbilt smearing, width (Ry)= 0.0100 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0046296 k( 2) = ( 0.0000000 0.0000000 0.0845874), wk = 0.0092593 k( 3) = ( 0.0000000 0.0000000 0.1691749), wk = 0.0092593 k( 4) = ( 0.0000000 0.0000000 -0.2537623), wk = 0.0046296 k( 5) = ( 0.0000000 0.1924501 0.0000000), wk = 0.0277778 k( 6) = ( 0.0000000 0.1924501 0.0845874), wk = 0.0277778 k( 7) = ( 0.0000000 0.1924501 0.1691749), wk = 0.0277778 k( 8) = ( 0.0000000 0.1924501 -0.2537623), wk = 0.0277778 k( 9) = ( 0.0000000 0.3849002 0.0000000), wk = 0.0277778 k( 10) = ( 0.0000000 0.3849002 0.0845874), wk = 0.0277778 k( 11) = ( 0.0000000 0.3849002 0.1691749), wk = 0.0277778 k( 12) = ( 0.0000000 0.3849002 -0.2537623), wk = 0.0277778 k( 13) = ( 0.0000000 -0.5773503 0.0000000), wk = 0.0138889 k( 14) = ( 0.0000000 -0.5773503 0.0845874), wk = 0.0277778 k( 15) = ( 0.0000000 -0.5773503 0.1691749), wk = 0.0277778 k( 16) = ( 0.0000000 -0.5773503 -0.2537623), wk = 0.0138889 k( 17) = ( 0.1666667 0.2886751 0.0000000), wk = 0.0277778 k( 18) = ( 0.1666667 0.2886751 0.0845874), wk = 0.0277778 k( 19) = ( 0.1666667 0.2886751 0.1691749), wk = 0.0277778 k( 20) = ( 0.1666667 0.2886751 -0.2537623), wk = 0.0277778 k( 21) = ( 0.1666667 0.4811252 0.0000000), wk = 0.0277778 k( 22) = ( 0.1666667 0.4811252 0.0845874), wk = 0.0277778 k( 23) = ( 0.1666667 0.4811252 0.1691749), wk = 0.0277778 k( 24) = ( 0.1666667 0.4811252 -0.2537623), wk = 0.0277778 k( 25) = ( 0.3333333 0.5773503 0.0000000), wk = 0.0092593 k( 26) = ( 0.3333333 0.5773503 0.0845874), wk = 0.0092593 k( 27) = ( 0.3333333 0.5773503 0.1691749), wk = 0.0092593 k( 28) = ( 0.3333333 0.5773503 -0.2537623), wk = 0.0092593 k( 29) = ( -0.1666667 0.0962250 0.0845874), wk = 0.0277778 k( 30) = ( -0.1666667 0.0962250 0.1691749), wk = 0.0277778 k( 31) = ( -0.3333333 0.1924501 0.0845874), wk = 0.0277778 k( 32) = ( -0.3333333 0.1924501 0.1691749), wk = 0.0277778 k( 33) = ( 0.1666667 -0.2886751 -0.0845874), wk = 0.0277778 k( 34) = ( 0.1666667 -0.2886751 -0.1691749), wk = 0.0277778 k( 35) = ( 0.1666667 -0.4811252 0.0000000), wk = 0.0277778 k( 36) = ( 0.1666667 -0.4811252 -0.0845874), wk = 0.0277778 k( 37) = ( -0.3333333 0.3849002 0.0845874), wk = 0.0277778 k( 38) = ( -0.1666667 0.4811252 -0.0845874), wk = 0.0277778 k( 39) = ( 0.1666667 -0.4811252 -0.1691749), wk = 0.0277778 k( 40) = ( -0.3333333 0.3849002 0.1691749), wk = 0.0277778 k( 41) = ( -0.1666667 0.4811252 -0.1691749), wk = 0.0277778 k( 42) = ( 0.1666667 -0.4811252 0.2537623), wk = 0.0277778 k( 43) = ( 0.3333333 -0.5773503 -0.0845874), wk = 0.0092593 k( 44) = ( 0.3333333 -0.5773503 -0.1691749), wk = 0.0092593 Dense grid: 279295 G-vectors FFT dimensions: ( 72, 72, 135) Estimated max dynamical RAM per process > 5709.07Mb Generating pointlists ... new r_m : 0.1612 (alat units) 1.8118 (a.u.) for type 1 new r_m : 0.1612 (alat units) 1.8118 (a.u.) for type 2 new r_m : 0.1612 (alat units) 1.8118 (a.u.) for type 3 Check: negative/imaginary core charge= -0.000001 0.000000 Initial potential from superposition of free atoms starting charge 167.99397, renormalised to 168.00000 Number of +U iterations with fixed ns = 0 Starting occupations: --- enter write_ns --- LDA+U parameters: U( 1) = 7.10000000 alpha( 1) = 0.00000000 U( 2) = 7.10000000 alpha( 2) = 0.00000000 atom 1 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 2 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 3 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 4 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 atom 5 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 atom 6 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 atom 7 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 8 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 9 Tr[ns(na)] (up, down, total) = 5.00000 1.00000 6.00000 spin 1 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 spin 2 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 atomic mag. moment = 4.000000 atom 10 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 atom 11 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 atom 12 Tr[ns(na)] (up, down, total) = 1.00000 5.00000 6.00000 spin 1 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 spin 2 eigenvalues: 1.000 1.000 1.000 1.000 1.000 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: 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 atomic mag. moment = -4.000000 N of occupied +U levels = 72.000000 --- exit write_ns --- Atomic wfc used for LDA+U Projector are NOT orthogonalized Starting wfc are 120 randomized atomic wfcs total cpu time spent up to now is 471.9 secs Self-consistent Calculation iteration # 1 ecut= 90.00 Ry beta=0.70 Davidson diagonalization with overlap