Dear CP2K developers, <br />    I am sorry to bother you for a technical detail regarding the treatment of atom mass in the mixing force_evals calculation.<br />    I intend to investigate the free energy difference between Na-Cl and K-Cl ion pairs in aqueous solutions through FPMD-based thermodynamic coupling constant integration, i.e., transitioning Na(+) gradually into K(+) using the &MIXED section:<div><img alt="CP1.png" width="112px" height="82px" src="cid:e65f9154-ae4f-4157-94bf-d682f8dc9766" /><br /></div><div>together with the following settings in the &SUBSYS section:<br /></div><div><img alt="CP2.png" width="126px" height="152px" src="cid:1ac8dc79-49db-4126-9ae0-7eb5a79b6906" /><br /></div><div>    During the simulation, energies and forces are mixed via lambda to solve the equation of motion. On the other hand, the masses of atoms are also required in the equation of motion. I am wondering how CP2K handles the mass of “Na” mentioned above in simulations with 0<lambda<1. Specifically, is it calculated as 0.5*(mass_Na+mass_K), lambda*mass_Na+(1- lambda)*mass_K, or utilizing another algorithm? Thanks in advance!<br /></div><div><br /></div><div>Sincerely,<br />Xu<br /></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/28957671-4290-4ca1-bccb-84e87f5826adn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/28957671-4290-4ca1-bccb-84e87f5826adn%40googlegroups.com</a>.<br />