<div dir="ltr">For clarity I'm attaching the output<br></div><br><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">albesantonocito@gmail.com</a>> ha scritto:<br></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">Hello qe users,<div>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>This is the input:<br></div><div>&control<br></div><div> calculation = 'relax',<br> verbosity = 'high',<br> restart_mode = 'from_scratch',<br> nstep=400,<br> etot_conv_thr = 1.d-8,<br> forc_conv_thr = 1.d-5,<br> tstress = .true.,<br> pseudo_dir = './PSEUDO_DIR/',<br> outdir = './TMP_DIR/'<br> /<br> &system<br> 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> ecutwfc = 80.0<br> ecutrho = 400.0<br> lspinorb=.true.<br> noncolin=.true.<br> / <br> &electrons<br> mixing_beta = 0.2<br> conv_thr = 1.0d-6<br> diagonalization='cg'<br> electron_maxstep=1000<br> /<br> &ions<br> ion_damping = 0.2,<br> ion_velocities = 'zero',<br> ion_nstepe = 10<br> trust_radius_min = 1.d-5<br> /<br>ATOMIC_SPECIES<br> Fe 55.84 Fe.pz-sp-hgh.UPF <br> C 12. C.pz-hgh.UPF<br> O 15.99 O.pz-hgh.UPF <br>ATOMIC_POSITIONS (crystal)<br>Fe -0.000083575 0.168404809 0.249866056<br>C -0.000047014 -0.084882975 0.249856908<br>O -0.000030546 -0.254341641 0.249868140<br>Fe 0.000083575 -0.168404809 0.750133944<br>C 0.000047014 0.084882975 0.750143092<br>O 0.000030546 0.254341641 0.750131860<br>Fe 0.499913517 0.668739566 0.249922673<br>C 0.499916311 0.415459437 0.249947904<br>O 0.499933946 0.246001224 0.249978986<br>Fe 0.500086483 0.331260434 0.750077327<br>C 0.500083689 0.584540563 0.750052096<br>O 0.500066054 0.753998776 0.750021014<br>C 0.076811157 0.298007544 0.407644795<br>C 0.125707127 0.162513245 0.189878936<br>O 0.128659623 0.383972117 0.514197502<br>O 0.208848172 0.154539483 0.149820949<br>C -0.076811157 -0.298007544 0.592355205<br>C -0.125707127 -0.162513245 0.810121064<br>O -0.128659623 -0.383972117 0.485802498<br>O -0.208848172 -0.154539483 0.850179051<br>C -0.077155467 0.297806630 0.092050345<br>C -0.125806766 0.162356585 0.309982931<br>O -0.129206806 0.383621342 -0.014488420<br>O -0.208866572 0.154312703 0.350201423<br>C 0.077155467 -0.297806630 0.907949655<br>C 0.125806766 -0.162356585 0.690017069<br>O 0.129206806 -0.383621342 1.014488420<br>O 0.208866572 -0.154312703 0.649798577<br>C 0.576872926 0.798305410 0.407693301<br>C 0.625677147 0.662669676 0.189890358<br>O 0.628796659 0.884237196 0.514228241<br>O 0.708789755 0.654527162 0.149783706<br>C 0.423127074 0.201694590 0.592306699<br>C 0.374322853 0.337330324 0.810109642<br>O 0.371203341 0.115762804 0.485771759<br>O 0.291210245 0.345472838 0.850216294<br>C 0.422938048 0.798247670 0.092123065<br>C 0.374177749 0.662699923 0.310017668<br>O 0.371003558 0.884146914 -0.014425910<br>O 0.291113571 0.654574932 0.350221762<br>C 0.577061952 0.201752330 0.907876935<br>C 0.625822251 0.337300077 0.689982332<br>O 0.628996442 0.115853086 1.014425910<br>O 0.708886429 0.345425068 0.649778238<br><br>K_POINTS {automatic}<br>1 1 1 0 0 0 <br></div><div><br></div><div>In the output I notice the following warning message:<br></div><div>"Total force = 0.000539 Total SCF correction = 0.000170</div><div> SCF correction compared to forces is large: reduce conv_thr to get better values"</div><div>Also not all the symmetry elements of my system are found:<br></div><div>"double point group C_i (-1) <br> there are 4 classes and 2 irreducible representations<br> the character table:<br><br> E -E i -i <br> <br>G_2+ 1.00 -1.00 1.00 -1.00<br>G_2- 1.00 -1.00 -1.00 1.00 "</div><div><br></div><div>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></div><div><br></div><div>Finally I don't understand if the calculation has finished correctly:<br></div><div><br></div><div>.....</div><div><br></div><div> End of self-consistent calculation<br><br> k = 0.0000 0.0000 0.0000 ( 58607 PWs) bands (ev):<br><br> -93.0586 -93.0586 -93.0499 -93.0499 -93.0496 -93.0496 -93.0495 -93.0495<br> -60.2668 -60.2668 -60.2666 -60.2666 -60.2664 -60.2664 -60.2663 -60.2663<br> -59.7960 -59.7960 -59.7960 -59.7960 -59.7958 -59.7958 -59.7957 -59.7957<br> -59.7492 -59.7492 -59.7491 -59.7491 -59.7487 -59.7487 -59.7485 -59.7485<br> -21.8852 -21.8852 -21.8665 -21.8665 -21.8635 -21.8635 -21.8604 -21.8604<br> -21.8538 -21.8538 -21.8495 -21.8495 -21.8483 -21.8483 -21.8442 -21.8442<br> -21.5947 -21.5947 -21.5916 -21.5916 -21.5793 -21.5793 -21.5792 -21.5792<br> -21.5351 -21.5351 -21.5135 -21.5135 -21.5089 -21.5089 -21.5078 -21.5078<br> -21.5028 -21.5028 -21.4967 -21.4967 -21.4845 -21.4845 -21.4798 -21.4798<br> -9.4894 -9.4894 -9.4712 -9.4712 -9.4699 -9.4699 -9.4678 -9.4678<br> -7.5924 -7.5924 -7.4724 -7.4724 -7.4722 -7.4722 -7.4715 -7.4715<br> -7.4141 -7.4141 -7.4018 -7.4018 -7.3890 -7.3890 -7.3542 -7.3542<br> -7.3512 -7.3512 -7.3299 -7.3299 -7.2802 -7.2802 -7.2168 -7.2168<br> -7.1215 -7.1215 -7.0898 -7.0898 -7.0224 -7.0224 -6.9905 -6.9905<br> -6.0635 -6.0635 -5.7884 -5.7884 -5.7565 -5.7565 -5.7434 -5.7434<br> -5.5964 -5.5964 -5.5476 -5.5476 -5.5187 -5.5187 -5.5159 -5.5159<br> -5.4978 -5.4978 -5.4942 -5.4942 -5.4869 -5.4869 -5.3868 -5.3868<br> -5.3547 -5.3547 -5.3463 -5.3463 -5.2984 -5.2984 -5.2748 -5.2748<br> -5.2549 -5.2549 -5.1991 -5.1991 -5.1533 -5.1533 -5.1514 -5.1514<br> -5.1322 -5.1322 -5.1126 -5.1126 -5.0616 -5.0616 -5.0080 -5.0080<br> -5.0035 -5.0035 -4.9605 -4.9605 -4.9483 -4.9483 -4.9064 -4.9064<br> -4.8033 -4.8033 -4.7701 -4.7701 -4.7300 -4.7300 -4.7200 -4.7200<br> -4.6509 -4.6509 -4.5755 -4.5755 -4.5438 -4.5438 -4.4974 -4.4974<br> -4.3646 -4.3646 -4.3629 -4.3629 -4.3117 -4.3117 -4.2716 -4.2716<br> -4.2712 -4.2712 -4.2576 -4.2576 -4.2547 -4.2547 -4.1542 -4.1542<br> -4.1521 -4.1521 -4.1473 -4.1473 -4.1334 -4.1334 -4.0947 -4.0947<br> -4.0746 -4.0746 -4.0531 -4.0531 -4.0458 -4.0458 -4.0119 -4.0119<br> -3.9876 -3.9876 -3.9804 -3.9804 -3.9614 -3.9614 -3.9458 -3.9458<br> -3.7926 -3.7926 -3.5548 -3.5548 -3.4339 -3.4339 -3.3468 -3.3468<br> -0.9636 -0.9636 -0.9096 -0.9096 -0.8934 -0.8934 -0.8879 -0.8879<br> -0.8574 -0.8574 -0.8399 -0.8399 -0.8197 -0.8197 -0.8051 -0.8051<br> 0.3615 0.3615 0.3761 0.3761 0.3855 0.3855 0.3879 0.3879<br> 0.3919 0.3919 0.3921 0.3921 0.4314 0.4314 0.4641 0.4641<br><br> occupation numbers <br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br> 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000<br><br> highest occupied level (ev): 0.4641<br><br>! total energy = -1837.03004023 Ry<br> Harris-Foulkes estimate = -1837.03004024 Ry<br> estimated scf accuracy < 0.00000001 Ry<br><br> The total energy is the sum of the following terms:<br><br> one-electron contribution = -1723.33773875 Ry<br> hartree contribution = 915.03055446 Ry<br> xc contribution = -301.88818855 Ry<br> ewald contribution = -726.83466739 Ry<br><br> convergence has been achieved in 14 iterations<br><br> Forces acting on atoms (Ry/au):<br><br> atom 1 type 1 force = -0.00008219 -0.00017955 0.00003643<br> atom 2 type 2 force = 0.00001135 -0.00006362 0.00001736<br> atom 3 type 3 force = 0.00000292 -0.00005114 0.00000176<br> atom 4 type 1 force = 0.00008219 0.00017955 -0.00003643<br> atom 5 type 2 force = -0.00001135 0.00006362 -0.00001736<br> atom 6 type 3 force = -0.00000292 0.00005114 -0.00000176<br> atom 7 type 1 force = 0.00002769 -0.00021894 0.00003150<br> atom 8 type 2 force = 0.00001118 -0.00006796 0.00000338<br> atom 9 type 3 force = 0.00000303 -0.00008140 -0.00000033<br> atom 10 type 1 force = -0.00002769 0.00021894 -0.00003150<br> atom 11 type 2 force = -0.00001118 0.00006796 -0.00000338<br> atom 12 type 3 force = -0.00000303 0.00008140 0.00000033<br> atom 13 type 2 force = 0.00000165 -0.00004705 0.00002128<br> atom 14 type 2 force = 0.00003679 -0.00003171 0.00001628<br> atom 15 type 3 force = -0.00002187 -0.00002144 0.00000888<br> atom 16 type 3 force = -0.00001651 0.00000359 0.00001682<br> atom 17 type 2 force = -0.00000165 0.00004705 -0.00002128<br> atom 18 type 2 force = -0.00003679 0.00003171 -0.00001628<br> atom 19 type 3 force = 0.00002187 0.00002144 -0.00000888<br> atom 20 type 3 force = 0.00001651 -0.00000359 -0.00001682<br> atom 21 type 2 force = 0.00001459 -0.00005592 0.00001950<br> atom 22 type 2 force = 0.00002824 -0.00002980 0.00000195<br> atom 23 type 3 force = 0.00001440 -0.00005234 0.00002546<br> atom 24 type 3 force = -0.00002875 0.00000424 -0.00001141<br> atom 25 type 2 force = -0.00001459 0.00005592 -0.00001950<br> atom 26 type 2 force = -0.00002824 0.00002980 -0.00000195<br> atom 27 type 3 force = -0.00001440 0.00005234 -0.00002546<br> atom 28 type 3 force = 0.00002875 -0.00000424 0.00001141<br> atom 29 type 2 force = -0.00000421 -0.00005722 0.00001115<br> atom 30 type 2 force = 0.00001879 -0.00003994 0.00002124<br> atom 31 type 3 force = -0.00003043 -0.00003201 -0.00000672<br> atom 32 type 3 force = 0.00002146 0.00000323 0.00002611<br> atom 33 type 2 force = 0.00000421 0.00005722 -0.00001115<br> atom 34 type 2 force = -0.00001879 0.00003994 -0.00002124<br> atom 35 type 3 force = 0.00003043 0.00003201 0.00000672<br> atom 36 type 3 force = -0.00002146 -0.00000323 -0.00002611<br> atom 37 type 2 force = 0.00001656 -0.00006534 0.00001467<br> atom 38 type 2 force = 0.00001128 -0.00004829 0.00000173<br> atom 39 type 3 force = 0.00002179 -0.00004646 0.00002699<br> atom 40 type 3 force = 0.00000256 -0.00000832 -0.00000418<br> atom 41 type 2 force = -0.00001656 0.00006534 -0.00001467<br> atom 42 type 2 force = -0.00001128 0.00004829 -0.00000173<br> atom 43 type 3 force = -0.00002179 0.00004646 -0.00002699<br> atom 44 type 3 force = -0.00000256 0.00000832 0.00000418<br> The non-local contrib. to forces<br> atom 1 type 1 force = -0.00000627 -0.00629403 0.00036022<br> atom 2 type 2 force = 0.00009600 -0.47733304 0.00004743<br> atom 3 type 3 force = -0.00037693 3.24340453 -0.00060572<br> atom 4 type 1 force = 0.00000627 0.00629403 -0.00036022<br> atom 5 type 2 force = -0.00009600 0.47733304 -0.00004743<br> atom 6 type 3 force = 0.00037693 -3.24340453 0.00060572<br> atom 7 type 1 force = 0.00010735 -0.00652077 0.00015855<br> atom 8 type 2 force = 0.00000184 -0.47728939 0.00010464<br> atom 9 type 3 force = -0.00031904 3.24344147 -0.00091287<br> atom 10 type 1 force = -0.00010735 0.00652077 -0.00015855<br> atom 11 type 2 force = -0.00000184 0.47728939 -0.00010464<br> atom 12 type 3 force = 0.00031904 -3.24344147 0.00091287<br> atom 13 type 2 force = 0.14002573 0.25641314 0.42533258<br> atom 14 type 2 force = 0.44434952 -0.02078299 -0.14427304<br> atom 15 type 3 force = -0.84464513 -1.62673322 -2.61122100<br> atom 16 type 3 force = -3.13583025 0.14830808 1.01950038<br> atom 17 type 2 force = -0.14002573 -0.25641314 -0.42533258<br> atom 18 type 2 force = -0.44434952 0.02078299 0.14427304<br> atom 19 type 3 force = 0.84464513 1.62673322 2.61122100<br> atom 20 type 3 force = 3.13583025 -0.14830808 -1.01950038<br> atom 21 type 2 force = -0.14092611 0.25623476 -0.42547810<br> atom 22 type 2 force = -0.44381899 -0.02092269 0.14475057<br> atom 23 type 3 force = 0.85165135 -1.62381035 2.61058715<br> atom 24 type 3 force = 3.13453384 0.14972090 -1.02375796<br> atom 25 type 2 force = 0.14092611 -0.25623476 0.42547810<br> atom 26 type 2 force = 0.44381899 0.02092269 -0.14475057<br> atom 27 type 3 force = -0.85165135 1.62381035 -2.61058715<br> atom 28 type 3 force = -3.13453384 -0.14972090 1.02375796<br> atom 29 type 2 force = 0.14040527 0.25637516 0.42539604<br> atom 30 type 2 force = 0.44414895 -0.02119419 -0.14443787<br> atom 31 type 3 force = -0.84738905 -1.62598838 -2.61064817<br> atom 32 type 3 force = -3.13522929 0.15173440 1.02089948<br> atom 33 type 2 force = -0.14040527 -0.25637516 -0.42539604<br> atom 34 type 2 force = -0.44414895 0.02119419 0.14443787<br> atom 35 type 3 force = 0.84738905 1.62598838 2.61064817<br> atom 36 type 3 force = 3.13522929 -0.15173440 -1.02089948<br> atom 37 type 2 force = -0.14046321 0.25640994 -0.42556324<br> atom 38 type 2 force = -0.44392624 -0.02112155 0.14468784<br> atom 39 type 3 force = 0.84759207 -1.62527888 2.61076778<br> atom 40 type 3 force = 3.13447742 0.15137810 -1.02334088<br> atom 41 type 2 force = 0.14046321 -0.25640994 0.42556324<br> atom 42 type 2 force = 0.44392624 0.02112155 -0.14468784<br> atom 43 type 3 force = -0.84759207 1.62527888 -2.61076778<br> atom 44 type 3 force = -3.13447742 -0.15137810 1.02334088<br> The ionic contribution to forces<br> atom 1 type 1 force = 0.01325092 0.33950651 0.01142842<br> atom 2 type 2 force = 0.00090117 -3.23575791 0.00456618<br> atom 3 type 3 force = -0.00145189 -11.93779430 0.00637171<br> atom 4 type 1 force = -0.01325092 -0.33950651 -0.01142842<br> atom 5 type 2 force = -0.00090117 3.23575791 -0.00456618<br> atom 6 type 3 force = 0.00145189 11.93779430 -0.00637171<br> atom 7 type 1 force = 0.00905180 0.32313388 0.01318593<br> atom 8 type 2 force = -0.00006212 -3.23544797 0.00367961<br> atom 9 type 3 force = 0.00200611 -11.92853126 0.00703592<br> atom 10 type 1 force = -0.00905180 -0.32313388 -0.01318593<br> atom 11 type 2 force = 0.00006212 3.23544797 -0.00367961<br> atom 12 type 3 force = -0.00200611 11.92853126 -0.00703592<br> atom 13 type 2 force = 1.10403563 1.82607712 2.45975076<br> atom 14 type 2 force = 3.33976300 0.29104673 -0.82075392<br> atom 15 type 3 force = 3.80265280 6.58609941 9.75859670<br> atom 16 type 3 force = 12.51132625 -0.17076722 -4.20439272<br> atom 17 type 2 force = -1.10403563 -1.82607712 -2.45975076<br> atom 18 type 2 force = -3.33976300 -0.29104673 0.82075392<br> atom 19 type 3 force = -3.80265280 -6.58609941 -9.75859670<br> atom 20 type 3 force = -12.51132625 0.17076722 4.20439272<br> atom 21 type 2 force = -1.08798238 1.82895036 -2.45502353<br> atom 22 type 2 force = -3.34212143 0.28287575 0.82855281<br> atom 23 type 3 force = -3.79634187 6.58990917 -9.74676391<br> atom 24 type 3 force = -12.51054256 -0.17989811 4.23449052<br> atom 25 type 2 force = 1.08798238 -1.82895036 2.45502353<br> atom 26 type 2 force = 3.34212143 -0.28287575 -0.82855281<br> atom 27 type 3 force = 3.79634187 -6.58990917 9.74676391<br> atom 28 type 3 force = 12.51054256 0.17989811 -4.23449052<br> atom 29 type 2 force = 1.10264271 1.81885527 2.46052242<br> atom 30 type 2 force = 3.34319951 0.28248311 -0.81964494<br> atom 31 type 3 force = 3.81034956 6.57120294 9.75391063<br> atom 32 type 3 force = 12.51471623 -0.19826128 -4.20710784<br> atom 33 type 2 force = -1.10264271 -1.81885527 -2.46052242<br> atom 34 type 2 force = -3.34319951 -0.28248311 0.81964494<br> atom 35 type 3 force = -3.81034956 -6.57120294 -9.75391063<br> atom 36 type 3 force = -12.51471623 0.19826128 4.20710784<br> atom 37 type 2 force = -1.09364900 1.82104545 -2.45560655<br> atom 38 type 2 force = -3.34123945 0.28124257 0.82899191<br> atom 39 type 3 force = -3.79163971 6.57779011 -9.75287422<br> atom 40 type 3 force = -12.51289492 -0.20135093 4.23804833<br> atom 41 type 2 force = 1.09364900 -1.82104545 2.45560655<br> atom 42 type 2 force = 3.34123945 -0.28124257 -0.82899191<br> atom 43 type 3 force = 3.79163971 -6.57779011 9.75287422<br> atom 44 type 3 force = 12.51289492 0.20135093 -4.23804833<br> The local contribution to forces<br> atom 1 type 1 force = -0.01332684 -0.33339203 -0.01175220<br> atom 2 type 2 force = -0.00098672 3.71305531 -0.00459235<br> atom 3 type 3 force = 0.00183076 8.69433704 -0.00575599<br> atom 4 type 1 force = 0.01332684 0.33339203 0.01175220<br> atom 5 type 2 force = 0.00098672 -3.71305531 0.00459235<br> atom 6 type 3 force = -0.00183076 -8.69433704 0.00575599<br> atom 7 type 1 force = -0.00913145 -0.31683206 -0.01331298<br> atom 8 type 2 force = 0.00007092 3.71268722 -0.00377871<br> atom 9 type 3 force = -0.00168518 8.68500300 -0.00611949<br> atom 10 type 1 force = 0.00913145 0.31683206 0.01331298<br> atom 11 type 2 force = -0.00007092 -3.71268722 0.00377871<br> atom 12 type 3 force = 0.00168518 -8.68500300 0.00611949<br> atom 13 type 2 force = -1.24407083 -2.08255937 -2.88509076<br> atom 14 type 2 force = -3.78410977 -0.27029645 0.96505564<br> atom 15 type 3 force = -2.95803108 -4.95939813 -7.14736903<br> atom 16 type 3 force = -9.37551516 0.02245792 3.18491378<br> atom 17 type 2 force = 1.24407083 2.08255937 2.88509076<br> atom 18 type 2 force = 3.78410977 0.27029645 -0.96505564<br> atom 19 type 3 force = 2.95803108 4.95939813 7.14736903<br> atom 20 type 3 force = 9.37551516 -0.02245792 -3.18491378<br> atom 21 type 2 force = 1.22893585 -2.08526812 2.88056169<br> atom 22 type 2 force = 3.78600757 -0.26198376 -0.97330890<br> atom 23 type 3 force = 2.94470964 -4.96616568 7.13621695<br> atom 24 type 3 force = 9.37599276 0.03017648 -3.21073948<br> atom 25 type 2 force = -1.22893585 2.08526812 -2.88056169<br> atom 26 type 2 force = -3.78600757 0.26198376 0.97330890<br> atom 27 type 3 force = -2.94470964 4.96616568 -7.13621695<br> atom 28 type 3 force = -9.37599276 -0.03017648 3.21073948<br> atom 29 type 2 force = -1.24306193 -2.07530904 -2.88593491<br> atom 30 type 2 force = -3.78735840 -0.26133038 0.96411419<br> atom 31 type 3 force = -2.96299256 -4.94526053 -7.14327525<br> atom 32 type 3 force = -9.37946750 0.04652380 3.18623710<br> atom 33 type 2 force = 1.24306193 2.07530904 2.88593491<br> atom 34 type 2 force = 3.78735840 0.26133038 -0.96411419<br> atom 35 type 3 force = 2.96299256 4.94526053 7.14327525<br> atom 36 type 3 force = 9.37946750 -0.04652380 -3.18623710<br> atom 37 type 2 force = 1.23413958 -2.07754255 2.88121601<br> atom 38 type 2 force = 3.78520861 -0.26017090 -0.97368363<br> atom 39 type 3 force = 2.94407248 -4.95256955 7.14214046<br> atom 40 type 3 force = 9.37842628 0.04995750 -3.21470680<br> atom 41 type 2 force = -1.23413958 2.07754255 -2.88121601<br> atom 42 type 2 force = -3.78520861 0.26017090 0.97368363<br> atom 43 type 3 force = -2.94407248 4.95256955 -7.14214046<br> atom 44 type 3 force = -9.37842628 -0.04995750 3.21470680<br> The core correction contribution to forces<br> atom 1 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 2 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 3 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 4 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 5 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 6 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 7 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 8 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 9 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 10 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 11 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 12 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 13 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 14 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 15 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 16 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 17 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 18 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 19 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 20 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 21 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 22 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 23 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 24 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 25 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 26 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 27 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 28 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 29 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 30 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 31 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 32 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 33 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 34 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 35 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 36 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 37 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 38 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 39 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 40 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 41 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 42 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 43 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 44 type 3 force = 0.00000000 0.00000000 0.00000000<br> The Hubbard contrib. to forces<br> atom 1 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 2 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 3 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 4 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 5 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 6 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 7 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 8 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 9 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 10 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 11 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 12 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 13 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 14 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 15 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 16 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 17 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 18 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 19 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 20 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 21 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 22 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 23 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 24 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 25 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 26 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 27 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 28 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 29 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 30 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 31 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 32 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 33 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 34 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 35 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 36 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 37 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 38 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 39 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 40 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 41 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 42 type 2 force = 0.00000000 0.00000000 0.00000000<br> atom 43 type 3 force = 0.00000000 0.00000000 0.00000000<br> atom 44 type 3 force = 0.00000000 0.00000000 0.00000000<br> The SCF correction term to forces<br> atom 1 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 2 type 2 force = 0.00000091 -0.00002797 -0.00000390<br> atom 3 type 3 force = 0.00000098 0.00000160 -0.00000823<br> atom 4 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 5 type 2 force = -0.00000091 0.00002797 0.00000390<br> atom 6 type 3 force = -0.00000098 -0.00000160 0.00000823<br> atom 7 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 8 type 2 force = 0.00000053 -0.00001782 -0.00000216<br> atom 9 type 3 force = 0.00000113 0.00000539 -0.00000389<br> atom 10 type 1 force = 0.00000000 0.00000000 0.00000000<br> atom 11 type 2 force = -0.00000053 0.00001782 0.00000216<br> atom 12 type 3 force = -0.00000113 -0.00000539 0.00000389<br> atom 13 type 2 force = 0.00001112 0.00002207 0.00002870<br> atom 14 type 2 force = 0.00003403 0.00000100 -0.00001240<br> atom 15 type 3 force = 0.00000154 0.00001049 0.00000221<br> atom 16 type 3 force = 0.00000265 0.00000481 -0.00000462<br> atom 17 type 2 force = -0.00001112 -0.00002207 -0.00002870<br> atom 18 type 2 force = -0.00003403 -0.00000100 0.00001240<br> atom 19 type 3 force = -0.00000154 -0.00001049 -0.00000221<br> atom 20 type 3 force = -0.00000265 -0.00000481 0.00000462<br> atom 21 type 2 force = -0.00001276 0.00002708 -0.00004056<br> atom 22 type 2 force = -0.00003890 0.00000089 0.00000747<br> atom 23 type 3 force = -0.00000472 0.00001452 -0.00001474<br> atom 24 type 3 force = -0.00001279 0.00000496 -0.00000448<br> atom 25 type 2 force = 0.00001276 -0.00002708 0.00004056<br> atom 26 type 2 force = 0.00003890 -0.00000089 -0.00000747<br> atom 27 type 3 force = 0.00000472 -0.00001452 0.00001474<br> atom 28 type 3 force = 0.00001279 -0.00000496 0.00000448<br> atom 29 type 2 force = 0.00000974 0.00002140 0.00002760<br> atom 30 type 2 force = 0.00002873 0.00000153 -0.00001014<br> atom 31 type 3 force = 0.00000162 0.00001396 0.00000607<br> atom 32 type 3 force = 0.00000201 0.00000631 -0.00000263<br> atom 33 type 2 force = -0.00000974 -0.00002140 -0.00002760<br> atom 34 type 2 force = -0.00002873 -0.00000153 0.00001014<br> atom 35 type 3 force = -0.00000162 -0.00001396 -0.00000607<br> atom 36 type 3 force = -0.00000201 -0.00000631 0.00000263<br> atom 37 type 2 force = -0.00001080 0.00002181 -0.00003154<br> atom 38 type 2 force = -0.00003164 0.00000159 0.00000561<br> atom 39 type 3 force = -0.00000306 0.00001186 -0.00000703<br> atom 40 type 3 force = -0.00000622 0.00000701 -0.00000483<br> atom 41 type 2 force = 0.00001080 -0.00002181 0.00003154<br> atom 42 type 2 force = 0.00003164 -0.00000159 -0.00000561<br> atom 43 type 3 force = 0.00000306 -0.00001186 0.00000703<br> atom 44 type 3 force = 0.00000622 -0.00000701 0.00000483<br><br> Total force = 0.000539 Total SCF correction = 0.000170<br> SCF correction compared to forces is large: reduce conv_thr to get better values<br><br><br> entering subroutine stress ...<br><br> total stress (Ry/bohr**3) (kbar) P=-1034.73<br> -0.00690874 0.00000013 -0.00008470 -1016.31 0.02 -12.46<br> 0.00000013 -0.00709583 0.00000022 0.02 -1043.83 0.03<br> -0.00008470 0.00000022 -0.00709737 -12.46 0.03 -1044.06<br><br> kinetic stress (kbar) 21980.66 0.01 174.00<br> 0.01 22476.70 -0.23<br> 174.00 -0.23 22461.37<br><br> local stress (kbar) -24888.17 4.00 2442.81<br> 4.00 -20335.62 -9.72<br> 2442.81 -9.72 -21900.02<br><br> nonloc. stress (kbar) 84.54 0.02 -55.75<br> 0.02 -69.24 0.06<br> -55.75 0.06 -68.81<br><br> hartree stress (kbar) 10002.52 -1.62 -912.08<br> -1.62 8528.48 4.04<br> -912.08 4.04 9227.00<br><br> exc-cor stress (kbar) -2851.55 0.00 0.00<br> 0.00 -2851.55 0.00<br> 0.00 0.00 -2851.55<br><br> corecor stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br> ewald stress (kbar) -5344.31 -2.38 -1661.45<br> -2.38 -8792.61 5.89<br> -1661.45 5.89 -7912.05<br><br> hubbard stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br> london stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br> XDM stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br> dft-nl stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br> TS-vdW stress (kbar) 0.00 0.00 0.00<br> 0.00 0.00 0.00<br> 0.00 0.00 0.00<br><br><br><br> number of scf cycles = 400<br> number of bfgs steps = 117<br><br> energy old = -1837.0300404593 Ry<br> energy new = -1837.0300402317 Ry<br><br> CASE: energy _new > energy _old<br><br> new trust radius = 0.0000243237 bohr<br> new conv_thr = 0.0000000128 Ry<br><br><br>ATOMIC_POSITIONS (crystal)<br>Fe -0.000083575 0.168404809 0.249866056<br>C -0.000047014 -0.084882975 0.249856908<br>O -0.000030546 -0.254341641 0.249868140<br>Fe 0.000083575 -0.168404809 0.750133944<br>C 0.000047014 0.084882975 0.750143092<br>O 0.000030546 0.254341641 0.750131860<br>Fe 0.499913517 0.668739566 0.249922673<br>C 0.499916311 0.415459437 0.249947904<br>O 0.499933946 0.246001224 0.249978986<br>Fe 0.500086483 0.331260434 0.750077327<br>C 0.500083689 0.584540563 0.750052096<br>O 0.500066054 0.753998776 0.750021014<br>C 0.076811157 0.298007544 0.407644795<br>C 0.125707127 0.162513245 0.189878936<br>O 0.128659623 0.383972117 0.514197502<br>O 0.208848172 0.154539483 0.149820949<br>C -0.076811157 -0.298007544 0.592355205<br>C -0.125707127 -0.162513245 0.810121064<br>O -0.128659623 -0.383972117 0.485802498<br>O -0.208848172 -0.154539483 0.850179051<br>C -0.077155467 0.297806630 0.092050345<br>C -0.125806766 0.162356585 0.309982931<br>O -0.129206806 0.383621342 -0.014488420<br>O -0.208866572 0.154312703 0.350201423<br>C 0.077155467 -0.297806630 0.907949655<br>C 0.125806766 -0.162356585 0.690017069<br>O 0.129206806 -0.383621342 1.014488420<br>O 0.208866572 -0.154312703 0.649798577<br>C 0.576872926 0.798305410 0.407693301<br>C 0.625677147 0.662669676 0.189890358<br>O 0.628796659 0.884237196 0.514228241<br>O 0.708789755 0.654527162 0.149783706<br>C 0.423127074 0.201694590 0.592306699<br>C 0.374322853 0.337330324 0.810109642<br>O 0.371203341 0.115762804 0.485771759<br>O 0.291210245 0.345472838 0.850216294<br>C 0.422938048 0.798247670 0.092123065<br>C 0.374177749 0.662699923 0.310017668<br>O 0.371003558 0.884146914 -0.014425910<br>O 0.291113571 0.654574932 0.350221762<br>C 0.577061952 0.201752330 0.907876935<br>C 0.625822251 0.337300077 0.689982332<br>O 0.628996442 0.115853086 1.014425910<br>O 0.708886429 0.345425068 0.649778238<br><br><br><br> Writing output data file pwscf.save<br> NEW-OLD atomic charge density approx. for the potential<br><br> Writing output data file pwscf.save<br> <br> init_run : 8.67s CPU 12.00s WALL ( 1 calls)<br> electrons : 97864.86s CPU 172396.43s WALL ( 400 calls)<br> update_pot : 2440.24s CPU 2869.61s WALL ( 400 calls)<br> forces : 1504.43s CPU 1676.43s WALL ( 400 calls)<br> stress : 2122.07s CPU 2503.75s WALL ( 400 calls)<br><br> Called by init_run:<br> wfcinit : 2.77s CPU 5.24s WALL ( 1 calls)<br> potinit : 3.30s CPU 3.65s WALL ( 1 calls)<br><br> Called by electrons:<br> c_bands : 92159.33s CPU 158168.32s WALL ( 5766 calls)<br> sum_band : 5076.73s CPU 12562.96s WALL ( 5766 calls)<br> v_of_rho : 91.17s CPU 161.96s WALL ( 5849 calls)<br> v_h : 35.55s CPU 82.77s WALL ( 5849 calls)<br> v_xc : 55.60s CPU 79.13s WALL ( 5849 calls)<br> mix_rho : 338.84s CPU 808.76s WALL ( 5766 calls)<br><br> Called by c_bands:<br> init_us_2 : 29.41s CPU 44.39s WALL ( 11533 calls)<br> ccgdiagg : 79579.15s CPU 133357.79s WALL ( 6098 calls)<br> wfcrot : 12558.97s CPU 24780.80s WALL ( 5381 calls)<br><br> Called by sum_band:<br><br> Called by *cgdiagg:<br> h_psi : 49287.89s CPU 95230.23s WALL ( 4647479 calls)<br> cdiaghg : 1226.89s CPU 1336.06s WALL ( 5381 calls)<br><br> Called by h_psi:<br> h_psi:vloc : 40289.45s CPU 77333.70s WALL ( 4647479 calls)<br> h_psi:vnl : 8861.47s CPU 17711.55s WALL ( 4647479 calls)<br> add_vuspsi : 3444.16s CPU 6516.57s WALL ( 4647479 calls)<br> h_1psi : 41477.62s CPU 76557.86s WALL ( 4642098 calls)<br><br> General routines<br> calbec : 10946.34s CPU 20460.58s WALL ( 9291577 calls)<br> fft : 565.20s CPU 2108.77s WALL ( 127532 calls)<br> ffts : 131.27s CPU 473.36s WALL ( 46460 calls)<br> fftw : 41280.27s CPU 81952.36s WALL (27295400 calls)<br> interpolate : 317.92s CPU 1137.22s WALL ( 46460 calls)<br> davcio : 0.05s CPU 4562.48s WALL ( 401 calls)<br> <br> Parallel routines<br> fft_scatter : 32987.94s CPU 66128.36s WALL (27469392 calls)<br> <br> PWSCF : 1d 5h 8m CPU 2d 3h38m WALL<br><br> <br> This run was terminated on: 6: 1: 5 13Jun2021 <br><br>=------------------------------------------------------------------------------=<br> JOB DONE.<br>=------------------------------------------------------------------------------=<br></div><div>Thanks!</div></div>
</blockquote></div>