Dear Marcella,<br />Thanks for your reply. I want to ask you for some clarification on the input part to be added:<br /><br />1) If I set DO_LANGEVIN_DEFAULT F and DO_LANGEVIN F I am associating my thermal region with an NVE molecular dynamics, right? I just want to "force" the titanium atoms to thermalise at 600 K, respecting a LANGEVIN. I must therefore choose .TRUE., right?<br /><br />2) If I exclude NOISY_GAMMA_REGION, I am associating the thermal region with the NOISY_GAMMA_REGION chosen on the total system, right?<br /><br />3) TEMP_TOL helps me with thermalisation, however doesn't this create 'non-physical' molecular dynamics, i.e. molecular dynamics from which I cannot extract dynamic properties such as the diffusion coefficient?  <br /><br />4) About LIST: I follow the implicit numbering of the atoms in the input file, right?<br /><br />Many thanks in advance,<br />Simone<br /><br />&THERMAL_REGION<br />       DO_LANGEVIN_DEFAULT  T<br />       &DEFINE_REGION<br />             DO_LANGEVIN T<br />             LIST  n1..n2<br />             TEMPERATURE     600<br />             #NOISY_GAMMA_REGION     2.1999999999999998E-04<br />             #TEMP_TOL 50<br />       &END DEFINE_REGION<br /> &END THERMAL_REGION<br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno mercoledì 8 febbraio 2023 alle 11:59:56 UTC+1 Marcella Iannuzzi ha scritto:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear Simone, <div><br></div><div>You may define the thermal regions and set tuned  noisy gamma for each region independently.</div><div><a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT.html&source=gmail&ust=1676077673960000&usg=AOvVaw0gcIyUv7tI1_1zb7GDpkAg">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION.html&source=gmail&ust=1676077673960000&usg=AOvVaw0Kt0nm87-bc-9wvr-HBKvT">MOTION</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD.html&source=gmail&ust=1676077673960000&usg=AOvVaw3yXhfvy-owbnDkZHhw2VSe">MD</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION.html&source=gmail&ust=1676077673960000&usg=AOvVaw1iY5iJ3zsOX8u40GFlntmm">THERMAL_REGION</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION/DEFINE_REGION.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION/DEFINE_REGION.html&source=gmail&ust=1676077673960000&usg=AOvVaw1HzUyMs3c8peZCzv8o3ROh">DEFINE_REGION</a></div><div><br>Regards</div><div>marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, February 6, 2023 at 5:25:48 PM UTC+1 <a href data-email-masked rel="nofollow">sim.ri...@gmail.com</a> 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>Dear cp2k-users,</div><br>I have thermalisation problems with my heterostructure, consisting of a GST + TiTe2 (633 atoms), where I consider van der Waals interactions (with DFT-D3). Initially I have a relaxed structure (CELL_OPT), which I want to thermalise at 600K. I am using the SGCPMD (&LANGEVIN). The problem is that the titanium atoms do not thermalise at 600K, but at a lower temperature, whereas the other atoms do so acceptably. I have also tried thermalising initially with an NVT and then using SGCPMD, but the problem persists. <br>I cannot use TEMP_TOL so as not to affect dynamic properties such as diffusion etc etc.<br><br>This also happens at 1300K (melting of GST). <br><br>In general, I have tried different gamma values. <br><br>I attached the input file.<br><br>Many thanks in advance,<br>Simone<br><br></blockquote></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/90d9faa6-de2e-4714-9d48-70ab56bf09ecn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/90d9faa6-de2e-4714-9d48-70ab56bf09ecn%40googlegroups.com</a>.<br />