Dear Marcella, <br />is the setting I proposed to "force" titanium atoms to thermalise at 600K (i.e. without giving a NOISY_GAMMA_REGION , keeping the LANGEVIN and without TEMP_TOL) equivalent to not defining a thermal region for titanium atoms? In other words: is making the temperature of 600K explicit for titanium atoms a non-trivial definition of a thermal region?<br /><div><br />Many thanks in advance,<br />Simone</div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno venerdì 10 febbraio 2023 alle 16:54:11 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;"><div><br></div><div><br></div>Dear Simone <div><br></div><div>It is all correct. The keyword TEMPERATURE within the thermal region is used for initialisation of the velocities</div><div><br></div><div>Kind regards</div><div>Marcella<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 10, 2023 at 2:48:21 AM UTC+1 Simone Ritarossi 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">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" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT.html&source=gmail&ust=1676161869305000&usg=AOvVaw0TbroAW8z4raoIm0RNOBQg">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION.html" rel="nofollow" target="_blank" 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=1676161869306000&usg=AOvVaw0s0ZoD99SkKO2fWTAW0wsv">MOTION</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD.html" rel="nofollow" target="_blank" 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=1676161869306000&usg=AOvVaw1Dy_ZNF4YPqR5CGtLV9_uv">MD</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION.html" rel="nofollow" target="_blank" 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=1676161869306000&usg=AOvVaw0uy5AFqJTPn8dUKP1SEpqA">THERMAL_REGION</a> / <a href="https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMAL_REGION/DEFINE_REGION.html" rel="nofollow" target="_blank" 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=1676161869306000&usg=AOvVaw3IMxZm8XSYeDXvkodFm4Fs">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 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></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/bdf5069f-0eaa-4578-ad01-0453f8b4eb85n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/bdf5069f-0eaa-4578-ad01-0453f8b4eb85n%40googlegroups.com</a>.<br />