<div dir="ltr">Dear Francesco,<div><br></div><div>I am interesting to obtain the hamiltonian of mixed system for the the harmiltonian of reactant and product using the eq. H(k) = k*H(k) + (1-k)H(k). I have know the second and third force_eval sections contain the coordinate of the reactant and product. But in the first force_eval, which one should be put into the &COORD ... &END COORD, reactant's or product's??</div><div><br></div><div>e.g.<br>for reaction:</div><div>H2SO4 --> HSO4- + H+</div><div><br></div><div>The second and third force_eval should put the COORD. of H2SO4 and that of HSO4-, respectively. What about in the first one?</div><div><br></div><div>Thanks,</div><div><br></div><div>Ning</div><br>在 2014年12月17日星期三 UTC+8下午11:14:57，Francesco Ambrosio写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Samuel<br><br>Thank you for your kind reply. I will test this asap.  <br><br>Kind Regards<br><br>Francesco<br><br>On Wednesday, December 17, 2014 4:03:17 PM UTC+1, Samuel Andermatt wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">You will need 3 force_eval sections. The first one has the keyword: FORCE_EVAL&MIXED&LINEAR&LAMBDA there you give your k parameter in.<br>The second and third force_eval sections contain the input for the charged and the neutral calculation.<br></div></blockquote></div></blockquote></div>