<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello,<div class=""><br class=""></div><div class="">In regards to the atomic positions, the crystal looks relaxed as the atomic forces are small. </div><div class=""><br class=""></div><div class="">Some remarks:</div><div class="">1. You have run the calculation at the gamma point (k points 1x1x1). It might happen that the structure is not fully relaxed if you take more k-points. As your cell is not that big (12x7x9 angs) you could try to use a larger grid, e.g. 2x4x4 or so, to see what happens.</div><div class="">2. You are relaxing the structure with spin-orbit (SOC) and non-collinear (N-coll) options. As far as I know, SOC forces are not implemented in QE (to be checked). If so, I would first relax without SOC nor N-coll. In any case you will save computer time. subsequently, when the structure is optimised, you can run SCF with SOC and N-coll.</div><div class="">3. The cell is far from relaxed: the pressure is large and negative, so your cell is too big. You may want to relax it with vc-relax. </div><div class=""><br class=""></div><div class="">HTH</div><div class="">Best regards,</div><div class=""><br class=""><div class="">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px;"><span class="" style="font-size: 16pt; font-family: Mistral;">Pascal Boulet</span></div><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px;"><font face="Lucida Handwriting" class=""><span class="" style="font-size: 21px;">—</span></font></div><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px;"><em class="" style="color: rgb(34, 187, 234); font-family: verdana; font-size: 11px; line-height: 22px;">Professor in computational materials chemistry - DEPARTMENT OF CHEMISTRY</em></div><div class=""><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span class="" style="font-size: 11px; font-family: verdana;">University of Aix-Marseille </span><span class="" style="font-size: 11px; font-family: verdana;">- Avenue Escadrille Normandie Niemen - F-13013 Marseille - FRANCE</span></div><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span class="" style="font-size: 11px; font-family: verdana;">Tél: +33(0)4 13 55 18 10 - Fax : +33(0)4 13 55 18 50</span></div><div class="" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; -webkit-text-stroke-width: 0px;"><font face="verdana" class=""><span class="" style="font-size: 12px; line-height: 22px;">Email : </span></font><font color="#22bbea" face="verdana" class=""><span style="font-size: 12px; line-height: 20px;" class=""><a href="mailto:pascal.boulet@univ-amu.fr" class="">pascal.boulet@univ-amu.fr</a></span></font></div><div class="" style="font-size: 14px;"><br class=""></div></div></div></div></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">Le 13 juin 2021 à 12:32, alberto santonocito <<a href="mailto:albesantonocito@gmail.com" class="">albesantonocito@gmail.com</a>> a écrit :</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class="">For clarity I'm attaching the output<br class=""></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno dom 13 giu 2021 alle ore 12:29 alberto santonocito <<a href="mailto:albesantonocito@gmail.com" class="">albesantonocito@gmail.com</a>> ha scritto:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr" class="">Hello qe users,<div class="">I'm doing a geometric optimization (calculation =' relax') of a supercell containing 4 molecules of Fe(CO)5 (C2/c, Z = 4, extracted from "Supplementary Information for Phase Diagram and Transformations of Iron Pentacarbonyl to nm Layered Hematite and Carbon-Oxygen Polymer under Pressure ").<br class="">This is the input:<br class=""></div><div class="">&control<br class=""></div><div class="">   calculation = 'relax',<br class="">   verbosity = 'high',<br class="">   restart_mode = 'from_scratch',<br class="">   nstep=400,<br class="">   etot_conv_thr = 1.d-8,<br class="">   forc_conv_thr = 1.d-5,<br class="">   tstress = .true.,<br class="">   pseudo_dir = './PSEUDO_DIR/',<br class="">   outdir = './TMP_DIR/'<br class=""> /<br class=""> &system<br class="">      ibrav = 14, A=11.807, B=6.821, C=9.367, cosAB=0.0, cosAC=-0.304365584, cosBC=0.0, nat= 44, ntyp= 3,<br class="">      ecutwfc = 80.0<br class="">      ecutrho = 400.0<br class="">      lspinorb=.true.<br class="">      noncolin=.true.<br class=""> / <br class=""> &electrons<br class="">    mixing_beta = 0.2<br class="">    conv_thr =  1.0d-6<br class="">    diagonalization='cg'<br class="">    electron_maxstep=1000<br class=""> /<br class=""> &ions<br class="">    ion_damping = 0.2,<br class="">    ion_velocities = 'zero',<br class="">    ion_nstepe = 10<br class="">    trust_radius_min = 1.d-5<br class="">    /<br class="">ATOMIC_SPECIES<br class=""> Fe  55.84  Fe.pz-sp-hgh.UPF   <br class=""> C   12.    C.pz-hgh.UPF<br class=""> O   15.99  O.pz-hgh.UPF   <br class="">ATOMIC_POSITIONS (crystal)<br class="">Fe      -0.000083575   0.168404809   0.249866056<br class="">C       -0.000047014  -0.084882975   0.249856908<br class="">O       -0.000030546  -0.254341641   0.249868140<br class="">Fe       0.000083575  -0.168404809   0.750133944<br class="">C        0.000047014   0.084882975   0.750143092<br class="">O        0.000030546   0.254341641   0.750131860<br class="">Fe       0.499913517   0.668739566   0.249922673<br class="">C        0.499916311   0.415459437   0.249947904<br class="">O        0.499933946   0.246001224   0.249978986<br class="">Fe       0.500086483   0.331260434   0.750077327<br class="">C        0.500083689   0.584540563   0.750052096<br class="">O        0.500066054   0.753998776   0.750021014<br class="">C        0.076811157   0.298007544   0.407644795<br class="">C        0.125707127   0.162513245   0.189878936<br class="">O        0.128659623   0.383972117   0.514197502<br class="">O        0.208848172   0.154539483   0.149820949<br class="">C       -0.076811157  -0.298007544   0.592355205<br class="">C       -0.125707127  -0.162513245   0.810121064<br class="">O       -0.128659623  -0.383972117   0.485802498<br class="">O       -0.208848172  -0.154539483   0.850179051<br class="">C       -0.077155467   0.297806630   0.092050345<br class="">C       -0.125806766   0.162356585   0.309982931<br class="">O       -0.129206806   0.383621342  -0.014488420<br class="">O       -0.208866572   0.154312703   0.350201423<br class="">C        0.077155467  -0.297806630   0.907949655<br class="">C        0.125806766  -0.162356585   0.690017069<br class="">O        0.129206806  -0.383621342   1.014488420<br class="">O        0.208866572  -0.154312703   0.649798577<br class="">C        0.576872926   0.798305410   0.407693301<br class="">C        0.625677147   0.662669676   0.189890358<br class="">O        0.628796659   0.884237196   0.514228241<br class="">O        0.708789755   0.654527162   0.149783706<br class="">C        0.423127074   0.201694590   0.592306699<br class="">C        0.374322853   0.337330324   0.810109642<br class="">O        0.371203341   0.115762804   0.485771759<br class="">O        0.291210245   0.345472838   0.850216294<br class="">C        0.422938048   0.798247670   0.092123065<br class="">C        0.374177749   0.662699923   0.310017668<br class="">O        0.371003558   0.884146914  -0.014425910<br class="">O        0.291113571   0.654574932   0.350221762<br class="">C        0.577061952   0.201752330   0.907876935<br class="">C        0.625822251   0.337300077   0.689982332<br class="">O        0.628996442   0.115853086   1.014425910<br class="">O        0.708886429   0.345425068   0.649778238<br class=""><br class="">K_POINTS {automatic}<br class="">1 1 1  0 0 0 <br class=""></div><div class=""><br class=""></div><div class="">In the output I notice the following warning message:<br class=""></div><div class="">"Total force =     0.000539     Total SCF correction =     0.000170</div><div class="">     SCF correction compared to forces is large: reduce conv_thr to get better values"</div><div class="">Also not all the symmetry elements of my system are found:<br class=""></div><div class="">"double point group C_i (-1)   <br class="">     there are  4 classes and  2 irreducible representations<br class="">     the character table:<br class=""><br class="">       E     -E    i     -i   <br class="">                              <br class="">G_2+   1.00 -1.00  1.00 -1.00<br class="">G_2-   1.00 -1.00 -1.00  1.00 "</div><div class=""><br class=""></div><div class="">The axis of rotation C2 and a reflection plane orthogonal to the main axis are missing (I don't know how serious the problem can be because the coordinates obtained seem to me consistent with the experimental data). <br class=""></div><div class=""><br class=""></div><div class="">Finally I don't understand if the calculation has finished correctly:<br class=""></div><div class=""><br class=""></div><div class="">.....</div><div class=""><br class=""></div><div class=""> End of self-consistent calculation<br class=""><br class="">          k = 0.0000 0.0000 0.0000 ( 58607 PWs)   bands (ev):<br class=""><br class="">   -93.0586 -93.0586 -93.0499 -93.0499 -93.0496 -93.0496 -93.0495 -93.0495<br class="">   -60.2668 -60.2668 -60.2666 -60.2666 -60.2664 -60.2664 -60.2663 -60.2663<br class="">   -59.7960 -59.7960 -59.7960 -59.7960 -59.7958 -59.7958 -59.7957 -59.7957<br class="">   -59.7492 -59.7492 -59.7491 -59.7491 -59.7487 -59.7487 -59.7485 -59.7485<br class="">   -21.8852 -21.8852 -21.8665 -21.8665 -21.8635 -21.8635 -21.8604 -21.8604<br class="">   -21.8538 -21.8538 -21.8495 -21.8495 -21.8483 -21.8483 -21.8442 -21.8442<br class="">   -21.5947 -21.5947 -21.5916 -21.5916 -21.5793 -21.5793 -21.5792 -21.5792<br class="">   -21.5351 -21.5351 -21.5135 -21.5135 -21.5089 -21.5089 -21.5078 -21.5078<br class="">   -21.5028 -21.5028 -21.4967 -21.4967 -21.4845 -21.4845 -21.4798 -21.4798<br class="">    -9.4894  -9.4894  -9.4712  -9.4712  -9.4699  -9.4699  -9.4678  -9.4678<br class="">    -7.5924  -7.5924  -7.4724  -7.4724  -7.4722  -7.4722  -7.4715  -7.4715<br class="">    -7.4141  -7.4141  -7.4018  -7.4018  -7.3890  -7.3890  -7.3542  -7.3542<br class="">    -7.3512  -7.3512  -7.3299  -7.3299  -7.2802  -7.2802  -7.2168  -7.2168<br class="">    -7.1215  -7.1215  -7.0898  -7.0898  -7.0224  -7.0224  -6.9905  -6.9905<br class="">    -6.0635  -6.0635  -5.7884  -5.7884  -5.7565  -5.7565  -5.7434  -5.7434<br class="">    -5.5964  -5.5964  -5.5476  -5.5476  -5.5187  -5.5187  -5.5159  -5.5159<br class="">    -5.4978  -5.4978  -5.4942  -5.4942  -5.4869  -5.4869  -5.3868  -5.3868<br class="">    -5.3547  -5.3547  -5.3463  -5.3463  -5.2984  -5.2984  -5.2748  -5.2748<br class="">    -5.2549  -5.2549  -5.1991  -5.1991  -5.1533  -5.1533  -5.1514  -5.1514<br class="">    -5.1322  -5.1322  -5.1126  -5.1126  -5.0616  -5.0616  -5.0080  -5.0080<br class="">    -5.0035  -5.0035  -4.9605  -4.9605  -4.9483  -4.9483  -4.9064  -4.9064<br class="">    -4.8033  -4.8033  -4.7701  -4.7701  -4.7300  -4.7300  -4.7200  -4.7200<br class="">    -4.6509  -4.6509  -4.5755  -4.5755  -4.5438  -4.5438  -4.4974  -4.4974<br class="">    -4.3646  -4.3646  -4.3629  -4.3629  -4.3117  -4.3117  -4.2716  -4.2716<br class="">    -4.2712  -4.2712  -4.2576  -4.2576  -4.2547  -4.2547  -4.1542  -4.1542<br class="">    -4.1521  -4.1521  -4.1473  -4.1473  -4.1334  -4.1334  -4.0947  -4.0947<br class="">    -4.0746  -4.0746  -4.0531  -4.0531  -4.0458  -4.0458  -4.0119  -4.0119<br class="">    -3.9876  -3.9876  -3.9804  -3.9804  -3.9614  -3.9614  -3.9458  -3.9458<br class="">    -3.7926  -3.7926  -3.5548  -3.5548  -3.4339  -3.4339  -3.3468  -3.3468<br class="">    -0.9636  -0.9636  -0.9096  -0.9096  -0.8934  -0.8934  -0.8879  -0.8879<br class="">    -0.8574  -0.8574  -0.8399  -0.8399  -0.8197  -0.8197  -0.8051  -0.8051<br class="">     0.3615   0.3615   0.3761   0.3761   0.3855   0.3855   0.3879   0.3879<br class="">     0.3919   0.3919   0.3921   0.3921   0.4314   0.4314   0.4641   0.4641<br class=""><br class="">     occupation numbers <br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class="">     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br class=""><br class="">     highest occupied level (ev):     0.4641<br class=""><br class="">!    total energy              =   -1837.03004023 Ry<br class="">     Harris-Foulkes estimate   =   -1837.03004024 Ry<br class="">     estimated scf accuracy    <       0.00000001 Ry<br class=""><br class="">     The total energy is the sum of the following terms:<br class=""><br class="">     one-electron contribution =   -1723.33773875 Ry<br class="">     hartree contribution      =     915.03055446 Ry<br class="">     xc contribution           =    -301.88818855 Ry<br class="">     ewald contribution        =    -726.83466739 Ry<br class=""><br class="">     convergence has been achieved in  14 iterations<br class=""><br class="">     Forces acting on atoms (Ry/au):<br class=""><br class="">     atom    1 type  1   force =    -0.00008219   -0.00017955    0.00003643<br class="">     atom    2 type  2   force =     0.00001135   -0.00006362    0.00001736<br class="">     atom    3 type  3   force =     0.00000292   -0.00005114    0.00000176<br class="">     atom    4 type  1   force =     0.00008219    0.00017955   -0.00003643<br class="">     atom    5 type  2   force =    -0.00001135    0.00006362   -0.00001736<br class="">     atom    6 type  3   force =    -0.00000292    0.00005114   -0.00000176<br class="">     atom    7 type  1   force =     0.00002769   -0.00021894    0.00003150<br class="">     atom    8 type  2   force =     0.00001118   -0.00006796    0.00000338<br class="">     atom    9 type  3   force =     0.00000303   -0.00008140   -0.00000033<br class="">     atom   10 type  1   force =    -0.00002769    0.00021894   -0.00003150<br class="">     atom   11 type  2   force =    -0.00001118    0.00006796   -0.00000338<br class="">     atom   12 type  3   force =    -0.00000303    0.00008140    0.00000033<br class="">     atom   13 type  2   force =     0.00000165   -0.00004705    0.00002128<br class="">     atom   14 type  2   force =     0.00003679   -0.00003171    0.00001628<br class="">     atom   15 type  3   force =    -0.00002187   -0.00002144    0.00000888<br class="">     atom   16 type  3   force =    -0.00001651    0.00000359    0.00001682<br class="">     atom   17 type  2   force =    -0.00000165    0.00004705   -0.00002128<br class="">     atom   18 type  2   force =    -0.00003679    0.00003171   -0.00001628<br class="">     atom   19 type  3   force =     0.00002187    0.00002144   -0.00000888<br class="">     atom   20 type  3   force =     0.00001651   -0.00000359   -0.00001682<br class="">     atom   21 type  2   force =     0.00001459   -0.00005592    0.00001950<br class="">     atom   22 type  2   force =     0.00002824   -0.00002980    0.00000195<br class="">     atom   23 type  3   force =     0.00001440   -0.00005234    0.00002546<br class="">     atom   24 type  3   force =    -0.00002875    0.00000424   -0.00001141<br class="">     atom   25 type  2   force =    -0.00001459    0.00005592   -0.00001950<br class="">     atom   26 type  2   force =    -0.00002824    0.00002980   -0.00000195<br class="">     atom   27 type  3   force =    -0.00001440    0.00005234   -0.00002546<br class="">     atom   28 type  3   force =     0.00002875   -0.00000424    0.00001141<br class="">     atom   29 type  2   force =    -0.00000421   -0.00005722    0.00001115<br class="">     atom   30 type  2   force =     0.00001879   -0.00003994    0.00002124<br class="">     atom   31 type  3   force =    -0.00003043   -0.00003201   -0.00000672<br class="">     atom   32 type  3   force =     0.00002146    0.00000323    0.00002611<br class="">     atom   33 type  2   force =     0.00000421    0.00005722   -0.00001115<br class="">     atom   34 type  2   force =    -0.00001879    0.00003994   -0.00002124<br class="">     atom   35 type  3   force =     0.00003043    0.00003201    0.00000672<br class="">     atom   36 type  3   force =    -0.00002146   -0.00000323   -0.00002611<br class="">     atom   37 type  2   force =     0.00001656   -0.00006534    0.00001467<br class="">     atom   38 type  2   force =     0.00001128   -0.00004829    0.00000173<br class="">     atom   39 type  3   force =     0.00002179   -0.00004646    0.00002699<br class="">     atom   40 type  3   force =     0.00000256   -0.00000832   -0.00000418<br class="">     atom   41 type  2   force =    -0.00001656    0.00006534   -0.00001467<br class="">     atom   42 type  2   force =    -0.00001128    0.00004829   -0.00000173<br class="">     atom   43 type  3   force =    -0.00002179    0.00004646   -0.00002699<br class="">     atom   44 type  3   force =    -0.00000256    0.00000832    0.00000418<br class="">     The non-local contrib.  to forces<br class="">     atom    1 type  1   force =    -0.00000627   -0.00629403    0.00036022<br class="">     atom    2 type  2   force =     0.00009600   -0.47733304    0.00004743<br class="">     atom    3 type  3   force =    -0.00037693    3.24340453   -0.00060572<br class="">     atom    4 type  1   force =     0.00000627    0.00629403   -0.00036022<br class="">     atom    5 type  2   force =    -0.00009600    0.47733304   -0.00004743<br class="">     atom    6 type  3   force =     0.00037693   -3.24340453    0.00060572<br class="">     atom    7 type  1   force =     0.00010735   -0.00652077    0.00015855<br class="">     atom    8 type  2   force =     0.00000184   -0.47728939    0.00010464<br class="">     atom    9 type  3   force =    -0.00031904    3.24344147   -0.00091287<br class="">     atom   10 type  1   force =    -0.00010735    0.00652077   -0.00015855<br class="">     atom   11 type  2   force =    -0.00000184    0.47728939   -0.00010464<br class="">     atom   12 type  3   force =     0.00031904   -3.24344147    0.00091287<br class="">     atom   13 type  2   force =     0.14002573    0.25641314    0.42533258<br class="">     atom   14 type  2   force =     0.44434952   -0.02078299   -0.14427304<br class="">     atom   15 type  3   force =    -0.84464513   -1.62673322   -2.61122100<br class="">     atom   16 type  3   force =    -3.13583025    0.14830808    1.01950038<br class="">     atom   17 type  2   force =    -0.14002573   -0.25641314   -0.42533258<br class="">     atom   18 type  2   force =    -0.44434952    0.02078299    0.14427304<br class="">     atom   19 type  3   force =     0.84464513    1.62673322    2.61122100<br class="">     atom   20 type  3   force =     3.13583025   -0.14830808   -1.01950038<br class="">     atom   21 type  2   force =    -0.14092611    0.25623476   -0.42547810<br class="">     atom   22 type  2   force =    -0.44381899   -0.02092269    0.14475057<br class="">     atom   23 type  3   force =     0.85165135   -1.62381035    2.61058715<br class="">     atom   24 type  3   force =     3.13453384    0.14972090   -1.02375796<br class="">     atom   25 type  2   force =     0.14092611   -0.25623476    0.42547810<br class="">     atom   26 type  2   force =     0.44381899    0.02092269   -0.14475057<br class="">     atom   27 type  3   force =    -0.85165135    1.62381035   -2.61058715<br class="">     atom   28 type  3   force =    -3.13453384   -0.14972090    1.02375796<br class="">     atom   29 type  2   force =     0.14040527    0.25637516    0.42539604<br class="">     atom   30 type  2   force =     0.44414895   -0.02119419   -0.14443787<br class="">     atom   31 type  3   force =    -0.84738905   -1.62598838   -2.61064817<br class="">     atom   32 type  3   force =    -3.13522929    0.15173440    1.02089948<br class="">     atom   33 type  2   force =    -0.14040527   -0.25637516   -0.42539604<br class="">     atom   34 type  2   force =    -0.44414895    0.02119419    0.14443787<br class="">     atom   35 type  3   force =     0.84738905    1.62598838    2.61064817<br class="">     atom   36 type  3   force =     3.13522929   -0.15173440   -1.02089948<br class="">     atom   37 type  2   force =    -0.14046321    0.25640994   -0.42556324<br class="">     atom   38 type  2   force =    -0.44392624   -0.02112155    0.14468784<br class="">     atom   39 type  3   force =     0.84759207   -1.62527888    2.61076778<br class="">     atom   40 type  3   force =     3.13447742    0.15137810   -1.02334088<br class="">     atom   41 type  2   force =     0.14046321   -0.25640994    0.42556324<br class="">     atom   42 type  2   force =     0.44392624    0.02112155   -0.14468784<br class="">     atom   43 type  3   force =    -0.84759207    1.62527888   -2.61076778<br class="">     atom   44 type  3   force =    -3.13447742   -0.15137810    1.02334088<br class="">     The ionic contribution  to forces<br class="">     atom    1 type  1   force =     0.01325092    0.33950651    0.01142842<br class="">     atom    2 type  2   force =     0.00090117   -3.23575791    0.00456618<br class="">     atom    3 type  3   force =    -0.00145189  -11.93779430    0.00637171<br class="">     atom    4 type  1   force =    -0.01325092   -0.33950651   -0.01142842<br class="">     atom    5 type  2   force =    -0.00090117    3.23575791   -0.00456618<br class="">     atom    6 type  3   force =     0.00145189   11.93779430   -0.00637171<br class="">     atom    7 type  1   force =     0.00905180    0.32313388    0.01318593<br class="">     atom    8 type  2   force =    -0.00006212   -3.23544797    0.00367961<br class="">     atom    9 type  3   force =     0.00200611  -11.92853126    0.00703592<br class="">     atom   10 type  1   force =    -0.00905180   -0.32313388   -0.01318593<br class="">     atom   11 type  2   force =     0.00006212    3.23544797   -0.00367961<br class="">     atom   12 type  3   force =    -0.00200611   11.92853126   -0.00703592<br class="">     atom   13 type  2   force =     1.10403563    1.82607712    2.45975076<br class="">     atom   14 type  2   force =     3.33976300    0.29104673   -0.82075392<br class="">     atom   15 type  3   force =     3.80265280    6.58609941    9.75859670<br class="">     atom   16 type  3   force =    12.51132625   -0.17076722   -4.20439272<br class="">     atom   17 type  2   force =    -1.10403563   -1.82607712   -2.45975076<br class="">     atom   18 type  2   force =    -3.33976300   -0.29104673    0.82075392<br class="">     atom   19 type  3   force =    -3.80265280   -6.58609941   -9.75859670<br class="">     atom   20 type  3   force =   -12.51132625    0.17076722    4.20439272<br class="">     atom   21 type  2   force =    -1.08798238    1.82895036   -2.45502353<br class="">     atom   22 type  2   force =    -3.34212143    0.28287575    0.82855281<br class="">     atom   23 type  3   force =    -3.79634187    6.58990917   -9.74676391<br class="">     atom   24 type  3   force =   -12.51054256   -0.17989811    4.23449052<br class="">     atom   25 type  2   force =     1.08798238   -1.82895036    2.45502353<br class="">     atom   26 type  2   force =     3.34212143   -0.28287575   -0.82855281<br class="">     atom   27 type  3   force =     3.79634187   -6.58990917    9.74676391<br class="">     atom   28 type  3   force =    12.51054256    0.17989811   -4.23449052<br class="">     atom   29 type  2   force =     1.10264271    1.81885527    2.46052242<br class="">     atom   30 type  2   force =     3.34319951    0.28248311   -0.81964494<br class="">     atom   31 type  3   force =     3.81034956    6.57120294    9.75391063<br class="">     atom   32 type  3   force =    12.51471623   -0.19826128   -4.20710784<br class="">     atom   33 type  2   force =    -1.10264271   -1.81885527   -2.46052242<br class="">     atom   34 type  2   force =    -3.34319951   -0.28248311    0.81964494<br class="">     atom   35 type  3   force =    -3.81034956   -6.57120294   -9.75391063<br class="">     atom   36 type  3   force =   -12.51471623    0.19826128    4.20710784<br class="">     atom   37 type  2   force =    -1.09364900    1.82104545   -2.45560655<br class="">     atom   38 type  2   force =    -3.34123945    0.28124257    0.82899191<br class="">     atom   39 type  3   force =    -3.79163971    6.57779011   -9.75287422<br class="">     atom   40 type  3   force =   -12.51289492   -0.20135093    4.23804833<br class="">     atom   41 type  2   force =     1.09364900   -1.82104545    2.45560655<br class="">     atom   42 type  2   force =     3.34123945   -0.28124257   -0.82899191<br class="">     atom   43 type  3   force =     3.79163971   -6.57779011    9.75287422<br class="">     atom   44 type  3   force =    12.51289492    0.20135093   -4.23804833<br class="">     The local contribution  to forces<br class="">     atom    1 type  1   force =    -0.01332684   -0.33339203   -0.01175220<br class="">     atom    2 type  2   force =    -0.00098672    3.71305531   -0.00459235<br class="">     atom    3 type  3   force =     0.00183076    8.69433704   -0.00575599<br class="">     atom    4 type  1   force =     0.01332684    0.33339203    0.01175220<br class="">     atom    5 type  2   force =     0.00098672   -3.71305531    0.00459235<br class="">     atom    6 type  3   force =    -0.00183076   -8.69433704    0.00575599<br class="">     atom    7 type  1   force =    -0.00913145   -0.31683206   -0.01331298<br class="">     atom    8 type  2   force =     0.00007092    3.71268722   -0.00377871<br class="">     atom    9 type  3   force =    -0.00168518    8.68500300   -0.00611949<br class="">     atom   10 type  1   force =     0.00913145    0.31683206    0.01331298<br class="">     atom   11 type  2   force =    -0.00007092   -3.71268722    0.00377871<br class="">     atom   12 type  3   force =     0.00168518   -8.68500300    0.00611949<br class="">     atom   13 type  2   force =    -1.24407083   -2.08255937   -2.88509076<br class="">     atom   14 type  2   force =    -3.78410977   -0.27029645    0.96505564<br class="">     atom   15 type  3   force =    -2.95803108   -4.95939813   -7.14736903<br class="">     atom   16 type  3   force =    -9.37551516    0.02245792    3.18491378<br class="">     atom   17 type  2   force =     1.24407083    2.08255937    2.88509076<br class="">     atom   18 type  2   force =     3.78410977    0.27029645   -0.96505564<br class="">     atom   19 type  3   force =     2.95803108    4.95939813    7.14736903<br class="">     atom   20 type  3   force =     9.37551516   -0.02245792   -3.18491378<br class="">     atom   21 type  2   force =     1.22893585   -2.08526812    2.88056169<br class="">     atom   22 type  2   force =     3.78600757   -0.26198376   -0.97330890<br class="">     atom   23 type  3   force =     2.94470964   -4.96616568    7.13621695<br class="">     atom   24 type  3   force =     9.37599276    0.03017648   -3.21073948<br class="">     atom   25 type  2   force =    -1.22893585    2.08526812   -2.88056169<br class="">     atom   26 type  2   force =    -3.78600757    0.26198376    0.97330890<br class="">     atom   27 type  3   force =    -2.94470964    4.96616568   -7.13621695<br class="">     atom   28 type  3   force =    -9.37599276   -0.03017648    3.21073948<br class="">     atom   29 type  2   force =    -1.24306193   -2.07530904   -2.88593491<br class="">     atom   30 type  2   force =    -3.78735840   -0.26133038    0.96411419<br class="">     atom   31 type  3   force =    -2.96299256   -4.94526053   -7.14327525<br class="">     atom   32 type  3   force =    -9.37946750    0.04652380    3.18623710<br class="">     atom   33 type  2   force =     1.24306193    2.07530904    2.88593491<br class="">     atom   34 type  2   force =     3.78735840    0.26133038   -0.96411419<br class="">     atom   35 type  3   force =     2.96299256    4.94526053    7.14327525<br class="">     atom   36 type  3   force =     9.37946750   -0.04652380   -3.18623710<br class="">     atom   37 type  2   force =     1.23413958   -2.07754255    2.88121601<br class="">     atom   38 type  2   force =     3.78520861   -0.26017090   -0.97368363<br class="">     atom   39 type  3   force =     2.94407248   -4.95256955    7.14214046<br class="">     atom   40 type  3   force =     9.37842628    0.04995750   -3.21470680<br class="">     atom   41 type  2   force =    -1.23413958    2.07754255   -2.88121601<br class="">     atom   42 type  2   force =    -3.78520861    0.26017090    0.97368363<br class="">     atom   43 type  3   force =    -2.94407248    4.95256955   -7.14214046<br class="">     atom   44 type  3   force =    -9.37842628   -0.04995750    3.21470680<br class="">     The core correction contribution to forces<br class="">     atom    1 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    2 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    3 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    4 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    5 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    6 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    7 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    8 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    9 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   10 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   11 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   12 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   13 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   14 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   15 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   16 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   17 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   18 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   19 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   20 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   21 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   22 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   23 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   24 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   25 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   26 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   27 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   28 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   29 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   30 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   31 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   32 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   33 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   34 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   35 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   36 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   37 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   38 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   39 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   40 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   41 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   42 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   43 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   44 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     The Hubbard contrib.    to forces<br class="">     atom    1 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    2 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    3 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    4 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    5 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    6 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    7 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    8 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    9 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   10 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   11 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   12 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   13 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   14 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   15 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   16 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   17 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   18 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   19 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   20 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   21 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   22 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   23 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   24 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   25 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   26 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   27 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   28 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   29 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   30 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   31 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   32 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   33 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   34 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   35 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   36 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   37 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   38 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   39 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   40 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   41 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   42 type  2   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   43 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   44 type  3   force =     0.00000000    0.00000000    0.00000000<br class="">     The SCF correction term to forces<br class="">     atom    1 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    2 type  2   force =     0.00000091   -0.00002797   -0.00000390<br class="">     atom    3 type  3   force =     0.00000098    0.00000160   -0.00000823<br class="">     atom    4 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    5 type  2   force =    -0.00000091    0.00002797    0.00000390<br class="">     atom    6 type  3   force =    -0.00000098   -0.00000160    0.00000823<br class="">     atom    7 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom    8 type  2   force =     0.00000053   -0.00001782   -0.00000216<br class="">     atom    9 type  3   force =     0.00000113    0.00000539   -0.00000389<br class="">     atom   10 type  1   force =     0.00000000    0.00000000    0.00000000<br class="">     atom   11 type  2   force =    -0.00000053    0.00001782    0.00000216<br class="">     atom   12 type  3   force =    -0.00000113   -0.00000539    0.00000389<br class="">     atom   13 type  2   force =     0.00001112    0.00002207    0.00002870<br class="">     atom   14 type  2   force =     0.00003403    0.00000100   -0.00001240<br class="">     atom   15 type  3   force =     0.00000154    0.00001049    0.00000221<br class="">     atom   16 type  3   force =     0.00000265    0.00000481   -0.00000462<br class="">     atom   17 type  2   force =    -0.00001112   -0.00002207   -0.00002870<br class="">     atom   18 type  2   force =    -0.00003403   -0.00000100    0.00001240<br class="">     atom   19 type  3   force =    -0.00000154   -0.00001049   -0.00000221<br class="">     atom   20 type  3   force =    -0.00000265   -0.00000481    0.00000462<br class="">     atom   21 type  2   force =    -0.00001276    0.00002708   -0.00004056<br class="">     atom   22 type  2   force =    -0.00003890    0.00000089    0.00000747<br class="">     atom   23 type  3   force =    -0.00000472    0.00001452   -0.00001474<br class="">     atom   24 type  3   force =    -0.00001279    0.00000496   -0.00000448<br class="">     atom   25 type  2   force =     0.00001276   -0.00002708    0.00004056<br class="">     atom   26 type  2   force =     0.00003890   -0.00000089   -0.00000747<br class="">     atom   27 type  3   force =     0.00000472   -0.00001452    0.00001474<br class="">     atom   28 type  3   force =     0.00001279   -0.00000496    0.00000448<br class="">     atom   29 type  2   force =     0.00000974    0.00002140    0.00002760<br class="">     atom   30 type  2   force =     0.00002873    0.00000153   -0.00001014<br class="">     atom   31 type  3   force =     0.00000162    0.00001396    0.00000607<br class="">     atom   32 type  3   force =     0.00000201    0.00000631   -0.00000263<br class="">     atom   33 type  2   force =    -0.00000974   -0.00002140   -0.00002760<br class="">     atom   34 type  2   force =    -0.00002873   -0.00000153    0.00001014<br class="">     atom   35 type  3   force =    -0.00000162   -0.00001396   -0.00000607<br class="">     atom   36 type  3   force =    -0.00000201   -0.00000631    0.00000263<br class="">     atom   37 type  2   force =    -0.00001080    0.00002181   -0.00003154<br class="">     atom   38 type  2   force =    -0.00003164    0.00000159    0.00000561<br class="">     atom   39 type  3   force =    -0.00000306    0.00001186   -0.00000703<br class="">     atom   40 type  3   force =    -0.00000622    0.00000701   -0.00000483<br class="">     atom   41 type  2   force =     0.00001080   -0.00002181    0.00003154<br class="">     atom   42 type  2   force =     0.00003164   -0.00000159   -0.00000561<br class="">     atom   43 type  3   force =     0.00000306   -0.00001186    0.00000703<br class="">     atom   44 type  3   force =     0.00000622   -0.00000701    0.00000483<br class=""><br class="">     Total force =     0.000539     Total SCF correction =     0.000170<br class="">     SCF correction compared to forces is large: reduce conv_thr to get better values<br class=""><br class=""><br class="">     entering subroutine stress ...<br class=""><br class="">          total   stress  (Ry/bohr**3)                   (kbar)     P=-1034.73<br class="">  -0.00690874   0.00000013  -0.00008470      -1016.31      0.02    -12.46<br class="">   0.00000013  -0.00709583   0.00000022          0.02  -1043.83      0.03<br class="">  -0.00008470   0.00000022  -0.00709737        -12.46      0.03  -1044.06<br class=""><br class="">     kinetic stress (kbar)  21980.66      0.01    174.00<br class="">                                0.01  22476.70     -0.23<br class="">                              174.00     -0.23  22461.37<br class=""><br class="">     local   stress (kbar) -24888.17      4.00   2442.81<br class="">                                4.00 -20335.62     -9.72<br class="">                             2442.81     -9.72 -21900.02<br class=""><br class="">     nonloc. stress (kbar)     84.54      0.02    -55.75<br class="">                                0.02    -69.24      0.06<br class="">                              -55.75      0.06    -68.81<br class=""><br class="">     hartree stress (kbar)  10002.52     -1.62   -912.08<br class="">                               -1.62   8528.48      4.04<br class="">                             -912.08      4.04   9227.00<br class=""><br class="">     exc-cor stress (kbar)  -2851.55      0.00      0.00<br class="">                                0.00  -2851.55      0.00<br class="">                                0.00      0.00  -2851.55<br class=""><br class="">     corecor stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class="">     ewald   stress (kbar)  -5344.31     -2.38  -1661.45<br class="">                               -2.38  -8792.61      5.89<br class="">                            -1661.45      5.89  -7912.05<br class=""><br class="">     hubbard stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class="">     london  stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class="">     XDM     stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class="">     dft-nl  stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class="">     TS-vdW  stress (kbar)      0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class="">                                0.00      0.00      0.00<br class=""><br class=""><br class=""><br class="">     number of scf cycles    = 400<br class="">     number of bfgs steps    = 117<br class=""><br class="">     energy   old            =   -1837.0300404593 Ry<br class="">     energy   new            =   -1837.0300402317 Ry<br class=""><br class="">     CASE: energy  _new > energy  _old<br class=""><br class="">     new trust radius        =       0.0000243237 bohr<br class="">     new conv_thr            =       0.0000000128 Ry<br class=""><br class=""><br class="">ATOMIC_POSITIONS (crystal)<br class="">Fe      -0.000083575   0.168404809   0.249866056<br class="">C       -0.000047014  -0.084882975   0.249856908<br class="">O       -0.000030546  -0.254341641   0.249868140<br class="">Fe       0.000083575  -0.168404809   0.750133944<br class="">C        0.000047014   0.084882975   0.750143092<br class="">O        0.000030546   0.254341641   0.750131860<br class="">Fe       0.499913517   0.668739566   0.249922673<br class="">C        0.499916311   0.415459437   0.249947904<br class="">O        0.499933946   0.246001224   0.249978986<br class="">Fe       0.500086483   0.331260434   0.750077327<br class="">C        0.500083689   0.584540563   0.750052096<br class="">O        0.500066054   0.753998776   0.750021014<br class="">C        0.076811157   0.298007544   0.407644795<br class="">C        0.125707127   0.162513245   0.189878936<br class="">O        0.128659623   0.383972117   0.514197502<br class="">O        0.208848172   0.154539483   0.149820949<br class="">C       -0.076811157  -0.298007544   0.592355205<br class="">C       -0.125707127  -0.162513245   0.810121064<br class="">O       -0.128659623  -0.383972117   0.485802498<br class="">O       -0.208848172  -0.154539483   0.850179051<br class="">C       -0.077155467   0.297806630   0.092050345<br class="">C       -0.125806766   0.162356585   0.309982931<br class="">O       -0.129206806   0.383621342  -0.014488420<br class="">O       -0.208866572   0.154312703   0.350201423<br class="">C        0.077155467  -0.297806630   0.907949655<br class="">C        0.125806766  -0.162356585   0.690017069<br class="">O        0.129206806  -0.383621342   1.014488420<br class="">O        0.208866572  -0.154312703   0.649798577<br class="">C        0.576872926   0.798305410   0.407693301<br class="">C        0.625677147   0.662669676   0.189890358<br class="">O        0.628796659   0.884237196   0.514228241<br class="">O        0.708789755   0.654527162   0.149783706<br class="">C        0.423127074   0.201694590   0.592306699<br class="">C        0.374322853   0.337330324   0.810109642<br class="">O        0.371203341   0.115762804   0.485771759<br class="">O        0.291210245   0.345472838   0.850216294<br class="">C        0.422938048   0.798247670   0.092123065<br class="">C        0.374177749   0.662699923   0.310017668<br class="">O        0.371003558   0.884146914  -0.014425910<br class="">O        0.291113571   0.654574932   0.350221762<br class="">C        0.577061952   0.201752330   0.907876935<br class="">C        0.625822251   0.337300077   0.689982332<br class="">O        0.628996442   0.115853086   1.014425910<br class="">O        0.708886429   0.345425068   0.649778238<br class=""><br class=""><br class=""><br class="">     Writing output data file pwscf.save<br class="">     NEW-OLD atomic charge density approx. for the potential<br class=""><br class="">     Writing output data file pwscf.save<br class=""> <br class="">     init_run     :      8.67s CPU     12.00s WALL (       1 calls)<br class="">     electrons    :  97864.86s CPU 172396.43s WALL (     400 calls)<br class="">     update_pot   :   2440.24s CPU   2869.61s WALL (     400 calls)<br class="">     forces       :   1504.43s CPU   1676.43s WALL (     400 calls)<br class="">     stress       :   2122.07s CPU   2503.75s WALL (     400 calls)<br class=""><br class="">     Called by init_run:<br class="">     wfcinit      :      2.77s CPU      5.24s WALL (       1 calls)<br class="">     potinit      :      3.30s CPU      3.65s WALL (       1 calls)<br class=""><br class="">     Called by electrons:<br class="">     c_bands      :  92159.33s CPU 158168.32s WALL (    5766 calls)<br class="">     sum_band     :   5076.73s CPU  12562.96s WALL (    5766 calls)<br class="">     v_of_rho     :     91.17s CPU    161.96s WALL (    5849 calls)<br class="">     v_h          :     35.55s CPU     82.77s WALL (    5849 calls)<br class="">     v_xc         :     55.60s CPU     79.13s WALL (    5849 calls)<br class="">     mix_rho      :    338.84s CPU    808.76s WALL (    5766 calls)<br class=""><br class="">     Called by c_bands:<br class="">     init_us_2    :     29.41s CPU     44.39s WALL (   11533 calls)<br class="">     ccgdiagg     :  79579.15s CPU 133357.79s WALL (    6098 calls)<br class="">     wfcrot       :  12558.97s CPU  24780.80s WALL (    5381 calls)<br class=""><br class="">     Called by sum_band:<br class=""><br class="">     Called by *cgdiagg:<br class="">     h_psi        :  49287.89s CPU  95230.23s WALL ( 4647479 calls)<br class="">     cdiaghg      :   1226.89s CPU   1336.06s WALL (    5381 calls)<br class=""><br class="">     Called by h_psi:<br class="">     h_psi:vloc   :  40289.45s CPU  77333.70s WALL ( 4647479 calls)<br class="">     h_psi:vnl    :   8861.47s CPU  17711.55s WALL ( 4647479 calls)<br class="">     add_vuspsi   :   3444.16s CPU   6516.57s WALL ( 4647479 calls)<br class="">     h_1psi       :  41477.62s CPU  76557.86s WALL ( 4642098 calls)<br class=""><br class="">     General routines<br class="">     calbec       :  10946.34s CPU  20460.58s WALL ( 9291577 calls)<br class="">     fft          :    565.20s CPU   2108.77s WALL (  127532 calls)<br class="">     ffts         :    131.27s CPU    473.36s WALL (   46460 calls)<br class="">     fftw         :  41280.27s CPU  81952.36s WALL (27295400 calls)<br class="">     interpolate  :    317.92s CPU   1137.22s WALL (   46460 calls)<br class="">     davcio       :      0.05s CPU   4562.48s WALL (     401 calls)<br class=""> <br class="">     Parallel routines<br class="">     fft_scatter  :  32987.94s CPU  66128.36s WALL (27469392 calls)<br class=""> <br class="">     PWSCF        :     1d    5h 8m CPU        2d    3h38m WALL<br class=""><br class=""> <br class="">   This run was terminated on:   6: 1: 5  13Jun2021            <br class=""><br class="">=------------------------------------------------------------------------------=<br class="">   JOB DONE.<br class="">=------------------------------------------------------------------------------=<br class=""></div><div class="">Thanks!</div></div>
</blockquote></div>
_______________________________________________<br class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" class="">www.max-centre.eu</a>)<br class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></blockquote></div><br class=""></div></body></html>