<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">I believe you should use grimme correction to get a converged binding distance, you can also plot the total energies ranging along 2-4 Å and you will get a parabolic shape, with equation fitting you can get the binding distance corresponding to the minimal energy.</div></div>