<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div style="line-height: normal;">Dear QE developers,</div><div style="line-height: normal;">The full input and output files are pasted as following:</div><div style="line-height: normal;"><br></div><div><div style="line-height: normal;"> &CONTROL</div><div style="line-height: normal;">                 calculation = 'relax' ,</div><div style="line-height: normal;">                restart_mode = 'from_scratch' ,</div><div style="line-height: normal;">                  pseudo_dir = '/home/wang/pseudo' ,</div><div style="line-height: normal;">                      prefix = 'pf' ,</div><div style="line-height: normal;">                       nstep = 300 ,</div><div style="line-height: normal;">                     tprnfor = .true. ,</div><div style="line-height: normal;"> /</div><div style="line-height: normal;"> &SYSTEM</div><div style="line-height: normal;">                       ibrav = 8,</div><div style="line-height: normal;">                   celldm(1) = 28.346,</div><div style="line-height: normal;">                   celldm(2) = 0.488,</div><div style="line-height: normal;">                   celldm(3) = 0.596,</div><div style="line-height: normal;">                         nat = 16,</div><div style="line-height: normal;">                        ntyp = 3,</div><div style="line-height: normal;">                     ecutwfc = 35 ,</div><div style="line-height: normal;">                     ecutrho = 260 ,</div><div style="line-height: normal;">                       nosym = .false. ,</div><div style="line-height: normal;">                 occupations = 'smearing' ,</div><div style="line-height: normal;">                     degauss = 0.02 ,</div><div style="line-height: normal;">                       nspin = 2 ,</div><div style="line-height: normal;">   starting_magnetization(3) = 0.5,</div><div style="line-height: normal;">            exxdiv_treatment = 'gygi-baldereschi' ,</div><div style="line-height: normal;">                  lda_plus_u = .true. ,</div><div style="line-height: normal;">             lda_plus_u_kind = 0 ,</div><div style="line-height: normal;">                Hubbard_U(3) = 4.5,</div><div style="line-height: normal;">               Hubbard_J0(3) = 0.5,</div><div style="line-height: normal;"> /</div><div style="line-height: normal;"> &ELECTRONS</div><div style="line-height: normal;">            electron_maxstep = 300,</div><div style="line-height: normal;">                 mixing_mode = 'local-TF' ,</div><div style="line-height: normal;">                 mixing_beta = 0.5 ,</div><div style="line-height: normal;">             diagonalization = 'david' ,</div><div style="line-height: normal;"> /</div><div style="line-height: normal;"> &IONS</div><div style="line-height: normal;">                ion_dynamics = 'bfgs' ,</div><div style="line-height: normal;">                      upscale= 1000 ,</div><div style="line-height: normal;"> /</div><div style="line-height: normal;">ATOMIC_SPECIES</div><div style="line-height: normal;">    S     32.06  S.pbe-n-kjpaw_psl.0.1.UPF</div><div style="line-height: normal;">    Ge    72.59  Ge.pbe-kjpaw.UPF</div><div><span style="line-height: normal;">    Co    58.93  Co.pbe-nd-rrkjus.UPF</span></div><div><span style="line-height: normal;">ATOMIC_POSITIONS crystal</span></div><div><span style="line-height: normal;">  S   0.2506999969482421   0.1243000030517581  -0.0000000000000000    1  1  1</span></div><div><span style="line-height: normal;">  S   0.1080000028014180   0.3743000030517582   0.2505500018596650    1  1  1</span></div><div><span style="line-height: normal;">  S   0.2506999969482421   0.6243000030517582   0.0000000000000000    1  1  1</span></div><div><span style="line-height: normal;">  S   0.1080000028014180   0.8743000030517582   0.2505500018596650    1  1  1</span></div><div><span style="line-height: normal;">  S   0.2506999969482421   0.1243000030517581   0.4999999999999998    1  1  1</span></div><div><span style="line-height: normal;">  S   0.1080000028014180   0.3743000030517582   0.7505500018596648    1  1  1</span></div><div><span style="line-height: normal;">  S   0.2506999969482421   0.6243000030517582   0.4999999999999998    1  1  1</span></div><div><span style="line-height: normal;">  S   0.1080000028014180   0.8743000030517582   0.7505500018596648    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.0939999967813492   0.1243000030517581   0.4340000003576278    1  1  1</span></div><div><span style="line-height: normal;"> Co   0.2646999955177311   0.3743000030517582   0.1834499984979630    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.0939999967813492   0.6243000030517582   0.4340000003576278    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.2646999955177311   0.8743000030517582   0.1834499984979631    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.0939999967813492   0.1243000030517581   0.9340000003576279    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.2646999955177311   0.3743000030517582   0.6834499984979631    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.0939999967813492   0.6243000030517582   0.9340000003576279    1  1  1</span></div><div><span style="line-height: normal;"> Ge   0.2646999955177311   0.8743000030517582   0.6834499984979631    1  1  1</span></div><div><span style="line-height: normal;">K_POINTS automatic</span></div><div><span style="line-height: normal;"> 1 3 3   0 0 0</span></div><div style="line-height: normal;"><br></div></div><br><div>The output file:</div><div><div><br></div><div>     Program PWSCF v.6.0 (svn rev. 13079) starts on  5May2017 at 15:16:24</div><div><br></div><div>     This program is part of the open-source Quantum ESPRESSO suite</div><div>     for quantum simulation of materials; please cite</div><div>         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);</div><div>          URL http://www.quantum-espresso.org",</div><div>     in publications or presentations arising from this work. More details at</div><div>     http://www.quantum-espresso.org/quote</div><div><br></div><div>     Parallel version (MPI), running on    32 processors</div><div>     R & G space division:  proc/nbgrp/npool/nimage =      32</div><div>     Waiting for input...</div><div>     Reading input from standard input</div><div><br></div><div>     Current dimensions of program PWSCF are:</div><div>     Max number of different atomic species (ntypx) = 10</div><div>     Max number of k-points (npk) =  40000</div><div>     Max angular momentum in pseudopotentials (lmaxx) =  3</div><div>               file S.pbe-n-kjpaw_psl.0.1.UPF: wavefunction(s)  3P renormalized</div><div>               file Ge.pbe-kjpaw.UPF: wavefunction(s)  4S renormalized</div><div>               file Co.pbe-nd-rrkjus.UPF: wavefunction(s)  4S 3D renormalized</div><div><br></div><div>     Subspace diagonalization in iterative solution of the eigenvalue problem:</div><div>     a serial algorithm will be used</div><div><br></div><div><br></div><div>     Parallelization info</div><div>     --------------------</div><div>     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW</div><div>     Min         253     136     36                14655     5785     798</div><div>     Max         254     138     38                14658     5787     801</div><div>     Sum        8107    4375   1169               468993   185157   25569</div><div><br></div><div><br></div><div><br></div><div>     bravais-lattice index     =            8</div><div>     lattice parameter (alat)  =      28.3460  a.u.</div><div>     unit-cell volume          =    6624.3219 (a.u.)^3</div><div>     number of atoms/cell      =           16</div><div>     number of atomic types    =            3</div></div><div><div>     number of electrons       =        85.00</div><div>     number of Kohn-Sham states=           52</div><div>     kinetic-energy cutoff     =      35.0000  Ry</div><div>     charge density cutoff     =     260.0000  Ry</div><div>     convergence threshold     =      1.0E-06</div><div>     mixing beta               =       0.5000</div><div>     number of iterations used =            8  local-TF  mixing</div><div>     Exchange-correlation      =  SLA  PW   PBE  PBE ( 1  4  3  4 0 0)</div><div>     nstep                     =          300</div><div><br></div><div><br></div><div>     celldm(1)=  28.346000  celldm(2)=   0.488000  celldm(3)=   0.596000</div><div>     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000</div><div><br></div><div>     crystal axes: (cart. coord. in units of alat)</div><div>               a(1) = (   1.000000   0.000000   0.000000 )</div><div>               a(2) = (   0.000000   0.488000   0.000000 )</div><div>               a(3) = (   0.000000   0.000000   0.596000 )</div><div><br></div><div>     reciprocal axes: (cart. coord. in units 2 pi/alat)</div><div>               b(1) = (  1.000000  0.000000  0.000000 )</div><div>               b(2) = (  0.000000  2.049180  0.000000 )</div><div>               b(3) = (  0.000000  0.000000  1.677852 )</div><div><br></div><div><br></div><div>     PseudoPot. # 1 for  S read from file:</div><div>     /home/wang/pseudo/S.pbe-n-kjpaw_psl.0.1.UPF</div><div>     MD5 check sum: 4d0aa7498161b95d6b336a838fdb998f</div><div>     Pseudo is Projector augmented-wave + core cor, Zval =  6.0</div><div>     Generated using "atomic" code by A. Dal Corso  v.5.0.2 svn rev. 9415</div><div>     Shape of augmentation charge: PSQ</div><div>     Using radial grid of 1151 points,  4 beta functions with:</div><div>                l(1) =   0</div><div>                l(2) =   0</div><div>                l(3) =   1</div><div>                l(4) =   1</div><div>     Q(r) pseudized with 0 coefficients</div></div><div>......................................................................</div><div><div>     one-electron contribution =    -910.78797779 Ry</div><div>     hartree contribution      =     484.08606418 Ry</div><div>     xc contribution           =    -201.99388284 Ry</div><div>     ewald contribution        =     230.35223553 Ry</div><div>     Hubbard energy            =       0.16968641 Ry</div><div>     one-center paw contrib.   =   -1355.33616563 Ry</div><div>     smearing contrib. (-TS)   =      -0.01567672 Ry</div><div><br></div><div>     total magnetization       =     2.31 Bohr mag/cell</div><div>     absolute magnetization    =     2.47 Bohr mag/cell</div><div><br></div><div>     convergence has been achieved in   6 iterations</div><div><br></div><div>     negative rho (up, down):  4.622E-02 4.591E-02</div><div><br></div><div>     Forces acting on atoms (Ry/au):</div><div><br></div><div>     atom    1 type  1   force =     0.00581021   -0.00263947    0.00678497</div><div>     atom    2 type  1   force =     0.01376296   -0.00466940    0.00777758</div><div>     atom    3 type  1   force =     0.00575713   -0.00674279    0.00683317</div><div>     atom    4 type  1   force =    -0.00437967   -0.00465536   -0.01043663</div><div>     atom    5 type  1   force =     0.00368567    0.00172134   -0.01762865</div><div>     atom    6 type  1   force =    -0.00169966   -0.00468960   -0.01594727</div><div>     atom    7 type  1   force =     0.00378233   -0.01108825   -0.01762898</div><div>     atom    8 type  1   force =    -0.00065837   -0.00465208   -0.01813616</div><div>     atom    9 type  2   force =    -0.00660559    0.00148699   -0.01948893</div><div>     atom   10 type  3   force =    -0.03560634    0.07015493    0.14605057</div><div>     atom   11 type  2   force =    -0.00677210   -0.01085718   -0.01934306</div><div>     atom   12 type  2   force =     0.00208925   -0.00465552   -0.03334974</div><div>     atom   13 type  2   force =     0.01013961   -0.00249984   -0.00484230</div><div>     atom   14 type  2   force =    -0.00181755   -0.00466404   -0.00422673</div><div>     atom   15 type  2   force =     0.01018410   -0.00687494   -0.00484869</div><div>     atom   16 type  2   force =     0.00232804   -0.00467477   -0.00156916</div><div><br></div><div>     Total force =     0.179000     Total SCF correction =     0.000018</div><div><br></div><div>     number of scf cycles    = 300</div><div>     number of bfgs steps    =   3</div><div><br></div><div>     energy   old            =   -1753.5257186124 Ry</div><div>     energy   new            =   -1753.5257168682 Ry</div></div><div><div><br></div><div>     CASE: energy  _new > energy  _old</div><div><br></div><div>     new trust radius        =       0.0001053782 bohr</div><div>     new conv_thr            =       0.0000000010 Ry</div><div><br></div><div><br></div><div>ATOMIC_POSITIONS (crystal)</div><div>S        0.250699886   0.124298313   0.000001713</div><div>S        0.108002516   0.374298238   0.250555374</div><div>S        0.250699717   0.624298263   0.000002169</div><div>S        0.107998206   0.874298339   0.250547990</div><div>S        0.250700259   0.124300230   0.499995739</div><div>S        0.107999360   0.374298188   0.750544959</div><div>S        0.250700414   0.624296176   0.499995743</div><div>S        0.107999190   0.874298396   0.750543956</div><div>Ge       0.093999561   0.124301391   0.433993435</div><div>Co       0.264696335   0.374326386   0.183488868</div><div>Ge       0.093999256   0.624295088   0.433993963</div><div>Ge       0.264700474   0.874298053   0.183439519</div><div>Ge       0.094002023   0.124298675   0.933999145</div><div>Ge       0.264699718   0.374298312   0.683448591</div><div>Ge       0.094002176   0.624297726   0.933999251</div><div>Ge       0.264700879   0.874298275   0.683449587</div><div><br></div><div><br></div><div><br></div><div>     Writing output data file pf.save</div><div>     Check: negative starting charge=(component1):   -0.085756</div><div>     NEW-OLD atomic charge density approx. for the potential</div><div><br></div><div>     Check: negative/imaginary core charge=   -0.000012    0.000000</div><div>     Check: negative starting charge=(component1):   -0.085756</div><div><br></div><div>     negative rho (up, down):  4.621E-02 4.592E-02</div><div> Atomic wfc used for LDA+U Projector are NOT orthogonalized</div><div><br></div><div>     Writing output data file pf.save</div></div><div><div>     init_run     :      9.60s CPU      9.67s WALL (       1 calls)</div><div>     electrons    :  98117.65s CPU  98674.16s WALL (     300 calls)</div><div>     update_pot   :   1008.23s CPU   1013.60s WALL (     300 calls)</div><div>     forces       :   1706.26s CPU   1713.28s WALL (     300 calls)</div><div><br></div><div>     Called by init_run:</div><div>     wfcinit      :      4.26s CPU      4.28s WALL (       1 calls)</div><div>     potinit      :      1.57s CPU      1.58s WALL (       1 calls)</div><div><br></div><div>     Called by electrons:</div><div>     c_bands      :  74036.29s CPU  74253.25s WALL (   10422 calls)</div><div>     sum_band     :  10557.00s CPU  10731.74s WALL (   10422 calls)</div><div>     v_of_rho     :   1705.61s CPU   1710.65s WALL (   10560 calls)</div><div>     newd         :   3309.18s CPU   3433.42s WALL (   10560 calls)</div><div>     PAW_pot      :   2338.64s CPU   2343.39s WALL (   10860 calls)</div><div>     mix_rho      :   6238.94s CPU   6264.71s WALL (   10422 calls)</div><div><br></div><div>     Called by c_bands:</div><div>     init_us_2    :    537.40s CPU    542.78s WALL (  220460 calls)</div><div>     cegterg      :  72223.84s CPU  72432.20s WALL (  104220 calls)</div><div><br></div><div>     Called by sum_band:</div><div>     sum_band:bec :     25.77s CPU     26.05s WALL (  107220 calls)</div><div>     addusdens    :   4240.46s CPU   4390.12s WALL (   10422 calls)</div><div><br></div><div>     Called by *egterg:</div><div>     h_psi        :  32733.76s CPU  32855.42s WALL (  453699 calls)</div><div>     s_psi        :   7559.41s CPU   7577.77s WALL (  459709 calls)</div><div>     g_psi        :     93.85s CPU     94.74s WALL (  349469 calls)</div><div>     cdiaghg      :   9307.35s CPU   9325.74s WALL (  449069 calls)</div><div><br></div><div>     Called by h_psi:</div><div>     h_psi:pot    :  29601.13s CPU  29715.55s WALL (  453699 calls)</div><div>     h_psi:calbec :   6257.47s CPU   6272.03s WALL (  453699 calls)</div><div>     vloc_psi     :  16212.46s CPU  16294.87s WALL (  453699 calls)</div><div>     add_vuspsi   :   7126.73s CPU   7143.92s WALL (  453699 calls)</div><div>     vhpsi        :   2930.81s CPU   2937.80s WALL (  453699 calls)</div><div><br></div><div>     General routines</div><div>     calbec       :  11990.58s CPU  12016.30s WALL ( 1571848 calls)</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Thank you very much for your kind help!</div><br><div style="position:relative;zoom:1">--<br><div>Best regards</div><div> </div><div>Q.J.Wang</div><div> </div><div>XiangTan University </div><div style="clear:both"></div></div></div><br><br><span title="neteasefooter"><p> </p></span>