Program PWSCF v.6.1 (svn rev. 13369) starts on 3Mar2017 at 0: 1:24 This program is part of the open-source Quantum ESPRESSO suite for quantum simulation of materials; please cite "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009); URL http://www.quantum-espresso.org", in publications or presentations arising from this work. More details at http://www.quantum-espresso.org/quote Serial multi-threaded version, running on 4 processor cores Reading input from vc-relax4.in Current dimensions of program PWSCF are: Max number of different atomic species (ntypx) = 10 Max number of k-points (npk) = 40000 Max angular momentum in pseudopotentials (lmaxx) = 3 G-vector sticks info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW 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/fs395/scratch/QE/espresso/test-suite/..//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 Methfessel-Paxton 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 > 3.79MB 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.2 secs per-process dynamical memory: 3.2 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.5 total cpu time spent up to now is 0.5 secs total energy = -25.43995304 Ry Harris-Foulkes estimate = -25.44370905 Ry estimated scf accuracy < 0.01555592 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.6 secs total energy = -25.44007840 Ry Harris-Foulkes estimate = -25.44026102 Ry estimated scf accuracy < 0.00088841 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.88E-06, avg # of iterations = 2.3 total cpu time spent up to now is 0.7 secs total energy = -25.44011434 Ry Harris-Foulkes estimate = -25.44011580 Ry estimated scf accuracy < 0.00000523 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.23E-08, avg # of iterations = 3.1 total cpu time spent up to now is 0.8 secs total energy = -25.44012214 Ry Harris-Foulkes estimate = -25.44012246 Ry estimated scf accuracy < 0.00000069 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.94E-09, avg # of iterations = 1.4 total cpu time spent up to now is 0.9 secs End of self-consistent calculation k = 0.0000 0.0000 0.1535 ( 531 PWs) bands (ev): -6.9960 4.5196 5.9667 5.9667 8.4360 11.0403 11.7601 11.7601 16.5645 k =-0.1436-0.2488 0.2558 ( 522 PWs) bands (ev): -5.9250 0.3917 5.3512 5.6501 9.2996 10.5303 11.7005 13.5632 15.7167 k = 0.2873 0.4976-0.0512 ( 520 PWs) bands (ev): -4.3490 -2.4704 4.7883 6.1554 7.8796 10.8149 12.5849 13.8261 17.7262 k = 0.1436 0.2488 0.0512 ( 525 PWs) bands (ev): -6.3695 1.3043 4.9860 7.1720 8.5435 10.8049 12.4702 13.9612 15.3511 k =-0.2873 0.0000 0.3581 ( 519 PWs) bands (ev): -5.5427 1.1264 3.5658 4.2978 7.5159 10.4217 13.7076 13.7746 16.9045 k = 0.1436 0.7464 0.0512 ( 510 PWs) bands (ev): -3.8393 -1.8099 2.3270 4.2466 8.0539 11.6204 13.3234 15.7202 17.3489 k = 0.0000 0.4976 0.1535 ( 521 PWs) bands (ev): -4.7124 -1.4722 3.0016 6.6926 7.7777 12.3034 13.0675 13.4304 16.0962 k = 0.5746 0.0000-0.2558 ( 510 PWs) bands (ev): -4.0542 -1.5061 3.7084 3.7296 6.0243 10.0593 15.9112 17.7151 18.4776 k = 0.0000 0.0000 0.4604 ( 522 PWs) bands (ev): -5.8586 0.8361 5.8840 5.8840 7.4114 10.0627 10.0627 12.1191 17.3944 k = 0.4309 0.7464 0.1535 ( 520 PWs) bands (ev): -4.8492 -0.0498 2.4338 4.7831 7.5088 11.6828 12.0642 14.4760 17.7700 the Fermi energy is 10.0033 ev ! total energy = -25.44012222 Ry Harris-Foulkes estimate = -25.44012223 Ry estimated scf accuracy < 0.00000002 Ry The total energy is the sum of the following terms: one-electron contribution = 7.72810355 Ry hartree contribution = 1.22165969 Ry xc contribution = -6.50440122 Ry ewald contribution = -27.88552965 Ry smearing contrib. (-TS) = 0.00004540 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.12659777 atom 2 type 1 force = -0.00000000 0.00000000 0.12659777 Total force = 0.179036 Total SCF correction = 0.000024 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 217.52 0.00172374 -0.00000000 -0.00000000 253.57 -0.00 -0.00 -0.00000000 0.00172374 0.00000000 -0.00 253.57 0.00 -0.00000000 0.00000000 0.00098853 -0.00 0.00 145.42 BFGS Geometry Optimization number of scf cycles = 1 number of bfgs steps = 0 enthalpy new = -24.6061248137 Ry new trust radius = 0.1887860850 bohr new conv_thr = 0.0000001000 Ry new unit-cell volume = 211.67521 a.u.^3 ( 31.36703 Ang^3 ) density = 7.93026 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.555833834 -0.000000000 0.765441019 -0.277916743 0.481366175 0.765441024 -0.277916743 -0.481366175 0.765441024 ATOMIC_POSITIONS (crystal) As 0.282619706 0.282619701 0.282619701 As -0.282619706 -0.282619701 -0.282619701 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 8.40823, renormalised to 10.00000 total cpu time spent up to now is 1.2 secs per-process dynamical memory: 3.5 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.1 total cpu time spent up to now is 1.5 secs total energy = -25.38370667 Ry Harris-Foulkes estimate = -24.32076633 Ry estimated scf accuracy < 0.02064018 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.06E-04, avg # of iterations = 2.6 total cpu time spent up to now is 1.6 secs total energy = -25.40138519 Ry Harris-Foulkes estimate = -25.40441072 Ry estimated scf accuracy < 0.00628332 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.28E-05, avg # of iterations = 1.1 total cpu time spent up to now is 1.7 secs total energy = -25.40159307 Ry Harris-Foulkes estimate = -25.40189699 Ry estimated scf accuracy < 0.00065384 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.54E-06, avg # of iterations = 2.0 total cpu time spent up to now is 1.8 secs total energy = -25.40163612 Ry Harris-Foulkes estimate = -25.40166440 Ry estimated scf accuracy < 0.00005044 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.04E-07, avg # of iterations = 2.8 total cpu time spent up to now is 2.0 secs total energy = -25.40166079 Ry Harris-Foulkes estimate = -25.40166153 Ry estimated scf accuracy < 0.00000331 Ry iteration # 6 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.31E-08, avg # of iterations = 1.3 total cpu time spent up to now is 2.1 secs total energy = -25.40165994 Ry Harris-Foulkes estimate = -25.40166086 Ry estimated scf accuracy < 0.00000165 Ry iteration # 7 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.65E-08, avg # of iterations = 2.0 total cpu time spent up to now is 2.2 secs End of self-consistent calculation k = 0.0000 0.0000 0.1633 ( 531 PWs) bands (ev): -5.9367 7.0817 8.2032 8.2032 11.2582 13.8350 14.3550 14.3550 18.4339 k =-0.1499-0.2597 0.2722 ( 522 PWs) bands (ev): -4.6714 2.2296 7.2109 8.1029 11.9028 13.1331 13.4310 16.5292 18.1830 k = 0.2998 0.5194-0.0544 ( 520 PWs) bands (ev): -2.8547 -0.9760 6.9923 7.9650 10.0196 13.7451 14.3129 16.3159 20.8095 k = 0.1499 0.2597 0.0544 ( 525 PWs) bands (ev): -5.2654 3.5949 7.0489 9.3439 10.6444 13.7211 15.1491 16.5920 18.0163 k =-0.2998 0.0000 0.3810 ( 519 PWs) bands (ev): -4.1322 3.3399 5.1342 6.0876 9.0188 12.9443 16.5259 17.0166 19.1509 k = 0.1499 0.7790 0.0544 ( 510 PWs) bands (ev): -2.0535 -0.0681 3.6590 5.6372 10.0896 14.2403 15.7481 18.9065 20.2716 k = 0.0000 0.5194 0.1633 ( 521 PWs) bands (ev): -3.3504 0.4109 4.7586 8.2224 10.1013 14.6107 15.9234 16.2643 18.6227 k = 0.5997 0.0000-0.2722 ( 510 PWs) bands (ev): -2.4288 0.5228 4.8732 5.4349 7.5396 12.5225 18.7270 20.6442 21.7479 k = 0.0000 0.0000 0.4899 ( 522 PWs) bands (ev): -4.3331 1.9783 8.1315 8.1315 10.5523 11.9553 11.9553 14.3778 20.4089 k = 0.4498 0.7790 0.1633 ( 520 PWs) bands (ev): -3.0367 1.3418 3.7838 6.7673 9.3545 14.2226 14.7924 17.7408 20.1087 the Fermi energy is 11.9332 ev ! total energy = -25.40166017 Ry Harris-Foulkes estimate = -25.40166018 Ry estimated scf accuracy < 3.5E-09 Ry The total energy is the sum of the following terms: one-electron contribution = 9.92264769 Ry hartree contribution = 0.89911181 Ry xc contribution = -6.68612417 Ry ewald contribution = -29.53716795 Ry smearing contrib. (-TS) = -0.00012755 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.13944956 atom 2 type 1 force = 0.00000000 0.00000000 0.13944956 Total force = 0.197211 Total SCF correction = 0.000043 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 444.96 0.00333265 0.00000000 -0.00000000 490.25 0.00 -0.00 0.00000000 0.00333265 0.00000000 0.00 490.25 0.00 -0.00000000 0.00000000 0.00240899 -0.00 0.00 354.37 number of scf cycles = 2 number of bfgs steps = 1 enthalpy old = -24.6061248137 Ry enthalpy new = -24.6821906498 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.3019728046 bohr new conv_thr = 0.0000001000 Ry new unit-cell volume = 183.03324 a.u.^3 ( 27.12273 Ang^3 ) density = 9.17123 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.539853984 0.000000000 0.701631417 -0.269926854 0.467527219 0.701631444 -0.269926854 -0.467527219 0.701631444 ATOMIC_POSITIONS (crystal) As 0.263861305 0.263861288 0.263861288 As -0.263861305 -0.263861288 -0.263861288 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 8.43521, renormalised to 10.00000 total cpu time spent up to now is 2.5 secs per-process dynamical memory: 3.5 Mb Self-consistent Calculation iteration # 1 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-06, avg # of iterations = 7.1 total cpu time spent up to now is 2.8 secs total energy = -25.34462710 Ry Harris-Foulkes estimate = -24.15521729 Ry estimated scf accuracy < 0.01769026 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.77E-04, avg # of iterations = 2.1 total cpu time spent up to now is 2.9 secs total energy = -25.35392007 Ry Harris-Foulkes estimate = -25.35521166 Ry estimated scf accuracy < 0.00290892 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.91E-05, avg # of iterations = 1.1 total cpu time spent up to now is 3.0 secs total energy = -25.35400719 Ry Harris-Foulkes estimate = -25.35412857 Ry estimated scf accuracy < 0.00027641 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.76E-06, avg # of iterations = 2.1 total cpu time spent up to now is 3.1 secs total energy = -25.35404740 Ry Harris-Foulkes estimate = -25.35405429 Ry estimated scf accuracy < 0.00001379 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.38E-07, avg # of iterations = 1.7 total cpu time spent up to now is 3.3 secs total energy = -25.35404840 Ry Harris-Foulkes estimate = -25.35404882 Ry estimated scf accuracy < 0.00000073 Ry iteration # 6 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.33E-09, avg # of iterations = 2.5 total cpu time spent up to now is 3.4 secs End of self-consistent calculation k = 0.0000-0.0000 0.1782 ( 531 PWs) bands (ev): -4.6205 9.9126 10.7890 10.7890 14.3254 17.6816 17.6816 17.9911 19.7716 k =-0.1544-0.2674 0.2969 ( 522 PWs) bands (ev): -3.1157 4.1034 9.8092 11.8603 14.1839 14.3307 16.3499 19.7370 21.1375 k = 0.3087 0.5347-0.0594 ( 520 PWs) bands (ev): -1.0934 0.6900 9.9550 10.4615 11.9453 15.7402 17.5257 19.6590 24.8247 k = 0.1544 0.2674 0.0594 ( 525 PWs) bands (ev): -3.9458 6.2430 9.4929 11.3902 13.3914 17.6838 18.5316 19.3202 20.5552 k =-0.3087-0.0000 0.4157 ( 519 PWs) bands (ev): -2.3094 5.6935 7.7049 8.0169 9.9729 16.1618 19.4720 20.3831 20.8735 k = 0.1544 0.8021 0.0594 ( 510 PWs) bands (ev): 0.3805 1.7730 5.1132 7.3127 12.5986 16.6469 18.2783 23.1595 23.6093 k = 0.0000 0.5347 0.1782 ( 521 PWs) bands (ev): -1.7598 2.5260 7.2643 9.5065 13.4608 16.6331 19.1892 19.8717 21.3138 k = 0.6175 0.0000-0.2969 ( 510 PWs) bands (ev): -0.4673 3.6384 5.4660 7.3446 9.0994 15.9645 21.7300 22.8442 25.5616 k = 0.0000 0.0000 0.5345 ( 522 PWs) bands (ev): -2.1711 3.2605 10.8670 10.8670 13.5264 13.5264 13.8805 15.9384 24.4486 k = 0.4631 0.8021 0.1782 ( 520 PWs) bands (ev): -0.1117 2.1284 5.6222 9.0777 10.9316 16.1662 18.1827 21.5569 22.8826 the Fermi energy is 14.1988 ev ! total energy = -25.35404859 Ry Harris-Foulkes estimate = -25.35404860 Ry estimated scf accuracy < 6.1E-09 Ry The total energy is the sum of the following terms: one-electron contribution = 12.37614057 Ry hartree contribution = 0.55960979 Ry xc contribution = -6.86700410 Ry ewald contribution = -31.42270908 Ry smearing contrib. (-TS) = -0.00008578 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.06408392 atom 2 type 1 force = 0.00000000 -0.00000000 0.06408392 Total force = 0.090628 Total SCF correction = 0.000068 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 681.16 0.00460007 0.00000000 -0.00000000 676.69 0.00 -0.00 0.00000000 0.00460007 -0.00000000 0.00 676.69 -0.00 -0.00000000 -0.00000000 0.00469112 -0.00 -0.00 690.09 number of scf cycles = 3 number of bfgs steps = 2 enthalpy old = -24.6821906498 Ry enthalpy new = -24.7319311401 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0715331011 bohr new conv_thr = 0.0000000641 Ry new unit-cell volume = 190.16871 a.u.^3 ( 28.18010 Ang^3 ) density = 8.82711 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.546565452 -0.000000000 0.711191207 -0.273282598 0.473339520 0.711191243 -0.273282598 -0.473339520 0.711191243 ATOMIC_POSITIONS (crystal) As 0.259013583 0.259013562 0.259013562 As -0.259013583 -0.259013562 -0.259013562 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 10.37520, renormalised to 10.00000 total cpu time spent up to now is 3.7 secs per-process dynamical memory: 3.5 Mb Self-consistent Calculation iteration # 1 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-06, avg # of iterations = 6.6 total cpu time spent up to now is 3.9 secs total energy = -25.38978062 Ry Harris-Foulkes estimate = -25.68123756 Ry estimated scf accuracy < 0.00086564 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 8.66E-06, avg # of iterations = 2.0 total cpu time spent up to now is 4.1 secs total energy = -25.39023431 Ry Harris-Foulkes estimate = -25.39031697 Ry estimated scf accuracy < 0.00017477 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.75E-06, avg # of iterations = 1.1 total cpu time spent up to now is 4.2 secs total energy = -25.39024960 Ry Harris-Foulkes estimate = -25.39025656 Ry estimated scf accuracy < 0.00001723 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.72E-07, avg # of iterations = 1.2 total cpu time spent up to now is 4.3 secs total energy = -25.39025124 Ry Harris-Foulkes estimate = -25.39025158 Ry estimated scf accuracy < 0.00000076 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.63E-09, avg # of iterations = 3.0 total cpu time spent up to now is 4.4 secs total energy = -25.39025163 Ry Harris-Foulkes estimate = -25.39025166 Ry estimated scf accuracy < 0.00000008 Ry iteration # 6 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.92E-10, avg # of iterations = 1.0 total cpu time spent up to now is 4.5 secs End of self-consistent calculation k = 0.0000 0.0000 0.1758 ( 531 PWs) bands (ev): -4.8198 9.0293 10.2830 10.2830 13.3627 17.0775 17.0775 17.4618 18.9632 k =-0.1525-0.2641 0.2929 ( 522 PWs) bands (ev): -3.3507 3.6000 9.2268 11.6229 13.1494 13.5911 15.4243 19.1739 20.0847 k = 0.3049 0.5282-0.0586 ( 520 PWs) bands (ev): -1.3712 0.3037 9.5724 9.9078 11.2118 14.8603 16.7332 18.6006 23.9921 k = 0.1525 0.2641 0.0586 ( 525 PWs) bands (ev): -4.1630 5.6772 9.0495 10.6816 12.7384 16.7394 17.9220 18.4116 19.3749 k =-0.3049 0.0000 0.4101 ( 519 PWs) bands (ev): -2.5650 4.9815 7.4063 7.6119 9.1986 15.7760 18.6847 19.3937 19.7776 k = 0.1525 0.7922 0.0586 ( 510 PWs) bands (ev): 0.0485 1.3379 4.7638 6.7971 11.8955 15.8351 17.7093 22.2404 22.5660 k = 0.0000 0.5282 0.1758 ( 521 PWs) bands (ev): -2.0359 2.0865 6.9289 8.8408 12.8263 15.6358 18.5853 18.9898 20.3177 k = 0.6099-0.0000-0.2929 ( 510 PWs) bands (ev): -0.8097 3.4110 4.7274 6.9658 8.4735 15.5451 20.7289 21.5690 24.3833 k = 0.0000 0.0000 0.5273 ( 522 PWs) bands (ev): -2.4035 2.6568 10.4489 10.4489 12.8312 12.8312 13.2317 14.9817 23.5667 k = 0.4574 0.7922 0.1758 ( 520 PWs) bands (ev): -0.3055 1.4323 5.3212 8.6726 10.2600 15.3862 17.3843 20.6460 22.2409 the Fermi energy is 13.1969 ev ! total energy = -25.39025163 Ry Harris-Foulkes estimate = -25.39025163 Ry estimated scf accuracy < 8.0E-09 Ry The total energy is the sum of the following terms: one-electron contribution = 11.89669985 Ry hartree contribution = 0.57917961 Ry xc contribution = -6.79963525 Ry ewald contribution = -31.06647876 Ry smearing contrib. (-TS) = -0.00001707 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.04260571 atom 2 type 1 force = 0.00000000 0.00000000 0.04260571 Total force = 0.060254 Total SCF correction = 0.000016 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 525.83 0.00339000 0.00000000 -0.00000000 498.69 0.00 -0.00 0.00000000 0.00339000 0.00000000 0.00 498.69 0.00 -0.00000000 0.00000000 0.00394347 -0.00 0.00 580.10 number of scf cycles = 4 number of bfgs steps = 3 enthalpy old = -24.7319311401 Ry enthalpy new = -24.7438811986 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.1595170423 bohr new conv_thr = 0.0000000426 Ry new unit-cell volume = 191.30781 a.u.^3 ( 28.34889 Ang^3 ) density = 8.77455 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.545422573 -0.000000000 0.718452599 -0.272711183 0.472349757 0.718452650 -0.272711183 -0.472349757 0.718452650 ATOMIC_POSITIONS (crystal) As 0.248348594 0.248348567 0.248348567 As -0.248348594 -0.248348567 -0.248348567 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 10.05954, renormalised to 10.00000 total cpu time spent up to now is 4.8 secs per-process dynamical memory: 3.5 Mb Self-consistent Calculation iteration # 1 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-06, avg # of iterations = 8.2 total cpu time spent up to now is 5.1 secs total energy = -25.40060177 Ry Harris-Foulkes estimate = -25.44671222 Ry estimated scf accuracy < 0.00124011 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.24E-05, avg # of iterations = 1.0 total cpu time spent up to now is 5.2 secs total energy = -25.40060421 Ry Harris-Foulkes estimate = -25.40061787 Ry estimated scf accuracy < 0.00009339 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.34E-07, avg # of iterations = 1.0 total cpu time spent up to now is 5.3 secs total energy = -25.40060605 Ry Harris-Foulkes estimate = -25.40060604 Ry estimated scf accuracy < 0.00000005 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.70E-10, avg # of iterations = 4.0 total cpu time spent up to now is 5.5 secs total energy = -25.40060672 Ry Harris-Foulkes estimate = -25.40060673 Ry estimated scf accuracy < 0.00000005 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.70E-10, avg # of iterations = 1.0 total cpu time spent up to now is 5.6 secs End of self-consistent calculation k = 0.0000 0.0000 0.1740 ( 531 PWs) bands (ev): -4.8323 8.7450 10.3299 10.3299 13.2295 17.0744 17.0744 17.7063 18.7930 k =-0.1528-0.2646 0.2900 ( 522 PWs) bands (ev): -3.3645 3.5628 9.0023 12.1614 12.4602 13.5133 15.3424 19.3063 19.7034 k = 0.3056 0.5293-0.0580 ( 520 PWs) bands (ev): -1.3483 0.2511 9.6430 9.7335 11.0817 14.5825 16.6495 18.1914 23.6342 k = 0.1528 0.2646 0.0580 ( 525 PWs) bands (ev): -4.1680 5.6533 9.0896 10.4995 12.5979 16.5035 18.1060 18.1527 18.7736 k =-0.3056 0.0000 0.4060 ( 519 PWs) bands (ev): -2.5861 4.7703 7.5474 7.6656 8.8881 15.8381 18.6425 19.3133 19.6537 k = 0.1528 0.7939 0.0580 ( 510 PWs) bands (ev): 0.0518 1.3039 4.7332 6.6123 11.7686 15.7585 17.7044 22.0343 22.5982 k = 0.0000 0.5293 0.1740 ( 521 PWs) bands (ev): -2.0359 2.0711 6.9582 8.6040 12.7322 15.2386 18.4923 18.9894 20.1208 k = 0.6111-0.0000-0.2900 ( 510 PWs) bands (ev): -0.8224 3.6601 4.2749 7.0146 8.3186 15.4954 20.5134 21.2593 24.2156 k = 0.0000 0.0000 0.5220 ( 522 PWs) bands (ev): -2.4222 2.3798 10.5756 10.5756 12.7333 12.7333 13.4269 14.8283 23.4466 k = 0.4584 0.7939 0.1740 ( 520 PWs) bands (ev): -0.1908 1.0657 5.3558 8.7323 10.2023 15.3102 17.4874 20.5948 22.2516 the Fermi energy is 12.7897 ev ! total energy = -25.40060672 Ry Harris-Foulkes estimate = -25.40060672 Ry estimated scf accuracy < 0.00000001 Ry The total energy is the sum of the following terms: one-electron contribution = 11.84274296 Ry hartree contribution = 0.57430078 Ry xc contribution = -6.78686844 Ry ewald contribution = -31.03085500 Ry smearing contrib. (-TS) = 0.00007299 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.00724301 atom 2 type 1 force = -0.00000000 -0.00000000 -0.00724301 Total force = 0.010243 Total SCF correction = 0.000021 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 494.44 0.00317558 -0.00000000 -0.00000000 467.14 -0.00 -0.00 0.00000000 0.00317558 0.00000000 0.00 467.14 0.00 -0.00000000 0.00000000 0.00373220 -0.00 0.00 549.02 number of scf cycles = 5 number of bfgs steps = 4 enthalpy old = -24.7438811986 Ry enthalpy new = -24.7503645951 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0303454225 bohr new conv_thr = 0.0000000072 Ry new unit-cell volume = 192.20705 a.u.^3 ( 28.48215 Ang^3 ) density = 8.73350 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.543728797 0.000000000 0.726333847 -0.271864296 0.470882903 0.726333898 -0.271864296 -0.470882903 0.726333898 ATOMIC_POSITIONS (crystal) As 0.249489770 0.249489744 0.249489744 As -0.249489770 -0.249489744 -0.249489744 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 10.04678, renormalised to 10.00000 total cpu time spent up to now is 5.9 secs per-process dynamical memory: 3.5 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 6.1 secs total energy = -25.40470684 Ry Harris-Foulkes estimate = -25.44069097 Ry estimated scf accuracy < 0.00001906 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.91E-07, avg # of iterations = 2.0 total cpu time spent up to now is 6.2 secs total energy = -25.40471158 Ry Harris-Foulkes estimate = -25.40471221 Ry estimated scf accuracy < 0.00000236 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.36E-08, avg # of iterations = 1.1 total cpu time spent up to now is 6.3 secs total energy = -25.40471165 Ry Harris-Foulkes estimate = -25.40471166 Ry estimated scf accuracy < 0.00000009 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.40E-10, avg # of iterations = 2.2 total cpu time spent up to now is 6.4 secs End of self-consistent calculation k = 0.0000 0.0000 0.1721 ( 531 PWs) bands (ev): -4.8684 8.5218 10.3471 10.3471 13.1538 16.9881 16.9881 17.6923 18.7833 k =-0.1533-0.2655 0.2868 ( 522 PWs) bands (ev): -3.4069 3.5408 8.7825 12.2042 12.3208 13.4821 15.2822 19.3381 19.4614 k = 0.3065 0.5309-0.0574 ( 520 PWs) bands (ev): -1.3612 0.2051 9.5378 9.6195 11.0434 14.4602 16.5245 17.8732 23.2228 k = 0.1533 0.2655 0.0574 ( 525 PWs) bands (ev): -4.1937 5.6149 9.0944 10.3726 12.4652 16.3434 17.8815 17.9764 18.6180 k =-0.3065 0.0000 0.4016 ( 519 PWs) bands (ev): -2.6456 4.6118 7.5214 7.7041 8.8034 15.7623 18.6183 19.3228 19.6757 k = 0.1533 0.7964 0.0574 ( 510 PWs) bands (ev): -0.0156 1.2829 4.7084 6.4701 11.6582 15.7371 17.6956 21.8549 22.5381 k = 0.0000 0.5309 0.1721 ( 521 PWs) bands (ev): -2.0589 2.0442 6.8979 8.4938 12.5607 15.0159 18.3718 18.9588 20.0346 k = 0.6131-0.0000-0.2868 ( 510 PWs) bands (ev): -0.8624 3.5819 4.1804 7.0437 8.2275 15.3333 20.3662 21.1607 24.0839 k = 0.0000 0.0000 0.5163 ( 522 PWs) bands (ev): -2.5081 2.1955 10.6186 10.6186 12.7373 12.7373 13.5006 14.8818 23.2834 k = 0.4598 0.7964 0.1721 ( 520 PWs) bands (ev): -0.2908 0.9360 5.3488 8.7523 10.2179 15.3101 17.5576 20.5942 22.0236 the Fermi energy is 12.7946 ev ! total energy = -25.40471166 Ry Harris-Foulkes estimate = -25.40471166 Ry estimated scf accuracy < 2.8E-09 Ry The total energy is the sum of the following terms: one-electron contribution = 11.77377374 Ry hartree contribution = 0.58046190 Ry xc contribution = -6.77983684 Ry ewald contribution = -30.97912861 Ry smearing contrib. (-TS) = 0.00001815 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.00217846 atom 2 type 1 force = -0.00000000 0.00000000 -0.00217846 Total force = 0.003081 Total SCF correction = 0.000004 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 480.35 0.00311301 -0.00000000 -0.00000000 457.94 -0.00 -0.00 0.00000000 0.00311301 -0.00000000 0.00 457.94 -0.00 -0.00000000 0.00000000 0.00356997 -0.00 0.00 525.16 number of scf cycles = 6 number of bfgs steps = 5 enthalpy old = -24.7503645951 Ry enthalpy new = -24.7514130786 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0451299230 bohr new conv_thr = 0.0000000022 Ry new unit-cell volume = 192.30794 a.u.^3 ( 28.49710 Ang^3 ) density = 8.72892 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.539498736 0.000000000 0.738155717 -0.269749272 0.467219562 0.738155770 -0.269749272 -0.467219562 0.738155770 ATOMIC_POSITIONS (crystal) As 0.250230244 0.250230218 0.250230218 As -0.250230244 -0.250230218 -0.250230218 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 10.00525, renormalised to 10.00000 total cpu time spent up to now is 6.7 secs per-process dynamical memory: 3.5 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 7.0 secs total energy = -25.40615028 Ry Harris-Foulkes estimate = -25.41017949 Ry estimated scf accuracy < 0.00001076 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.08E-07, avg # of iterations = 1.0 total cpu time spent up to now is 7.1 secs total energy = -25.40615049 Ry Harris-Foulkes estimate = -25.40615045 Ry estimated scf accuracy < 0.00000091 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 9.06E-09, avg # of iterations = 1.0 total cpu time spent up to now is 7.2 secs End of self-consistent calculation k = 0.0000-0.0000 0.1693 ( 531 PWs) bands (ev): -4.8845 8.2744 10.4947 10.4947 13.2148 16.9902 16.9902 17.8216 18.9005 k =-0.1545-0.2675 0.2822 ( 522 PWs) bands (ev): -3.4250 3.6024 8.5226 12.2806 12.3366 13.5717 15.3496 19.2306 19.5376 k = 0.3089 0.5351-0.0564 ( 520 PWs) bands (ev): -1.3175 0.2074 9.3179 9.6824 11.1041 14.4010 16.4804 17.5058 22.6964 k = 0.1545 0.2675 0.0564 ( 525 PWs) bands (ev): -4.1881 5.6729 9.2086 10.2617 12.3700 16.2290 17.5992 17.8446 18.5714 k =-0.3089-0.0000 0.3951 ( 519 PWs) bands (ev): -2.6869 4.4588 7.5629 7.8634 8.7963 15.7402 18.7361 19.5191 19.9072 k = 0.1545 0.8026 0.0564 ( 510 PWs) bands (ev): -0.0481 1.3360 4.7393 6.3245 11.5988 15.8452 17.8154 21.7161 22.6164 k = 0.0000 0.5351 0.1693 ( 521 PWs) bands (ev): -2.0351 2.0940 6.8807 8.4199 12.3839 14.7946 18.3278 19.0670 20.0763 k = 0.6179 0.0000-0.2822 ( 510 PWs) bands (ev): -0.8532 3.5335 4.1230 7.1823 8.1787 15.1753 20.2845 21.1800 24.0390 k = 0.0000 0.0000 0.5080 ( 522 PWs) bands (ev): -2.5952 1.9839 10.8005 10.8005 12.8741 12.8741 13.7766 15.1431 23.1789 k = 0.4634 0.8026 0.1693 ( 520 PWs) bands (ev): -0.3963 0.8217 5.4146 8.8845 10.3686 15.4528 17.8460 20.7223 21.8295 the Fermi energy is 12.9314 ev ! total energy = -25.40615051 Ry Harris-Foulkes estimate = -25.40615051 Ry estimated scf accuracy < 8.2E-10 Ry The total energy is the sum of the following terms: one-electron contribution = 11.75888329 Ry hartree contribution = 0.58344315 Ry xc contribution = -6.77985140 Ry ewald contribution = -30.96864367 Ry smearing contrib. (-TS) = 0.00001812 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.00093854 atom 2 type 1 force = 0.00000000 -0.00000000 0.00093854 Total force = 0.001327 Total SCF correction = 0.000002 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 478.07 0.00316738 0.00000000 -0.00000000 465.94 0.00 -0.00 0.00000000 0.00316738 0.00000000 0.00 465.94 0.00 -0.00000000 0.00000000 0.00341472 -0.00 0.00 502.32 number of scf cycles = 7 number of bfgs steps = 6 enthalpy old = -24.7514130786 Ry enthalpy new = -24.7525090096 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0451322681 bohr new conv_thr = 0.0000000010 Ry new unit-cell volume = 191.31213 a.u.^3 ( 28.34953 Ang^3 ) density = 8.77435 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.533774669 0.000000000 0.750167473 -0.266887248 0.462262374 0.750167529 -0.266887248 -0.462262374 0.750167529 ATOMIC_POSITIONS (crystal) As 0.250427008 0.250426982 0.250426982 As -0.250427008 -0.250426982 -0.250426982 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 9.94795, renormalised to 10.00000 total cpu time spent up to now is 7.5 secs per-process dynamical memory: 3.5 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.8 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 6.52E-08, avg # of iterations = 1.0 total cpu time spent up to now is 7.8 secs total energy = -25.40319975 Ry Harris-Foulkes estimate = -25.36315409 Ry estimated scf accuracy < 0.00000653 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.53E-08, avg # of iterations = 2.3 total cpu time spent up to now is 7.9 secs total energy = -25.40320497 Ry Harris-Foulkes estimate = -25.40320548 Ry estimated scf accuracy < 0.00000197 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.97E-08, avg # of iterations = 1.5 total cpu time spent up to now is 8.0 secs total energy = -25.40320497 Ry Harris-Foulkes estimate = -25.40320503 Ry estimated scf accuracy < 0.00000028 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.75E-09, avg # of iterations = 1.8 total cpu time spent up to now is 8.1 secs total energy = -25.40320499 Ry Harris-Foulkes estimate = -25.40320500 Ry estimated scf accuracy < 7.1E-09 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.05E-11, avg # of iterations = 2.5 total cpu time spent up to now is 8.3 secs End of self-consistent calculation k = 0.0000-0.0000 0.1666 ( 531 PWs) bands (ev): -4.8648 8.0791 10.7614 10.7614 13.4382 17.1183 17.1183 18.1032 18.8110 k =-0.1561-0.2704 0.2777 ( 522 PWs) bands (ev): -3.3999 3.7480 8.3246 12.4254 12.5044 13.7855 15.5462 19.1058 19.8990 k = 0.3122 0.5408-0.0555 ( 520 PWs) bands (ev): -1.2129 0.2756 9.1631 9.8428 11.2660 14.4495 16.5654 17.2214 22.2448 k = 0.1561 0.2704 0.0555 ( 525 PWs) bands (ev): -4.1408 5.8389 9.4270 10.2041 12.3694 16.1906 17.3900 17.8302 18.6646 k =-0.3122-0.0000 0.3888 ( 519 PWs) bands (ev): -2.6828 4.3672 7.6884 8.1219 8.8992 15.8104 18.9944 19.8753 20.3088 k = 0.1561 0.8112 0.0555 ( 510 PWs) bands (ev): -0.0151 1.4671 4.8345 6.2378 11.6279 16.0704 18.0698 21.6664 22.6764 k = 0.0000 0.5408 0.1666 ( 521 PWs) bands (ev): -1.9558 2.2290 6.9341 8.4242 12.2676 14.6625 18.4131 19.3123 20.2687 k = 0.6245 0.0000-0.2777 ( 510 PWs) bands (ev): -0.7781 3.5648 4.1202 7.4133 8.2051 15.1057 20.3147 21.3265 24.0848 k = 0.0000 0.0000 0.4999 ( 522 PWs) bands (ev): -2.6427 1.8223 11.0990 11.0990 13.1312 13.1312 14.2173 15.5656 23.1898 k = 0.4684 0.8112 0.1666 ( 520 PWs) bands (ev): -0.4525 0.7648 5.5583 9.1166 10.6356 15.7228 18.3036 20.8722 21.8699 the Fermi energy is 13.1885 ev ! total energy = -25.40320499 Ry Harris-Foulkes estimate = -25.40320499 Ry estimated scf accuracy < 9.4E-10 Ry The total energy is the sum of the following terms: one-electron contribution = 11.82566107 Ry hartree contribution = 0.57945328 Ry xc contribution = -6.78865716 Ry ewald contribution = -31.01968031 Ry smearing contrib. (-TS) = 0.00001812 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.00184042 atom 2 type 1 force = 0.00000000 0.00000000 0.00184042 Total force = 0.002603 Total SCF correction = 0.000002 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 491.79 0.00336290 -0.00000000 -0.00000000 494.70 -0.00 -0.00 0.00000000 0.00336290 0.00000000 0.00 494.70 0.00 -0.00000000 0.00000000 0.00330352 -0.00 0.00 485.96 number of scf cycles = 8 number of bfgs steps = 7 enthalpy old = -24.7525090096 Ry enthalpy new = -24.7529481677 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0098897105 bohr new conv_thr = 0.0000000018 Ry new unit-cell volume = 190.99432 a.u.^3 ( 28.30244 Ang^3 ) density = 8.78895 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.534287031 0.000000000 0.747485585 -0.267143429 0.462706093 0.747485642 -0.267143429 -0.462706093 0.747485642 ATOMIC_POSITIONS (crystal) As 0.250142892 0.250142866 0.250142866 As -0.250142892 -0.250142866 -0.250142866 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 9.98336, renormalised to 10.00000 total cpu time spent up to now is 8.6 secs per-process dynamical memory: 3.5 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.1 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 1.64E-08, avg # of iterations = 1.4 total cpu time spent up to now is 8.8 secs total energy = -25.40217353 Ry Harris-Foulkes estimate = -25.38935431 Ry estimated scf accuracy < 0.00000165 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.65E-08, avg # of iterations = 2.0 total cpu time spent up to now is 9.0 secs total energy = -25.40217418 Ry Harris-Foulkes estimate = -25.40217424 Ry estimated scf accuracy < 0.00000022 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 2.16E-09, avg # of iterations = 1.0 total cpu time spent up to now is 9.1 secs total energy = -25.40217418 Ry Harris-Foulkes estimate = -25.40217418 Ry estimated scf accuracy < 0.00000002 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.60E-10, avg # of iterations = 1.4 total cpu time spent up to now is 9.2 secs End of self-consistent calculation k = 0.0000-0.0000 0.1672 ( 531 PWs) bands (ev): -4.8525 8.1616 10.7561 10.7561 13.4668 17.1473 17.1473 18.1149 18.8774 k =-0.1560-0.2701 0.2787 ( 522 PWs) bands (ev): -3.3857 3.7597 8.3907 12.4669 12.4934 13.7962 15.5767 19.1863 19.8837 k = 0.3119 0.5403-0.0557 ( 520 PWs) bands (ev): -1.2075 0.2904 9.2246 9.8511 11.2821 14.4864 16.6081 17.3273 22.3742 k = 0.1560 0.2701 0.0557 ( 525 PWs) bands (ev): -4.1319 5.8526 9.4254 10.2578 12.4110 16.2650 17.4679 17.8849 18.7063 k =-0.3119-0.0000 0.3902 ( 519 PWs) bands (ev): -2.6633 4.4236 7.6998 8.1100 8.9186 15.8332 19.0068 19.8771 20.3132 k = 0.1560 0.8104 0.0557 ( 510 PWs) bands (ev): 0.0074 1.4761 4.8424 6.2799 11.6691 16.0869 18.0660 21.7441 22.7623 k = 0.0000 0.5403 0.1672 ( 521 PWs) bands (ev): -1.9477 2.2387 6.9547 8.4567 12.3356 14.7332 18.4479 19.3294 20.2978 k = 0.6239 0.0000-0.2787 ( 510 PWs) bands (ev): -0.7644 3.5991 4.1419 7.4043 8.2385 15.1500 20.3719 21.3713 24.1598 k = 0.0000 0.0000 0.5017 ( 522 PWs) bands (ev): -2.6146 1.8796 11.0874 11.0874 13.1292 13.1292 14.1982 15.5487 23.2516 k = 0.4679 0.8104 0.1672 ( 520 PWs) bands (ev): -0.4174 0.8038 5.5593 9.1109 10.6328 15.7275 18.2904 20.9121 21.9142 the Fermi energy is 13.4095 ev ! total energy = -25.40217418 Ry Harris-Foulkes estimate = -25.40217418 Ry estimated scf accuracy < 2.6E-10 Ry The total energy is the sum of the following terms: one-electron contribution = 11.84843316 Ry hartree contribution = 0.57791914 Ry xc contribution = -6.79138614 Ry ewald contribution = -31.03714940 Ry smearing contrib. (-TS) = 0.00000906 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.00062490 atom 2 type 1 force = 0.00000000 -0.00000000 0.00062490 Total force = 0.000884 Total SCF correction = 0.000001 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 496.88 0.00337773 0.00000000 -0.00000000 496.88 0.00 -0.00 0.00000000 0.00337773 0.00000000 0.00 496.88 0.00 -0.00000000 -0.00000000 0.00337777 -0.00 -0.00 496.89 number of scf cycles = 9 number of bfgs steps = 8 enthalpy old = -24.7529481677 Ry enthalpy new = -24.7529975804 Ry CASE: enthalpy_new < enthalpy_old new trust radius = 0.0021564219 bohr new conv_thr = 0.0000000010 Ry new unit-cell volume = 190.79974 a.u.^3 ( 28.27360 Ang^3 ) density = 8.79791 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.534112779 0.000000000 0.747211384 -0.267056304 0.462555186 0.747211442 -0.267056304 -0.462555186 0.747211442 ATOMIC_POSITIONS (crystal) As 0.250005719 0.250005692 0.250005692 As -0.250005719 -0.250005692 -0.250005692 Writing output data file pwscf.save NEW-OLD atomic charge density approx. for the potential extrapolated charge 9.98980, renormalised to 10.00000 total cpu time spent up to now is 9.5 secs per-process dynamical memory: 3.5 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.0 Threshold (ethr) on eigenvalues was too large: Diagonalizing with lowered threshold Davidson diagonalization with overlap ethr = 3.74E-09, avg # of iterations = 2.3 total cpu time spent up to now is 9.7 secs total energy = -25.40151612 Ry Harris-Foulkes estimate = -25.39365474 Ry estimated scf accuracy < 0.00000036 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 3.55E-09, avg # of iterations = 2.0 total cpu time spent up to now is 9.8 secs total energy = -25.40151634 Ry Harris-Foulkes estimate = -25.40151636 Ry estimated scf accuracy < 0.00000007 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 6.87E-10, avg # of iterations = 1.0 total cpu time spent up to now is 9.9 secs total energy = -25.40151634 Ry Harris-Foulkes estimate = -25.40151634 Ry estimated scf accuracy < 7.7E-09 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 7.70E-11, avg # of iterations = 1.8 total cpu time spent up to now is 10.1 secs End of self-consistent calculation k = 0.0000-0.0000 0.1673 ( 531 PWs) bands (ev): -4.8461 8.1811 10.7721 10.7721 13.4935 17.1679 17.1679 18.1382 18.8949 k =-0.1560-0.2702 0.2788 ( 522 PWs) bands (ev): -3.3782 3.7741 8.4047 12.4946 12.4996 13.8151 15.6016 19.2104 19.9027 k = 0.3120 0.5405-0.0558 ( 520 PWs) bands (ev): -1.1982 0.3011 9.2393 9.8655 11.2994 14.5050 16.6312 17.3510 22.3967 k = 0.1560 0.2702 0.0558 ( 525 PWs) bands (ev): -4.1250 5.8694 9.4397 10.2749 12.4280 16.2907 17.4857 17.9064 18.7298 k =-0.3120-0.0000 0.3903 ( 519 PWs) bands (ev): -2.6550 4.4400 7.7129 8.1227 8.9342 15.8475 19.0291 19.9017 20.3414 k = 0.1560 0.8107 0.0558 ( 510 PWs) bands (ev): 0.0189 1.4886 4.8522 6.2915 11.6875 16.1082 18.0823 21.7709 22.7868 k = 0.0000 0.5405 0.1673 ( 521 PWs) bands (ev): -1.9389 2.2518 6.9667 8.4705 12.3546 14.7531 18.4681 19.3526 20.3226 k = 0.6241 0.0000-0.2788 ( 510 PWs) bands (ev): -0.7537 3.6147 4.1510 7.4163 8.2538 15.1640 20.3966 21.3985 24.1908 k = 0.0000 0.0000 0.5019 ( 522 PWs) bands (ev): -2.6064 1.8921 11.1029 11.1029 13.1460 13.1460 14.2202 15.5704 23.2768 k = 0.4681 0.8107 0.1673 ( 520 PWs) bands (ev): -0.4069 0.8159 5.5696 9.1243 10.6498 15.7474 18.3154 20.9338 21.9393 the Fermi energy is 13.4363 ev ! total energy = -25.40151634 Ry Harris-Foulkes estimate = -25.40151634 Ry estimated scf accuracy < 1.3E-10 Ry The total energy is the sum of the following terms: one-electron contribution = 11.86226070 Ry hartree contribution = 0.57698784 Ry xc contribution = -6.79306467 Ry ewald contribution = -31.04770927 Ry smearing contrib. (-TS) = 0.00000906 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.00002533 atom 2 type 1 force = 0.00000000 -0.00000000 0.00002533 Total force = 0.000036 Total SCF correction = 0.000000 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 499.90 0.00339784 -0.00000000 -0.00000000 499.84 -0.00 -0.00 -0.00000000 0.00339784 0.00000000 -0.00 499.84 0.00 -0.00000000 0.00000000 0.00339904 -0.00 0.00 500.02 bfgs converged in 10 scf cycles and 9 bfgs steps (criteria: energy < 1.0E-04 Ry, force < 1.0E-03Ry/Bohr, cell < 5.0E-01kbar) End of BFGS Geometry Optimization Final enthalpy = -24.7530010968 Ry Begin final coordinates new unit-cell volume = 190.79974 a.u.^3 ( 28.27360 Ang^3 ) density = 8.79791 g/cm^3 CELL_PARAMETERS (alat= 7.01033623) 0.534112779 0.000000000 0.747211384 -0.267056304 0.462555186 0.747211442 -0.267056304 -0.462555186 0.747211442 ATOMIC_POSITIONS (crystal) As 0.250005719 0.250005692 0.250005692 As -0.250005719 -0.250005692 -0.250005692 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. G-vector sticks info -------------------- sticks: dense smooth PW G-vecs: dense smooth PW Sum 287 287 93 3221 3221 633 bravais-lattice index = 0 lattice parameter (alat) = 7.0103 a.u. unit-cell volume = 190.7997 (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.534113 0.000000 0.747211 ) a(2) = ( -0.267056 0.462555 0.747211 ) a(3) = ( -0.267056 -0.462555 0.747211 ) reciprocal axes: (cart. coord. in units 2 pi/alat) b(1) = ( 1.248176 0.000000 0.446103 ) b(2) = ( -0.624088 1.080952 0.446103 ) b(3) = ( -0.624088 -1.080952 0.446103 ) PseudoPot. # 1 for As read from file: /home/fs395/scratch/QE/espresso/test-suite/..//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.5604213 ) 2 As tau( 2) = ( -0.0000001 0.0000000 -0.5604213 ) number of k points= 10 Methfessel-Paxton smearing, width (Ry)= 0.0050 cart. coord. in units 2pi/alat k( 1) = ( 0.0000000 -0.0000000 0.1672887), wk = 0.0625000 k( 2) = ( -0.1560220 -0.2702380 0.2788145), wk = 0.1875000 k( 3) = ( 0.3120440 0.5404761 -0.0557629), wk = 0.1875000 k( 4) = ( 0.1560220 0.2702380 0.0557629), wk = 0.1875000 k( 5) = ( -0.3120440 -0.0000000 0.3903402), wk = 0.1875000 k( 6) = ( 0.1560220 0.8107141 0.0557629), wk = 0.3750000 k( 7) = ( 0.0000000 0.5404761 0.1672887), wk = 0.3750000 k( 8) = ( 0.6240880 0.0000000 -0.2788145), wk = 0.1875000 k( 9) = ( 0.0000000 0.0000000 0.5018660), wk = 0.0625000 k( 10) = ( 0.4680660 0.8107141 0.1672886), wk = 0.1875000 Dense grid: 3221 G-vectors FFT dimensions: ( 24, 24, 24) Estimated max dynamical RAM per process > 3.16MB 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 10.4 secs per-process dynamical memory: 2.9 Mb Self-consistent Calculation iteration # 1 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.00E-06, avg # of iterations = 10.9 total cpu time spent up to now is 10.7 secs total energy = -25.39786349 Ry Harris-Foulkes estimate = -25.40003088 Ry estimated scf accuracy < 0.01490712 Ry iteration # 2 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.49E-04, avg # of iterations = 1.0 total cpu time spent up to now is 10.8 secs total energy = -25.39784199 Ry Harris-Foulkes estimate = -25.39797660 Ry estimated scf accuracy < 0.00110532 Ry iteration # 3 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.11E-05, avg # of iterations = 1.0 total cpu time spent up to now is 10.9 secs total energy = -25.39785212 Ry Harris-Foulkes estimate = -25.39785260 Ry estimated scf accuracy < 0.00000171 Ry iteration # 4 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.71E-08, avg # of iterations = 3.0 total cpu time spent up to now is 11.0 secs total energy = -25.39785383 Ry Harris-Foulkes estimate = -25.39785388 Ry estimated scf accuracy < 0.00000012 Ry iteration # 5 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 1.18E-09, avg # of iterations = 1.5 total cpu time spent up to now is 11.1 secs total energy = -25.39785384 Ry Harris-Foulkes estimate = -25.39785384 Ry estimated scf accuracy < 5.7E-09 Ry iteration # 6 ecut= 25.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 5.68E-11, avg # of iterations = 2.0 total cpu time spent up to now is 11.2 secs End of self-consistent calculation k = 0.0000-0.0000 0.1673 ( 396 PWs) bands (ev): -4.8415 8.1940 10.7721 10.7721 13.5173 17.1701 17.1701 18.1436 18.8986 k =-0.1560-0.2702 0.2788 ( 397 PWs) bands (ev): -3.3731 3.7836 8.4063 12.4956 12.5008 13.8201 15.6150 19.2148 19.9070 k = 0.3120 0.5405-0.0558 ( 401 PWs) bands (ev): -1.1925 0.3079 9.2414 9.8659 11.3046 14.5063 16.6392 17.3575 22.4006 k = 0.1560 0.2702 0.0558 ( 396 PWs) bands (ev): -4.1203 5.8812 9.4403 10.2794 12.4305 16.3064 17.4908 17.9082 18.7352 k =-0.3120-0.0000 0.3903 ( 407 PWs) bands (ev): -2.6500 4.4492 7.7150 8.1229 8.9379 15.8499 19.0313 19.9060 20.3503 k = 0.1560 0.8107 0.0558 ( 402 PWs) bands (ev): 0.0252 1.4955 4.8533 6.2926 11.6907 16.1116 18.0851 21.7742 22.7876 k = 0.0000 0.5405 0.1673 ( 405 PWs) bands (ev): -1.9336 2.2602 6.9681 8.4725 12.3583 14.7563 18.4723 19.3595 20.3261 k = 0.6241 0.0000-0.2788 ( 410 PWs) bands (ev): -0.7483 3.6230 4.1519 7.4165 8.2564 15.1665 20.3980 21.4016 24.1952 k = 0.0000 0.0000 0.5019 ( 407 PWs) bands (ev): -2.6011 1.8994 11.1036 11.1036 13.1472 13.1472 14.2308 15.5784 23.2830 k = 0.4681 0.8107 0.1673 ( 403 PWs) bands (ev): -0.4008 0.8235 5.5709 9.1248 10.6549 15.7489 18.3234 20.9368 21.9430 the Fermi energy is 13.4600 ev ! total energy = -25.39785384 Ry Harris-Foulkes estimate = -25.39785384 Ry estimated scf accuracy < 4.1E-11 Ry The total energy is the sum of the following terms: one-electron contribution = 11.86640263 Ry hartree contribution = 0.57607634 Ry xc contribution = -6.79263256 Ry ewald contribution = -31.04770931 Ry smearing contrib. (-TS) = 0.00000906 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.00002481 atom 2 type 1 force = 0.00000000 0.00000000 0.00002481 Total force = 0.000035 Total SCF correction = 0.000000 Computing stress (Cartesian axis) and pressure total stress (Ry/bohr**3) (kbar) P= 501.85 0.00341036 0.00000000 -0.00000000 501.68 0.00 -0.00 0.00000000 0.00341036 0.00000000 0.00 501.68 0.00 -0.00000000 -0.00000000 0.00341382 -0.00 -0.00 502.19 Writing output data file pwscf.save init_run : 0.53s CPU 0.29s WALL ( 2 calls) electrons : 31.31s CPU 8.00s WALL ( 11 calls) update_pot : 1.04s CPU 1.05s WALL ( 9 calls) forces : 0.64s CPU 0.45s WALL ( 11 calls) stress : 1.68s CPU 1.04s WALL ( 11 calls) Called by init_run: wfcinit : 0.29s CPU 0.07s WALL ( 2 calls) potinit : 0.05s CPU 0.05s WALL ( 2 calls) Called by electrons: c_bands : 26.64s CPU 6.79s WALL ( 59 calls) sum_band : 4.24s CPU 1.06s WALL ( 59 calls) v_of_rho : 0.19s CPU 0.05s WALL ( 66 calls) mix_rho : 0.11s CPU 0.03s WALL ( 59 calls) Called by c_bands: init_us_2 : 0.45s CPU 0.11s WALL ( 1420 calls) cegterg : 25.67s CPU 6.49s WALL ( 590 calls) Called by sum_band: Called by *egterg: h_psi : 21.32s CPU 5.37s WALL ( 2052 calls) g_psi : 0.15s CPU 0.04s WALL ( 1442 calls) cdiaghg : 1.12s CPU 0.28s WALL ( 1902 calls) Called by h_psi: h_psi:pot : 21.24s CPU 5.34s WALL ( 2052 calls) h_psi:calbec : 0.44s CPU 0.11s WALL ( 2052 calls) vloc_psi : 20.31s CPU 5.11s WALL ( 2052 calls) add_vuspsi : 0.48s CPU 0.12s WALL ( 2052 calls) General routines calbec : 0.57s CPU 0.14s WALL ( 2602 calls) fft : 0.25s CPU 0.08s WALL ( 325 calls) fftw : 21.47s CPU 5.40s WALL ( 34672 calls) davcio : 0.00s CPU 0.00s WALL ( 20 calls) PWSCF : 35.46s CPU 11.30s WALL This run was terminated on: 0: 1:35 3Mar2017 =------------------------------------------------------------------------------= JOB DONE. =------------------------------------------------------------------------------=