[Pw_forum] an abnormal error

nisihara225 at gmail.com nisihara225 at gmail.com
Mon Nov 14 15:04:20 CET 2016


Dear Paolo Giannozzi and Shaofeng Wang,

I think this error occurs in subroutine dprojdtau_gamma (PW/src/force_hub.f90).
At the line 557, DGEMM is called:
      CALL DGEMM('T','N',ldim, nbnd, 2*npw, 2.0_dp,  &
                  dwfc, 2*npwx, spsi, 2*npwx, 0.0_dp,&
                  dproj0, ldim)
However, parameters ‘dwfc’ and ‘spsi’ is not real number but complex number.
It seems this is the reason of the error.

Regards,
Satomichi Nishihara


差出人: Paolo Giannozzi
送信日時: 2016年11月14日 2:52
宛先: Shaofeng Wang; PWSCF Forum
件名: Re: [Pw_forum] an abnormal error

If the error message is printed between "convergence has been achieved in  29 iterations" and "Forces acting on atoms", it may come only from the calls to DGEMM in PW/src/force_us.f90 or (more likely) PW/src/force_hub.f90. Locate those DGEMM, print argument n.10. There might be something funny in your pseudopotentials.
Paolo

On Fri, Nov 11, 2016 at 9:15 AM, Shaofeng Wang <wangshaofeng at iae.ac.cn> wrote:
Dear expert,

I met an error "MKL ERROR: Parameter 10 was incorrect on entry to DGEMM ."
when I trying to optimize a lepidocrocite structure. Even though, the
program did not stop and still run correctly, but this error showed at every
scf cycle. The following is my input and output file. Could anyone help me
to see what has happened?

&CONTROL
                 calculation = 'relax' ,
                      outdir = './tmp' ,
                  pseudo_dir = '../pseudo/ncpp' ,
                      prefix = 'vc' ,
                 etot_conv_thr = 1.0D-4 ,
                 forc_conv_thr = 1.0d-3 ,
                     tstress = .true. ,
                     tprnfor = .true. ,
                        nstep = 150 ,
/
&SYSTEM
                       ibrav = 0,
                    celldm(1) = 1.889726,
                         nat = 94,
                        ntyp = 5,
                     ecutwfc = 70 ,
!                    ecutrho = 400 ,
                        input_dft = pw91 ,
                        vdw_corr = 'DFT-D' ,
                  occupations = 'smearing',
                     smearing = 'mp',
                      degauss = 0.02,
                          nspin = 2 ,
                       starting_magnetization(1) = 0.5 ,
                       starting_magnetization(2) = -0.5 ,
!                         tot_magnetization = 0 ,
              lda_plus_u=.true.  Hubbard_U(1)=5, Hubbard_U(2)=5,
/
&ELECTRONS
            electron_maxstep = 100,
                    conv_thr = 1.0d-8 ,
                 mixing_beta = 0.3 ,
/
&IONS
                ion_dynamics = 'bfgs' ,
/
&CELL
               cell_dynamics = 'bfgs' ,
                   cell_dofree = xyz ,
/
ATOMIC_SPECIES
    Fe   55.845     Fe_00LDA_OP.ncpp  ! Fe_00PBE_OP.ncpp  !
    Fe1  55.845     Fe_00LDA_OP.ncpp  ! Fe_00PBE_OP.ncpp  !
    O   15.99940    O_00PBE.ncpp  !    O_00LDA_OP.ncpp   !
    H   1.00794     H.pbe-hgh.UPF  !     H_00LDA_OP.ncpp  !  H.pz-vbc.UPF  !
H.pz-hgh.UPF    !
   As   74.92160    As.pz-hgh.UPF  !     As.pbe-hgh.UPF  !  As_lda.ncpp  !
CELL_PARAMETERS (alat=  1.88972600)
      12.368930000000004       0.000000000000000       0.000000000000000
       0.000000000000000       7.823999999999999       0.000000000000000
       0.000000000000000       0.000000000000000       9.262440000000003

ATOMIC_POSITIONS (crystal)
  H   0.0151640000000000   0.4738560000000008   0.1665000000000002
  H   0.4840120000000009   0.0265040000000001   0.0002243333333333
  H   0.9845150000000018   0.2240600000000000   0.1666423333333336
  H   0.0151640000000000   0.9738560000000009   0.1665000000000002
  H   0.4840120000000009   0.5265040000000001   0.0002243333333333
  H   0.5158370000000003   0.7766940000000006   0.3330460000000002
  H   0.9845150000000018   0.7240600000000000   0.1666423333333336
  H   0.0151640000000000   0.4738560000000008   0.4998333333333336
  H   0.4840120000000009   0.0265040000000001   0.3335576666666668
  H   0.9845150000000018   0.2240600000000000   0.4999756666666670
  H   0.0151640000000000   0.9738560000000009   0.4998333333333336
  H   0.4840120000000009   0.5265040000000001   0.3335576666666668
  H   0.5158370000000003   0.7766940000000006   0.6663793333333337
  H   0.9845150000000018   0.7240600000000000   0.4999756666666670
  H   0.0151640000000000   0.4738560000000008   0.8331666666666672
  H   0.4840120000000009   0.0265040000000001   0.6668910000000002
  H   0.5158370000000003   0.2766940000000007   0.9997126666666671
  H   0.9845150000000018   0.2240600000000000   0.8333090000000004
  H   0.0151640000000000   0.9738560000000009   0.8331666666666672
  H   0.4840120000000009   0.5265040000000001   0.6668910000000002
  H   0.5158370000000003   0.7766940000000006   0.9997126666666671
  H   0.9845150000000018   0.7240600000000000   0.8333090000000004
  O   0.2097940000000006   0.1255450000000000   0.0000370000000000
  O   0.4317220000000011   0.1261805000000004   0.0001056666666666
  O   0.7112330000000000   0.1245665000000006   0.1664996666666668
  O   0.9325770000000020   0.1236255000000004   0.1666353333333334
  O   0.7902350000000022   0.3752430000000010   0.3331590000000004
  O   0.6108676311362189   0.3923812510369004   0.3664960629839341
  O   0.2887220000000009   0.3748540000000010   0.1666983333333335
  O   0.0671390000000002   0.3735600000000011   0.1665803333333334
  O   0.2097940000000006   0.6255450000000000   0.0000370000000000
  O   0.4317220000000011   0.6261805000000006   0.0001056666666666
  O   0.7112330000000000   0.6245665000000008   0.1664996666666668
  O   0.9325770000000020   0.6236255000000006   0.1666353333333334
  O   0.7902350000000022   0.8752430000000011   0.3331590000000004
  O   0.5681570000000016   0.8762365000000006   0.3331353333333336
  O   0.2887220000000009   0.8748540000000011   0.1666983333333335
  O   0.0671390000000002   0.8735600000000012   0.1665803333333334
  O   0.2097940000000006   0.1255450000000000   0.3333703333333334
  O   0.4317220000000011   0.1261805000000004   0.3334390000000001
  O   0.7112330000000000   0.1245665000000006   0.4998330000000003
  O   0.9325770000000020   0.1236255000000004   0.4999686666666670
  O   0.7902350000000022   0.3752430000000010   0.6664923333333338
  O   0.6066537996602838   0.3679112890129118   0.6208237138383496
  O   0.2887220000000009   0.3748540000000010   0.5000316666666670
  O   0.0671390000000002   0.3735600000000011   0.4999136666666668
  O   0.2097940000000006   0.6255450000000000   0.3333703333333334
  O   0.4317220000000011   0.6261805000000006   0.3334390000000001
  O   0.7112330000000000   0.6245665000000008   0.4998330000000003
  O   0.9325770000000020   0.6236255000000006   0.4999686666666670
  O   0.7902350000000022   0.8752430000000011   0.6664923333333338
  O   0.5681570000000016   0.8762365000000006   0.6664686666666670
  O   0.2887220000000009   0.8748540000000011   0.5000316666666670
  O   0.0671390000000002   0.8735600000000012   0.4999136666666668
...............
Fe   0.6746790000000003   0.8745660000000006   0.4998173333333336
Fe1   0.8268580000000024   0.1250805000000003   0.9997880000000007
Fe   0.3257760000000003   0.1250645000000000   0.8334123333333336
Fe1   0.1735250000000001   0.3755685000000005   0.9999753333333339
Fe   0.6746790000000003   0.3745660000000005   0.8331506666666670
Fe1   0.8268580000000024   0.6250805000000003   0.9997880000000007
Fe   0.3257760000000003   0.6250645000000000   0.8334123333333336
Fe1   0.1735250000000001   0.8755685000000006   0.9999753333333339
Fe   0.6746790000000003   0.8745660000000006   0.8331506666666670
As   0.7077822966962698   0.4210551096941716   0.4959948415593733
K_POINTS gamma

The output is
     Program PWSCF v.6.0 (svn rev. 13079) starts on 11Nov2016 at 15:51:19

     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    12 processors
     R & G space division:  proc/nbgrp/npool/nimage =      12
     Waiting for input...
     Reading input from standard input
Warning: card  &CELL ignored
Warning: card                CELL_DYNAMICS = 'BFGS' , ignored
Warning: card                    CELL_DOFREE = XYZ , 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

     IMPORTANT: XC functional enforced from input :
     Exchange-correlation      = PW91 ( 1  4  2  2 0 0)
     Any further DFT definition will be discarded
     Please, verify this is what you really want


     -------------------------------------
     Parameters for Dispersion Correction:
     -------------------------------------
       atom      VdW radius       C_6

        Fe         2.952        374.666
        Fe1        2.952        374.666
        O          2.536         24.284
        H          1.892          4.857
        As         3.326        567.896

     gamma-point specific algorithms are used

     Subspace diagonalization in iterative solution of the eigenvalue
problem:
     one sub-group per band group will be used
     scalapack distributed-memory algorithm (size of sub-group:  2*  2
procs)


     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min         641     641    160                39886    39886    4992
     Max         644     644    162                39896    39896    4998
     Sum        7711    7711   1925               478687   478687   59919



     bravais-lattice index     =            0
     lattice parameter (alat)  =       1.8897  a.u.
     unit-cell volume          =    6048.9904 (a.u.)^3
     number of atoms/cell      =           94
     number of atomic types    =            5
     number of electrons       =       683.00
     number of Kohn-Sham states=          410
     kinetic-energy cutoff     =      70.0000  Ry
     charge density cutoff     =     280.0000  Ry
     convergence threshold     =      1.0E-08
     mixing beta               =       0.3000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      = PW91 ( 1  4  2  2 0 0)
     nstep                     =          150


     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) = (  12.368930   0.000000   0.000000 )
               a(2) = (   0.000000   7.824000   0.000000 )
               a(3) = (   0.000000   0.000000   9.262440 )

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  0.080848  0.000000  0.000000 )
               b(2) = (  0.000000  0.127812  0.000000 )
               b(3) = (  0.000000  0.000000  0.107963 )



     atomic species   valence    mass     pseudopotential
        Fe            16.00    55.84500     Fe( 1.00)
        Fe1           16.00    55.84500     Fe( 1.00)
        O              6.00    15.99940     O ( 1.00)
        H              1.00     1.00794      H( 1.00)
        As             5.00    74.92160     As( 1.00)

     Starting magnetic structure
     atomic species   magnetization
        Fe           0.500
        Fe1         -0.500
        O            0.000
        H            0.000
        As           0.000


     Simplified LDA+U calculation (l_max = 2) with parameters (eV):
     atomic species    L          U    alpha       J0     beta
        Fe             2     5.0000   0.0000   0.0000   0.0000
        Fe1            2     5.0000   0.0000   0.0000   0.0000



     No symmetry found



   Cartesian axes

     site n.     atom                  positions (alat units)
         1           H   tau(   1) = (   0.1875625   3.7074493
 1.5421963  )
         2           H   tau(   2) = (   5.9867105   0.2073673
 0.0020779  )
         3           H   tau(   3) = (  12.1773971   1.7530454
 1.5435146  )
         4           H   tau(   4) = (   0.1875625   7.6194493
 1.5421963  )
         5           H   tau(   5) = (   5.9867105   4.1193673
 0.0020779  )
         6           H   tau(   6) = (   6.3803517   6.0768539
 3.0848186  )
         7           H   tau(   7) = (  12.1773971   5.6650454
 1.5435146  )
         8           H   tau(   8) = (   0.1875625   3.7074493
 4.6296763  )
         9           H   tau(   9) = (   5.9867105   0.2073673
 3.0895579  )
        10           H   tau(  10) = (  12.1773971   1.7530454
 4.6309946  )
        11           H   tau(  11) = (   0.1875625   7.6194493
 4.6296763  )
        12           H   tau(  12) = (   5.9867105   4.1193673
 3.0895579  )
        13           H   tau(  13) = (   6.3803517   6.0768539
 6.1722986  )
        14           H   tau(  14) = (  12.1773971   5.6650454
 4.6309946  )
        15           H   tau(  15) = (   0.1875625   3.7074493
 7.7171563  )
        16           H   tau(  16) = (   5.9867105   0.2073673
 6.1770379  )
        17           H   tau(  17) = (   6.3803517   2.1648539
 9.2597786  )
        18           H   tau(  18) = (  12.1773971   1.7530454
 7.7184746  )
        19           H   tau(  19) = (   0.1875625   7.6194493
 7.7171563  )
        20           H   tau(  20) = (   5.9867105   4.1193673
 6.1770379  )
        21           H   tau(  21) = (   6.3803517   6.0768539
 9.2597786  )
        22           H   tau(  22) = (  12.1773971   5.6650454
 7.7184746  )
        23           O   tau(  23) = (   2.5949273   0.9822641
 0.0003427  )
        24           O   tau(  24) = (   5.3399392   0.9872362
 0.0009787  )
        25           O   tau(  25) = (   8.7971912   0.9746083
 1.5421932  )
        26           O   tau(  26) = (  11.5349796   0.9672459
 1.5434498  )
        27           O   tau(  27) = (   9.7743614   2.9359012
 3.0858652  )
        28           O   tau(  28) = (   7.5557790   3.0699909
 3.3946478  )
        29           O   tau(  29) = (   3.5711822   2.9328577
 1.5440333  )
        30           O   tau(  30) = (   0.8304376   2.9227334
 1.5429403  )
        31           O   tau(  31) = (   2.5949273   4.8942641
 0.0003427  )
        32           O   tau(  32) = (   5.3399392   4.8992362
 0.0009787  )
        33           O   tau(  33) = (   8.7971912   4.8866083
 1.5421932  )
        34           O   tau(  34) = (  11.5349796   4.8792459
 1.5434498  )
        35           O   tau(  35) = (   9.7743614   6.8479012
 3.0858652  )
        36           O   tau(  36) = (   7.0274942   6.8556744
 3.0856460  )
        37           O   tau(  37) = (   3.5711822   6.8448577
 1.5440333  )
        38           O   tau(  38) = (   0.8304376   6.8347334
 1.5429403  )
        39           O   tau(  39) = (   2.5949273   0.9822641
 3.0878227  )
        40           O   tau(  40) = (   5.3399392   0.9872362
 3.0884587  )
        41           O   tau(  41) = (   8.7971912   0.9746083
 4.6296732  )
        42           O   tau(  42) = (  11.5349796   0.9672459
 4.6309298  )
        43           O   tau(  43) = (   9.7743614   2.9359012
 6.1733452  )
        44           O   tau(  44) = (   7.5036584   2.8785379
 5.7503424  )
        45           O   tau(  45) = (   3.5711822   2.9328577
 4.6315133  )
        46           O   tau(  46) = (   0.8304376   2.9227334
 4.6304203  )
        47           O   tau(  47) = (   2.5949273   4.8942641
 3.0878227  )
        48           O   tau(  48) = (   5.3399392   4.8992362
 3.0884587  )
        49           O   tau(  49) = (   8.7971912   4.8866083
 4.6296732  )
        50           O   tau(  50) = (  11.5349796   4.8792459
 4.6309298  )
        51           O   tau(  51) = (   9.7743614   6.8479012
 6.1733452  )
        52           O   tau(  52) = (   7.0274942   6.8556744
 6.1731260  )
        53           O   tau(  53) = (   3.5711822   6.8448577
 4.6315133  )
        54           O   tau(  54) = (   0.8304376   6.8347334
 4.6304203  )
        55           O   tau(  55) = (   2.5949273   0.9822641
 6.1753027  )
        56           O   tau(  56) = (   5.3399392   0.9872362
 6.1759387  )
        57           O   tau(  57) = (   8.7971912   0.9746083
 7.7171532  )
        58           O   tau(  58) = (  11.5349796   0.9672459
 7.7184098  )
        59           O   tau(  59) = (   9.7743614   2.9359012
 9.2608252  )
        60           O   tau(  60) = (   7.0274942   2.9436744
 9.2606060  )
        61           O   tau(  61) = (   3.5711822   2.9328577
 7.7189933  )
        62           O   tau(  62) = (   0.8304376   2.9227334
 7.7179003  )
        63           O   tau(  63) = (   2.5949273   4.8942641
 6.1753027  )
        64           O   tau(  64) = (   5.3399392   4.8992362
 6.1759387  )
        65           O   tau(  65) = (   8.7971912   4.8866083
 7.7171532  )
        66           O   tau(  66) = (  11.5349796   4.8792459
 7.7184098  )
        67           O   tau(  67) = (   9.7743614   6.8479012
 9.2608252  )
        68           O   tau(  68) = (   7.0274942   6.8556744
 9.2606060  )
        69           O   tau(  69) = (   3.5711822   6.8448577
 7.7189933  )
        70           O   tau(  70) = (   0.8304376   6.8347334
 7.7179003  )
        71           Fe  tau(  71) = (  10.2273487   0.9786298
 3.0855164  )
        72           Fe1 tau(  72) = (   4.0295005   0.9785046
 1.5444717  )
        73           Fe  tau(  73) = (   2.1463186   2.9384479
 3.0872515  )
        74           Fe1 tau(  74) = (   8.3450573   2.9306044
 1.5420481  )
        75           Fe  tau(  75) = (  10.2273487   4.8906298
 3.0855164  )
        76           Fe1 tau(  76) = (   4.0295005   4.8905046
 1.5444717  )
        77           Fe  tau(  77) = (   2.1463186   6.8504479
 3.0872515  )
        78           Fe1 tau(  78) = (   8.3450573   6.8426044
 1.5420481  )
        79           Fe  tau(  79) = (  10.2273487   0.9786298
 6.1729964  )
        80           Fe1 tau(  80) = (   4.0295005   0.9785046
 4.6319517  )
        81           Fe  tau(  81) = (   2.1463186   2.9384479
 6.1747315  )
        82           Fe1 tau(  82) = (  10.2273487   4.8906298
 6.1729964  )
        83           Fe  tau(  83) = (   4.0295005   4.8905046
 4.6319517  )
        84           Fe1 tau(  84) = (   2.1463186   6.8504479
 6.1747315  )
        85           Fe  tau(  85) = (   8.3450573   6.8426044
 4.6295281  )
        86           Fe1 tau(  86) = (  10.2273487   0.9786298
 9.2604764  )
        87           Fe  tau(  87) = (   4.0295005   0.9785046
 7.7194317  )
        88           Fe1 tau(  88) = (   2.1463186   2.9384479
 9.2622115  )
        89           Fe  tau(  89) = (   8.3450573   2.9306044
 7.7170081  )
        90           Fe1 tau(  90) = (  10.2273487   4.8906298
 9.2604764  )
        91           Fe  tau(  91) = (   4.0295005   4.8905046
 7.7194317  )
        92           Fe1 tau(  92) = (   2.1463186   6.8504479
 9.2622115  )
        93           Fe  tau(  93) = (   8.3450573   6.8426044
 7.7170081  )
        94           As  tau(  94) = (   8.7545097   3.2943352
 4.5941225  )

     number of k points=     1  Methfessel-Paxton smearing, width (Ry)=
0.0200
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   1.0000000

     Dense  grid:   239344 G-vectors     FFT dimensions: ( 125,  80,  96)

     Estimated max dynamical RAM per process >     161.68Mb

     Estimated total allocated dynamical RAM >    1940.12Mb
     Generating pointlists ...
     new r_m :   0.4125 (alat units)  0.7795 (a.u.) for type    1
     new r_m :   0.4125 (alat units)  0.7795 (a.u.) for type    2
     new r_m :   0.4125 (alat units)  0.7795 (a.u.) for type    3
     new r_m :   0.4125 (alat units)  0.7795 (a.u.) for type    4
     new r_m :   0.4125 (alat units)  0.7795 (a.u.) for type    5

     Initial potential from superposition of free atoms

     starting charge  635.99972, renormalised to  683.00000
     Number of +U iterations with fixed ns =  0
     Starting occupations:
--- enter write_ns ---
LDA+U parameters:
U( 1)     =  5.00000000
alpha( 1) =  0.00000000
U( 2)     =  5.00000000
....................



0.943 -0.009  0.001 -0.004  0.001
-0.009  0.958  0.001 -0.009  0.000
  0.001  0.001  0.946 -0.002  0.000
-0.004 -0.009 -0.002  0.932  0.000
  0.001  0.000  0.000  0.000  0.953
   spin  2
    eigenvalues:
  0.026  0.027  0.030  0.085  0.086
    eigenvectors:
  0.559  0.018  0.078  0.346  0.000
  0.000  0.000  0.000  0.000  0.999
  0.017  0.969  0.013  0.000  0.000
  0.291  0.009  0.045  0.654  0.000
  0.133  0.004  0.864  0.000  0.000
    occupations:
  0.047  0.000  0.000  0.028  0.001
  0.000  0.086  0.001  0.000  0.000
  0.000  0.001  0.027  0.000  0.000
  0.028  0.000  0.000  0.065  0.000
  0.001  0.000  0.000  0.000  0.029
atomic mag. moment =     4.477510
N of occupied +U levels =  112.611038
--- exit write_ns ---

     total cpu time spent up to now is       34.0 secs

     total energy              =   -7068.43154138 Ry
     Harris-Foulkes estimate   =   -7115.00280712 Ry
     estimated scf accuracy    <     121.70212063 Ry

     total magnetization       =     1.70 Bohr mag/cell
     absolute magnetization    =   107.49 Bohr mag/cell

     iteration #  2     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  4.5

     total cpu time spent up to now is       55.3 secs

     total energy              =   -7072.71471679 Ry
     Harris-Foulkes estimate   =   -7085.37420668 Ry
     estimated scf accuracy    <      47.91542590 Ry

     total magnetization       =    -2.76 Bohr mag/cell
     absolute magnetization    =   104.21 Bohr mag/cell

     iteration #  3     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  7.02E-03,  avg # of iterations =  4.5

     total cpu time spent up to now is       68.6 secs

     total energy              =   -7066.24788223 Ry
     Harris-Foulkes estimate   =   -7092.17327634 Ry
     estimated scf accuracy    <     239.18012931 Ry

     total magnetization       =    -4.18 Bohr mag/cell
     absolute magnetization    =    79.89 Bohr mag/cell

     iteration #  4     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  7.02E-03,  avg # of iterations =  3.5

     total cpu time spent up to now is       80.1 secs

     total energy              =   -7079.10076803 Ry
     Harris-Foulkes estimate   =   -7086.60938132 Ry
     estimated scf accuracy    <      71.74959017 Ry

     total magnetization       =     0.12 Bohr mag/cell
     absolute magnetization    =   100.81 Bohr mag/cell

     iteration #  5     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  7.02E-03,  avg # of iterations =  3.0

     total cpu time spent up to now is       90.7 secs

     total energy              =   -7083.86123332 Ry
     Harris-Foulkes estimate   =   -7085.01672264 Ry
     estimated scf accuracy    <       6.72771422 Ry

     total magnetization       =     4.84 Bohr mag/cell
     absolute magnetization    =   110.44 Bohr mag/cell

     iteration #  6     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  9.85E-04,  avg # of iterations =  6.0

     total cpu time spent up to now is      107.2 secs

     total energy              =   -7084.50665628 Ry
     Harris-Foulkes estimate   =   -7084.47878355 Ry
     estimated scf accuracy    <       1.99977621 Ry

     total magnetization       =     4.67 Bohr mag/cell
     absolute magnetization    =   107.59 Bohr mag/cell

     iteration #  7     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  2.93E-04,  avg # of iterations =  3.5

     total cpu time spent up to now is      119.1 secs

     total energy              =   -7084.59409710 Ry
     Harris-Foulkes estimate   =   -7084.92314947 Ry
     estimated scf accuracy    <       4.54490680 Ry

     total magnetization       =     1.73 Bohr mag/cell
     absolute magnetization    =   105.64 Bohr mag/cell

     iteration #  8     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  2.93E-04,  avg # of iterations =  3.0

     total cpu time spent up to now is      129.7 secs

     total energy              =   -7084.99390279 Ry
     Harris-Foulkes estimate   =   -7084.95654173 Ry
     estimated scf accuracy    <       0.55976201 Ry

     total magnetization       =     5.09 Bohr mag/cell
     absolute magnetization    =   107.02 Bohr mag/cell

     iteration #  9     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  8.20E-05,  avg # of iterations =  3.5

     total cpu time spent up to now is      141.0 secs

     total energy              =   -7085.10524491 Ry
     Harris-Foulkes estimate   =   -7085.01995111 Ry
     estimated scf accuracy    <       0.23858404 Ry

     total magnetization       =     5.02 Bohr mag/cell
     absolute magnetization    =   106.77 Bohr mag/cell

     iteration # 10     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.49E-05,  avg # of iterations =  4.0

     total cpu time spent up to now is      152.7 secs

     total energy              =   -7085.14413492 Ry
     Harris-Foulkes estimate   =   -7085.13041918 Ry
     estimated scf accuracy    <       0.06394377 Ry

     total magnetization       =     4.91 Bohr mag/cell
     absolute magnetization    =   106.75 Bohr mag/cell

     iteration # 11     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  9.36E-06,  avg # of iterations =  8.0

     total cpu time spent up to now is      168.7 secs

     total energy              =   -7085.15604888 Ry
     Harris-Foulkes estimate   =   -7085.15083256 Ry
     estimated scf accuracy    <       0.02424419 Ry

     total magnetization       =     4.94 Bohr mag/cell
     absolute magnetization    =   106.77 Bohr mag/cell

     iteration # 12     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.55E-06,  avg # of iterations =  4.0

     total cpu time spent up to now is      180.5 secs

     total energy              =   -7085.16032753 Ry
     Harris-Foulkes estimate   =   -7085.15916698 Ry
     estimated scf accuracy    <       0.01152866 Ry

     total magnetization       =     4.91 Bohr mag/cell
     absolute magnetization    =   106.76 Bohr mag/cell

     iteration # 13     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.69E-06,  avg # of iterations =  4.0

     total cpu time spent up to now is      192.2 secs

     total energy              =   -7085.16255976 Ry
     Harris-Foulkes estimate   =   -7085.16212065 Ry
     estimated scf accuracy    <       0.00539513 Ry

     total magnetization       =     4.94 Bohr mag/cell
     absolute magnetization    =   106.83 Bohr mag/cell

     iteration # 14     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  7.90E-07,  avg # of iterations =  3.0

     total cpu time spent up to now is      202.8 secs

     total energy              =   -7085.16339835 Ry
     Harris-Foulkes estimate   =   -7085.16332415 Ry
     estimated scf accuracy    <       0.00198389 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.79 Bohr mag/cell

     iteration # 15     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  2.90E-07,  avg # of iterations =  2.0

     total cpu time spent up to now is      212.5 secs

     total energy              =   -7085.16387928 Ry
     Harris-Foulkes estimate   =   -7085.16363161 Ry
     estimated scf accuracy    <       0.00055735 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.79 Bohr mag/cell

     iteration # 16     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  8.16E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is      222.4 secs

     total energy              =   -7085.16401376 Ry
     Harris-Foulkes estimate   =   -7085.16399335 Ry
     estimated scf accuracy    <       0.00025799 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.79 Bohr mag/cell

     iteration # 17     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.78E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is      232.3 secs

     total energy              =   -7085.16396494 Ry
     Harris-Foulkes estimate   =   -7085.16404675 Ry
     estimated scf accuracy    <       0.00008559 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.79 Bohr mag/cell

     iteration # 18     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.25E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is      243.9 secs

     total energy              =   -7085.16397259 Ry
     Harris-Foulkes estimate   =   -7085.16399875 Ry
     estimated scf accuracy    <       0.00002734 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 19     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  4.00E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is      253.6 secs

     total energy              =   -7085.16396086 Ry
     Harris-Foulkes estimate   =   -7085.16397631 Ry
     estimated scf accuracy    <       0.00001122 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 20     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.64E-09,  avg # of iterations =  2.5

     total cpu time spent up to now is      264.1 secs

     total energy              =   -7085.16396766 Ry
     Harris-Foulkes estimate   =   -7085.16396297 Ry
     estimated scf accuracy    <       0.00000553 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 21     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  8.09E-10,  avg # of iterations =  2.5

     total cpu time spent up to now is      274.2 secs

     total energy              =   -7085.16397199 Ry
     Harris-Foulkes estimate   =   -7085.16396894 Ry
     estimated scf accuracy    <       0.00000240 Ry

     total magnetization       =     4.92 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 22     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.51E-10,  avg # of iterations =  3.0

     total cpu time spent up to now is      287.1 secs

     total energy              =   -7085.16397997 Ry
     Harris-Foulkes estimate   =   -7085.16397315 Ry
     estimated scf accuracy    <       0.00000074 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 23     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.09E-10,  avg # of iterations =  2.0

     total cpu time spent up to now is      296.8 secs

     total energy              =   -7085.16398438 Ry
     Harris-Foulkes estimate   =   -7085.16398014 Ry
     estimated scf accuracy    <       0.00000048 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 24     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  7.08E-11,  avg # of iterations =  2.5

     total cpu time spent up to now is      308.2 secs

     total energy              =   -7085.16398573 Ry
     Harris-Foulkes estimate   =   -7085.16398460 Ry
     estimated scf accuracy    <       0.00000015 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 25     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  2.22E-11,  avg # of iterations =  2.5

     total cpu time spent up to now is      318.9 secs

     total energy              =   -7085.16398673 Ry
     Harris-Foulkes estimate   =   -7085.16398578 Ry
     estimated scf accuracy    <       0.00000010 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 26     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.41E-11,  avg # of iterations =  3.0

     total cpu time spent up to now is      329.8 secs

     total energy              =   -7085.16398657 Ry
     Harris-Foulkes estimate   =   -7085.16398676 Ry
     estimated scf accuracy    <       0.00000004 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 27     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  6.40E-12,  avg # of iterations =  2.5

     total cpu time spent up to now is      340.8 secs

     total energy              =   -7085.16398619 Ry
     Harris-Foulkes estimate   =   -7085.16398658 Ry
     estimated scf accuracy    <       0.00000002 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 28     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  3.03E-12,  avg # of iterations =  2.5

     total cpu time spent up to now is      351.0 secs

     total energy              =   -7085.16398585 Ry
     Harris-Foulkes estimate   =   -7085.16398620 Ry
     estimated scf accuracy    <       0.00000001 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     iteration # 29     ecut=    70.00 Ry     beta=0.30
     Davidson diagonalization with overlap
     ethr =  1.86E-12,  avg # of iterations =  2.5

     Magnetic moment per site:
     atom:    1    charge:    0.3840    magn:    0.0000    constr:    0.0000
     atom:    2    charge:    0.3834    magn:    0.0001    constr:    0.0000
     atom:    3    charge:    0.3816    magn:    0.0001    constr:    0.0000
     atom:    4    charge:    0.3832    magn:    0.0001    constr:    0.0000
     atom:    5    charge:    0.3831    magn:    0.0000    constr:    0.0000
     atom:    6    charge:    0.3833    magn:    0.0001    constr:    0.0000
     atom:    7    charge:    0.3806    magn:    0.0001    constr:    0.0000
     atom:    8    charge:    0.3848    magn:    0.0027    constr:    0.0000
     atom:    9    charge:    0.3849    magn:   -0.0032    constr:    0.0000
     atom:   10    charge:    0.3818    magn:    0.0030    constr:    0.0000
     atom:   11    charge:    0.3833    magn:    0.0001    constr:    0.0000
     atom:   12    charge:    0.3764    magn:    0.0000    constr:    0.0000
     atom:   13    charge:    0.3834    magn:    0.0031    constr:    0.0000
     atom:   14    charge:    0.3785    magn:    0.0001    constr:    0.0000
     atom:   15    charge:    0.3839    magn:   -0.0002    constr:    0.0000
     atom:   16    charge:    0.3850    magn:    0.0001    constr:    0.0000
     atom:   17    charge:    0.3815    magn:    0.0000    constr:    0.0000
     atom:   18    charge:    0.3816    magn:    0.0000    constr:    0.0000
     atom:   19    charge:    0.3831    magn:   -0.0027    constr:    0.0000
     atom:   20    charge:    0.3768    magn:    0.0029    constr:    0.0000
     atom:   21    charge:    0.3829    magn:    0.0001    constr:    0.0000
     atom:   22    charge:    0.3807    magn:   -0.0029    constr:    0.0000
     atom:   23    charge:    1.9148    magn:   -0.0759    constr:    0.0000
     atom:   24    charge:    1.9339    magn:   -0.0015    constr:    0.0000
     atom:   25    charge:    1.9137    magn:   -0.0677    constr:    0.0000
     atom:   26    charge:    1.9307    magn:    0.0007    constr:    0.0000
     atom:   27    charge:    1.9297    magn:    0.0499    constr:    0.0000
     atom:   28    charge:    1.9193    magn:   -0.0558    constr:    0.0000
     atom:   29    charge:    1.9144    magn:   -0.0717    constr:    0.0000
     atom:   30    charge:    1.9310    magn:    0.0006    constr:    0.0000
     atom:   31    charge:    1.9147    magn:   -0.0731    constr:    0.0000
     atom:   32    charge:    1.9341    magn:    0.0007    constr:    0.0000
     atom:   33    charge:    1.9144    magn:   -0.0612    constr:    0.0000
     atom:   34    charge:    1.9304    magn:    0.0078    constr:    0.0000
     atom:   35    charge:    1.9147    magn:    0.0836    constr:    0.0000
     atom:   36    charge:    1.9327    magn:    0.0006    constr:    0.0000
     atom:   37    charge:    1.9145    magn:   -0.0722    constr:    0.0000
     atom:   38    charge:    1.9310    magn:   -0.0012    constr:    0.0000
     atom:   39    charge:    1.9137    magn:    0.0118    constr:    0.0000
     atom:   40    charge:    1.9249    magn:   -0.0837    constr:    0.0000
     atom:   41    charge:    1.9030    magn:    0.1599    constr:    0.0000
     atom:   42    charge:    1.9307    magn:    0.0674    constr:    0.0000
     atom:   43    charge:    1.9243    magn:    0.0217    constr:    0.0000
     atom:   44    charge:    1.9078    magn:    0.0728    constr:    0.0000
     atom:   45    charge:    1.9125    magn:    0.0858    constr:    0.0000
     atom:   46    charge:    1.9309    magn:    0.0617    constr:    0.0000
     atom:   47    charge:    1.9141    magn:    0.0685    constr:    0.0000
     atom:   48    charge:    1.9306    magn:   -0.0009    constr:    0.0000
     atom:   49    charge:    1.9609    magn:    0.0204    constr:    0.0000
     atom:   50    charge:    1.9310    magn:    0.0029    constr:    0.0000
     atom:   51    charge:    1.9133    magn:    0.0653    constr:    0.0000
     atom:   52    charge:    1.9327    magn:    0.0701    constr:    0.0000
     atom:   53    charge:    1.9129    magn:    0.0113    constr:    0.0000
     atom:   54    charge:    1.9309    magn:    0.0031    constr:    0.0000
     atom:   55    charge:    1.9132    magn:    0.0101    constr:    0.0000
     atom:   56    charge:    1.9244    magn:    0.0074    constr:    0.0000
     atom:   57    charge:    1.9134    magn:    0.0642    constr:    0.0000
     atom:   58    charge:    1.9307    magn:   -0.0013    constr:    0.0000
     atom:   59    charge:    1.9159    magn:   -0.0740    constr:    0.0000
     atom:   60    charge:    1.9326    magn:   -0.0003    constr:    0.0000
     atom:   61    charge:    1.9137    magn:    0.0645    constr:    0.0000
     atom:   62    charge:    1.9311    magn:   -0.0027    constr:    0.0000
     atom:   63    charge:    1.9130    magn:    0.0746    constr:    0.0000
     atom:   64    charge:    1.9302    magn:    0.0610    constr:    0.0000
     atom:   65    charge:    1.9135    magn:   -0.0094    constr:    0.0000
     atom:   66    charge:    1.9307    magn:   -0.0707    constr:    0.0000
     atom:   67    charge:    1.9154    magn:   -0.0733    constr:    0.0000
     atom:   68    charge:    1.9335    magn:    0.0006    constr:    0.0000
     atom:   69    charge:    1.9135    magn:   -0.0018    constr:    0.0000
     atom:   70    charge:    1.9312    magn:   -0.0620    constr:    0.0000
     atom:   71    charge:    6.3246    magn:    1.4509    constr:    0.0000
     atom:   72    charge:    6.3305    magn:   -1.4671    constr:    0.0000
     atom:   73    charge:    6.3261    magn:    1.4604    constr:    0.0000
     atom:   74    charge:    6.3283    magn:   -1.4426    constr:    0.0000
     atom:   75    charge:    6.3262    magn:    1.4412    constr:    0.0000
     atom:   76    charge:    6.3306    magn:   -1.4654    constr:    0.0000
     atom:   77    charge:    6.3261    magn:    1.4588    constr:    0.0000
     atom:   78    charge:    6.3264    magn:   -1.4652    constr:    0.0000
     atom:   79    charge:    6.3247    magn:    1.4407    constr:    0.0000
     atom:   80    charge:    6.3306    magn:   -1.4527    constr:    0.0000
     atom:   81    charge:    6.3276    magn:    1.4446    constr:    0.0000
     atom:   82    charge:    6.3292    magn:   -1.4308    constr:    0.0000
     atom:   83    charge:    6.3302    magn:    1.4459    constr:    0.0000
     atom:   84    charge:    6.3295    magn:   -1.4458    constr:    0.0000
     atom:   85    charge:    6.3256    magn:    1.4429    constr:    0.0000
     atom:   86    charge:    6.3255    magn:   -1.4635    constr:    0.0000
     atom:   87    charge:    6.3282    magn:    1.4584    constr:    0.0000
     atom:   88    charge:    6.3284    magn:   -1.4656    constr:    0.0000
     atom:   89    charge:    6.3258    magn:    1.4309    constr:    0.0000
     atom:   90    charge:    6.3254    magn:   -1.4653    constr:    0.0000
     atom:   91    charge:    6.3282    magn:    1.4605    constr:    0.0000
     atom:   92    charge:    6.3284    magn:   -1.4669    constr:    0.0000
     atom:   93    charge:    6.3242    magn:    1.4626    constr:    0.0000
     atom:   94    charge:    0.0556    magn:    0.0002    constr:    0.0000

     total cpu time spent up to now is      361.5 secs

     End of self-consistent calculation
--- enter write_ns ---
LDA+U parameters:
U( 1)     =  5.00000000
alpha( 1) =  0.00000000
U( 2)     =  5.00000000
alpha( 2) =  0.00000000
atom   71   Tr[ns(na)] (up, down, total) =   4.99245  0.92760  5.92004
   spin  1
    eigenvalues:
  0.997  0.998  0.999  0.999  0.999
    eigenvectors:
  0.238  0.002  0.017  0.082  0.661
  0.089  0.837  0.000  0.006  0.067
  0.009  0.046  0.004  0.804  0.137
  0.662  0.112  0.021  0.088  0.117
  0.002  0.003  0.958  0.020  0.017
.............

-0.003  0.010  0.078 -0.003 -0.004
  0.118 -0.002 -0.003  0.251 -0.001
  0.007 -0.003 -0.004 -0.001  0.090
atomic mag. moment =     4.083643
N of occupied +U levels =  136.091883
--- exit write_ns ---

------ SPIN UP ------------


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

   -81.4543 -81.0779 -81.0480 -81.0085 -80.8793 -80.8037 -80.7250 -80.7191
   -80.6397 -80.5329 -80.4853 -80.4845 -76.7942 -76.5235 -76.4942 -76.4087
   -76.2364 -76.2035 -76.1943 -76.1836 -76.0191 -75.9536 -75.6732 -49.9857
   -49.8737 -49.8127 -49.5679 -49.5463 -49.5101 -49.4876 -49.4754 -49.4589
   -49.4439 -49.4367 -49.3688 -49.3364 -49.3276 -49.2453 -49.2368 -49.2308
   -49.2143 -49.1982 -49.1676 -49.1616 -49.1599 -49.1502 -49.1454 -49.1428
   -49.1367 -49.0539 -48.9973 -48.9934 -48.9601 -48.9311 -48.9290 -48.8443
   -48.8406 -48.8368 -48.8234 -44.9497 -44.8467 -44.8347 -44.6129 -44.5823
   -44.5744 -44.5671 -44.5469 -44.5432 -44.4922 -44.4639 -44.4613 -44.3808
   -44.3211 -44.2958 -44.2865 -44.2779 -44.2684 -44.2618 -44.2592 -44.2515
   -44.2378 -44.2305 -44.1606 -44.1150 -44.0864 -44.0637 -44.0332 -44.0245
   -43.9853 -43.7521 -43.7410 -43.7178 -14.8811 -12.5668 -12.4905 -12.3138
   -12.2000 -12.1474 -12.0411 -12.0215 -11.9161 -11.8586 -11.6997 -11.5937
   -11.5304 -11.5233 -11.4384 -11.4328 -11.3117 -11.2979 -11.2165 -11.1071
   -10.8894 -10.7927 -10.7461 -10.6870 -10.6162 -10.5592 -10.4673 -10.3865
   -10.3544 -10.2252 -10.1497 -10.1461 -10.1000  -9.9609  -9.9015  -9.8786
    -9.8328  -9.8182  -9.7998  -9.7711  -9.6990  -9.6128  -9.5585  -9.5547
    -9.5112  -9.4961  -9.2076  -8.9177  -3.7449  -1.9273  -1.7489  -1.6449
    -1.5622  -1.4843  -1.4395  -1.3368  -1.3330  -1.2777  -1.1960  -1.1599
    -1.1322  -1.1100  -1.0709  -1.0303  -1.0167  -0.9895  -0.9558  -0.9350
    -0.9192  -0.8908  -0.8780  -0.8427  -0.8171  -0.8157  -0.7648  -0.7363
    -0.7230  -0.6988  -0.6798  -0.6759  -0.6387  -0.6095  -0.5990  -0.5939
    -0.5452  -0.4894  -0.4780  -0.4651  -0.4165  -0.3999  -0.3757  -0.3552
    -0.3278  -0.3141  -0.3036  -0.2776  -0.2373  -0.2153  -0.1992  -0.1701
    -0.1592  -0.1270  -0.0768  -0.0241   0.0291   0.0777   0.1224   0.2133
     0.2646   0.3232   0.3706   0.4179   0.5609   0.6071   0.7128   0.7379
     0.7856   0.8223   0.9119   0.9409   0.9826   1.0449   1.1265   1.1599
     1.1989   1.2694   1.3165   1.3461   1.4794   1.5609   1.6011   1.7642
     1.8470   1.8773   2.0559   2.0782   2.2092   2.2624   2.3101   2.3581
     2.4917   2.5094   2.5538   2.6130   2.6551   2.6851   2.8155   2.8543
     2.8811   2.9493   2.9836   3.0274   3.0714   3.1172   3.1369   3.1863
     3.2784   3.3376   3.3977   3.4807   3.5201   3.5658   3.6088   3.6489
     3.6757   3.7085   3.7236   3.7579   3.7949   3.8156   3.8381   3.8510
     3.9011   3.9766   4.0220   4.0818   4.1086   4.1211   4.1703   4.2044
     4.2428   4.2661   4.3258   4.3743   4.4167   4.4210   4.4560   4.4848
     4.5126   4.5269   4.5714   4.6305   4.6668   4.6930   4.7028   4.7261
     4.7849   4.7888   4.8215   4.8472   4.9219   4.9441   4.9522   4.9753
     4.9830   5.0512   5.0559   5.0746   5.1349   5.1548   5.1832   5.1926
     5.2440   5.2719   5.2799   5.2944   5.3383   5.3785   5.3951   5.4601
     5.4784   5.5499   5.5698   5.6336   5.7287   5.7427   5.7576   5.8488
     5.8696   5.9098   5.9815   6.0476   6.1373   6.1410   6.1722   6.2594
     6.3359   6.3895   6.4303   6.4599   6.5275   6.5968   6.6482   6.7306
     6.7627   6.8170   6.8512   7.0056   7.1295   7.2631   7.2824   7.3407
     8.9620   9.0499   9.0927   9.1760   9.2360   9.3012   9.3171   9.3380
     9.3717   9.3818   9.4591   9.4684   9.5455   9.5577   9.5792   9.6359
     9.6585   9.6710   9.6828   9.6928   9.6976   9.7108   9.7305   9.7591
     9.7924   9.8021   9.8229   9.8306   9.8538   9.8830   9.9465   9.9753
     9.9758   9.9944  10.0182  10.0709  10.0849  10.0932  10.1016  10.1200
    10.1542  10.1657  10.1731  10.2178  10.2335  10.3391  10.3543  10.4194
    10.4727  10.4769  10.4852  10.5218  10.5443  10.6358  10.6527  10.7756
    11.3671  12.6386  12.8873  13.5230  13.7098  13.9410  14.0815  14.2446
    14.3789  14.4333

------ SPIN DOWN ----------


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

   -81.2546 -81.0955 -81.0639 -80.9751 -80.8023 -80.7771 -80.7526 -80.6908
   -80.5258 -80.5254 -80.1990 -76.9605 -76.5789 -76.4845 -76.4791 -76.3818
   -76.2510 -76.1754 -76.1737 -76.1644 -76.0254 -75.9790 -75.9387 -49.8113
   -49.6712 -49.5866 -49.5378 -49.5151 -49.5068 -49.4904 -49.4893 -49.4777
   -49.4117 -49.4019 -49.3847 -49.2496 -49.2430 -49.2247 -49.2144 -49.2077
   -49.1968 -49.1957 -49.1918 -49.1831 -49.1739 -49.1713 -49.0224 -48.9888
   -48.9671 -48.9605 -48.9040 -48.8943 -48.8687 -48.6525 -48.6161 -48.5989
   -45.1187 -45.0142 -44.9975 -44.7024 -44.6318 -44.6247 -44.6052 -44.5687
   -44.5511 -44.5421 -44.5396 -44.5027 -44.4994 -44.4585 -44.4161 -44.3527
   -44.3370 -44.3137 -44.3016 -44.2657 -44.2525 -44.2469 -44.2345 -44.2294
   -44.2234 -44.2203 -44.1504 -44.1352 -44.1098 -44.0836 -44.0544 -44.0489
   -44.0280 -44.0206 -44.0177 -43.9735 -14.7863 -12.4814 -12.3632 -12.2968
   -12.2362 -12.1390 -11.9497 -11.9018 -11.7917 -11.6962 -11.6690 -11.6157
   -11.5695 -11.5064 -11.4289 -11.3667 -11.3255 -11.2607 -11.1319 -11.0119
   -10.8606 -10.8043 -10.7461 -10.7080 -10.5182 -10.4848 -10.4063 -10.3213
   -10.3037 -10.2565 -10.2481 -10.2078 -10.1885 -10.0998 -10.0116  -9.9762
    -9.9145  -9.8999  -9.7541  -9.7102  -9.6979  -9.6286  -9.5527  -9.5217
    -9.3980  -9.2516  -9.0621  -8.7388  -3.3111  -1.7582  -1.4925  -1.4594
    -1.3783  -1.2769  -1.2592  -1.2142  -1.2075  -1.1929  -1.1142  -1.1052
    -1.0549  -1.0225  -0.9979  -0.9681  -0.9447  -0.9243  -0.9142  -0.8901
    -0.8754  -0.8201  -0.7978  -0.7908  -0.7425  -0.7141  -0.6996  -0.6784
    -0.6255  -0.6089  -0.6065  -0.5646  -0.5227  -0.5139  -0.4732  -0.4320
    -0.4147  -0.4010  -0.3735  -0.3584  -0.3234  -0.2970  -0.2814  -0.2619
    -0.2213  -0.1663  -0.1486  -0.1125  -0.0409  -0.0214   0.0384   0.1047
     0.1606   0.1833   0.2145   0.2493   0.3029   0.3602   0.4197   0.4940
     0.5353   0.6083   0.6378   0.7139   0.8105   0.8490   0.8736   0.8813
     0.9375   0.9828   1.0642   1.0837   1.1108   1.1846   1.2702   1.4219
     1.4933   1.6210   1.6692   1.8182   1.8880   1.9713   2.0445   2.2055
     2.2334   2.2937   2.3835   2.4733   2.5266   2.5585   2.5642   2.6196
     2.6758   2.7500   2.8096   2.8977   2.9439   2.9846   3.0234   3.0543
     3.0950   3.1647   3.2059   3.2283   3.2907   3.3790   3.3902   3.4619
     3.5323   3.5625   3.5899   3.6141   3.6554   3.6747   3.7124   3.7583
     3.7915   3.8168   3.8210   3.9022   3.9144   3.9836   4.0052   4.0329
     4.0689   4.1075   4.1489   4.2071   4.2262   4.2727   4.2931   4.3406
     4.3876   4.3976   4.4684   4.4868   4.5326   4.5744   4.6115   4.6243
     4.6402   4.6911   4.7150   4.7742   4.8096   4.8391   4.8631   4.8687
     4.8814   4.9331   4.9386   4.9653   5.0309   5.0517   5.0724   5.0947
     5.1171   5.1479   5.1726   5.1894   5.2267   5.2723   5.3123   5.3230
     5.3538   5.4044   5.4096   5.4389   5.4875   5.5237   5.5770   5.6330
     5.6565   5.6976   5.7882   5.8142   5.8847   5.9113   6.0419   6.1105
     6.1821   6.2411   6.2799   6.3414   6.3594   6.4339   6.4606   6.5376
     6.6057   6.6840   6.7166   6.8100   6.8842   6.9375   7.0230   7.1696
     7.2477   7.2712   7.3129   8.7311   8.7889   8.8100   8.9740   9.0612
     9.0763   9.1908   9.2321   9.2894   9.3090   9.3723   9.4370   9.4567
     9.4728   9.5079   9.5200   9.5298   9.5381   9.5721   9.6007   9.6131
     9.6426   9.6692   9.6763   9.6917   9.7057   9.7385   9.7490   9.7531
     9.7866   9.7969   9.8082   9.8258   9.8528   9.8735   9.8827   9.9166
     9.9343   9.9429   9.9605   9.9913  10.0115  10.0346  10.0498  10.0604
    10.0965  10.1194  10.1268  10.1535  10.1814  10.1899  10.2318  10.2636
    10.2876  10.3396  10.3945  10.4323  10.4980  10.5861  10.6295  10.6800
    11.4856  12.7291  12.9304  13.6379  13.8050  14.0260  14.1594  14.3288
    14.4686  14.6194

     the Fermi energy is     7.5012 ev

!    total energy              =   -7085.16398555 Ry
     Harris-Foulkes estimate   =   -7085.16398586 Ry
     estimated scf accuracy    <          9.3E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =   -5053.67689612 Ry
     hartree contribution      =    2784.75257815 Ry
     xc contribution           =    -949.37013145 Ry
     ewald contribution        =   -3868.73345437 Ry
     Dispersion Correction     =      -0.99120707 Ry
     Hubbard energy            =       2.84880267 Ry
     smearing contrib. (-TS)   =       0.00632264 Ry

     total magnetization       =     4.93 Bohr mag/cell
     absolute magnetization    =   106.78 Bohr mag/cell

     convergence has been achieved in  29 iterations












MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .
MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .

MKL ERROR: Parameter 10 was incorrect on entry to DGEMM .



Best,

Shaofeng



--------------------------------------
Shaofeng Wang, Ph.D of Geochemistry
Environmental Molecular Science Group
Institute of Applied Ecology, Chinese Academy of Sciences
Shenyang, 110016, China
wangshaofeng at iae.ac.cn
www.iae.cas.cn


_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://pwscf.org/mailman/listinfo/pw_forum



-- 
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
Phone +39-0432-558216, fax +39-0432-558222

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20161114/326dd454/attachment.html>


More information about the users mailing list