Program PWSCF v.6.1 (svn rev. 13369) starts on 2Mar2017 at 23:50:12 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 Serial multi-threaded version, running on 4 processor cores Reading input from lda+U_force.in 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 O.pz-rrkjus.UPF: wavefunction(s) 2S renormalized file Fe.pz-nd-rrkjus.UPF: wavefunction(s) 4S renormalized file Fe.pz-nd-rrkjus.UPF: wavefunction(s) 4S renormalized G-vector sticks info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Sum 1061 539 163 17255 6111 1081 bravais-lattice index = 0 lattice parameter (alat) = 8.1900 a.u. unit-cell volume = 274.6766 (a.u.)^3 number of atoms/cell = 4 number of atomic types = 3 number of electrons = 28.00 number of Kohn-Sham states= 20 kinetic-energy cutoff = 30.0000 Ry charge density cutoff = 240.0000 Ry convergence threshold = 1.0E-09 mixing beta = 0.2000 number of iterations used = 8 plain mixing Exchange-correlation = SLA PZ NOGX NOGC ( 1 1 0 0 0 0) celldm(1)= 8.190000 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) = ( 0.500000 0.500000 1.000000 ) a(2) = ( 0.500000 1.000000 0.500000 ) a(3) = ( 1.000000 0.500000 0.500000 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( -0.500000 -0.500000 1.500000 ) b(2) = ( -0.500000 1.500000 -0.500000 ) b(3) = ( 1.500000 -0.500000 -0.500000 ) PseudoPot. # 1 for O read from file: /home/fs395/scratch/QE/espresso/test-suite/..//pseudo/O.pz-rrkjus.UPF MD5 check sum: 24fb942a68ef5d262e498166c462ef4a Pseudo is Ultrasoft, Zval = 6.0 Generated by new atomic code, or converted to UPF format Using radial grid of 1269 points, 4 beta functions with: l(1) = 0 l(2) = 0 l(3) = 1 l(4) = 1 Q(r) pseudized with 0 coefficients PseudoPot. # 2 for Fe read from file: /home/fs395/scratch/QE/espresso/test-suite/..//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 Fe read from file: /home/fs395/scratch/QE/espresso/test-suite/..//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 atomic species valence mass pseudopotential O1 6.00 1.00000 O ( 1.00) Fe1 8.00 1.00000 Fe( 1.00) Fe2 8.00 1.00000 Fe( 1.00) Starting magnetic structure atomic species magnetization O1 0.000 Fe1 0.500 Fe2 -0.500 Simplified LDA+U calculation (l_max = 2) with parameters (eV): atomic species L U alpha J0 beta Fe1 2 4.3000 0.0000 0.0000 0.0000 Fe2 2 4.3000 0.0000 0.0000 0.0000 6 Sym. Ops. (no inversion) found Cartesian axes site n. atom positions (alat units) 1 O1 tau( 1) = ( 0.5000000 0.5000000 0.5000000 ) 2 O1 tau( 2) = ( 1.5000000 1.5000000 1.5000000 ) 3 Fe1 tau( 3) = ( 0.0500000 0.0500000 0.0500000 ) 4 Fe2 tau( 4) = ( 0.9500000 0.9500000 0.9500000 ) number of k points= 4 gaussian smearing, width (Ry)= 0.0100 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.1250000 k( 2) = ( -0.7500000 0.2500000 0.2500000), wk = 0.3750000 k( 3) = ( -0.5000000 -0.5000000 0.5000000), wk = 0.3750000 k( 4) = ( -0.2500000 -0.2500000 -0.2500000), wk = 0.1250000 Dense grid: 17255 G-vectors FFT dimensions: ( 50, 50, 50) Smooth grid: 6111 G-vectors FFT dimensions: ( 36, 36, 36) Estimated max dynamical RAM per process > 30.86MB Generating pointlists ... new r_m : 0.1879 (alat units) 1.5389 (a.u.) for type 1 new r_m : 0.1879 (alat units) 1.5389 (a.u.) for type 2 new r_m : 0.1879 (alat units) 1.5389 (a.u.) for type 3 Check: negative/imaginary core charge= -0.000004 0.000000 Initial potential from superposition of free atoms starting charge 27.99905, renormalised to 28.00000 Number of +U iterations with fixed ns = 0 Starting occupations: --- enter write_ns --- LDA+U parameters: U( 2) = 4.30000000 alpha( 2) = 0.00000000 U( 3) = 4.30000000 alpha( 3) = 0.00000000 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 N of occupied +U levels = 12.000000 --- exit write_ns --- Atomic wfc used for LDA+U Projector are NOT orthogonalized Starting wfc are 20 randomized atomic wfcs total cpu time spent up to now is 2.3 secs per-process dynamical memory: 30.8 Mb Self-consistent Calculation iteration # 1 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 1.00E-02, avg # of iterations = 2.8 --- enter write_ns --- LDA+U parameters: U( 2) = 4.30000000 alpha( 2) = 0.00000000 U( 3) = 4.30000000 alpha( 3) = 0.00000000 atom 3 Tr[ns(na)] (up, down, total) = 4.95724 1.20131 6.15855 spin 1 eigenvalues: 0.986 0.986 0.992 0.992 1.000 eigenvectors: 0.169 0.051 0.169 0.611 0.000 0.000 0.520 0.146 0.000 0.333 0.381 0.139 0.042 0.105 0.333 0.051 0.169 0.611 0.169 0.000 0.399 0.121 0.032 0.115 0.333 occupations: 0.991 0.001 0.001 -0.000 0.002 0.001 0.992 0.004 0.002 -0.004 0.001 0.004 0.992 -0.002 -0.004 -0.000 0.002 -0.002 0.991 -0.000 0.002 -0.004 -0.004 -0.000 0.992 spin 2 eigenvalues: 0.164 0.164 0.212 0.330 0.330 eigenvectors: 0.014 0.975 0.000 0.011 0.000 0.006 0.001 0.333 0.157 0.503 0.005 0.003 0.333 0.173 0.486 0.975 0.014 0.000 0.000 0.011 0.000 0.007 0.333 0.659 0.000 occupations: 0.166 0.007 0.007 -0.000 0.014 0.007 0.290 -0.039 0.012 0.039 0.007 -0.039 0.290 -0.012 0.039 -0.000 0.012 -0.012 0.166 -0.000 0.014 0.039 0.039 -0.000 0.290 atomic mag. moment = 3.755924 atom 4 Tr[ns(na)] (up, down, total) = 1.20125 4.95726 6.15851 spin 1 eigenvalues: 0.164 0.164 0.212 0.330 0.330 eigenvectors: 0.054 0.935 0.000 0.011 0.000 0.007 0.001 0.333 0.138 0.521 0.004 0.003 0.333 0.193 0.466 0.935 0.054 0.000 0.000 0.011 0.000 0.007 0.333 0.658 0.002 occupations: 0.166 0.007 0.007 -0.000 0.014 0.007 0.290 -0.039 0.012 0.039 0.007 -0.039 0.290 -0.012 0.039 -0.000 0.012 -0.012 0.166 -0.000 0.014 0.039 0.039 -0.000 0.290 spin 2 eigenvalues: 0.986 0.986 0.992 0.992 1.000 eigenvectors: 0.157 0.060 0.161 0.622 0.000 0.001 0.521 0.145 0.000 0.333 0.405 0.116 0.043 0.102 0.333 0.060 0.157 0.622 0.161 0.000 0.377 0.145 0.030 0.115 0.333 occupations: 0.991 0.001 0.001 -0.000 0.002 0.001 0.992 0.004 0.002 -0.004 0.001 0.004 0.992 -0.002 -0.004 -0.000 0.002 -0.002 0.991 -0.000 0.002 -0.004 -0.004 -0.000 0.992 atomic mag. moment = -3.756016 N of occupied +U levels = 12.317063 --- exit write_ns --- Modify starting ns matrices according to input values --- enter write_ns --- LDA+U parameters: U( 2) = 4.30000000 alpha( 2) = 0.00000000 U( 3) = 4.30000000 alpha( 3) = 0.00000000 atom 3 Tr[ns(na)] (up, down, total) = 4.95724 1.98965 6.94689 spin 1 eigenvalues: 0.986 0.986 0.992 0.992 1.000 eigenvectors: 0.164 0.056 0.150 0.630 0.000 0.000 0.520 0.146 0.001 0.333 0.392 0.128 0.046 0.101 0.333 0.056 0.164 0.630 0.150 0.000 0.388 0.132 0.028 0.118 0.333 occupations: 0.991 0.001 0.001 -0.000 0.002 0.001 0.992 0.004 0.002 -0.004 0.001 0.004 0.992 -0.002 -0.004 -0.000 0.002 -0.002 0.991 -0.000 0.002 -0.004 -0.004 -0.000 0.992 spin 2 eigenvalues: 0.164 0.164 0.330 0.330 1.000 eigenvectors: 0.018 0.971 0.011 0.000 0.000 0.006 0.001 0.080 0.579 0.333 0.005 0.003 0.268 0.391 0.333 0.971 0.018 0.000 0.011 0.000 0.000 0.007 0.641 0.018 0.333 occupations: 0.166 0.007 0.007 -0.000 0.014 0.007 0.552 0.224 0.012 -0.224 0.007 0.224 0.552 -0.012 -0.224 -0.000 0.012 -0.012 0.166 -0.000 0.014 -0.224 -0.224 -0.000 0.552 atomic mag. moment = 2.967585 atom 4 Tr[ns(na)] (up, down, total) = 1.98961 4.95726 6.94687 spin 1 eigenvalues: 0.164 0.164 0.330 0.330 1.000 eigenvectors: 0.030 0.959 0.011 0.000 0.000 0.006 0.001 0.109 0.551 0.333 0.004 0.003 0.228 0.431 0.333 0.959 0.030 0.000 0.011 0.000 0.000 0.007 0.652 0.007 0.333 occupations: 0.166 0.007 0.007 -0.000 0.014 0.007 0.552 0.224 0.012 -0.224 0.007 0.224 0.552 -0.012 -0.224 -0.000 0.012 -0.012 0.166 -0.000 0.014 -0.224 -0.224 -0.000 0.552 spin 2 eigenvalues: 0.986 0.986 0.992 0.992 1.000 eigenvectors: 0.151 0.067 0.090 0.692 0.000 0.002 0.520 0.140 0.004 0.333 0.418 0.103 0.060 0.085 0.333 0.067 0.151 0.692 0.090 0.000 0.362 0.160 0.017 0.128 0.333 occupations: 0.991 0.001 0.001 -0.000 0.002 0.001 0.992 0.004 0.002 -0.004 0.001 0.004 0.992 -0.002 -0.004 -0.000 0.002 -0.002 0.991 -0.000 0.002 -0.004 -0.004 -0.000 0.992 atomic mag. moment = -2.967658 N of occupied +U levels = 13.893761 --- exit write_ns --- total cpu time spent up to now is 3.6 secs total energy = -174.18604535 Ry Harris-Foulkes estimate = -174.79123784 Ry estimated scf accuracy < 2.56513173 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 8.24 Bohr mag/cell iteration # 2 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 9.16E-03, avg # of iterations = 1.5 total cpu time spent up to now is 4.7 secs total energy = -174.30056999 Ry Harris-Foulkes estimate = -174.33044072 Ry estimated scf accuracy < 0.27058696 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 7.06 Bohr mag/cell iteration # 3 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 9.66E-04, avg # of iterations = 2.8 total cpu time spent up to now is 5.9 secs total energy = -174.33844310 Ry Harris-Foulkes estimate = -174.33504497 Ry estimated scf accuracy < 0.11449333 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 6.86 Bohr mag/cell iteration # 4 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 4.09E-04, avg # of iterations = 1.2 total cpu time spent up to now is 7.0 secs total energy = -174.34928424 Ry Harris-Foulkes estimate = -174.34636263 Ry estimated scf accuracy < 0.05868962 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 6.95 Bohr mag/cell iteration # 5 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 2.10E-04, avg # of iterations = 1.0 total cpu time spent up to now is 8.1 secs total energy = -174.35359647 Ry Harris-Foulkes estimate = -174.35156451 Ry estimated scf accuracy < 0.01541826 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 6.90 Bohr mag/cell iteration # 6 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 5.51E-05, avg # of iterations = 1.8 total cpu time spent up to now is 9.3 secs total energy = -174.35621602 Ry Harris-Foulkes estimate = -174.35650315 Ry estimated scf accuracy < 0.00308054 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.83 Bohr mag/cell iteration # 7 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 1.10E-05, avg # of iterations = 2.4 total cpu time spent up to now is 10.4 secs total energy = -174.35684609 Ry Harris-Foulkes estimate = -174.35703636 Ry estimated scf accuracy < 0.00138500 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.86 Bohr mag/cell iteration # 8 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 4.95E-06, avg # of iterations = 2.0 total cpu time spent up to now is 11.6 secs total energy = -174.35704157 Ry Harris-Foulkes estimate = -174.35701713 Ry estimated scf accuracy < 0.00009676 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.85 Bohr mag/cell iteration # 9 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 3.46E-07, avg # of iterations = 4.0 total cpu time spent up to now is 13.0 secs total energy = -174.35709784 Ry Harris-Foulkes estimate = -174.35708862 Ry estimated scf accuracy < 0.00002283 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 10 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 8.15E-08, avg # of iterations = 2.2 total cpu time spent up to now is 14.1 secs total energy = -174.35710342 Ry Harris-Foulkes estimate = -174.35710021 Ry estimated scf accuracy < 0.00001290 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 11 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 4.61E-08, avg # of iterations = 1.1 total cpu time spent up to now is 15.2 secs total energy = -174.35710646 Ry Harris-Foulkes estimate = -174.35710423 Ry estimated scf accuracy < 0.00000204 Ry total magnetization = -0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 12 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 7.28E-09, avg # of iterations = 2.6 total cpu time spent up to now is 16.4 secs total energy = -174.35710831 Ry Harris-Foulkes estimate = -174.35710685 Ry estimated scf accuracy < 0.00000031 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 13 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 1.09E-09, avg # of iterations = 3.1 total cpu time spent up to now is 17.7 secs total energy = -174.35710932 Ry Harris-Foulkes estimate = -174.35710839 Ry estimated scf accuracy < 0.00000005 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 14 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 1.96E-10, avg # of iterations = 2.5 total cpu time spent up to now is 18.8 secs total energy = -174.35710993 Ry Harris-Foulkes estimate = -174.35710933 Ry estimated scf accuracy < 0.00000004 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 15 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 1.26E-10, avg # of iterations = 1.6 total cpu time spent up to now is 19.9 secs total energy = -174.35711031 Ry Harris-Foulkes estimate = -174.35710993 Ry estimated scf accuracy < 1.9E-09 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell iteration # 16 ecut= 30.00 Ry beta=0.20 Davidson diagonalization with overlap ethr = 6.73E-12, avg # of iterations = 4.0 Magnetic moment per site: atom: 1 charge: 5.2453 magn: 0.0000 constr: 0.0000 atom: 2 charge: 5.5752 magn: -0.0000 constr: 0.0000 atom: 3 charge: 5.4036 magn: 3.0174 constr: 0.0000 atom: 4 charge: 5.4036 magn: -3.0174 constr: 0.0000 total cpu time spent up to now is 21.6 secs End of self-consistent calculation --- enter write_ns --- LDA+U parameters: U( 2) = 4.30000000 alpha( 2) = 0.00000000 U( 3) = 4.30000000 alpha( 3) = 0.00000000 atom 3 Tr[ns(na)] (up, down, total) = 4.95171 1.94192 6.89363 spin 1 eigenvalues: 0.987 0.987 0.990 0.994 0.994 eigenvectors: 0.136 0.564 0.000 0.152 0.148 0.062 0.138 0.333 0.434 0.033 0.199 0.001 0.333 0.030 0.437 0.564 0.136 0.000 0.148 0.152 0.039 0.161 0.333 0.237 0.230 occupations: 0.989 0.001 0.001 -0.000 0.002 0.001 0.991 -0.001 0.002 0.001 0.001 -0.001 0.991 -0.002 0.001 -0.000 0.002 -0.002 0.989 -0.000 0.002 0.001 0.001 -0.000 0.991 spin 2 eigenvalues: 0.130 0.130 0.345 0.345 0.992 eigenvectors: 0.044 0.002 0.070 0.884 0.000 0.056 0.580 0.029 0.002 0.333 0.293 0.343 0.015 0.016 0.333 0.002 0.044 0.884 0.070 0.000 0.605 0.031 0.002 0.029 0.333 occupations: 0.335 0.018 0.018 -0.000 0.037 0.018 0.424 0.284 0.032 -0.284 0.018 0.284 0.424 -0.032 -0.284 -0.000 0.032 -0.032 0.335 -0.000 0.037 -0.284 -0.284 -0.000 0.424 atomic mag. moment = 3.009787 atom 4 Tr[ns(na)] (up, down, total) = 1.94192 4.95171 6.89363 spin 1 eigenvalues: 0.130 0.130 0.345 0.345 0.992 eigenvectors: 0.045 0.001 0.100 0.854 0.000 0.079 0.556 0.030 0.001 0.333 0.255 0.381 0.013 0.018 0.333 0.001 0.045 0.854 0.100 0.000 0.619 0.017 0.003 0.028 0.333 occupations: 0.335 0.018 0.018 -0.000 0.037 0.018 0.424 0.284 0.032 -0.284 0.018 0.284 0.424 -0.032 -0.284 -0.000 0.032 -0.032 0.335 -0.000 0.037 -0.284 -0.284 -0.000 0.424 spin 2 eigenvalues: 0.987 0.987 0.990 0.994 0.994 eigenvectors: 0.120 0.580 0.000 0.147 0.153 0.068 0.133 0.333 0.438 0.029 0.198 0.002 0.333 0.034 0.433 0.580 0.120 0.000 0.153 0.147 0.034 0.166 0.333 0.229 0.238 occupations: 0.989 0.001 0.001 -0.000 0.002 0.001 0.991 -0.001 0.002 0.001 0.001 -0.001 0.991 -0.002 0.001 -0.000 0.002 -0.002 0.989 -0.000 0.002 0.001 0.001 -0.000 0.991 atomic mag. moment = -3.009787 N of occupied +U levels = 13.787264 --- exit write_ns --- ------ SPIN UP ------------ k = 0.0000 0.0000 0.0000 ( 731 PWs) bands (ev): -11.1248 -6.1210 0.3144 1.4045 1.4045 4.6819 4.6819 7.1317 7.2022 7.2022 7.8852 9.3139 9.7408 9.7408 10.8372 11.8145 11.8145 12.8076 12.8076 15.9708 k =-0.7500 0.2500 0.2500 ( 764 PWs) bands (ev): -10.7534 -5.6560 0.4964 1.9148 2.7428 4.7384 4.9826 5.1443 5.8428 6.7548 7.0426 8.4920 9.1241 9.5992 11.6066 11.7491 12.3100 12.8910 17.0318 17.2642 k =-0.5000-0.5000 0.5000 ( 780 PWs) bands (ev): -10.7395 -5.6254 0.2854 1.9206 2.7931 4.8778 4.9808 5.1444 5.7796 6.8321 7.2543 8.5125 9.3543 9.3601 11.6445 11.8631 12.5792 12.9855 15.5103 16.6922 k =-0.2500-0.2500-0.2500 ( 748 PWs) bands (ev): -11.1643 -6.2050 1.1182 1.4278 1.4278 4.6921 4.6921 6.0449 6.9851 6.9851 7.3866 9.1313 10.0954 10.0954 11.7970 11.7970 12.5813 12.5813 13.1209 14.7930 ------ SPIN DOWN ---------- k = 0.0000 0.0000 0.0000 ( 731 PWs) bands (ev): -11.1248 -6.1210 0.3144 1.4045 1.4045 4.6819 4.6819 7.1317 7.2022 7.2022 7.8852 9.3139 9.7408 9.7408 10.8372 11.8145 11.8145 12.8076 12.8076 15.9708 k =-0.7500 0.2500 0.2500 ( 764 PWs) bands (ev): -10.7534 -5.6560 0.4964 1.9148 2.7428 4.7384 4.9826 5.1443 5.8428 6.7548 7.0426 8.4920 9.1241 9.5992 11.6066 11.7491 12.3100 12.8910 17.0318 17.2642 k =-0.5000-0.5000 0.5000 ( 780 PWs) bands (ev): -10.7395 -5.6254 0.2854 1.9206 2.7931 4.8778 4.9808 5.1444 5.7796 6.8321 7.2543 8.5125 9.3543 9.3601 11.6445 11.8631 12.5792 12.9855 15.5103 16.6922 k =-0.2500-0.2500-0.2500 ( 748 PWs) bands (ev): -11.1643 -6.2049 1.1182 1.4278 1.4278 4.6921 4.6921 6.0449 6.9851 6.9851 7.3866 9.1313 10.0954 10.0954 11.7970 11.7970 12.5813 12.5813 13.1209 14.7930 the Fermi energy is 10.4744 ev ! total energy = -174.35711099 Ry Harris-Foulkes estimate = -174.35711031 Ry estimated scf accuracy < 7.1E-10 Ry The total energy is the sum of the following terms: one-electron contribution = -11.12517498 Ry hartree contribution = 33.36687165 Ry xc contribution = -65.91454662 Ry ewald contribution = -130.91621185 Ry Hubbard energy = 0.23195199 Ry smearing contrib. (-TS) = -0.00000118 Ry total magnetization = 0.00 Bohr mag/cell absolute magnetization = 6.84 Bohr mag/cell convergence has been achieved in 16 iterations Forces acting on atoms (cartesian axes, Ry/au): atom 1 type 1 force = 0.00000040 0.00000040 0.00000040 atom 2 type 1 force = -0.00000101 -0.00000101 -0.00000101 atom 3 type 2 force = -0.14189797 -0.14189797 -0.14189797 atom 4 type 3 force = 0.14189858 0.14189858 0.14189858 Total force = 0.347578 Total SCF correction = 0.000021 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 44.52 0.00030264 -0.00079819 -0.00079819 44.52 -117.42 -117.42 -0.00079819 0.00030264 -0.00079819 -117.42 44.52 -117.42 -0.00079819 -0.00079819 0.00030264 -117.42 -117.42 44.52 Writing output data file pwscf.save init_run : 3.42s CPU 2.20s WALL ( 1 calls) electrons : 62.67s CPU 19.25s WALL ( 1 calls) forces : 0.95s CPU 0.45s WALL ( 1 calls) stress : 3.97s CPU 1.38s WALL ( 1 calls) Called by init_run: wfcinit : 1.26s CPU 0.37s WALL ( 1 calls) potinit : 0.14s CPU 0.08s WALL ( 1 calls) Called by electrons: c_bands : 37.05s CPU 11.15s WALL ( 16 calls) sum_band : 18.67s CPU 5.63s WALL ( 16 calls) v_of_rho : 1.25s CPU 0.34s WALL ( 17 calls) newd : 5.74s CPU 1.83s WALL ( 17 calls) mix_rho : 0.25s CPU 0.25s WALL ( 16 calls) Called by c_bands: init_us_2 : 0.78s CPU 0.20s WALL ( 368 calls) cegterg : 36.18s CPU 10.88s WALL ( 128 calls) Called by sum_band: sum_band:bec : 0.03s CPU 0.01s WALL ( 128 calls) addusdens : 9.22s CPU 2.97s WALL ( 16 calls) Called by *egterg: h_psi : 35.21s CPU 9.47s WALL ( 429 calls) s_psi : 0.61s CPU 0.53s WALL ( 517 calls) g_psi : 0.02s CPU 0.02s WALL ( 293 calls) cdiaghg : 0.51s CPU 0.33s WALL ( 421 calls) Called by h_psi: h_psi:pot : 34.99s CPU 9.29s WALL ( 429 calls) h_psi:calbec : 1.71s CPU 0.44s WALL ( 429 calls) vloc_psi : 32.69s CPU 8.43s WALL ( 429 calls) add_vuspsi : 0.58s CPU 0.42s WALL ( 429 calls) vhpsi : 0.21s CPU 0.16s WALL ( 429 calls) General routines calbec : 3.76s CPU 1.20s WALL ( 2546 calls) fft : 1.12s CPU 0.87s WALL ( 289 calls) ffts : 0.11s CPU 0.06s WALL ( 66 calls) fftw : 34.66s CPU 8.88s WALL ( 14640 calls) interpolate : 0.35s CPU 0.28s WALL ( 66 calls) davcio : 0.00s CPU 0.00s WALL ( 16 calls) Hubbard U routines new_ns : 0.04s CPU 0.04s WALL ( 16 calls) vhpsi : 0.21s CPU 0.16s WALL ( 429 calls) force_hub : 0.18s CPU 0.11s WALL ( 1 calls) stres_hub : 2.04s CPU 0.72s WALL ( 1 calls) PWSCF : 1m11.15s CPU 0m23.47s WALL This run was terminated on: 23:50:36 2Mar2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=