[CP2K-user] [CP2K:18807] energy from GEO_OPT job does not equal ENERGY job?

[Aaron Ledray]

Hello! 

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: 

One thing I am wondering is if it is necessary to perform an *energy*  job 
after the *geo_opt* job, or if I can use the total force_eval energy from 
the final optimized structure of *geo_opt* as my energy evaluation? 

I am finding that the energy I am getting at the end of *geo_opt *jobs are 
consistently different by about ~0.05 a.u. from *energy* 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?

I've attached my input files. Energy.inp runs from the restart file made by 
geo_opt.inp.

Thank you kindly,
Aaron


-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/21d053ca-b333-4f4d-a89f-d8494af5dc15n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230512/752f612e/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: geo_opt.inp
Type: chemical/x-gamess-input
Size: 4552 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230512/752f612e/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: energy.inp
Type: chemical/x-gamess-input
Size: 2957 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230512/752f612e/attachment-0001.inp>