Hello! <div><br /></div><div>I'm trying to run a DFT calculation for a porphyrin bound to a small molecule and compare energy to that of the unbound porphyrin state and the small molecule alone: <div><div><br /></div><div>One thing I am wondering is if it is necessary to perform an <i>energy</i>  job after the <i>geo_opt</i> job, or if I can use the total force_eval energy from the final optimized structure of <i>geo_opt</i> as my energy evaluation? </div><div><br /></div><div>I am finding that the energy I am getting at the end of <i>geo_opt </i>jobs are consistently different by about ~0.05 a.u. from <i>energy</i> jobs of the same molecule (the small-molecule bound state), with the geo_opt being lower in energy. Is this expected, is there something I'm not considering?</div><div><br /></div><div>I've attached my input files. Energy.inp runs from the restart file made by geo_opt.inp.</div><div><br /></div><div>Thank you kindly,</div><div>Aaron</div><div><br /></div><div><br /></div></div></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/21d053ca-b333-4f4d-a89f-d8494af5dc15n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/21d053ca-b333-4f4d-a89f-d8494af5dc15n%40googlegroups.com</a>.<br />