Dear Laura,<div><br></div><div>Difficult to guess what is going wrong with graph.popt. Maybe there is some normalisation factor that is wrongly set for the calculation of a multi-dimensional FES.</div><div>Maybe more information would help.</div><div>The python script you are indicating, though, it is not generating a FES from a metadynamics  potential. It builds a probability histogram from a standard molecular dynamics sampling.</div><div>This cannot be directly applied to an enhanced sampling, because you are not taking into account the bias potential in collecting the probability distribution. </div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, June 3, 2022 at 6:51:17 PM UTC+2 laura.b...@edu.unito.it 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;">Hi everyone!<div>I am using CP2K for running metadynamics simulations. I used the graph.popt tool for obtaining the free energy of my system, but actually I can't understand how it is calculated. I read that the final energy results are in hartree but if I convert them in eV they do not make any sense in terms of order of magnitude, even if the obtained curve is in agreement with what I expected. I know that the free energy should be calculated considering the probability distribution of the CVs but if I try to reconstruct the curve in this way, I don't obtain the same results obtained by exploiting the graph.popt tool. I also tried to use the python script found at  <a href="https://www.cp2k.org/exercises:2021_uzh_acpc2:ex03" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://www.cp2k.org/exercises:2021_uzh_acpc2:ex03&source=gmail&ust=1654586951797000&usg=AOvVaw0_svH9WGzVXrUIPbLgWpBt">https://www.cp2k.org/exercises:2021_uzh_acpc2:ex03</a> and the range of free energy is reasonable in this case, but completely different respect to what obtained by using graph.popt, so I think that I'm not understanding what is printed by this program.</div><div>Can you help me to solve this question?</div><div><br></div><div>Kind regards,</div><div>Laura</div><div><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/5e483eb7-80c7-4a54-a9b6-1a3c588a9459n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5e483eb7-80c7-4a54-a9b6-1a3c588a9459n%40googlegroups.com</a>.<br />