Hi,<br><br>I have manage to (wrongly) run a force matching run based on a  ab initio trajectory I got from siesta. I have however a doubt about the energy that should be included in the<br><table class="default"><tbody><tr><td class="l"><br></td>
                  <td class="r"><big class="uctt">REF_FORCE_FILE_NAME FRC.XYZ</big></td></tr></tbody></table><br>the lines containing the energy have the format<br> i =        0, time =        0.000, E =        -0.0053069965<br>I was wondering what energy what energy should be put there from my DFT MD run? I suppose that I could use the Kohn-Sham energy and convert it in Hartree units but where should I put the 0?<br><br>Are the energies needed if I use <font size="2"><span style="font-family: arial,sans-serif;"><a href="http://manual.cp2k.org/trunk/CP2K_INPUT/OPTIMIZE_INPUT/FORCE_MATCHING/COMPARE_FORCES.html">COMPARE_FORCES</a> instead of</span></font><a href="http://manual.cp2k.org/trunk/index.html#CP2K_INPUT/OPTIMIZE_INPUT/FORCE_MATCHING/COMPARE_ENERGIES.html"> COMPARE_ENERGIES</a>? <font size="2"><span style="font-family: arial,sans-serif;"> <font size="2"><font size="2">W</font>hat should I put in EACH <font size="2">for co<font size="2">m<font size="2">paring every <font size="2">n ste<font size="2">p of <font size="2">my refere<font size="2">nce<font size="2"> force<font size="2">/</font>trajectory file?</font></font></font></font></font></font></font></font></font><br></span></font><br><br><br><br>I suppose that they should be <br><br>Le lundi 24 septembre 2012 13:30:46 UTC+2, Chuong a écrit :<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Sorry,<br><br> I forgot to run argon-ref.inp before the other scripts.<br><br>Le vendredi 27 juillet 2012 10:52:26 UTC+2, Chuong a écrit :<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br><br>There is a force matching section in the manual of cp2k <a href="http://manual.cp2k.org/trunk/CP2K_INPUT/OPTIMIZE_INPUT/FORCE_MATCHING.html" target="_blank">http://manual.cp2k.org/trunk/<wbr>CP2K_INPUT/OPTIMIZE_INPUT/<wbr>FORCE_MATCHING.html</a> but it is not very clear to me how can I define which parameters I would like to do fit or does it try to fit all the parameters? It also says that <i>the filename of the reference forces, should also contain the energy</i>  but in which order?, units? Is there an example available for the force matching procedure with cp2k?<br><br>Yours,<br><br><pre cols="72">Nguyen Huu Chuong
Institut für Festkörpertheorie
Westfälische Wilhelms-Universität
Wilhelm-Klemm-Straße 10, 48149 Münster</pre><br></blockquote></blockquote>