<div dir="ltr"><br><br>On Wednesday, October 15, 2014 1:19:59 PM UTC-4, Noam Bernstein wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi all - I'm trying to do a calculation for a molecule in vacuum, with
<br>the nonperiodic MT solver, and I'd like to get the most accurate
<br>forces I can, and by that I mean with minimal egg-box effects, and
<br>consistent with the energy.
<br>
<br>It's an organic molecule, which I'm doing with BLYP and DZVP basis
<br>sets.  What do I need to converge with respect to, and do these values
<br>seem reasonable?
<br>
<br>CUTOFF 400 eV
<br>CUTOFF_REL 60 eV
<br>
<br>EPS_SCF 1.0e-8 (inner and outer)
<br>XC_DERIV NN10_SMOOTH
<br>XC_SMOOTH_RHO NN10
<br>
<br>I can't use USE_FINER_GRID in &XC_GRID because it conflicts with the
<br>MT solver I'm using.
<br><br></blockquote><div><br></div><div>I just noticed that the original message said eV for the cutoff, which is wrong, since I'm using the default units which are Ry. </div><div><br></div><div>In any case, does no one have any suggestions for how to get egg-box forces below 0.01 eV/A? I can relax the molecule and get the forces lower, but if I rotate it (to align with the cartesian axes) the forces jump back to 0.01 eV/A, and I'd like them to stay about 10 times smaller.</div><div><br></div><div>      thanks,</div><div>      Noam</div></div>