<div>Added note:</div><div>I just saw another conversation which is similar: https://groups.google.com/g/cp2k/c/EMxxtX0sFK4</div><div>I missed it initially as it is too new, and does not come up on google. <br></div><div><br></div><div>It is stated that correct Gamma would result in stable energy partition. For determining the parameters, I was simulating the metal system at 1000K as a simple LJ liquid.</div><div>I observed in output file that K.E. is stable, however P.E. was drifting in the time scale that I simulated it for (512 atoms, 2 ns), which might be due to melting lattice. Is that ok?<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, January 2, 2023 at 4:41:04 PM UTC-6 Amit Gupta 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;"><div>Following the conversation of this thread: <a href="https://groups.google.com/g/cp2k/c/fmZX3_2kivc" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/g/cp2k/c/fmZX3_2kivc&source=gmail&ust=1672786552533000&usg=AOvVaw0RfAX0e7DY3U3fhOUC4WAe">https://groups.google.com/g/cp2k/c/fmZX3_2kivc</a></div><div><br></div><div>I did a parameter sweep of my system using FIST dynamics, and determined that<br></div><div>GAMMA = 0.35</div><div>NOISY_GAMMA=1e-5 </div><div><br></div><div>yields best thermostat stability in my system. The paper mentioned in above conversation however uses slightly different terminology of Gamma_L (Langevin friction) and Gamma_D(intrinsic friction). And their values are rather low (1e-8 for Gamma_D and 0 for Gamma_L ). I would like to understand</div><div><br></div><div>1. How are GAMMA, NOISY_GAMMA, and SHADOW_GAMMA are related to Gamma_L and Gamma_D?</div><div>2. Does the value I got (0.3 and 1e-5) pass the initial smell test? Or do they look non-physical (I mostly use Nose thermostat but needs to use Langevin for this study for replicating a paper)? </div><div><br></div><div>Thank you<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/64318746-dc51-4537-a690-c4d15d8f44e9n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/64318746-dc51-4537-a690-c4d15d8f44e9n%40googlegroups.com</a>.<br />