<div dir="ltr">Dear all, <br><br>I&#39;m trying to do cell relaxation on my Mg-Boron
systems. For some of the structures, I got extremely strange results.
The total energy diverges starting from some ionic step. For example,
here are the total energy in each ionic step for one of the problematic
run: <br>
<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -24.47168660 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -35.02062718 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -38.37822321 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -40.20237747 Ry<br>
!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -40.96387007 Ry<br>
!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.20956341 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.41389316 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.46640700 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.49930067 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.51263104 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.53987139 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.59139496 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.55162726 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.64537650 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.67440874 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -41.73712135 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -51.60413428 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -688.94041793 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -817.40812316 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -914.65691139 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; -980.98769592 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1145.13424734 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1079.87331043 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1148.61890203 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1147.67898328 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1147.50188067 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1148.72218340 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1149.83861471 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1690.39538237 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -2822.35642442 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -2718.63827111 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -2790.48198795 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -2823.74590405 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -2878.31317290 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -4478.65409327 Ry<br>

!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -5464.15854938 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =-14818.93675687 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =-17158.30103958 Ry<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =-18287.28161766 Ry<br><br>
and the forces:<br>
<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 39.327209&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000071<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 11.517813&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000029<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 8.564152&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000058<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 3.973671&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000035<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 1.056249&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000106<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.617165&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000092<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.326250&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000046<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.416678&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000033<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.227223&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000033<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.501912&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000042<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.269395&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000104<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.369080&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000049<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.334142&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000026<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.148318&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000034<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.145030&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000012<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.143331&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000064<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 116.856723&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000084<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 65.417648&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000066<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 171.381517&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000070<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 155.577114&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000070<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 184.877111&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000057<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 109.156982&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000037<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 142.602012&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000080<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 88.522593&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000117<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 127.562744&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000058<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 92.956163&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000108<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 88.462372&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000117<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp; 84.376621&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000056<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 134.526750&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000106<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 282.688334&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000132<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 600.035535&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000078<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 224.477965&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000122<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 398.869729&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000123<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 396.502768&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000149<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 420.885999&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000136<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp; 443.182054&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000031<br>

&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp; 3819.452849&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000261<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp; 1265.093632&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000207<br>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp; 1027.278240&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000338<br>

<br>For this problematic system, this behavior is very stable. No
matter I use &#39;cg&#39; or &#39;david&#39;, change the mixing parameter from 0.9 to
0.1, I still get similar problem. <br><br>I cannot find anything
relevant in the archive. I wonder if anyone knows what&#39;s going on here
and can help me get rid of this annoying problem. <br><br>btw, the input, output files and pseudo-potentials can be found at: <a href="http://volga.eng.yale.edu/hui/MgonBS.tar.gz">http://volga.eng.yale.edu/hui/MgonBS.tar.gz</a><br>

<br><br><br><br>Thank you very much,<br>---------------------------------------------<br>Hui Tang<br>Graduate Student<br>Department of Applied Physics<br>Yale University<br>Email: <a href="mailto:hui.tang@yale.edu">hui.tang@yale.edu</a><br>
Phone: 1-203-432-7763<br><br><br></div>