Dear Massimo<div><br></div><div>In the MD the reference cell is used, which is quite different from the cell of the single point calculation.</div><div>Kind regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, September 22, 2020 at 10:55:14 AM UTC+2 mas...@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;">Dear CP2K users/developers,<div><br></div><div>I am extracting some snapshots from an umbrella sampling simulation (with PLUMED as dependency of CP2K 5.1) and I would like to recalculate their energy at a higher level of theory. So to start I performed a single point calculation at the same level of theory of the US simulation (PBE-D3) and, surprisingly for me, I got very different results in the energy. In particular, the energy of the MD step is:</div><div>ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -3520.987951048081868<br></div><div>While for the single point calculation:</div><div>ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):            -3520.903314334958850<br></div><div>Which is an enormous difference. Attached the two input scripts, the main differences being the RUN_TIPE, the cell parameters (for the single point they are just taken from the md output at the selected step) and the lack of the MOTION section.</div><div>Does anybody have an idea on where this large energy difference could come from?</div><div><br></div><div>Thank you in advance, </div><div>Massimo Bocus</div></blockquote></div>