si.out 25.5 KB
Newer Older
psd's avatar
psd committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667

     Program PWSCF v.6.3MaX starts on 16May2019 at  9:36:37 

     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);
         "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);
          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 & OpenMP), running on      12 processor cores
     Number of MPI processes:                12
     Threads/MPI process:                     1

     MPI processes distributed on     1 nodes
     R & G space division:  proc/nbgrp/npool/nimage =      12
     Waiting for input...
     Reading input from standard input
Warning: card  &IONS ignored
Warning: card  / ignored

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
               file Si.pbe-mt_fhi.UPF: wavefunction(s)  4f renormalized

     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)

     Found symmetry operation: I + (  0.0000 -0.5000 -0.5000)
     This is a supercell, fractional translations are disabled
 
     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min          83      83     24                 1989     1989     324
     Max          85      85     25                 1990     1990     329
     Sum        1005    1005    293                23871    23871    3911
 


     bravais-lattice index     =            1
     lattice parameter (alat)  =      10.2612  a.u.
     unit-cell volume          =    1080.4299 (a.u.)^3
     number of atoms/cell      =            8
     number of atomic types    =            1
     number of electrons       =        32.00
     number of Kohn-Sham states=           20
     kinetic-energy cutoff     =      30.0000  Ry
     charge density cutoff     =     120.0000  Ry
     convergence threshold     =      1.0E-08
     mixing beta               =       0.5000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      = PBE ( 1  4  3  4 0 0)

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

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   1.000000   0.000000   0.000000 )  
               a(2) = (   0.000000   1.000000   0.000000 )  
               a(3) = (   0.000000   0.000000   1.000000 )  

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


     PseudoPot. # 1 for Si read from file:
     ./Si.pbe-mt_fhi.UPF
     MD5 check sum: 259057c2adf9952f9ced37ae3d9e311d
     Pseudo is Norm-conserving, Zval =  4.0
     Generated using FHI98PP, converted with fhi2upf.x v.5.0.2
     Using radial grid of  495 points,  3 beta functions with: 
                l(1) =   0
                l(2) =   1
                l(3) =   3

     atomic species   valence    mass     pseudopotential
        Si             4.00    28.08600     Si( 1.00)

     24 Sym. Ops. (no inversion) found


                                    s                        frac. trans.

      isym =  1     identity                                     

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

 cart.    s( 1) = (  1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  1.0000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


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

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

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


     point group T_d (-43m) 
     there are  5 classes
     the character table:

       E     8C3   3C2   6S4   6s_d 
A_1    1.00  1.00  1.00  1.00  1.00
A_2    1.00  1.00  1.00 -1.00 -1.00
E      2.00 -1.00  2.00  0.00  0.00
T_1    3.00  0.00 -1.00  1.00 -1.00
T_2    3.00  0.00 -1.00 -1.00  1.00

     the symmetry operations in each class and the name of the first element:

     E        1
          identity                                               
     8C3      5    7    8    6   12    9   10   11
          120 deg rotation - cart. axis [-1,-1,-1]               
     3C2      2    3    4
          180 deg rotation - cart. axis [0,0,1]                  
     6S4     15   16   20   19   23   24
          inv.  90 deg rotation - cart. axis [0,0,-1]            
     6s_d    13   14   17   18   21   22
          inv. 180 deg rotation - cart. axis [1,1,0]             

   Cartesian axes

     site n.     atom                  positions (alat units)
         1           Si  tau(   1) = (   0.0000000   0.0000000   0.0000000  )
         2           Si  tau(   2) = (   0.0000000   0.5000000   0.5000000  )
         3           Si  tau(   3) = (   0.5000000   0.0000000   0.5000000  )
         4           Si  tau(   4) = (   0.5000000   0.5000000   0.0000000  )
         5           Si  tau(   5) = (   0.7500000   0.2500000   0.7500000  )
         6           Si  tau(   6) = (   0.2500000   0.2500000   0.2500000  )
         7           Si  tau(   7) = (   0.2500000   0.7500000   0.7500000  )
         8           Si  tau(   8) = (   0.7500000   0.7500000   0.2500000  )

   Crystallographic axes

     site n.     atom                  positions (cryst. coord.)
         1           Si  tau(   1) = (  0.0000000  0.0000000  0.0000000  )
         2           Si  tau(   2) = (  0.0000000  0.5000000  0.5000000  )
         3           Si  tau(   3) = (  0.5000000  0.0000000  0.5000000  )
         4           Si  tau(   4) = (  0.5000000  0.5000000  0.0000000  )
         5           Si  tau(   5) = (  0.7500000  0.2500000  0.7500000  )
         6           Si  tau(   6) = (  0.2500000  0.2500000  0.2500000  )
         7           Si  tau(   7) = (  0.2500000  0.7500000  0.7500000  )
         8           Si  tau(   8) = (  0.7500000  0.7500000  0.2500000  )

     number of k points=     4  gaussian smearing, width (Ry)=  0.0010
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.2500000
        k(    2) = (   0.0000000   0.0000000  -0.5000000), wk =   0.7500000
        k(    3) = (   0.0000000  -0.5000000  -0.5000000), wk =   0.7500000
        k(    4) = (  -0.5000000  -0.5000000  -0.5000000), wk =   0.2500000

                       cryst. coord.
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.2500000
        k(    2) = (   0.0000000   0.0000000  -0.5000000), wk =   0.7500000
        k(    3) = (   0.0000000  -0.5000000  -0.5000000), wk =   0.7500000
        k(    4) = (  -0.5000000  -0.5000000  -0.5000000), wk =   0.2500000

     Dense  grid:    23871 G-vectors     FFT dimensions: (  36,  36,  36)

     Dynamical RAM for                 wfc:       0.08 MB

     Dynamical RAM for     wfc (w. buffer):       0.38 MB

     Dynamical RAM for           str. fact:       0.03 MB

     Dynamical RAM for           local pot:       0.00 MB

     Dynamical RAM for          nlocal pot:       0.33 MB

     Dynamical RAM for                qrad:       0.57 MB

     Dynamical RAM for          rho,v,vnew:       0.18 MB

     Dynamical RAM for               rhoin:       0.06 MB

     Dynamical RAM for            rho*nmix:       0.49 MB

     Dynamical RAM for           G-vectors:       0.13 MB

     Dynamical RAM for          h,s,v(r/c):       0.07 MB

     Dynamical RAM for          <psi|beta>:       0.03 MB

     Dynamical RAM for                 psi:       0.30 MB

     Dynamical RAM for                hpsi:       0.30 MB

     Dynamical RAM for      wfcinit/wfcrot:       1.35 MB

     Estimated static dynamical RAM per process >       2.38 MB

     Estimated max dynamical RAM per process >       3.73 MB

     Estimated total dynamical RAM >      44.78 MB

     Initial potential from superposition of free atoms

     starting charge   31.99436, renormalised to   32.00000
     Starting wfcs are  128 randomized atomic wfcs

     total cpu time spent up to now is        0.3 secs

     Self-consistent Calculation

     iteration #  1     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  2.0

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

     Davidson diagonalization with overlap
     ethr =  7.77E-04,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.4 secs

     total energy              =     -62.97275731 Ry
     Harris-Foulkes estimate   =     -63.02198623 Ry
     estimated scf accuracy    <       0.25419698 Ry

     iteration #  2     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  7.94E-04,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.5 secs

     total energy              =     -62.97154068 Ry
     Harris-Foulkes estimate   =     -62.97827111 Ry
     estimated scf accuracy    <       0.04696407 Ry

     iteration #  3     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  1.47E-04,  avg # of iterations =  1.0

     total cpu time spent up to now is        0.5 secs

     total energy              =     -62.97422892 Ry
     Harris-Foulkes estimate   =     -62.97424394 Ry
     estimated scf accuracy    <       0.00074031 Ry

     iteration #  4     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  2.31E-06,  avg # of iterations =  2.5

     total cpu time spent up to now is        0.6 secs

     total energy              =     -62.97450557 Ry
     Harris-Foulkes estimate   =     -62.97450606 Ry
     estimated scf accuracy    <       0.00000792 Ry

     iteration #  5     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  2.48E-08,  avg # of iterations =  3.2

     total cpu time spent up to now is        0.7 secs

     total energy              =     -62.97450966 Ry
     Harris-Foulkes estimate   =     -62.97450939 Ry
     estimated scf accuracy    <       0.00000086 Ry

     iteration #  6     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  2.68E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is        0.7 secs

     total energy              =     -62.97450979 Ry
     Harris-Foulkes estimate   =     -62.97450984 Ry
     estimated scf accuracy    <       0.00000015 Ry

     iteration #  7     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  4.67E-10,  avg # of iterations =  1.5

     total cpu time spent up to now is        0.8 secs

     total energy              =     -62.97450979 Ry
     Harris-Foulkes estimate   =     -62.97450980 Ry
     estimated scf accuracy    <       0.00000002 Ry

     iteration #  8     ecut=    30.00 Ry     beta= 0.50
     Davidson diagonalization with overlap
     ethr =  5.63E-11,  avg # of iterations =  2.2

     total cpu time spent up to now is        0.8 secs

     End of self-consistent calculation

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

    -5.6975  -1.5426  -1.5426  -1.5426  -1.5426  -1.5426  -1.5426   3.4134
     3.4134   3.4134   3.4134   3.4134   3.4134   6.2765   6.2765   6.2765
     6.9781   6.9781   6.9781   6.9781

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.0000-0.5000 (  2946 PWs)   bands (ev):

    -4.6114  -4.6114  -1.3772  -1.3772  -1.3772  -1.3772   2.4107   2.4107
     2.4107   2.4107   2.7616   2.7616   4.3789   4.3789   4.3790   4.3790
     7.3985   7.3985   9.3968   9.3968

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     0.0000   0.0000   0.0000   0.0000

          k = 0.0000-0.5000-0.5000 (  3016 PWs)   bands (ev):

    -3.6476  -3.6476  -3.6476  -3.6476   0.7945   0.7945   0.7945   0.7945
     2.5122   2.5122   2.5122   2.5122   4.9432   4.9432   4.9432   4.9432
     8.7128   8.7128   8.7128   8.7128

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     0.0000   0.0000   0.0000   0.0000

          k =-0.5000-0.5000-0.5000 (  3016 PWs)   bands (ev):

    -3.3564  -3.3564  -3.3564  -3.3564  -0.7047  -0.7047  -0.7047  -0.7047
     5.0729   5.0729   5.0729   5.0729   5.0729   5.0729   5.0729   5.0729
     7.7952   7.7952   7.7952   7.7952

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     0.0000   0.0000   0.0000   0.0000

     the Fermi energy is     6.5836 ev

!    total energy              =     -62.97450980 Ry
     Harris-Foulkes estimate   =     -62.97450980 Ry
     estimated scf accuracy    <          5.2E-10 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =      19.05885882 Ry
     hartree contribution      =       4.48586405 Ry
     xc contribution           =     -19.32348450 Ry
     ewald contribution        =     -67.19574816 Ry
     smearing contrib. (-TS)   =      -0.00000000 Ry

     convergence has been achieved in   8 iterations

     Writing output data file SI.save/
 
     init_run     :      0.14s CPU      0.18s WALL (       1 calls)
     electrons    :      0.45s CPU      0.50s WALL (       1 calls)

     Called by init_run:
     wfcinit      :      0.12s CPU      0.14s WALL (       1 calls)
     wfcinit:atom :      0.00s CPU      0.00s WALL (       4 calls)
     wfcinit:wfcr :      0.11s CPU      0.14s WALL (       4 calls)
     potinit      :      0.00s CPU      0.00s WALL (       1 calls)
     hinit0       :      0.00s CPU      0.02s WALL (       1 calls)

     Called by electrons:
     c_bands      :      0.37s CPU      0.41s WALL (       9 calls)
     sum_band     :      0.06s CPU      0.06s WALL (       9 calls)
     v_of_rho     :      0.02s CPU      0.02s WALL (       9 calls)
     v_h          :      0.00s CPU      0.00s WALL (       9 calls)
     v_xc         :      0.02s CPU      0.02s WALL (       9 calls)
     mix_rho      :      0.00s CPU      0.00s WALL (       9 calls)

     Called by c_bands:
     init_us_2    :      0.00s CPU      0.01s WALL (      76 calls)
     cegterg      :      0.32s CPU      0.37s WALL (      36 calls)

     Called by sum_band:

     Called by *egterg:
     h_psi        :      0.28s CPU      0.34s WALL (     106 calls)
     g_psi        :      0.00s CPU      0.00s WALL (      66 calls)
     cdiaghg      :      0.12s CPU      0.13s WALL (      98 calls)
     cegterg:over :      0.00s CPU      0.01s WALL (      66 calls)
     cegterg:upda :      0.01s CPU      0.01s WALL (      66 calls)
     cegterg:last :      0.01s CPU      0.00s WALL (      36 calls)
     cdiaghg:chol :      0.01s CPU      0.01s WALL (      98 calls)
     cdiaghg:inve :      0.00s CPU      0.00s WALL (      98 calls)
     cdiaghg:para :      0.00s CPU      0.01s WALL (     196 calls)

     Called by h_psi:
     h_psi:pot    :      0.28s CPU      0.34s WALL (     106 calls)
     h_psi:calbec :      0.02s CPU      0.03s WALL (     106 calls)
     vloc_psi     :      0.25s CPU      0.30s WALL (     106 calls)
     add_vuspsi   :      0.01s CPU      0.01s WALL (     106 calls)

     General routines
     calbec       :      0.02s CPU      0.03s WALL (     106 calls)
     fft          :      0.01s CPU      0.03s WALL (      91 calls)
     ffts         :      0.00s CPU      0.00s WALL (       9 calls)
     fftw         :      0.24s CPU      0.32s WALL (    5138 calls)
 
     Parallel routines
     fft_scatt_xy :      0.02s CPU      0.03s WALL (    5238 calls)
     fft_scatt_yz :      0.17s CPU      0.18s WALL (    5238 calls)
 
     PWSCF        :     0.67s CPU         0.92s WALL

 
   This run was terminated on:   9:36:38  16May2019            

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