<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>