
     Program PWSCF v.6.0 (svn rev. 13188M) starts on  9Dec2016 at 16: 4:29 

     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     2 processors
     R & G space division:  proc/nbgrp/npool/nimage =       2
     Waiting for input...
     Reading input from standard input
     Message from routine read_cards :
     DEPRECATED: no units specified in CELL_PARAMETERS card

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3

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


     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min         174     174     60                 2079     2079     416
     Max         175     175     61                 2080     2080     417
     Sum         349     349    121                 4159     4159     833



     bravais-lattice index     =            0
     lattice parameter (alat)  =       7.0103  a.u.
     unit-cell volume          =     245.3705 (a.u.)^3
     number of atoms/cell      =            2
     number of atomic types    =            1
     number of electrons       =        10.00
     number of Kohn-Sham states=            9
     kinetic-energy cutoff     =      25.0000  Ry
     charge density cutoff     =     100.0000  Ry
     convergence threshold     =      1.0E-07
     mixing beta               =       0.7000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      =  SLA  PZ   NOGX NOGC ( 1  1  0  0 0 0)
     nstep                     =           50


     celldm(1)=   7.010336  celldm(2)=   0.000000  celldm(3)=   0.000000
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   0.580130   0.000000   0.814524 )  
               a(2) = (  -0.290065   0.502407   0.814524 )  
               a(3) = (  -0.290065  -0.502407   0.814524 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  1.149169  0.000000  0.409237 )  
               b(2) = ( -0.574584  0.995209  0.409237 )  
               b(3) = ( -0.574584 -0.995209  0.409237 )  


     PseudoPot. # 1 for As read from file:
     /home/pietro/espresso-svn/pseudo/As.pz-bhs.UPF
     MD5 check sum: 451cd3365afcfc94d28b1934951c34a8
     Pseudo is Norm-conserving, Zval =  5.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  525 points,  2 beta functions with: 
                l(1) =   0
                l(2) =   1

     atomic species   valence    mass     pseudopotential
        As             5.00    74.90000     As( 1.00)

     12 Sym. Ops., with inversion, found



   Cartesian axes

     site n.     atom                  positions (alat units)
         1           As  tau(   1) = (   0.0000001   0.0000000   0.7086605  )
         2           As  tau(   2) = (  -0.0000001   0.0000000  -0.7086605  )

     number of k points=    10  Marzari-Vanderbilt smearing, width (Ry)=  0.0050
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.1534638), wk =   0.0625000
        k(    2) = (  -0.1436461  -0.2488023   0.2557731), wk =   0.1875000
        k(    3) = (   0.2872922   0.4976046  -0.0511547), wk =   0.1875000
        k(    4) = (   0.1436461   0.2488023   0.0511546), wk =   0.1875000
        k(    5) = (  -0.2872922   0.0000000   0.3580823), wk =   0.1875000
        k(    6) = (   0.1436461   0.7464070   0.0511546), wk =   0.3750000
        k(    7) = (   0.0000000   0.4976046   0.1534638), wk =   0.3750000
        k(    8) = (   0.5745844   0.0000000  -0.2557731), wk =   0.1875000
        k(    9) = (   0.0000000   0.0000000   0.4603915), wk =   0.0625000
        k(   10) = (   0.4309383   0.7464070   0.1534638), wk =   0.1875000

     Dense  grid:     4159 G-vectors     FFT dimensions: (  24,  24,  24)

     Estimated max dynamical RAM per process >       1.91Mb

     Estimated total allocated dynamical RAM >       3.83Mb

     Initial potential from superposition of free atoms

     starting charge    9.99960, renormalised to   10.00000
     Starting wfc are    8 randomized atomic wfcs +    1 random wfc

     total cpu time spent up to now is        0.1 secs

     per-process dynamical memory:     5.0 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  4.3

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

     Davidson diagonalization with overlap
     ethr =  1.54E-04,  avg # of iterations =  1.3

     total cpu time spent up to now is        0.2 secs

     total energy              =     -25.43995582 Ry
     Harris-Foulkes estimate   =     -25.44371120 Ry
     estimated scf accuracy    <       0.01555792 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.56E-04,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.3 secs

     total energy              =     -25.44012469 Ry
     Harris-Foulkes estimate   =     -25.44030751 Ry
     estimated scf accuracy    <       0.00088879 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.89E-06,  avg # of iterations =  1.6

     total cpu time spent up to now is        0.3 secs

     total energy              =     -25.44015903 Ry
     Harris-Foulkes estimate   =     -25.44016035 Ry
     estimated scf accuracy    <       0.00000501 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.01E-08,  avg # of iterations =  3.1

     total cpu time spent up to now is        0.4 secs

     total energy              =     -25.44016731 Ry
     Harris-Foulkes estimate   =     -25.44016767 Ry
     estimated scf accuracy    <       0.00000077 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  7.74E-09,  avg # of iterations =  1.4

     total cpu time spent up to now is        0.4 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1535 (   531 PWs)   bands (ev):

    -6.9960   4.5198   5.9668   5.9668   8.4358  11.0403  11.7602  11.7603
    16.5646

          k =-0.1436-0.2488 0.2558 (   522 PWs)   bands (ev):

    -5.9250   0.3918   5.3512   5.6502   9.2996  10.5304  11.7007  13.5632
    15.7167

          k = 0.2873 0.4976-0.0512 (   520 PWs)   bands (ev):

    -4.3490  -2.4704   4.7884   6.1554   7.8796  10.8149  12.5850  13.8262
    17.7262

          k = 0.1436 0.2488 0.0512 (   525 PWs)   bands (ev):

    -6.3695   1.3043   4.9860   7.1721   8.5435  10.8048  12.4703  13.9613
    15.3510

          k =-0.2873 0.0000 0.3581 (   519 PWs)   bands (ev):

    -5.5427   1.1265   3.5658   4.2978   7.5160  10.4216  13.7076  13.7747
    16.9047

          k = 0.1436 0.7464 0.0512 (   510 PWs)   bands (ev):

    -3.8393  -1.8099   2.3271   4.2466   8.0539  11.6204  13.3234  15.7203
    17.3490

          k = 0.0000 0.4976 0.1535 (   521 PWs)   bands (ev):

    -4.7124  -1.4722   3.0016   6.6926   7.7778  12.3034  13.0675  13.4305
    16.0962

          k = 0.5746 0.0000-0.2558 (   510 PWs)   bands (ev):

    -4.0542  -1.5061   3.7085   3.7297   6.0243  10.0592  15.9113  17.7151
    18.4776

          k = 0.0000 0.0000 0.4604 (   522 PWs)   bands (ev):

    -5.8586   0.8362   5.8840   5.8840   7.4111  10.0628  10.0628  12.1193
    17.3945

          k = 0.4309 0.7464 0.1535 (   520 PWs)   bands (ev):

    -4.8493  -0.0497   2.4337   4.7831   7.5089  11.6828  12.0643  14.4760
    17.7702

     the Fermi energy is     9.7636 ev

!    total energy              =     -25.44016740 Ry
     Harris-Foulkes estimate   =     -25.44016741 Ry
     estimated scf accuracy    <       0.00000002 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       7.72807894 Ry
     hartree contribution      =       1.22170646 Ry
     xc contribution           =      -6.50442316 Ry
     ewald contribution        =     -27.88552965 Ry
     smearing contrib. (-TS)   =      -0.00000000 Ry

     convergence has been achieved in   5 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000   -0.12672086
     atom    2 type  1   force =    -0.00000000    0.00000000    0.12672086

     Total force =     0.179210     Total SCF correction =     0.000017


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=  217.55
   0.00172500  -0.00000000  -0.00000000        253.76     -0.00     -0.00
  -0.00000000   0.00172500  -0.00000000         -0.00    253.76     -0.00
  -0.00000000   0.00000000   0.00098671         -0.00      0.00    145.15


     BFGS Geometry Optimization

     number of scf cycles    =   1
     number of bfgs steps    =   0

     enthalpy new            =     -25.4401674000 Ry

     new trust radius        =       0.1351044206 bohr
     new conv_thr            =       0.0000001000 Ry

     new unit-cell volume =    273.57599 a.u.^3 (    40.53978 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.605147602  -0.000000000   0.834616652
  -0.302573611   0.524073148   0.834616658
  -0.302573611  -0.524073148   0.834616658

ATOMIC_POSITIONS (crystal)
As       0.282612520   0.282612515   0.282612515
As      -0.282612520  -0.282612515  -0.282612515



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   11.03095, renormalised to   10.00000

     total cpu time spent up to now is        0.6 secs

     per-process dynamical memory:     5.8 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  5.9

     total cpu time spent up to now is        0.7 secs

     total energy              =     -25.47728472 Ry
     Harris-Foulkes estimate   =     -26.08217532 Ry
     estimated scf accuracy    <       0.00508025 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.08E-05,  avg # of iterations =  3.1

     total cpu time spent up to now is        0.8 secs

     total energy              =     -25.48676619 Ry
     Harris-Foulkes estimate   =     -25.48848155 Ry
     estimated scf accuracy    <       0.00433134 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.33E-05,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.8 secs

     total energy              =     -25.48660142 Ry
     Harris-Foulkes estimate   =     -25.48695092 Ry
     estimated scf accuracy    <       0.00097272 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  9.73E-06,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.9 secs

     total energy              =     -25.48650510 Ry
     Harris-Foulkes estimate   =     -25.48664514 Ry
     estimated scf accuracy    <       0.00026347 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.63E-06,  avg # of iterations =  2.8

     total cpu time spent up to now is        1.0 secs

     total energy              =     -25.48655857 Ry
     Harris-Foulkes estimate   =     -25.48655931 Ry
     estimated scf accuracy    <       0.00000251 Ry

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.51E-08,  avg # of iterations =  2.2

     total cpu time spent up to now is        1.0 secs

     total energy              =     -25.48655903 Ry
     Harris-Foulkes estimate   =     -25.48655946 Ry
     estimated scf accuracy    <       0.00000108 Ry

     iteration #  7     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.08E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is        1.1 secs

     total energy              =     -25.48655912 Ry
     Harris-Foulkes estimate   =     -25.48655917 Ry
     estimated scf accuracy    <       0.00000011 Ry

     iteration #  8     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.13E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        1.1 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1498 (   531 PWs)   bands (ev):

    -7.2730   2.2868   5.0421   5.0421   6.9072   9.5375  10.3837  10.3837
    14.6597

          k =-0.1377-0.2385 0.2496 (   522 PWs)   bands (ev):

    -6.2823  -0.8521   4.3357   4.9766   7.6864   8.5615   9.6220  11.9330
    13.9339

          k = 0.2754 0.4770-0.0499 (   520 PWs)   bands (ev):

    -4.8738  -3.2673   4.1656   5.0068   6.3102   9.2341  10.4184  11.6360
    16.1380

          k = 0.1377 0.2385 0.0499 (   525 PWs)   bands (ev):

    -6.7201   0.0194   4.2085   5.6907   7.0523   9.7000  10.1591  11.9480
    13.8704

          k =-0.2754 0.0000 0.3495 (   519 PWs)   bands (ev):

    -5.9136  -0.3515   2.7965   3.5200   5.8875   9.3588  11.6622  11.9512
    14.3956

          k = 0.1377 0.7155 0.0499 (   510 PWs)   bands (ev):

    -4.4002  -2.7153   1.7325   3.1867   6.4968   9.7746  12.0134  13.7560
    14.8120

          k = 0.0000 0.4770 0.1498 (   521 PWs)   bands (ev):

    -5.2248  -2.3973   2.4901   5.2347   6.4410  10.2387  11.1192  12.0033
    14.0429

          k = 0.5508 0.0000-0.2496 (   510 PWs)   bands (ev):

    -4.6715  -2.2059   2.3724   3.0121   4.6417   9.2661  13.5073  15.1141
    15.5041

          k = 0.0000 0.0000 0.4493 (   522 PWs)   bands (ev):

    -6.1278  -0.6897   5.0548   5.0548   5.9417   8.4854   8.4854   9.8230
    15.7278

          k = 0.4131 0.7155 0.1498 (   520 PWs)   bands (ev):

    -5.1967  -1.5133   1.8547   4.0502   5.9453   9.8231  10.1482  12.7569
    15.7093

     the Fermi energy is     8.1330 ev

!    total energy              =     -25.48655913 Ry
     Harris-Foulkes estimate   =     -25.48655913 Ry
     estimated scf accuracy    <          6.8E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.69378678 Ry
     hartree contribution      =       1.27320853 Ry
     xc contribution           =      -6.33961648 Ry
     ewald contribution        =     -27.11393796 Ry
     smearing contrib. (-TS)   =       0.00000000 Ry

     convergence has been achieved in   8 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000   -0.00000000   -0.05672083
     atom    2 type  1   force =     0.00000000    0.00000000    0.05672083

     Total force =     0.080215     Total SCF correction =     0.000003


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=   35.25
   0.00014812   0.00000000  -0.00000000         21.79      0.00     -0.00
  -0.00000000   0.00014812  -0.00000000         -0.00     21.79     -0.00
  -0.00000000   0.00000000   0.00042269         -0.00      0.00     62.18


     number of scf cycles    =   2
     number of bfgs steps    =   1

     enthalpy old            =     -25.4401674000 Ry
     enthalpy new            =     -25.4865591321 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0911513437 bohr
     new conv_thr            =       0.0000001000 Ry

     new unit-cell volume =    282.71828 a.u.^3 (    41.89453 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.609967692  -0.000000000   0.848930076
  -0.304983656   0.528247467   0.848930083
  -0.304983656  -0.528247467   0.848930083

ATOMIC_POSITIONS (crystal)
As       0.277419551   0.277419542   0.277419542
As      -0.277419551  -0.277419542  -0.277419542



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.32336, renormalised to   10.00000

     total cpu time spent up to now is        1.3 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  5.0

     total cpu time spent up to now is        1.4 secs

     total energy              =     -25.49417391 Ry
     Harris-Foulkes estimate   =     -25.67858569 Ry
     estimated scf accuracy    <       0.00064463 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.45E-06,  avg # of iterations =  3.0

     total cpu time spent up to now is        1.5 secs

     total energy              =     -25.49502312 Ry
     Harris-Foulkes estimate   =     -25.49516132 Ry
     estimated scf accuracy    <       0.00032922 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  3.29E-06,  avg # of iterations =  1.0

     total cpu time spent up to now is        1.5 secs

     total energy              =     -25.49501488 Ry
     Harris-Foulkes estimate   =     -25.49504039 Ry
     estimated scf accuracy    <       0.00005574 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.57E-07,  avg # of iterations =  1.2

     total cpu time spent up to now is        1.6 secs

     total energy              =     -25.49501747 Ry
     Harris-Foulkes estimate   =     -25.49501978 Ry
     estimated scf accuracy    <       0.00000426 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.26E-08,  avg # of iterations =  3.0

     total cpu time spent up to now is        1.6 secs

     total energy              =     -25.49501993 Ry
     Harris-Foulkes estimate   =     -25.49501994 Ry
     estimated scf accuracy    <       0.00000022 Ry

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.19E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        1.7 secs

     total energy              =     -25.49501985 Ry
     Harris-Foulkes estimate   =     -25.49501993 Ry
     estimated scf accuracy    <       0.00000017 Ry

     iteration #  7     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.68E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        1.7 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1472 (   531 PWs)   bands (ev):

    -7.3116   1.5356   4.9273   4.9273   6.3974   9.2622  10.0287  10.0287
    14.0911

          k =-0.1366-0.2366 0.2454 (   522 PWs)   bands (ev):

    -6.3360  -1.2276   3.9743   4.9826   7.3558   7.9426   8.8718  11.4928
    13.4184

          k = 0.2732 0.4733-0.0491 (   520 PWs)   bands (ev):

    -4.9466  -3.4944   4.1521   4.5889   5.8650   8.8110   9.5957  10.6957
    15.7190

          k = 0.1366 0.2366 0.0491 (   525 PWs)   bands (ev):

    -6.7753  -0.3334   4.1253   5.2259   6.5546   9.2907   9.6089  11.2565
    13.3180

          k =-0.2732 0.0000 0.3436 (   519 PWs)   bands (ev):

    -5.9690  -0.8520   2.6380   3.4344   5.2739   9.3156  11.1848  11.4551
    13.5091

          k = 0.1366 0.7099 0.0491 (   510 PWs)   bands (ev):

    -4.4954  -2.9514   1.5961   2.8176   5.9893   9.3173  11.7785  13.1755
    13.9178

          k = 0.0000 0.4733 0.1472 (   521 PWs)   bands (ev):

    -5.3170  -2.6220   2.4376   4.7080   6.0029   9.4631  10.6173  11.5928
    13.4045

          k = 0.5465 0.0000-0.2454 (   510 PWs)   bands (ev):

    -4.7987  -2.3493   1.8772   2.9436   4.1371   9.2137  12.7084  14.1861
    14.5614

          k = 0.0000 0.0000 0.4417 (   522 PWs)   bands (ev):

    -6.1409  -1.3815   5.0163   5.0163   5.8703   8.0227   8.0227   9.1000
    15.2896

          k = 0.4099 0.7099 0.1472 (   520 PWs)   bands (ev):

    -5.2096  -2.0879   1.7602   3.9972   5.5060   9.4665   9.5765  12.3601
    15.0545

     the Fermi energy is     7.4184 ev

!    total energy              =     -25.49501986 Ry
     Harris-Foulkes estimate   =     -25.49501986 Ry
     estimated scf accuracy    <          4.5E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.44090272 Ry
     hartree contribution      =       1.27007810 Ry
     xc contribution           =      -6.28731981 Ry
     ewald contribution        =     -26.91875732 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   7 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000    0.00000000   -0.02258260
     atom    2 type  1   force =     0.00000000   -0.00000000    0.02258260

     Total force =     0.031937     Total SCF correction =     0.000036


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=  -12.19
  -0.00021023   0.00000000   0.00000000        -30.93      0.00      0.00
   0.00000000  -0.00021023   0.00000000          0.00    -30.93      0.00
   0.00000000  -0.00000000   0.00017189          0.00     -0.00     25.29


     number of scf cycles    =   3
     number of bfgs steps    =   2

     enthalpy old            =     -25.4865591321 Ry
     enthalpy new            =     -25.4950198607 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0616076975 bohr
     new conv_thr            =       0.0000001000 Ry

     new unit-cell volume =    283.63570 a.u.^3 (    42.03047 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.607526414  -0.000000000   0.858543415
  -0.303763017   0.526133258   0.858543421
  -0.303763017  -0.526133258   0.858543421

ATOMIC_POSITIONS (crystal)
As       0.273968886   0.273968875   0.273968875
As      -0.273968886  -0.273968875  -0.273968875



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.03234, renormalised to   10.00000

     total cpu time spent up to now is        1.9 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  4.7

     total cpu time spent up to now is        2.0 secs

     total energy              =     -25.49749568 Ry
     Harris-Foulkes estimate   =     -25.51601799 Ry
     estimated scf accuracy    <       0.00016539 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.65E-06,  avg # of iterations =  2.0

     total cpu time spent up to now is        2.0 secs

     total energy              =     -25.49752956 Ry
     Harris-Foulkes estimate   =     -25.49753059 Ry
     estimated scf accuracy    <       0.00000402 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.02E-08,  avg # of iterations =  2.4

     total cpu time spent up to now is        2.1 secs

     total energy              =     -25.49753051 Ry
     Harris-Foulkes estimate   =     -25.49753044 Ry
     estimated scf accuracy    <       0.00000022 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.19E-09,  avg # of iterations =  1.2

     total cpu time spent up to now is        2.1 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1456 (   531 PWs)   bands (ev):

    -7.2663   1.3584   5.0705   5.0705   6.2472   9.4010  10.0441  10.0441
    14.0072

          k =-0.1372-0.2376 0.2427 (   522 PWs)   bands (ev):

    -6.2825  -1.2639   3.8765   5.1916   7.4313   7.8271   8.6372  11.4348
    13.4015

          k = 0.2743 0.4752-0.0485 (   520 PWs)   bands (ev):

    -4.8651  -3.4957   4.3224   4.4619   5.8116   8.8119   9.3187  10.3214
    15.5531

          k = 0.1372 0.2376 0.0485 (   525 PWs)   bands (ev):

    -6.7267  -0.3216   4.2625   5.1090   6.4466   9.0724   9.7258  11.0756
    13.1442

          k =-0.2743 0.0000 0.3397 (   519 PWs)   bands (ev):

    -5.9115  -0.9639   2.6938   3.5661   5.0630   9.5687  11.2167  11.4382
    13.2486

          k = 0.1372 0.7127 0.0485 (   510 PWs)   bands (ev):

    -4.4229  -2.9361   1.6195   2.7149   5.8821   9.3090  11.8980  13.1446
    13.6482

          k = 0.0000 0.4752 0.1456 (   521 PWs)   bands (ev):

    -5.2603  -2.5852   2.5322   4.5448   5.9112   9.1761  10.5968  11.6147
    13.2885

          k = 0.5487-0.0000-0.2427 (   510 PWs)   bands (ev):

    -4.7486  -2.2742   1.7075   3.0717   3.9675   9.3839  12.5133  13.8649
    14.3773

          k = 0.0000-0.0000 0.4368 (   522 PWs)   bands (ev):

    -6.0660  -1.6657   5.2146   5.2146   6.1495   7.9725   7.9725   8.9718
    15.2699

          k = 0.4115 0.7127 0.1456 (   520 PWs)   bands (ev):

    -5.1064  -2.2952   1.8312   4.1516   5.4583   9.4963   9.5958  12.4430
    14.8915

     the Fermi energy is     7.4940 ev

!    total energy              =     -25.49753053 Ry
     Harris-Foulkes estimate   =     -25.49753053 Ry
     estimated scf accuracy    <          4.1E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.47830256 Ry
     hartree contribution      =       1.24382992 Ry
     xc contribution           =      -6.27599985 Ry
     ewald contribution        =     -26.94373961 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   4 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000    0.00000000   -0.00372975
     atom    2 type  1   force =     0.00000000   -0.00000000    0.00372975

     Total force =     0.005275     Total SCF correction =     0.000048


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=  -24.52
  -0.00026668   0.00000000   0.00000000        -39.23      0.00      0.00
   0.00000000  -0.00026668  -0.00000000          0.00    -39.23     -0.00
   0.00000000  -0.00000000   0.00003322          0.00     -0.00      4.89


     number of scf cycles    =   4
     number of bfgs steps    =   3

     enthalpy old            =     -25.4950198607 Ry
     enthalpy new            =     -25.4975305264 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0333195651 bohr
     new conv_thr            =       0.0000000373 Ry

     new unit-cell volume =    280.10915 a.u.^3 (    41.50789 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.601896405  -0.000000000   0.863804576
  -0.300948012   0.521257528   0.863804580
  -0.300948012  -0.521257528   0.863804580

ATOMIC_POSITIONS (crystal)
As       0.272123544   0.272123533   0.272123533
As      -0.272123544  -0.272123533  -0.272123533



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge    9.87411, renormalised to   10.00000

     total cpu time spent up to now is        2.3 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  4.2

     total cpu time spent up to now is        2.4 secs

     total energy              =     -25.49857806 Ry
     Harris-Foulkes estimate   =     -25.42608085 Ry
     estimated scf accuracy    <       0.00015869 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.59E-06,  avg # of iterations =  2.7

     total cpu time spent up to now is        2.5 secs

     total energy              =     -25.49873831 Ry
     Harris-Foulkes estimate   =     -25.49877781 Ry
     estimated scf accuracy    <       0.00008845 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.85E-07,  avg # of iterations =  1.0

     total cpu time spent up to now is        2.5 secs

     total energy              =     -25.49873964 Ry
     Harris-Foulkes estimate   =     -25.49874522 Ry
     estimated scf accuracy    <       0.00001113 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.11E-07,  avg # of iterations =  1.7

     total cpu time spent up to now is        2.6 secs

     total energy              =     -25.49874087 Ry
     Harris-Foulkes estimate   =     -25.49874102 Ry
     estimated scf accuracy    <       0.00000029 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.87E-09,  avg # of iterations =  2.9

     total cpu time spent up to now is        2.6 secs

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

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  7.39E-10,  avg # of iterations =  1.0

     total cpu time spent up to now is        2.7 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1447 (   531 PWs)   bands (ev):

    -7.1840   1.4886   5.3261   5.3261   6.3254   9.7045  10.2567  10.2567
    14.1785

          k =-0.1385-0.2398 0.2412 (   522 PWs)   bands (ev):

    -6.1817  -1.1114   3.9219   5.4641   7.6940   7.9906   8.7147  11.6018
    13.6213

          k = 0.2769 0.4796-0.0482 (   520 PWs)   bands (ev):

    -4.7161  -3.3736   4.5010   4.5570   5.9660   9.0280   9.3692  10.2900
    15.5625

          k = 0.1385 0.2398 0.0482 (   525 PWs)   bands (ev):

    -6.6306  -0.1146   4.4926   5.1762   6.5474   9.1480   9.9938  11.1896
    13.2169

          k =-0.2769 0.0000 0.3377 (   519 PWs)   bands (ev):

    -5.8051  -0.8528   2.8396   3.7904   5.1031   9.9074  11.5116  11.6753
    13.3606

          k = 0.1385 0.7194 0.0482 (   510 PWs)   bands (ev):

    -4.2751  -2.7843   1.7265   2.7583   5.9914   9.5352  12.1777  13.3908
    13.7359

          k = 0.0000 0.4796 0.1447 (   521 PWs)   bands (ev):

    -5.1333  -2.4141   2.6861   4.5961   6.0004   9.1941  10.8213  11.8682
    13.4375

          k = 0.5538-0.0000-0.2412 (   510 PWs)   bands (ev):

    -4.6122  -2.0901   1.7234   3.2820   3.9973   9.6225  12.6525  13.9399
    14.5695

          k = 0.0000-0.0000 0.4341 (   522 PWs)   bands (ev):

    -5.9619  -1.6922   5.5010   5.5010   6.5553   8.1506   8.1506   9.1765
    15.4603

          k = 0.4154 0.7194 0.1447 (   520 PWs)   bands (ev):

    -4.9670  -2.2709   1.9745   4.3879   5.6249   9.6852   9.9456  12.7674
    14.9915

     the Fermi energy is     7.7567 ev

!    total energy              =     -25.49874122 Ry
     Harris-Foulkes estimate   =     -25.49874124 Ry
     estimated scf accuracy    <       0.00000003 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.65775231 Ry
     hartree contribution      =       1.21059515 Ry
     xc contribution           =      -6.28758593 Ry
     ewald contribution        =     -27.07957920 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   6 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =     0.00000000   -0.00000000    0.00387824
     atom    2 type  1   force =    -0.00000000    0.00000000   -0.00387824

     Total force =     0.005485     Total SCF correction =     0.000141


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=  -19.26
  -0.00018924   0.00000000   0.00000000        -27.84      0.00      0.00
   0.00000000  -0.00018924   0.00000000          0.00    -27.84      0.00
   0.00000000   0.00000000  -0.00001438          0.00      0.00     -2.12


     number of scf cycles    =   5
     number of bfgs steps    =   4

     enthalpy old            =     -25.4975305264 Ry
     enthalpy new            =     -25.4987412187 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0367520472 bohr
     new conv_thr            =       0.0000000388 Ry

     new unit-cell volume =    274.78818 a.u.^3 (    40.71941 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.595287803  -0.000000000   0.866314900
  -0.297643712   0.515534312   0.866314903
  -0.297643712  -0.515534312   0.866314903

ATOMIC_POSITIONS (crystal)
As       0.271447979   0.271447967   0.271447967
As      -0.271447979  -0.271447967  -0.271447967



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge    9.80637, renormalised to   10.00000

     total cpu time spent up to now is        2.9 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  4.0

     total cpu time spent up to now is        3.0 secs

     total energy              =     -25.49906744 Ry
     Harris-Foulkes estimate   =     -25.38589553 Ry
     estimated scf accuracy    <       0.00020362 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.04E-06,  avg # of iterations =  3.0

     total cpu time spent up to now is        3.0 secs

     total energy              =     -25.49936865 Ry
     Harris-Foulkes estimate   =     -25.49943642 Ry
     estimated scf accuracy    <       0.00016130 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.61E-06,  avg # of iterations =  1.0

     total cpu time spent up to now is        3.1 secs

     total energy              =     -25.49937210 Ry
     Harris-Foulkes estimate   =     -25.49938006 Ry
     estimated scf accuracy    <       0.00002492 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.49E-07,  avg # of iterations =  1.0

     total cpu time spent up to now is        3.1 secs

     total energy              =     -25.49936867 Ry
     Harris-Foulkes estimate   =     -25.49937312 Ry
     estimated scf accuracy    <       0.00000807 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.07E-08,  avg # of iterations =  2.4

     total cpu time spent up to now is        3.2 secs

     total energy              =     -25.49937053 Ry
     Harris-Foulkes estimate   =     -25.49937053 Ry
     estimated scf accuracy    <       0.00000006 Ry

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.75E-10,  avg # of iterations =  1.9

     total cpu time spent up to now is        3.2 secs

     total energy              =     -25.49937052 Ry
     Harris-Foulkes estimate   =     -25.49937054 Ry
     estimated scf accuracy    <       0.00000004 Ry

     iteration #  7     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.10E-10,  avg # of iterations =  2.0

     total cpu time spent up to now is        3.3 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1443 (   531 PWs)   bands (ev):

    -7.0925   1.7596   5.6164   5.6164   6.5317  10.0417  10.5547  10.5547
    14.4677

          k =-0.1400-0.2425 0.2405 (   522 PWs)   bands (ev):

    -6.0685  -0.8707   4.0313   5.7325   8.0281   8.2882   8.9636  11.8859
    13.9280

          k = 0.2800 0.4849-0.0481 (   520 PWs)   bands (ev):

    -4.5489  -3.1973   4.6198   4.7965   6.2181   9.3317   9.5942  10.4366
    15.6894

          k = 0.1400 0.2425 0.0481 (   525 PWs)   bands (ev):

    -6.5205   0.1765   4.7482   5.3331   6.7377   9.3855  10.2926  11.4517
    13.4318

          k =-0.2800 0.0000 0.3367 (   519 PWs)   bands (ev):

    -5.6876  -0.6327   3.0075   4.0425   5.2838  10.2367  11.9271  12.0268
    13.6726

          k = 0.1400 0.7274 0.0481 (   510 PWs)   bands (ev):

    -4.1073  -2.5713   1.8695   2.8715   6.2004   9.8731  12.5103  13.7411
    14.0069

          k = 0.0000 0.4849 0.1443 (   521 PWs)   bands (ev):

    -4.9829  -2.1887   2.8517   4.7604   6.1650   9.3800  11.1586  12.2146
    13.7029

          k = 0.5600 0.0000-0.2405 (   510 PWs)   bands (ev):

    -4.4455  -1.8705   1.8381   3.5154   4.1326   9.8579  12.9534  14.2355
    14.9297

          k = 0.0000 0.0000 0.4329 (   522 PWs)   bands (ev):

    -5.8607  -1.5793   5.8026   5.8026   6.9926   8.4394   8.4394   9.5425
    15.7385

          k = 0.4200 0.7274 0.1443 (   520 PWs)   bands (ev):

    -4.8349  -2.1169   2.1391   4.6429   5.8961  10.0037  10.3744  13.1867
    15.2207

     the Fermi energy is     8.0907 ev

!    total energy              =     -25.49937052 Ry
     Harris-Foulkes estimate   =     -25.49937053 Ry
     estimated scf accuracy    <          9.9E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.88665965 Ry
     hartree contribution      =       1.18019809 Ry
     xc contribution           =      -6.31054790 Ry
     ewald contribution        =     -27.25575681 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   7 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000    0.00403082
     atom    2 type  1   force =    -0.00000000   -0.00000000   -0.00403082

     Total force =     0.005700     Total SCF correction =     0.000084


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=   -4.58
  -0.00004771   0.00000000   0.00000000         -7.02      0.00      0.00
   0.00000000  -0.00004771   0.00000000          0.00     -7.02      0.00
   0.00000000   0.00000000   0.00000204          0.00      0.00      0.30


     number of scf cycles    =   6
     number of bfgs steps    =   5

     enthalpy old            =     -25.4987412187 Ry
     enthalpy new            =     -25.4993705209 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0082525329 bohr
     new conv_thr            =       0.0000000403 Ry

     new unit-cell volume =    273.44166 a.u.^3 (    40.51987 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593806018   0.000000000   0.866377544
  -0.296902819   0.514251049   0.866377547
  -0.296902819  -0.514251049   0.866377547

ATOMIC_POSITIONS (crystal)
As       0.271764507   0.271764495   0.271764495
As      -0.271764507  -0.271764495  -0.271764495



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge    9.95076, renormalised to   10.00000

     total cpu time spent up to now is        3.5 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  3.0

     total cpu time spent up to now is        3.5 secs

     total energy              =     -25.49941606 Ry
     Harris-Foulkes estimate   =     -25.47047576 Ry
     estimated scf accuracy    <       0.00001415 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.41E-07,  avg # of iterations =  3.0

     total cpu time spent up to now is        3.6 secs

     total energy              =     -25.49943468 Ry
     Harris-Foulkes estimate   =     -25.49943827 Ry
     estimated scf accuracy    <       0.00000822 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.22E-08,  avg # of iterations =  1.0

     total cpu time spent up to now is        3.7 secs

     total energy              =     -25.49943486 Ry
     Harris-Foulkes estimate   =     -25.49943529 Ry
     estimated scf accuracy    <       0.00000114 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.14E-08,  avg # of iterations =  1.0

     total cpu time spent up to now is        3.7 secs

     total energy              =     -25.49943480 Ry
     Harris-Foulkes estimate   =     -25.49943493 Ry
     estimated scf accuracy    <       0.00000025 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.51E-09,  avg # of iterations =  2.9

     total cpu time spent up to now is        3.8 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1443 (   531 PWs)   bands (ev):

    -7.0793   1.8448   5.6696   5.6696   6.6007  10.0929  10.6241  10.6242
    14.5523

          k =-0.1403-0.2431 0.2405 (   522 PWs)   bands (ev):

    -6.0516  -0.8063   4.0610   5.7642   8.1050   8.3722   9.0617  11.9695
    13.9926

          k = 0.2807 0.4861-0.0481 (   520 PWs)   bands (ev):

    -4.5215  -3.1551   4.6567   4.8287   6.2884   9.4010   9.6859  10.5086
    15.7442

          k = 0.1403 0.2431 0.0481 (   525 PWs)   bands (ev):

    -6.5023   0.2452   4.7926   5.3837   6.7865   9.4778  10.3376  11.5369
    13.5117

          k =-0.2807-0.0000 0.3367 (   519 PWs)   bands (ev):

    -5.6716  -0.5643   3.0320   4.0892   5.3581  10.2820  12.0300  12.1162
    13.7874

          k = 0.1403 0.7292 0.0481 (   510 PWs)   bands (ev):

    -4.0816  -2.5207   1.9016   2.9074   6.2600   9.9663  12.5758  13.8253
    14.0995

          k = 0.0000 0.4861 0.1443 (   521 PWs)   bands (ev):

    -4.9557  -2.1401   2.8759   4.8217   6.2048   9.4635  11.2503  12.2881
    13.7754

          k = 0.5614 0.0000-0.2405 (   510 PWs)   bands (ev):

    -4.4147  -1.8320   1.8882   3.5571   4.1851   9.8867  13.0491  14.3533
    15.0387

          k = 0.0000 0.0000 0.4328 (   522 PWs)   bands (ev):

    -5.8540  -1.5186   5.8506   5.8506   7.0717   8.5219   8.5219   9.6514
    15.8064

          k = 0.4210 0.7292 0.1443 (   520 PWs)   bands (ev):

    -4.8268  -2.0487   2.1671   4.6856   5.9738  10.0966  10.4647  13.2780
    15.2970

     the Fermi energy is     8.1677 ev

!    total energy              =     -25.49943486 Ry
     Harris-Foulkes estimate   =     -25.49943486 Ry
     estimated scf accuracy    <          1.4E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.93129596 Ry
     hartree contribution      =       1.17870479 Ry
     xc contribution           =      -6.31795917 Ry
     ewald contribution        =     -27.29155289 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   5 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000    0.00139570
     atom    2 type  1   force =    -0.00000000    0.00000000   -0.00139570

     Total force =     0.001974     Total SCF correction =     0.000003


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    1.36
  -0.00000006  -0.00000000   0.00000000         -0.01     -0.00      0.00
  -0.00000000  -0.00000006  -0.00000000         -0.00     -0.01     -0.00
   0.00000000  -0.00000000   0.00002780          0.00     -0.00      4.09


     number of scf cycles    =   7
     number of bfgs steps    =   6

     enthalpy old            =     -25.4993705209 Ry
     enthalpy new            =     -25.4994348609 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0034073691 bohr
     new conv_thr            =       0.0000000064 Ry

     new unit-cell volume =    273.54534 a.u.^3 (    40.53524 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593734386   0.000000000   0.866915186
  -0.296867003   0.514189013   0.866915189
  -0.296867003  -0.514189013   0.866915189

ATOMIC_POSITIONS (crystal)
As       0.271951512   0.271951500   0.271951500
As      -0.271951512  -0.271951500  -0.271951500



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.00379, renormalised to   10.00000

     total cpu time spent up to now is        3.9 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  2.4

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

     Davidson diagonalization with overlap
     ethr =  3.82E-09,  avg # of iterations =  2.1

     total cpu time spent up to now is        4.0 secs

     total energy              =     -25.49944583 Ry
     Harris-Foulkes estimate   =     -25.50167335 Ry
     estimated scf accuracy    <       0.00000034 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  3.39E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        4.1 secs

     total energy              =     -25.49944596 Ry
     Harris-Foulkes estimate   =     -25.49944603 Ry
     estimated scf accuracy    <       0.00000014 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.37E-09,  avg # of iterations =  1.5

     total cpu time spent up to now is        4.2 secs

     total energy              =     -25.49944598 Ry
     Harris-Foulkes estimate   =     -25.49944598 Ry
     estimated scf accuracy    <          9.0E-09 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.97E-11,  avg # of iterations =  2.5

     total cpu time spent up to now is        4.2 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1442 (   531 PWs)   bands (ev):

    -7.0863   1.8420   5.6612   5.6612   6.5973  10.0735  10.6162  10.6162
    14.5564

          k =-0.1404-0.2431 0.2403 (   522 PWs)   bands (ev):

    -6.0597  -0.8085   4.0525   5.7462   8.1006   8.3687   9.0734  11.9644
    13.9845

          k = 0.2807 0.4862-0.0481 (   520 PWs)   bands (ev):

    -4.5301  -3.1589   4.6494   4.8144   6.2871   9.3906   9.6944  10.5063
    15.7360

          k = 0.1404 0.2431 0.0481 (   525 PWs)   bands (ev):

    -6.5086   0.2389   4.7835   5.3789   6.7771   9.4818  10.3184  11.5356
    13.5142

          k =-0.2807-0.0000 0.3364 (   519 PWs)   bands (ev):

    -5.6815  -0.5658   3.0203   4.0833   5.3653  10.2652  12.0253  12.1119
    13.7982

          k = 0.1404 0.7293 0.0481 (   510 PWs)   bands (ev):

    -4.0931  -2.5249   1.8979   2.9044   6.2559   9.9668  12.5659  13.8144
    14.0986

          k = 0.0000 0.4862 0.1442 (   521 PWs)   bands (ev):

    -4.9626  -2.1467   2.8651   4.8246   6.1924   9.4682  11.2462  12.2767
    13.7715

          k = 0.5614 0.0000-0.2403 (   510 PWs)   bands (ev):

    -4.4230  -1.8441   1.8929   3.5507   4.1858   9.8686  13.0473  14.3623
    15.0342

          k = 0.0000 0.0000 0.4326 (   522 PWs)   bands (ev):

    -5.8677  -1.5132   5.8394   5.8394   7.0610   8.5255   8.5255   9.6566
    15.7967

          k = 0.4211 0.7293 0.1442 (   520 PWs)   bands (ev):

    -4.8438  -2.0417   2.1598   4.6765   5.9766  10.0998  10.4558  13.2681
    15.2958

     the Fermi energy is     8.1633 ev

!    total energy              =     -25.49944598 Ry
     Harris-Foulkes estimate   =     -25.49944598 Ry
     estimated scf accuracy    <          1.6E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.91965466 Ry
     hartree contribution      =       1.18244407 Ry
     xc contribution           =      -6.31820696 Ry
     ewald contribution        =     -27.28341420 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   4 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =     0.00000000   -0.00000000    0.00048903
     atom    2 type  1   force =    -0.00000000    0.00000000   -0.00048903

     Total force =     0.000692     Total SCF correction =     0.000036


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    1.91
   0.00000490  -0.00000000  -0.00000000          0.72     -0.00     -0.00
  -0.00000000   0.00000490   0.00000000         -0.00      0.72      0.00
  -0.00000000   0.00000000   0.00002918         -0.00      0.00      4.29


     number of scf cycles    =   8
     number of bfgs steps    =   7

     enthalpy old            =     -25.4994348609 Ry
     enthalpy new            =     -25.4994459783 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0037504320 bohr
     new conv_thr            =       0.0000000011 Ry

     new unit-cell volume =    273.83828 a.u.^3 (    40.57865 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593715742  -0.000000000   0.867898095
  -0.296857682   0.514172867   0.867898098
  -0.296857682  -0.514172867   0.867898098

ATOMIC_POSITIONS (crystal)
As       0.272093399   0.272093387   0.272093387
As      -0.272093399  -0.272093387  -0.272093387



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.01070, renormalised to   10.00000

     total cpu time spent up to now is        4.4 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  2.2

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

     Davidson diagonalization with overlap
     ethr =  5.01E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        4.5 secs

     total energy              =     -25.49945467 Ry
     Harris-Foulkes estimate   =     -25.50573181 Ry
     estimated scf accuracy    <       0.00000046 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.61E-09,  avg # of iterations =  3.0

     total cpu time spent up to now is        4.6 secs

     total energy              =     -25.49945550 Ry
     Harris-Foulkes estimate   =     -25.49945572 Ry
     estimated scf accuracy    <       0.00000058 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.61E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        4.6 secs

     total energy              =     -25.49945546 Ry
     Harris-Foulkes estimate   =     -25.49945553 Ry
     estimated scf accuracy    <       0.00000013 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.27E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        4.7 secs

     total energy              =     -25.49945548 Ry
     Harris-Foulkes estimate   =     -25.49945548 Ry
     estimated scf accuracy    <          1.7E-09 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.68E-11,  avg # of iterations =  3.0

     total cpu time spent up to now is        4.7 secs

     End of self-consistent calculation

          k = 0.0000-0.0000 0.1440 (   531 PWs)   bands (ev):

    -7.0960   1.8250   5.6483   5.6483   6.5818  10.0468  10.5988  10.5988
    14.5523

          k =-0.1404-0.2431 0.2400 (   522 PWs)   bands (ev):

    -6.0712  -0.8189   4.0358   5.7242   8.0882   8.3537   9.0739  11.9455
    13.9697

          k = 0.2807 0.4862-0.0480 (   520 PWs)   bands (ev):

    -4.5425  -3.1683   4.6329   4.7962   6.2770   9.3700   9.6908  10.4870
    15.7121

          k = 0.1404 0.2431 0.0480 (   525 PWs)   bands (ev):

    -6.5179   0.2241   4.7703   5.3636   6.7577   9.4710  10.2921  11.5211
    13.5048

          k =-0.2807-0.0000 0.3361 (   519 PWs)   bands (ev):

    -5.6952  -0.5779   3.0044   4.0742   5.3629  10.2438  12.0091  12.0968
    13.7959

          k = 0.1404 0.7293 0.0480 (   510 PWs)   bands (ev):

    -4.1094  -2.5354   1.8903   2.8946   6.2423   9.9569  12.5490  13.7902
    14.0829

          k = 0.0000 0.4862 0.1440 (   521 PWs)   bands (ev):

    -4.9732  -2.1596   2.8503   4.8187   6.1698   9.4584  11.2292  12.2559
    13.7577

          k = 0.5614 0.0000-0.2400 (   510 PWs)   bands (ev):

    -4.4360  -1.8624   1.8903   3.5411   4.1775   9.8451  13.0305  14.3549
    15.0118

          k = 0.0000-0.0000 0.4321 (   522 PWs)   bands (ev):

    -5.8851  -1.5186   5.8246   5.8246   7.0448   8.5213   8.5213   9.6509
    15.7759

          k = 0.4211 0.7293 0.1440 (   520 PWs)   bands (ev):

    -4.8650  -2.0444   2.1492   4.6639   5.9719  10.0933  10.4379  13.2479
    15.2823

     the Fermi energy is     8.1509 ev

!    total energy              =     -25.49945548 Ry
     Harris-Foulkes estimate   =     -25.49945548 Ry
     estimated scf accuracy    <          3.2E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.90061087 Ry
     hartree contribution      =       1.18697204 Ry
     xc contribution           =      -6.31753370 Ry
     ewald contribution        =     -27.26958114 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   5 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000    0.00000000   -0.00001565
     atom    2 type  1   force =     0.00000000   -0.00000000    0.00001565

     Total force =     0.000022     Total SCF correction =     0.000015
     SCF correction compared to forces is large: reduce conv_thr to get better values


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    1.77
   0.00000578   0.00000000  -0.00000000          0.85      0.00     -0.00
   0.00000000   0.00000578  -0.00000000          0.00      0.85     -0.00
  -0.00000000  -0.00000000   0.00002445         -0.00     -0.00      3.60


     number of scf cycles    =   9
     number of bfgs steps    =   8

     enthalpy old            =     -25.4994459783 Ry
     enthalpy new            =     -25.4994554792 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0056222850 bohr
     new conv_thr            =       0.0000000010 Ry

     new unit-cell volume =    274.28938 a.u.^3 (    40.64549 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593700494  -0.000000000   0.869372458
  -0.296850057   0.514159662   0.869372461
  -0.296850057  -0.514159662   0.869372461

ATOMIC_POSITIONS (crystal)
As       0.272214989   0.272214977   0.272214977
As      -0.272214989  -0.272214977  -0.272214977



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.01645, renormalised to   10.00000

     total cpu time spent up to now is        4.9 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  2.5

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

     Davidson diagonalization with overlap
     ethr =  9.38E-09,  avg # of iterations =  1.4

     total cpu time spent up to now is        5.0 secs

     total energy              =     -25.49946189 Ry
     Harris-Foulkes estimate   =     -25.50909192 Ry
     estimated scf accuracy    <       0.00000091 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  9.14E-09,  avg # of iterations =  3.0

     total cpu time spent up to now is        5.1 secs

     total energy              =     -25.49946384 Ry
     Harris-Foulkes estimate   =     -25.49946431 Ry
     estimated scf accuracy    <       0.00000125 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  9.14E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        5.1 secs

     total energy              =     -25.49946379 Ry
     Harris-Foulkes estimate   =     -25.49946391 Ry
     estimated scf accuracy    <       0.00000029 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.90E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        5.2 secs

     total energy              =     -25.49946378 Ry
     Harris-Foulkes estimate   =     -25.49946381 Ry
     estimated scf accuracy    <       0.00000006 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.71E-10,  avg # of iterations =  2.6

     total cpu time spent up to now is        5.2 secs

     total energy              =     -25.49946379 Ry
     Harris-Foulkes estimate   =     -25.49946379 Ry
     estimated scf accuracy    <          1.8E-09 Ry

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.76E-11,  avg # of iterations =  1.5

     total cpu time spent up to now is        5.3 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1438 (   531 PWs)   bands (ev):

    -7.1086   1.7962   5.6318   5.6318   6.5561  10.0130  10.5737  10.5737
    14.5429

          k =-0.1404-0.2431 0.2396 (   522 PWs)   bands (ev):

    -6.0860  -0.8357   4.0112   5.6977   8.0700   8.3293   9.0671  11.9149
    13.9495

          k = 0.2807 0.4862-0.0479 (   520 PWs)   bands (ev):

    -4.5584  -3.1822   4.6078   4.7740   6.2602   9.3403   9.6784  10.4528
    15.6736

          k = 0.1404 0.2431 0.0479 (   525 PWs)   bands (ev):

    -6.5299   0.2023   4.7539   5.3392   6.7290   9.4486  10.2586  11.4959
    13.4862

          k =-0.2807 0.0000 0.3355 (   519 PWs)   bands (ev):

    -5.7129  -0.5986   2.9841   4.0627   5.3537  10.2179  11.9842  12.0733
    13.7848

          k = 0.1404 0.7293 0.0479 (   510 PWs)   bands (ev):

    -4.1305  -2.5508   1.8795   2.8788   6.2207   9.9396  12.5263  13.7542
    14.0551

          k = 0.0000 0.4862 0.1438 (   521 PWs)   bands (ev):

    -4.9872  -2.1776   2.8316   4.8059   6.1371   9.4373  11.2019  12.2273
    13.7357

          k = 0.5615 0.0000-0.2396 (   510 PWs)   bands (ev):

    -4.4534  -1.8866   1.8825   3.5292   4.1619   9.8157  13.0016  14.3354
    14.9744

          k = 0.0000 0.0000 0.4313 (   522 PWs)   bands (ev):

    -5.9072  -1.5324   5.8069   5.8069   7.0251   8.5122   8.5122   9.6381
    15.7452

          k = 0.4211 0.7293 0.1438 (   520 PWs)   bands (ev):

    -4.8915  -2.0542   2.1358   4.6484   5.9621  10.0801  10.4132  13.2195
    15.2589

     the Fermi energy is     8.1326 ev

!    total energy              =     -25.49946379 Ry
     Harris-Foulkes estimate   =     -25.49946379 Ry
     estimated scf accuracy    <          4.4E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.87456826 Ry
     hartree contribution      =       1.19252137 Ry
     xc contribution           =      -6.31617112 Ry
     ewald contribution        =     -27.25045874 Ry
     smearing contrib. (-TS)   =       0.00007645 Ry

     convergence has been achieved in   6 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000   -0.00000000   -0.00027060
     atom    2 type  1   force =     0.00000000    0.00000000    0.00027060

     Total force =     0.000383     Total SCF correction =     0.000017


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    1.11
   0.00000426  -0.00000000  -0.00000000          0.63     -0.00     -0.00
  -0.00000000   0.00000426   0.00000000         -0.00      0.63      0.00
  -0.00000000   0.00000000   0.00001411         -0.00      0.00      2.08


     number of scf cycles    =  10
     number of bfgs steps    =   9

     enthalpy old            =     -25.4994554792 Ry
     enthalpy new            =     -25.4994637914 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0055664058 bohr
     new conv_thr            =       0.0000000010 Ry

     new unit-cell volume =    274.73647 a.u.^3 (    40.71174 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593685521  -0.000000000   0.870833445
  -0.296842571   0.514146695   0.870833447
  -0.296842571  -0.514146695   0.870833447

ATOMIC_POSITIONS (crystal)
As       0.272274708   0.272274696   0.272274696
As      -0.272274708  -0.272274696  -0.272274696



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.01627, renormalised to   10.00000

     total cpu time spent up to now is        5.5 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  2.4

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

     Davidson diagonalization with overlap
     ethr =  9.05E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        5.6 secs

     total energy              =     -25.49946516 Ry
     Harris-Foulkes estimate   =     -25.50897545 Ry
     estimated scf accuracy    <       0.00000084 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.38E-09,  avg # of iterations =  3.0

     total cpu time spent up to now is        5.6 secs

     total energy              =     -25.49946707 Ry
     Harris-Foulkes estimate   =     -25.49946750 Ry
     estimated scf accuracy    <       0.00000116 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.38E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        5.7 secs

     total energy              =     -25.49946703 Ry
     Harris-Foulkes estimate   =     -25.49946712 Ry
     estimated scf accuracy    <       0.00000027 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.72E-09,  avg # of iterations =  1.0

     total cpu time spent up to now is        5.8 secs

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

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  7.04E-10,  avg # of iterations =  2.5

     total cpu time spent up to now is        5.8 secs

     total energy              =     -25.49946702 Ry
     Harris-Foulkes estimate   =     -25.49946702 Ry
     estimated scf accuracy    <          1.9E-09 Ry

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.91E-11,  avg # of iterations =  1.8

     total cpu time spent up to now is        5.9 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1435 (   531 PWs)   bands (ev):

    -7.1196   1.7659   5.6180   5.6180   6.5293   9.9840  10.5498  10.5498
    14.5319

          k =-0.1404-0.2431 0.2392 (   522 PWs)   bands (ev):

    -6.0990  -0.8527   3.9875   5.6762   8.0526   8.3043   9.0560  11.8836
    13.9313

          k = 0.2807 0.4862-0.0478 (   520 PWs)   bands (ev):

    -4.5722  -3.1956   4.5829   4.7557   6.2427   9.3122   9.6620  10.4158
    15.6341

          k = 0.1404 0.2431 0.0478 (   525 PWs)   bands (ev):

    -6.5404   0.1813   4.7400   5.3142   6.7012   9.4226  10.2299  11.4690
    13.4652

          k =-0.2807 0.0000 0.3349 (   519 PWs)   bands (ev):

    -5.7284  -0.6204   2.9667   4.0533   5.3413  10.1968  11.9596  12.0499
    13.7692

          k = 0.1404 0.7294 0.0478 (   510 PWs)   bands (ev):

    -4.1491  -2.5656   1.8694   2.8626   6.1988   9.9211  12.5060  13.7195
    14.0253

          k = 0.0000 0.4862 0.1435 (   521 PWs)   bands (ev):

    -4.9996  -2.1943   2.8155   4.7911   6.1059   9.4122  11.1741  12.2011
    13.7137

          k = 0.5615 0.0000-0.2392 (   510 PWs)   bands (ev):

    -4.4690  -1.9082   1.8722   3.5194   4.1446   9.7908  12.9709  14.3109
    14.9354

          k = 0.0000 0.0000 0.4306 (   522 PWs)   bands (ev):

    -5.9261  -1.5496   5.7926   5.7926   7.0089   8.5016   8.5016   9.6229
    15.7156

          k = 0.4211 0.7294 0.1435 (   520 PWs)   bands (ev):

    -4.9139  -2.0674   2.1244   4.6357   5.9510  10.0652  10.3907  13.1932
    15.2339

     the Fermi energy is     8.1152 ev

!    total energy              =     -25.49946702 Ry
     Harris-Foulkes estimate   =     -25.49946702 Ry
     estimated scf accuracy    <          4.4E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.85079404 Ry
     hartree contribution      =       1.19716323 Ry
     xc contribution           =      -6.31462782 Ry
     ewald contribution        =     -27.23287292 Ry
     smearing contrib. (-TS)   =       0.00007644 Ry

     convergence has been achieved in   6 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000   -0.00000000   -0.00018966
     atom    2 type  1   force =     0.00000000    0.00000000    0.00018966

     Total force =     0.000268     Total SCF correction =     0.000017


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.21
   0.00000116  -0.00000000  -0.00000000          0.17     -0.00     -0.00
  -0.00000000   0.00000116   0.00000000         -0.00      0.17      0.00
  -0.00000000   0.00000000   0.00000199         -0.00      0.00      0.29


     number of scf cycles    =  11
     number of bfgs steps    =  10

     enthalpy old            =     -25.4994637914 Ry
     enthalpy new            =     -25.4994670164 Ry

     CASE: enthalpy_new < enthalpy_old

     new trust radius        =       0.0004346788 bohr
     new conv_thr            =       0.0000000010 Ry

     new unit-cell volume =    274.77212 a.u.^3 (    40.71703 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593685121  -0.000000000   0.870947631
  -0.296842371   0.514146348   0.870947634
  -0.296842371  -0.514146348   0.870947634

ATOMIC_POSITIONS (crystal)
As       0.272257059   0.272257047   0.272257047
As      -0.272257059  -0.272257047  -0.272257047



     Writing output data file pwscf.save
     NEW-OLD atomic charge density approx. for the potential
     extrapolated charge   10.00130, renormalised to   10.00000

     total cpu time spent up to now is        6.2 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations =  1.0

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

     Davidson diagonalization with overlap
     ethr =  1.74E-10,  avg # of iterations =  3.6

     total cpu time spent up to now is        6.3 secs

     total energy              =     -25.49946713 Ry
     Harris-Foulkes estimate   =     -25.50022515 Ry
     estimated scf accuracy    <          8.6E-09 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.56E-11,  avg # of iterations =  3.0

     total cpu time spent up to now is        6.4 secs

     total energy              =     -25.49946714 Ry
     Harris-Foulkes estimate   =     -25.49946714 Ry
     estimated scf accuracy    <          4.6E-09 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  4.59E-11,  avg # of iterations =  1.0

     total cpu time spent up to now is        6.5 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1435 (   531 PWs)   bands (ev):

    -7.1200   1.7627   5.6176   5.6176   6.5266   9.9832  10.5481  10.5482
    14.5303

          k =-0.1404-0.2431 0.2392 (   522 PWs)   bands (ev):

    -6.0995  -0.8543   3.9858   5.6761   8.0513   8.3020   9.0533  11.8806
    13.9305

          k = 0.2807 0.4862-0.0478 (   520 PWs)   bands (ev):

    -4.5727  -3.1966   4.5809   4.7555   6.2410   9.3104   9.6591  10.4116
    15.6304

          k = 0.1404 0.2431 0.0478 (   525 PWs)   bands (ev):

    -6.5408   0.1798   4.7397   5.3119   6.6991   9.4191  10.2291  11.4660
    13.4625

          k =-0.2807 0.0000 0.3349 (   519 PWs)   bands (ev):

    -5.7289  -0.6228   2.9661   4.0531   5.3390  10.1966  11.9576  12.0479
    13.7662

          k = 0.1404 0.7294 0.0478 (   510 PWs)   bands (ev):

    -4.1498  -2.5668   1.8688   2.8611   6.1968   9.9191  12.5051  13.7169
    14.0220

          k = 0.0000 0.4862 0.1435 (   521 PWs)   bands (ev):

    -5.0002  -2.1953   2.8151   4.7891   6.1038   9.4087  11.1714  12.1997
    13.7118

          k = 0.5615 0.0000-0.2392 (   510 PWs)   bands (ev):

    -4.4697  -1.9091   1.8704   3.5193   4.1425   9.7903  12.9676  14.3070
    14.9314

          k = 0.0000 0.0000 0.4306 (   522 PWs)   bands (ev):

    -5.9266  -1.5524   5.7926   5.7926   7.0088   8.5001   8.5001   9.6208
    15.7135

          k = 0.4211 0.7294 0.1435 (   520 PWs)   bands (ev):

    -4.9144  -2.0698   2.1241   4.6355   5.9495  10.0632  10.3895  13.1916
    15.2312

     the Fermi energy is     8.1140 ev

!    total energy              =     -25.49946714 Ry
     Harris-Foulkes estimate   =     -25.49946714 Ry
     estimated scf accuracy    <          7.1E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.84964540 Ry
     hartree contribution      =       1.19722925 Ry
     xc contribution           =      -6.31443859 Ry
     ewald contribution        =     -27.23197964 Ry
     smearing contrib. (-TS)   =       0.00007644 Ry

     convergence has been achieved in   3 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000    0.00000000   -0.00005862
     atom    2 type  1   force =     0.00000000    0.00000000    0.00005862

     Total force =     0.000083     Total SCF correction =     0.000010
     SCF correction compared to forces is large: reduce conv_thr to get better values


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.07
   0.00000043  -0.00000000  -0.00000000          0.06     -0.00     -0.00
  -0.00000000   0.00000043   0.00000000         -0.00      0.06      0.00
  -0.00000000   0.00000000   0.00000058         -0.00      0.00      0.09


     bfgs converged in  12 scf cycles and  11 bfgs steps
     (criteria: energy <  1.0E-05, force <  1.0E-04, cell <  5.0E-01)

     End of BFGS Geometry Optimization

     Final enthalpy =     -25.4994671418 Ry
Begin final coordinates
     new unit-cell volume =    274.77212 a.u.^3 (    40.71703 Ang^3 )

CELL_PARAMETERS (alat=  7.01033623)
   0.593685121  -0.000000000   0.870947631
  -0.296842371   0.514146348   0.870947634
  -0.296842371  -0.514146348   0.870947634

ATOMIC_POSITIONS (crystal)
As       0.272257059   0.272257047   0.272257047
As      -0.272257059  -0.272257047  -0.272257047
End final coordinates



     A final scf calculation at the relaxed structure.
     The G-vectors are recalculated for the final unit cell
     Results may differ from those at the preceding step.

     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min         185     185     60                 2337     2337     446
     Max         186     186     61                 2338     2338     447
     Sum         371     371    121                 4675     4675     893



     bravais-lattice index     =            0
     lattice parameter (alat)  =       7.0103  a.u.
     unit-cell volume          =     274.7721 (a.u.)^3
     number of atoms/cell      =            2
     number of atomic types    =            1
     number of electrons       =        10.00
     number of Kohn-Sham states=            9
     kinetic-energy cutoff     =      25.0000  Ry
     charge density cutoff     =     100.0000  Ry
     convergence threshold     =      1.0E-09
     mixing beta               =       0.7000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      =  SLA  PZ   NOGX NOGC ( 1  1  0  0 0 0)

     celldm(1)=   7.010336  celldm(2)=   0.000000  celldm(3)=   0.000000
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   0.593685  -0.000000   0.870948 )  
               a(2) = (  -0.296842   0.514146   0.870948 )  
               a(3) = (  -0.296842  -0.514146   0.870948 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  1.122930 -0.000000  0.382725 )  
               b(2) = ( -0.561465  0.972486  0.382725 )  
               b(3) = ( -0.561465 -0.972486  0.382725 )  


     PseudoPot. # 1 for As read from file:
     /home/pietro/espresso-svn/pseudo/As.pz-bhs.UPF
     MD5 check sum: 451cd3365afcfc94d28b1934951c34a8
     Pseudo is Norm-conserving, Zval =  5.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  525 points,  2 beta functions with: 
                l(1) =   0
                l(2) =   1

     atomic species   valence    mass     pseudopotential
        As             5.00    74.90000     As( 1.00)

     12 Sym. Ops., with inversion, found



   Cartesian axes

     site n.     atom                  positions (alat units)
         1           As  tau(   1) = (   0.0000001   0.0000000   0.7113649  )
         2           As  tau(   2) = (  -0.0000001   0.0000000  -0.7113649  )

     number of k points=    10  Marzari-Vanderbilt smearing, width (Ry)=  0.0050
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.1435218), wk =   0.0625000
        k(    2) = (  -0.1403662  -0.2431214   0.2392031), wk =   0.1875000
        k(    3) = (   0.2807325   0.4862429  -0.0478407), wk =   0.1875000
        k(    4) = (   0.1403662   0.2431214   0.0478406), wk =   0.1875000
        k(    5) = (  -0.2807325   0.0000000   0.3348843), wk =   0.1875000
        k(    6) = (   0.1403662   0.7293643   0.0478406), wk =   0.3750000
        k(    7) = (   0.0000000   0.4862429   0.1435218), wk =   0.3750000
        k(    8) = (   0.5614650   0.0000000  -0.2392031), wk =   0.1875000
        k(    9) = (   0.0000000   0.0000000   0.4305655), wk =   0.0625000
        k(   10) = (   0.4210987   0.7293643   0.1435218), wk =   0.1875000

     Dense  grid:     4675 G-vectors     FFT dimensions: (  24,  24,  24)

     Estimated max dynamical RAM per process >       2.08Mb

     Estimated total allocated dynamical RAM >       4.17Mb

     Initial potential from superposition of free atoms

     starting charge    9.99960, renormalised to   10.00000
     Starting wfc are    8 randomized atomic wfcs +    1 random wfc

     Writing output data file pwscf.save

     total cpu time spent up to now is        6.8 secs

     per-process dynamical memory:     6.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-06,  avg # of iterations = 11.1

     total cpu time spent up to now is        7.1 secs

     total energy              =     -25.50068169 Ry
     Harris-Foulkes estimate   =     -25.50280761 Ry
     estimated scf accuracy    <       0.01040035 Ry

     iteration #  2     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.04E-04,  avg # of iterations =  1.0

     total cpu time spent up to now is        7.1 secs

     total energy              =     -25.50080524 Ry
     Harris-Foulkes estimate   =     -25.50099151 Ry
     estimated scf accuracy    <       0.00082266 Ry

     iteration #  3     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  8.23E-06,  avg # of iterations =  1.0

     total cpu time spent up to now is        7.2 secs

     total energy              =     -25.50084094 Ry
     Harris-Foulkes estimate   =     -25.50086041 Ry
     estimated scf accuracy    <       0.00003689 Ry

     iteration #  4     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  3.69E-07,  avg # of iterations =  2.1

     total cpu time spent up to now is        7.3 secs

     total energy              =     -25.50085018 Ry
     Harris-Foulkes estimate   =     -25.50085022 Ry
     estimated scf accuracy    <       0.00000029 Ry

     iteration #  5     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.94E-09,  avg # of iterations =  2.2

     total cpu time spent up to now is        7.4 secs

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

     iteration #  6     ecut=    25.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.28E-10,  avg # of iterations =  2.2

     total cpu time spent up to now is        7.5 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.1435 (   567 PWs)   bands (ev):

    -7.1197   1.7623   5.6174   5.6174   6.5241   9.9823  10.5476  10.5476
    14.5301

          k =-0.1404-0.2431 0.2392 (   579 PWs)   bands (ev):

    -6.0997  -0.8551   3.9849   5.6754   8.0488   8.2988   9.0528  11.8793
    13.9293

          k = 0.2807 0.4862-0.0478 (   582 PWs)   bands (ev):

    -4.5731  -3.1971   4.5793   4.7544   6.2384   9.3084   9.6580  10.4097
    15.6293

          k = 0.1404 0.2431 0.0478 (   581 PWs)   bands (ev):

    -6.5408   0.1781   4.7389   5.3105   6.6970   9.4175  10.2277  11.4650
    13.4618

          k =-0.2807 0.0000 0.3349 (   579 PWs)   bands (ev):

    -5.7293  -0.6238   2.9653   4.0521   5.3378  10.1952  11.9549  12.0456
    13.7636

          k = 0.1404 0.7294 0.0478 (   579 PWs)   bands (ev):

    -4.1504  -2.5677   1.8677   2.8599   6.1949   9.9170  12.5032  13.7126
    14.0184

          k = 0.0000 0.4862 0.1435 (   574 PWs)   bands (ev):

    -5.0004  -2.1954   2.8143   4.7880   6.1021   9.4073  11.1691  12.1977
    13.7103

          k = 0.5615 0.0000-0.2392 (   585 PWs)   bands (ev):

    -4.4705  -1.9106   1.8691   3.5181   4.1407   9.7886  12.9650  14.3028
    14.9256

          k = 0.0000 0.0000 0.4306 (   592 PWs)   bands (ev):

    -5.9271  -1.5533   5.7910   5.7910   7.0052   8.4994   8.4994   9.6156
    15.7102

          k = 0.4211 0.7294 0.1435 (   583 PWs)   bands (ev):

    -4.9148  -2.0712   2.1231   4.6345   5.9481  10.0620  10.3867  13.1878
    15.2303

     the Fermi energy is     8.1114 ev

!    total energy              =     -25.50085031 Ry
     Harris-Foulkes estimate   =     -25.50085031 Ry
     estimated scf accuracy    <          7.9E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =       6.84794426 Ry
     hartree contribution      =       1.19788936 Ry
     xc contribution           =      -6.31478085 Ry
     ewald contribution        =     -27.23197952 Ry
     smearing contrib. (-TS)   =       0.00007644 Ry

     convergence has been achieved in   6 iterations

     Forces acting on atoms (Cartesian axes, Ry/au):

     atom    1 type  1   force =    -0.00000000   -0.00000000   -0.00005381
     atom    2 type  1   force =     0.00000000    0.00000000    0.00005381

     Total force =     0.000076     Total SCF correction =     0.000012
     SCF correction compared to forces is large: reduce conv_thr to get better values


     Computing stress (Cartesian axis) and pressure

          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.06
  -0.00000207  -0.00000000   0.00000000         -0.30     -0.00      0.00
  -0.00000000  -0.00000207  -0.00000000         -0.00     -0.30     -0.00
   0.00000000  -0.00000000   0.00000288          0.00     -0.00      0.42


     Writing output data file pwscf.save

     init_run     :      0.18s CPU      0.19s WALL (       2 calls)
     electrons    :      4.98s CPU      5.05s WALL (      13 calls)
     update_pot   :      0.72s CPU      0.71s WALL (      11 calls)
     forces       :      0.32s CPU      0.33s WALL (      13 calls)
     stress       :      0.67s CPU      0.68s WALL (      13 calls)

     Called by init_run:
     wfcinit      :      0.04s CPU      0.05s WALL (       2 calls)
     potinit      :      0.04s CPU      0.04s WALL (       2 calls)

     Called by electrons:
     c_bands      :      4.23s CPU      4.28s WALL (      78 calls)
     sum_band     :      0.64s CPU      0.66s WALL (      78 calls)
     v_of_rho     :      0.06s CPU      0.07s WALL (      85 calls)
     mix_rho      :      0.02s CPU      0.02s WALL (      78 calls)

     Called by c_bands:
     init_us_2    :      0.08s CPU      0.10s WALL (    1840 calls)
     cegterg      :      4.16s CPU      4.20s WALL (     780 calls)

     Called by sum_band:

     Called by *egterg:
     h_psi        :      3.23s CPU      3.26s WALL (    2565 calls)
     g_psi        :      0.01s CPU      0.01s WALL (    1765 calls)
     cdiaghg      :      0.32s CPU      0.34s WALL (    2375 calls)

     Called by h_psi:
     h_psi:pot    :      3.22s CPU      3.25s WALL (    2565 calls)
     h_psi:calbec :      0.09s CPU      0.09s WALL (    2565 calls)
     vloc_psi     :      3.06s CPU      3.08s WALL (    2565 calls)
     add_vuspsi   :      0.07s CPU      0.07s WALL (    2565 calls)

     General routines
     calbec       :      0.12s CPU      0.13s WALL (    3215 calls)
     fft          :      0.06s CPU      0.06s WALL (     415 calls)
     fftw         :      3.37s CPU      3.37s WALL (   44052 calls)
     davcio       :      0.00s CPU      0.00s WALL (      20 calls)

     Parallel routines
     fft_scatter  :      0.66s CPU      0.59s WALL (   44467 calls)

     PWSCF        :     7.51s CPU         7.61s WALL


   This run was terminated on:  16: 4:37   9Dec2016            

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