<div dir="ltr"><div dir="ltr"><div dir="ltr">HI.</div><div>To me the mixing_beta is slightly too high.</div><div dir="ltr">Try with 0.02 </div><div dir="ltr"><br></div><div>Regards</div><div>ff</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 6, 2019 at 10:03 AM Gui Wei <<a href="mailto:201707021019@cqu.edu.cn">201707021019@cqu.edu.cn</a>> wrote:<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>Dear QE users,</div><div>I am trying to optimize the adsorption system that contain Fe(110) substrate and single-layer graphene connected to an oxygen atom.I got the total force and the energy like this.At first, everything was fine,but after 0.006846 the total force began to increase ,causing the total force not to converge ,it is a extremely strange result.Is there any way to solve this problem?</div><div><br></div><div>Any suggestions would be greatly appreciated.</div><div>Thank you.</div><div><br></div><div><br></div><div>Total force = 0.316790 energy new = -4107.2627566975 Ry</div><div>Total force = 0.246608 energy new = -4107.3151393348 Ry</div><div>Total force = 0.156503 energy new = -4107.3635865834 Ry</div><div>Total force = 0.119628 energy new = -4107.3895957292 Ry</div><div>Total force = 0.078444 energy new = -4107.3969303863 Ry</div><div>Total force = 0.052366 energy new = -4107.4017754787 Ry</div><div>Total force = 0.040309 energy new = -4107.4042249914 Ry</div><div>Total force = 0.033754 energy new = -4107.4058726129 Ry</div><div>Total force = 0.028026 energy new = -4107.4071437156 Ry</div><div>Total force = 0.025319 energy new = -4107.4081667548 Ry</div><div>Total force = 0.024093 energy new = -4107.4089392392 Ry</div><div>Total force = 0.022277 energy new = -4107.4096467799 Ry</div><div>Total force = 0.019561 energy new = -4107.4102126980 Ry</div><div>Total force = 0.016523 energy new = -4107.4106348116 Ry</div><div>Total force = 0.014292 energy new = -4107.4109663982 Ry</div><div>Total force = 0.013121 energy new = -4107.4112149183 Ry</div><div>Total force = 0.011934 energy new = -4107.4114021825 Ry</div><div>Total force = 0.011261 energy new = -4107.4115692713 Ry</div><div>Total force = 0.010521 energy new = -4107.4117253881 Ry</div><div>Total force = 0.009179 energy new = -4107.4118535327 Ry</div><div>Total force = 0.008233 energy new = -4107.4119532022 Ry</div><div>Total force = 0.007949 energy new = -4107.4120378076 Ry</div><div>Total force = 0.007614 energy new = -4107.4121212783 Ry</div><div>Total force = 0.007100 energy new = -4107.4121953505 Ry</div><div>Total force = 0.006846 energy new = -4107.4122538477 Ry</div><div>Total force = 0.006967 energy new = -4107.4123143593 Ry</div><div>Total force = 0.007559 energy new = -4107.4123893617 Ry</div><div>Total force = 0.007912 energy new = -4107.4124672492 Ry</div><div>Total force = 0.007587 energy new = -4107.4125386594 Ry</div><div>Total force = 0.007532 energy new = -4107.4126143150 Ry</div><div>Total force = 0.008046 energy new = -4107.4126949028 Ry</div><div>Total force = 0.008145 energy new = -4107.4127766088 Ry</div><div>Total force = 0.008142 energy new = -4107.4128651944 Ry</div><div>Total force = 0.009356 energy new = -4107.4129774552 Ry</div><div>Total force = 0.011161 energy new = -4107.4131286975 Ry</div><div>Total force = 0.012473 energy new = -4107.4133389405 Ry</div><div>Total force = 0.013542 energy new = -4107.4136064842 Ry</div><div>Total force = 0.014769 energy new = -4107.4138982197 Ry</div><div>Total force = 0.015571 energy new = -4107.4142168735 Ry</div><div>Total force = 0.015344 energy new = -4107.4145642683 Ry</div><div>Total force = 0.015134 energy new = -4107.4149145704 Ry</div><div>Total force = 0.016511 energy new = -4107.4152819587 Ry</div><div>Total force = 0.016788 energy new = -4107.4156041418 Ry</div><div><br></div><div><br></div><div> &control</div><div> calculation = 'relax',</div><div> prefix='epoxy1',</div><div> restart_mode='from_scratch' </div><div> pseudo_dir = '/public/home/duan1/GuiWei/pseudo/',</div><div> outdir='/public/home/duan1/GuiWei/tempdir/'</div><div> tstress = .true.</div><div> tprnfor = .true.</div><div> etot_conv_thr=1.0D-4</div><div> forc_conv_thr=1.0D-3</div><div>/</div><div> &system </div><div> ibrav= 0, </div><div> nat=125, </div><div> ntyp=3,</div><div> nspin = 2 </div><div> starting_magnetization(1) =0.4 </div><div> occupations='smearing',</div><div> smearing='mp',</div><div> degauss=0.02,</div><div> ecutwfc =37.0</div><div> ecutrho =300.0 </div><div> vdw_corr='DFT-D' </div><div> london_s6 = 0.75</div><div>/</div><div>&electrons</div><div> conv_thr = 1.0d-8</div><div> mixing_beta = 0.2</div><div> mixing_mode = 'local-TF'</div><div> mixing_ndim=16</div><div> electron_maxstep=200</div><div>/</div><div>&IONS</div><div> ion_dynamics='bfgs'</div><div>/</div><div>CELL_PARAMETERS angstrom</div><div>19.68000000000000 0.000000000000011 0.000000000000000</div><div>0.000000000000000 8.521689973238880 0.000000000000000</div><div>0.000000000000000 0.000000000000000 25.00000000000000</div><div>ATOMIC_SPECIES</div><div>Fe 55.847 Fe.pbe-nd-rrkjus.UPF</div><div>C 12.0107 C.pbe-rrkjus.UPF</div><div>O 16.00 O.pbe-rrkjus.UPF</div><div>ATOMIC_POSITIONS crystal</div><div>......</div><div>K_POINTS automatic</div><div>2 2 1 0 0 0</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div>Wei gui</div><div>School of Mechanical Engineering，Chongqing University, China</div></div><div><br></div><div><br></div>_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><span style="background-color:transparent;font-family:Arial;font-size:10pt;font-weight:700;white-space:pre-wrap">Federico IORI</span><br></div><div><span><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:10pt;font-family:Arial;background-color:transparent;font-style:italic;vertical-align:baseline;white-space:pre-wrap">Computational material scientist</span></p><p style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="background-color:transparent;color:rgb(61,133,198);font-family:Arial;font-size:11pt;font-style:italic;font-weight:700;white-space:pre-wrap">Paris-Saclay Research Center</span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:9pt;font-family:Arial;background-color:transparent;vertical-align:baseline;white-space:pre-wrap"><a href="https://www.google.com/maps/place/Air+Liquide/@48.8297381,2.2016685,11.75z/data=!4m5!3m4!1s0x47e67e61a4fbbdc7:0xca3bea9e80059880!8m2!3d48.7630415!4d2.1333045" style="color:rgb(17,85,204);font-family:arial,sans-serif;font-size:12.8px;white-space:normal" target="_blank">1 chemin de la Porte des Loges</a><br style="font-family:arial,sans-serif;font-size:12.8px;white-space:normal"><span style="font-family:arial,sans-serif;font-size:12.8px;white-space:normal">Les Loges en Josas – 78354 Jouy en Josas cedex</span><br style="font-family:arial,sans-serif;font-size:12.8px;white-space:normal"><span style="font-family:arial,sans-serif;font-size:12.8px;white-space:normal">Mail</span><span style="font-family:arial,sans-serif;font-size:12.8px;white-space:normal">: <a href="mailto:federico.iori@airliquide.com" target="_blank">federico.iori@airliquide.com</a></span><br></span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:9pt;font-family:Arial;background-color:transparent;vertical-align:baseline;white-space:pre-wrap">Phone: +33 7 621 605 15</span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><span style="font-size:9pt;font-family:Arial;background-color:transparent;vertical-align:baseline;white-space:pre-wrap"><img src="https://lh5.googleusercontent.com/jWge08q7BtEKOlxuBF_cj804ywzhJmkjw9rEb5dBZLKyUL_HQkVwhcwNi9vtzVTj-EiY-wQhE18VPSACNFaAFCAWvxp6bEwkTDR6-IjGBqp-Y3Afzx-OxN8OHxXYPLS0Dz9tP8fr" width="165" height="40" style="border: none;"></span></p><p dir="ltr" style="line-height:1.38;margin-top:0pt;margin-bottom:0pt"><br></p></span></div></div></div></div></div></div></div></div></div></div></div></div></div>