Program PWSCF v.5.1.1 starts on 28Oct2015 at 19: 1:56 

     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    48 processors
     R & G space division:  proc/nbgrp/npool/nimage =      48
     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 Nb.rel-pbe-spn-rrkjus_psl.1.0.0.UPF: wavefunction(s)  4S 4P 4P renormalized
               file Cu.rel-pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s)  3D renormalized

     Subspace diagonalization in iterative solution of the eigenvalue problem:
     a serial algorithm will be used

 
     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min          47      34      9                 1720     1040     155
     Max          49      35     10                 1724     1061     160
     Sum        2305    1649    465                82663    50541    7497
 


     bravais-lattice index     =           14
     lattice parameter (alat)  =      10.3935  a.u.
     unit-cell volume          =    1122.7539 (a.u.)^3
     number of atoms/cell      =            8
     number of atomic types    =            3
     number of electrons       =        70.00
     number of Kohn-Sham states=           84
     kinetic-energy cutoff     =      48.0000  Ry
     charge density cutoff     =     267.0000  Ry
     convergence threshold     =      1.0E-08
     mixing beta               =       0.3000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      = PBE ( 1  4  3  4 0 0)
     Non magnetic calculation with spin-orbit


     celldm(1)=  10.393493  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  S read from file:
     /home/autes/Pseudo/S.rel-pbe-n-rrkjus_psl.1.0.0.UPF
     MD5 check sum: 95f1324244b600346090892298cb4451
     Pseudo is Ultrasoft + core correction, Zval =  6.0
     Generated using "atomic" code by A. Dal Corso  v.5.0.2 svn rev. 9656
     Using radial grid of 1151 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 Nb read from file:
     /home/autes/Pseudo/Nb.rel-pbe-spn-rrkjus_psl.1.0.0.UPF
     MD5 check sum: c1e153ab096482f446ae869583c6ae55
     Pseudo is Ultrasoft + core correction, Zval = 13.0
     Generated using "atomic" code by A. Dal Corso  v.5.0.2 svn rev. 9656
     Using radial grid of 1227 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. # 3 for Cu read from file:
     /home/autes/Pseudo/Cu.rel-pbe-dn-rrkjus_psl.1.0.0.UPF
     MD5 check sum: bfb9b73be46f5385f72b338ff57911ad
     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 1199 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
        S              6.00    32.06500      S( 1.00)
        Nb            13.00    92.90640     Nb( 1.00)
        Cu            11.00    63.54600     Cu( 1.00)

     24 Sym. Ops. (no 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     120 deg rotation - cart. axis [-1,-1,-1]     

 cryst.   s( 5) = (     0          0          1      )
                  (     1          0          0      )
                  (     0          1          0      )

 cart.    s( 5) = (  0.0000000  1.0000000  0.0000000 )
                  (  0.0000000 -0.0000000  1.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )


      isym =  6     120 deg rotation - cart. axis [-1,1,1]       

 cryst.   s( 6) = (     0          0         -1      )
                  (    -1          0          0      )
                  (     0          1          0      )

 cart.    s( 6) = ( -0.0000000 -1.0000000  0.0000000 )
                  ( -0.0000000  0.0000000  1.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )


      isym =  7     120 deg rotation - cart. axis [1,1,-1]       

 cryst.   s( 7) = (     0          0         -1      )
                  (     1          0          0      )
                  (     0         -1          0      )

 cart.    s( 7) = ( -0.0000000  1.0000000 -0.0000000 )
                  ( -0.0000000  0.0000000 -1.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )


      isym =  8     120 deg rotation - cart. axis [1,-1,1]       

 cryst.   s( 8) = (     0          0          1      )
                  (    -1          0          0      )
                  (     0         -1          0      )

 cart.    s( 8) = (  0.0000000 -1.0000000 -0.0000000 )
                  (  0.0000000 -0.0000000 -1.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )


      isym =  9     120 deg rotation - cart. axis [1,1,1]        

 cryst.   s( 9) = (     0          1          0      )
                  (     0          0          1      )
                  (     1          0          0      )

 cart.    s( 9) = (  0.0000000  0.0000000  1.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )
                  (  0.0000000  1.0000000 -0.0000000 )


      isym = 10     120 deg rotation - cart. axis [-1,1,-1]      

 cryst.   s(10) = (     0         -1          0      )
                  (     0          0         -1      )
                  (     1          0          0      )

 cart.    s(10) = ( -0.0000000 -0.0000000  1.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )
                  (  0.0000000 -1.0000000  0.0000000 )


      isym = 11     120 deg rotation - cart. axis [1,-1,-1]      

 cryst.   s(11) = (     0         -1          0      )
                  (     0          0          1      )
                  (    -1          0          0      )

 cart.    s(11) = ( -0.0000000  0.0000000 -1.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  1.0000000 -0.0000000 )


      isym = 12     120 deg rotation - cart. axis [-1,-1,1]      

 cryst.   s(12) = (     0          1          0      )
                  (     0          0         -1      )
                  (    -1          0          0      )

 cart.    s(12) = (  0.0000000 -0.0000000 -1.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )
                  (  0.0000000 -1.0000000  0.0000000 )


      isym = 13     inv. 180 deg rotation - cart. axis [1,1,0]   

 cryst.   s(13) = (     0         -1          0      )
                  (    -1          0          0      )
                  (     0          0          1      )

 cart.    s(13) = ( -0.0000000 -1.0000000  0.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


      isym = 14     inv. 180 deg rotation - cart. axis [1,-1,0]  

 cryst.   s(14) = (     0          1          0      )
                  (     1          0          0      )
                  (     0          0          1      )

 cart.    s(14) = (  0.0000000  1.0000000  0.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


      isym = 15     inv.  90 deg rotation - cart. axis [0,0,-1]  

 cryst.   s(15) = (     0          1          0      )
                  (    -1          0          0      )
                  (     0          0         -1      )

 cart.    s(15) = (  0.0000000 -1.0000000 -0.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )
                  (  0.0000000  0.0000000 -1.0000000 )


      isym = 16     inv.  90 deg rotation - cart. axis [0,0,1]   

 cryst.   s(16) = (     0         -1          0      )
                  (     1          0          0      )
                  (     0          0         -1      )

 cart.    s(16) = ( -0.0000000  1.0000000 -0.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  0.0000000 -1.0000000 )


      isym = 17     inv. 180 deg rotation - cart. axis [1,0,1]   

 cryst.   s(17) = (     0          0         -1      )
                  (     0          1          0      )
                  (    -1          0          0      )

 cart.    s(17) = ( -0.0000000  0.0000000 -1.0000000 )
                  ( -0.0000000  1.0000000 -0.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )


      isym = 18     inv. 180 deg rotation - cart. axis [-1,0,1]  

 cryst.   s(18) = (     0          0          1      )
                  (     0          1          0      )
                  (     1          0          0      )

 cart.    s(18) = (  0.0000000  0.0000000  1.0000000 )
                  (  0.0000000  1.0000000 -0.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )


      isym = 19     inv.  90 deg rotation - cart. axis [0,1,0]   

 cryst.   s(19) = (     0          0          1      )
                  (     0         -1          0      )
                  (    -1          0          0      )

 cart.    s(19) = (  0.0000000 -0.0000000 -1.0000000 )
                  (  0.0000000 -1.0000000  0.0000000 )
                  (  1.0000000 -0.0000000 -0.0000000 )


      isym = 20     inv.  90 deg rotation - cart. axis [0,-1,0]  

 cryst.   s(20) = (     0          0         -1      )
                  (     0         -1          0      )
                  (     1          0          0      )

 cart.    s(20) = ( -0.0000000 -0.0000000  1.0000000 )
                  ( -0.0000000 -1.0000000  0.0000000 )
                  ( -1.0000000  0.0000000  0.0000000 )


      isym = 21     inv. 180 deg rotation - cart. axis [0,1,1]   

 cryst.   s(21) = (     1          0          0      )
                  (     0          0         -1      )
                  (     0         -1          0      )

 cart.    s(21) = (  1.0000000 -0.0000000 -0.0000000 )
                  (  0.0000000 -0.0000000 -1.0000000 )
                  (  0.0000000 -1.0000000  0.0000000 )


      isym = 22     inv. 180 deg rotation - cart. axis [0,1,-1]  

 cryst.   s(22) = (     1          0          0      )
                  (     0          0          1      )
                  (     0          1          0      )

 cart.    s(22) = (  1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )
                  (  0.0000000  1.0000000 -0.0000000 )


      isym = 23     inv.  90 deg rotation - cart. axis [-1,0,0]  

 cryst.   s(23) = (    -1          0          0      )
                  (     0          0          1      )
                  (     0         -1          0      )

 cart.    s(23) = ( -1.0000000  0.0000000 -0.0000000 )
                  (  0.0000000  0.0000000 -1.0000000 )
                  (  0.0000000  1.0000000 -0.0000000 )


      isym = 24     inv.  90 deg rotation - cart. axis [1,0,0]   

 cryst.   s(24) = (    -1          0          0      )
                  (     0          0         -1      )
                  (     0          1          0      )

 cart.    s(24) = ( -1.0000000  0.0000000  0.0000000 )
                  (  0.0000000 -0.0000000  1.0000000 )
                  (  0.0000000 -1.0000000  0.0000000 )


     double point group T_d (-43m) 
     there are  8 classes and  3 irreducible representations
     the character table:

       E     -E    8C3   -8C3   3C2  6s_d  6S4   -6S4 
                               -3C2  -6s_d            
G_6    2.00 -2.00  1.00 -1.00  0.00  0.00  1.41 -1.41
G_7    2.00 -2.00  1.00 -1.00  0.00  0.00 -1.41  1.41
G_8    4.00 -4.00 -1.00  1.00  0.00  0.00  0.00  0.00

     the symmetry operations in each class:
     E             1
      3C2 -3C2     2   -2    3   -3    4   -4
     8C3           5    7    8    6   12    9   10   11
     6s_d -6s_d   13  -13   14  -14   17  -17   18  -18   21  -21   22  -22
     6S4          15   16   20   19   23   24
     -E           -1
     -8C3         -5   -7   -8   -6  -12   -9  -10  -11
     -6S4        -15  -16  -20  -19  -23  -24

   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:    82663 G-vectors     FFT dimensions: (  60,  60,  60)

     Smooth grid:    50541 G-vectors     FFT dimensions: (  45,  45,  45)

     Largest allocated arrays     est. size (Mb)     dimensions
        Kohn-Sham Wavefunctions         0.34 Mb     (     266,   84)
        NL pseudopotentials             0.55 Mb     (     133,  272)
        Each V/rho on FFT grid          0.11 Mb     (    7200)
        Each G-vector array             0.01 Mb     (    1722)
        G-vector shells                 0.00 Mb     (     427)
     Largest temporary arrays     est. size (Mb)     dimensions
        Auxiliary wavefunctions         1.36 Mb     (     266,  336)
        Each subspace H/S matrix        1.72 Mb     (     336,  336)
        Each <psi_i|beta_j> matrix      0.70 Mb     (     272,   2,   84)
        Arrays for rho mixing           0.88 Mb     (    7200,    8)

     Check: negative/imaginary core charge=   -0.000005    0.000000

     Initial potential from superposition of free atoms

     starting charge   69.99785, renormalised to   70.00000
     Starting wfc are  112 randomized atomic wfcs

     total cpu time spent up to now is       46.3 secs

     per-process dynamical memory:    51.4 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  3.0

     Threshold (ethr) on eigenvalues was too large:
     Diagonalizing with lowered threshold

     Davidson diagonalization with overlap
     ethr =  4.22E-04,  avg # of iterations =  2.5

     total cpu time spent up to now is       61.8 secs

     total energy              =    -570.29634349 Ry
     Harris-Foulkes estimate   =    -570.43239333 Ry
     estimated scf accuracy    <       0.34278103 Ry

     iteration #  2     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  4.90E-04,  avg # of iterations =  2.0

     total cpu time spent up to now is       68.9 secs

     total energy              =    -570.32820122 Ry
     Harris-Foulkes estimate   =    -570.38940184 Ry
     estimated scf accuracy    <       0.23051692 Ry

     iteration #  3     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.29E-04,  avg # of iterations =  1.5

     total cpu time spent up to now is       74.4 secs

     total energy              =    -570.34713163 Ry
     Harris-Foulkes estimate   =    -570.35647002 Ry
     estimated scf accuracy    <       0.03137102 Ry

     iteration #  4     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  4.48E-05,  avg # of iterations =  4.0

     total cpu time spent up to now is       84.5 secs

     total energy              =    -570.35177192 Ry
     Harris-Foulkes estimate   =    -570.35306387 Ry
     estimated scf accuracy    <       0.00782993 Ry

     iteration #  5     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.12E-05,  avg # of iterations =  2.0

     total cpu time spent up to now is       90.3 secs

     total energy              =    -570.35209950 Ry
     Harris-Foulkes estimate   =    -570.35240960 Ry
     estimated scf accuracy    <       0.00225736 Ry

     iteration #  6     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.22E-06,  avg # of iterations =  2.2

     total cpu time spent up to now is       96.7 secs

     total energy              =    -570.35230942 Ry
     Harris-Foulkes estimate   =    -570.35251957 Ry
     estimated scf accuracy    <       0.00101179 Ry

     iteration #  7     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.45E-06,  avg # of iterations =  1.3

     total cpu time spent up to now is      102.0 secs

     total energy              =    -570.35238598 Ry
     Harris-Foulkes estimate   =    -570.35240242 Ry
     estimated scf accuracy    <       0.00004649 Ry

     iteration #  8     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  6.64E-08,  avg # of iterations =  3.2

     total cpu time spent up to now is      110.3 secs

     total energy              =    -570.35240224 Ry
     Harris-Foulkes estimate   =    -570.35240559 Ry
     estimated scf accuracy    <       0.00000911 Ry

     iteration #  9     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.30E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is      117.1 secs

     total energy              =    -570.35240197 Ry
     Harris-Foulkes estimate   =    -570.35240409 Ry
     estimated scf accuracy    <       0.00000665 Ry

     iteration # 10     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  9.51E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is      123.0 secs

     total energy              =    -570.35240305 Ry
     Harris-Foulkes estimate   =    -570.35240307 Ry
     estimated scf accuracy    <       0.00000016 Ry

     iteration # 11     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  2.30E-10,  avg # of iterations =  3.7

     total cpu time spent up to now is      131.3 secs

     total energy              =    -570.35240313 Ry
     Harris-Foulkes estimate   =    -570.35240313 Ry
     estimated scf accuracy    <       0.00000007 Ry

     iteration # 12     ecut=    48.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  9.32E-11,  avg # of iterations =  2.0

     total cpu time spent up to now is      137.6 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.0000 (  6355 PWs)   bands (ev):

   -48.7714 -48.7714 -26.3894 -26.3894 -24.3060 -24.3060 -24.3060 -24.3060
    -7.0259  -7.0259  -6.0223  -6.0223  -5.9918  -5.9918  -5.9918  -5.9918
     0.7035   0.7035   0.7035   0.7035   1.6008   1.6008   2.3126   2.3126
     2.3126   2.3126   2.3473   2.3473   2.9724   2.9724   3.0106   3.0106
     3.0106   3.0106   3.4705   3.4705   3.4913   3.4913   3.4913   3.4913
     3.5520   3.5520   3.5520   3.5520   4.2955   4.2955   4.3715   4.3715
     4.3715   4.3715   4.5010   4.5010   4.5010   4.5010   4.5292   4.5292
     4.7935   4.7935   4.9518   4.9518   4.9518   4.9518   4.9704   4.9704
     6.0024   6.0024   6.0024   6.0024   6.0892   6.0892   8.1803   8.1803
     8.2022   8.2022   8.2022   8.2022   8.2685   8.2685   8.2685   8.2685
     8.4583   8.4583  12.0682  12.0682

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.0000 0.2000 (  6327 PWs)   bands (ev):

   -48.7713 -48.7713 -26.3897 -26.3897 -24.3069 -24.3069 -24.3060 -24.3060
    -6.9438  -6.9438  -6.1015  -6.1015  -6.0091  -6.0091  -5.9920  -5.9920
     0.7890   0.7890   0.8672   0.8672   1.7041   1.7041   2.1650   2.1650
     2.3536   2.3536   2.3755   2.3755   2.7990   2.7990   2.9591   2.9591
     2.9828   2.9828   3.2204   3.2204   3.4877   3.4877   3.5057   3.5057
     3.6047   3.6047   3.7784   3.7784   4.1905   4.1905   4.3155   4.3155
     4.3308   4.3308   4.4747   4.4747   4.5017   4.5017   4.5580   4.5580
     4.8062   4.8062   4.8906   4.8906   4.9521   4.9521   5.0045   5.0045
     5.7412   5.7412   5.9374   5.9374   5.9869   5.9869   7.9544   7.9544
     8.2953   8.2953   8.3523   8.3523   8.3904   8.3904   8.4091   8.4091
     9.1454   9.1454  11.8464  11.8464

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.0000 0.4000 (  6311 PWs)   bands (ev):

   -48.7711 -48.7711 -26.3902 -26.3902 -24.3082 -24.3082 -24.3060 -24.3060
    -6.7433  -6.7433  -6.3246  -6.3246  -6.0095  -6.0095  -5.9923  -5.9923
     0.9608   0.9608   1.2657   1.2657   1.7669   1.7669   2.0393   2.0393
     2.4265   2.4265   2.4432   2.4432   2.4757   2.4757   2.7555   2.7555
     2.9122   2.9122   2.9362   2.9362   3.5089   3.5089   3.5389   3.5389
     3.6638   3.6638   4.0594   4.0594   4.0857   4.0857   4.2131   4.2131
     4.2654   4.2654   4.3975   4.3975   4.4705   4.4705   4.5003   4.5003
     4.8023   4.8023   4.8448   4.8448   5.0142   5.0142   5.0622   5.0622
     5.3532   5.3532   5.8547   5.8547   5.8876   5.8876   7.8003   7.8003
     8.3553   8.3553   8.5723   8.5723   8.6029   8.6029   8.6272   8.6272
    10.1606  10.1606  11.6239  11.6240

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.2000 0.2000 (  6330 PWs)   bands (ev):

   -48.7712 -48.7712 -26.3900 -26.3900 -24.3073 -24.3073 -24.3064 -24.3064
    -6.8663  -6.8660  -6.1055  -6.0982  -6.0805  -6.0771  -6.0115  -6.0034
     0.9015   0.9066   0.9666   0.9700   1.7894   1.7920   2.0719   2.0803
     2.3334   2.3528   2.3985   2.3991   2.7324   2.7340   2.8105   2.8411
     2.9340   2.9626   3.2088   3.2390   3.3017   3.3226   3.4661   3.4890
     3.7024   3.7117   3.8726   3.8807   4.0985   4.1815   4.2055   4.2298
     4.3145   4.3626   4.4013   4.4092   4.4975   4.5227   4.6457   4.6495
     4.7126   4.8059   4.8271   4.8351   4.9479   4.9535   5.1724   5.2260
     5.4858   5.4939   5.8228   5.8312   5.9884   6.0010   8.0692   8.1124
     8.1475   8.1572   8.3800   8.3888   8.4720   8.5420   8.5773   8.5866
     9.5422   9.5458  11.9614  11.9699

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.2000 0.4000 (  6284 PWs)   bands (ev):

   -48.7710 -48.7710 -26.3905 -26.3905 -24.3085 -24.3085 -24.3066 -24.3066
    -6.6770  -6.6767  -6.2992  -6.2926  -6.0718  -6.0628  -6.0317  -6.0223
     1.0587   1.0650   1.3000   1.3013   1.7323   1.7350   2.0797   2.0859
     2.3922   2.4073   2.4446   2.4473   2.4772   2.4869   2.6784   2.6961
     2.8867   2.8921   2.9114   2.9229   3.3432   3.3571   3.4120   3.4390
     3.7443   3.7486   4.0254   4.0750   4.0877   4.1078   4.1596   4.1780
     4.2447   4.3112   4.3641   4.3956   4.4985   4.5344   4.6160   4.6506
     4.6991   4.7236   4.8074   4.8247   4.9671   5.0000   5.1901   5.2306
     5.3466   5.4036   5.7014   5.7373   5.8387   5.8551   7.9574   7.9844
     8.3377   8.3874   8.4708   8.4723   8.5784   8.6271   8.8071   8.8567
    10.3559  10.3614  11.7321  11.7467

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.4000 0.4000 (  6308 PWs)   bands (ev):

   -48.7710 -48.7710 -26.3910 -26.3910 -24.3093 -24.3093 -24.3072 -24.3072
    -6.5157  -6.5152  -6.2422  -6.2346  -6.2161  -6.2120  -6.1007  -6.0897
     1.2327   1.2411   1.3949   1.3976   1.5848   1.5868   2.2828   2.2886
     2.4020   2.4138   2.4371   2.4480   2.5260   2.5320   2.5603   2.5698
     2.7503   2.7686   2.8886   2.9011   3.1367   3.1537   3.2126   3.2308
     3.8049   3.8089   4.0704   4.1296   4.1334   4.1602   4.1668   4.1711
     4.2883   4.2914   4.3058   4.3284   4.4020   4.5098   4.6748   4.6904
     4.7853   4.8049   4.8406   4.8574   4.8797   4.8940   5.0319   5.0843
     5.4874   5.5207   5.5338   5.5458   5.9310   5.9410   8.2430   8.2551
     8.2838   8.3298   8.3676   8.3732   8.8698   8.9286   8.9435   8.9439
    10.8475  10.8539  11.5338  11.5409

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.2000 0.2000 0.2000 (  6322 PWs)   bands (ev):

   -48.7712 -48.7712 -26.3903 -26.3903 -24.3073 -24.3073 -24.3073 -24.3073
    -6.7935  -6.7935  -6.1078  -6.1078  -6.0830  -6.0798  -6.0727  -6.0727
     1.0220   1.0223   1.0223   1.0315   1.8461   1.8461   1.9803   1.9803
     2.3598   2.3877   2.3920   2.3920   2.7199   2.7199   2.7240   2.7258
     2.8897   2.8897   3.1819   3.1819   3.3193   3.3193   3.3200   3.3451
     3.8355   3.8355   3.8745   3.8815   4.1008   4.1515   4.1515   4.2813
     4.2813   4.2874   4.3803   4.3803   4.5380   4.5380   4.5864   4.6248
     4.7872   4.8084   4.8084   4.9065   4.9065   4.9083   5.2211   5.2211
     5.3824   5.3824   5.8701   5.9054   5.9437   5.9437   8.1743   8.1743
     8.1829   8.2189   8.2640   8.2640   8.6833   8.6940   8.6940   8.7248
     9.7076   9.7076  11.9134  11.9221

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.2000 0.2000 0.4000 (  6315 PWs)   bands (ev):

   -48.7711 -48.7711 -26.3908 -26.3908 -24.3087 -24.3087 -24.3073 -24.3073
    -6.6165  -6.6164  -6.2785  -6.2698  -6.0976  -6.0890  -6.0751  -6.0711
     1.1290   1.1366   1.3096   1.3116   1.6708   1.6709   2.0597   2.0675
     2.3968   2.4163   2.4312   2.4335   2.5150   2.5244   2.6863   2.6887
     2.8429   2.8623   2.8887   2.8909   3.2867   3.3091   3.3222   3.3273
     3.7959   3.7979   4.0304   4.0505   4.0946   4.1221   4.1612   4.2435
     4.2494   4.3329   4.3435   4.3443   4.5482   4.5727   4.6000   4.6503
     4.6837   4.7032   4.7674   4.8194   4.9146   4.9421   5.0742   5.1919
     5.5220   5.5529   5.7621   5.7889   5.8405   5.8506   8.0658   8.0768
     8.2403   8.2724   8.4949   8.5096   8.7679   8.7732   8.9168   8.9501
    10.2989  10.3001  11.5310  11.5400

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.2000 0.4000 0.4000 (  6280 PWs)   bands (ev):

   -48.7709 -48.7709 -26.3913 -26.3913 -24.3093 -24.3093 -24.3080 -24.3080
    -6.4707  -6.4707  -6.2383  -6.2315  -6.2124  -6.2077  -6.1302  -6.1202
     1.2462   1.2551   1.3482   1.3526   1.5038   1.5050   2.1046   2.1141
     2.4313   2.4427   2.4529   2.4551   2.5838   2.5905   2.6036   2.6132
     2.8304   2.8407   2.8565   2.8692   3.1355   3.1520   3.1929   3.2031
     3.8648   3.8727   4.0538   4.0944   4.1348   4.1508   4.1880   4.2136
     4.2754   4.2979   4.3377   4.3489   4.4175   4.5098   4.6072   4.6371
     4.6936   4.7344   4.7958   4.8029   4.8728   4.8830   4.9808   4.9920
     5.6935   5.7143   5.7531   5.7568   5.9778   5.9841   8.1194   8.1308
     8.2305   8.2352   8.6520   8.6801   8.9832   9.0109   9.0122   9.0405
    10.5182  10.5188  11.1621  11.1662

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.4000 0.4000 0.4000 (  6299 PWs)   bands (ev):

   -48.7709 -48.7709 -26.3918 -26.3918 -24.3094 -24.3094 -24.3094 -24.3094
    -6.3650  -6.3650  -6.2359  -6.2359  -6.2132  -6.2113  -6.1916  -6.1916
     1.2630   1.2630   1.2803   1.2900   1.3587   1.3587   1.9279   1.9279
     2.4973   2.5097   2.5097   2.5126   2.6718   2.6731   2.6898   2.6898
     2.7612   2.7612   2.9771   2.9771   3.1041   3.1166   3.1166   3.1224
     3.9599   3.9599   4.0692   4.1235   4.1694   4.1694   4.2313   4.2313
     4.2579   4.3268   4.3268   4.3535   4.4701   4.5027   4.5027   4.5637
     4.5776   4.5776   4.7067   4.7067   4.8032   4.8320   4.8320   4.8637
     6.0272   6.0335   6.0518   6.0518   6.1000   6.1000   8.1501   8.1619
     8.1619   8.1759   8.9018   8.9018   9.0883   9.0984   9.0984   9.1161
    10.3243  10.3243  10.7203  10.7329

     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
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   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   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000

     the Fermi energy is     6.3866 ev

!    total energy              =    -570.35240314 Ry
     Harris-Foulkes estimate   =    -570.35240314 Ry
     estimated scf accuracy    <          5.7E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =    -308.33327403 Ry
     hartree contribution      =     210.64594463 Ry
     xc contribution           =    -172.71282275 Ry
     ewald contribution        =    -299.95225099 Ry
     smearing contrib. (-TS)   =      -0.00000000 Ry

     convergence has been achieved in  12 iterations

     Writing output data file NbCu3S4.save
 
     init_run     :     11.01s CPU     23.72s WALL (       1 calls)
     electrons    :     88.02s CPU     91.51s WALL (       1 calls)

     Called by init_run:
     wfcinit      :      2.24s CPU      3.47s WALL (       1 calls)
     potinit      :      0.21s CPU      1.46s WALL (       1 calls)

     Called by electrons:
     c_bands      :     68.30s CPU     69.62s WALL (      13 calls)
     sum_band     :     12.32s CPU     12.72s WALL (      13 calls)
     v_of_rho     :      0.23s CPU      1.03s WALL (      13 calls)
     v_h          :      0.06s CPU      0.06s WALL (      13 calls)
     v_xc         :      0.16s CPU      0.74s WALL (      13 calls)
     newd         :      7.45s CPU      7.63s WALL (      13 calls)
     mix_rho      :      0.22s CPU      1.14s WALL (      13 calls)

     Called by c_bands:
     init_us_2    :      0.14s CPU      0.22s WALL (     270 calls)
     cegterg      :     64.70s CPU     65.86s WALL (     130 calls)

     Called by sum_band:
     sum_band:bec :      1.71s CPU      1.72s WALL (     130 calls)
     addusdens    :      2.68s CPU      2.68s WALL (      13 calls)

     Called by *egterg:
     h_psi        :     32.13s CPU     33.24s WALL (     454 calls)
     s_psi        :      6.52s CPU      6.63s WALL (     454 calls)
     g_psi        :      0.05s CPU      0.06s WALL (     314 calls)
     cdiaghg      :     15.33s CPU     15.29s WALL (     434 calls)
     cegterg:over :      5.58s CPU      5.51s WALL (     314 calls)
     cegterg:upda :      1.16s CPU      1.24s WALL (     314 calls)
     cegterg:last :      0.69s CPU      0.73s WALL (     130 calls)

     Called by h_psi:
     h_psi:vloc   :     20.31s CPU     20.58s WALL (     454 calls)
     h_psi:vnl    :     11.76s CPU     12.56s WALL (     454 calls)
     add_vuspsi   :      4.89s CPU      5.04s WALL (     454 calls)

     General routines
     calbec       :      9.56s CPU     10.18s WALL (     584 calls)
     fft          :      0.69s CPU      1.92s WALL (     397 calls)
     ffts         :      0.08s CPU      0.09s WALL (     104 calls)
     fftw         :     23.32s CPU     23.40s WALL (  131300 calls)
     interpolate  :      0.16s CPU      0.19s WALL (     104 calls)
 
     Parallel routines
     fft_scatter  :     16.24s CPU     16.17s WALL (  131801 calls)
 
     PWSCF        :  1m47.98s CPU     2m24.11s WALL

 
   This run was terminated on:  19: 4:19  28Oct2015            

=------------------------------------------------------------------------------=
   JOB DONE.
=------------------------------------------------------------------------------=