Hi Alex,<div><br></div><div>The quality of the forces you calculate will depend on the accuracy to which the SCF problem is solved. The tolerance for the latter is given by EPS_SCF - a typical value for EPS_SCF during a geometry optimisation would be around 1e-7 Ha. Another important parameter for the SCF tolerance is EPS_DEFAULT, this should be decreased from the default value of 1e-10 to 1e-14 to complement the value of 1e-7 for EPS_SCF. The error in the forces with respect to these parameters should be sqrt(EPS_SCF) < error on force; and sqrt(EPS_DEFAULT)  < EPS_SCF.</div><div><br></div><div>The default units in CP2K are always atomic units, so the force is in Hartree/Bohr. If you want to change the units you put in the input, this is possible but I have forgotten the exact syntax. I think it is something like:</div><div><br></div><div>KEYWORD [unit] 0.025</div><div><br></div><div>and the "unit" is replaced by one of the units in this list:</div><div><br></div><div>https://manual.cp2k.org/trunk/units.html<br></div><div><br></div><div>Best wishes,</div><div>Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, 25 April 2022 at 21:56:44 UTC+1 alexvd...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi,</div><div><br></div><div>I'm curious what the precision on forces in a cp2k energy_forces job. In particular, how does it depend on various input variables, such as EPS_SCF? <br></div><div><br></div><div>Also, in the forces section, the printout says the units are in atomic units. Is that also the default unit for all force dependent convergence criteria, such as in a NEB job?</div><div><br></div><div>Thanks,</div><div>Alex<br></div></blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0dfdc809-93cf-41af-976f-68e7bf7a56b9n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/0dfdc809-93cf-41af-976f-68e7bf7a56b9n%40googlegroups.com</a>.<br />