Hi Matt, <div><br></div><div>You should be able to use different gammas by activating thermal regions addressing specific subsystems. </div><div>You need a LANGEVIN section with GAMMA and a THERMAL_REGION section with DEFINE_REGION subsections, </div><div>each of them containing a NOISY_GAMMA_REGION keyword and a LIST of particles. </div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, July 14, 2021 at 5:41:18 PM UTC+2 Matt D 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;">Hi all,<div><br></div><div>I'm trying to run second generation Car-Parrinello MD (i.e. Thomas Kühne's method in Kuhne, TD; Krack, M; Mohamed, FR; Parrinello, M., PHYSICAL REVIEW LETTERS, 98 (6), 066401 (2007) ) in CP2K but using different values of NOISY_GAMMA for different species.</div><div><br></div><div>In Kühne's paper it suggests this is possible, but I can't seem to find the relevant keyword in the manual. I've also seen different values of NOISY_GAMMA being used in recent papers using SGCPMD, for example in work from Marcella Ianuzzi and Jun Cheng's groups (JACS Au 2021, 1, 569−577) so I believe it is possible.</div><div><br></div><div>So far I have tried:</div><div><br></div><div>NOISY_GAMMA {value1} {value2} {value 2}</div><div>and </div><div>GAMMA_D {value1} {value2} {value3}</div><div><br></div><div>but CP2K version 7.1 doesn't recognise GAMMA_D as a parameter (I informally heard it used to used in earlier implementations) and NOISY_GAMMA is expected to be a single value, rather than multiple values.</div><div><br></div><div>If anybody can share/suggest what might be the correct syntax I'd greatly appreciate it.</div><div><br></div><div>Thanks.</div><div><br></div><div>Best wishes,</div><div>Matt</div></blockquote></div>