<div dir="ltr">Hi Linda,<div><br></div><div>Thanks for your reply.</div><div><br></div><div><i>But if you use MULTIPLE_FORCE_EVALS,you can draw your initial model for product and reactant in one system, then arrange the coordinate in accordance with the same type of atom.</i></div><div><br></div><div>I don't understand you. Would you give me an example about this?</div><div><br></div><div>Looking forward to your responses.</div><div><br></div><div>Happy new year.</div><div><br></div><div>Ning<br><br>在 2015年6月24日星期三 UTC+8上午12:07:06，Linda Sundarti写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hai Ning,<div><br></div><div>Maybe you can see this exercise link to solve your problem:</div><div><br></div><div><a href="http://wiki.wpi.edu/deskinsgroup/Learn_CP2K" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwiki.wpi.edu%2Fdeskinsgroup%2FLearn_CP2K\46sa\75D\46sntz\0751\46usg\75AFQjCNFf66nRbIynRUN66vaZauOkUyBi9g';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwiki.wpi.edu%2Fdeskinsgroup%2FLearn_CP2K\46sa\75D\46sntz\0751\46usg\75AFQjCNFf66nRbIynRUN66vaZauOkUyBi9g';return true;">http://wiki.wpi.edu/<wbr>deskinsgroup/Learn_CP2K</a><br></div><div><br></div><div>May I simplify your problem, I suggest you to entry the coordinate of H2SO4 firstly and run FORCE_EVALS, then you can entry the next coordinate molecule to get next FORCE_EVALS and so on, in the end you can compare the different of total energy to get gap energy. But if you use MULTIPLE_FORCE_EVALS,you can draw your initial model for product and reactant in one system, then arrange the coordinate in accordance with the same type of atom.</div><div><br></div><div>Good luck</div><div><br></div><div>Cheers,</div><div>Linda</div><div><br></div><div><br></div><div><br></div><div>On Saturday, June 20, 2015 at 10:58:54 AM UTC+7, <a>ning....@gmail.com</a> wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi cp2k expert,</div><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 known that it need to use the MULTIPLE_FORCE_EVALS, in which 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.</div><div>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? Furthermore, How can get the average energy gap of reaction from the Hamiltonian??</div><div>Any suggestions are appreciate.</div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Ning</div></div></blockquote></div></div></blockquote></div></div>