I am trying to calculate an H-center in CaF2, and I am having trouble with the geometry optimization. <br>I have tried adjusting various parameters such as ion dynamics, trust radius, convergence tresholds, but it seems that convergence cannot be improved. <br>
<br>Using bfgs I always get the following error: <br><br>%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     from bfgs : error #         1<br>     bfgs history already reset at previous step<br> %%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>Alternatively the electonic convergence gets extremely slow. <br>
<br>In addition, no matter what parameters I adjust the total forces in the system never converge, e.g.<br>     Total force =     0.026394     Total SCF correction =     0.000314<br>     Total force =     0.018379     Total SCF correction =     0.000282<br>
     Total force =     0.018824     Total SCF correction =     0.000120<br>     Total force =     0.018857     Total SCF correction =     0.000155<br>     Total force =     0.020215     Total SCF correction =     0.000193<br>
     Total force =     0.058765     Total SCF correction =     0.000241<br>     Total force =     0.016998     Total SCF correction =     0.000205<br>     Total force =     0.042967     Total SCF correction =     0.000238<br>
     Total force =     0.106226     Total SCF correction =     0.000073<br>     Total force =     0.137142     Total SCF correction =     0.000071<br>     Total force =     0.109549     Total SCF correction =     0.000108<br>
     Total force =     0.068723     Total SCF correction =     0.000103<br>     Total force =     0.076599     Total SCF correction =     0.000168<br>     Total force =     0.033233     Total SCF correction =     0.000342<br>
     Total force =     0.026883     Total SCF correction =     0.000248<br>     Total force =     0.025846     Total SCF correction =     0.000378<br>     Total force =     0.025596     Total SCF correction =     0.000225<br>
     Total force =     0.025482     Total SCF correction =     0.000275<br>     Total force =     0.025508     Total SCF correction =     0.000128<br><br>I am thankful for any suggestions. <br><br>Stephan Rix <br>Institute for Inorganic und Analytical Chemistry <br>
Johannes Gutenberg-University Mainz, Germany <br>