Dear Marcella,<div><br /></div><div>Before SGCPMD, I equilibrated the system with an NVT.<br />MAX_SCF_HIST i commented out in this file, but in the run it was not commented on (was set to 6). </div><div>In generale, you recommend a preconditioner such as FULL_ALL? <br /><div>"EPS_SCF  is rather tight to be a SGCPMD", but I stop the scf after 6 OT steps.</div><div><br /></div><div>Many thanks in advance,<br />Simone<br /></div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno venerdì 17 febbraio 2023 alle 15:46:24 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;"><br><br><div>Dear Simone, </div><div><br></div><div>Did you run BOMD with this system? </div><div>It is always better to start sampling with SGCPMD from a stable BOMD.</div><div><br></div><div><br></div><div>I don't know how you have determined the cutoff. One way to do it would be to calculate the BO forces with increasing cutoff to check the deviation. </div><div>Since the parameters have been tested for the "almost identical system", I understand that for that system under the same conditions the COM did not drift. Is it so?</div><div><br></div><div>Concerning the setting:</div><div>The time step is very large. 0.5 would be safer.</div><div>EPS_SCF  is rather tight to be a SGCPMD</div><div>The preconditions is poor</div><div>MAX_SCF_HIST is commented out</div><div><br></div><div>However it is difficult to judge without knowing how the wfn's convergence is behaving.</div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 17, 2023 at 2:58:55 PM 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><br>I chose the cutoff as the minimum possible. Without COMVEL_TOL there is drifting (along x) that must not be there. I would like to solve the drifting without the forcing of COMVEL_TOL (e.g. by increasing the cutoff so that the forces can be calculated better? )<br>The noisygamma was chosen from previous simulations on almost identical systems, as advised by the supervisor.<br>Why do you say that "

Your settings for SGCPMD look strange." ?<br><br>Many thanks in advance,<br>Simone<br><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno venerdì 17 febbraio 2023 alle 10:12:09 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>The cutoff looks indeed small. </div><div>Does the COM drift by BOMD? Have you tried without COMVEL_TOL? </div><div>Your settings for SGCPMD look strange. </div><div>Have you run the needed tests for the nosy-gamma to find the right parameters?</div><div><br></div><div>Regards</div><div>Marcella<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, February 17, 2023 at 2:25:00 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 cp2k-users,<br><br>I have a problem with COM (center of mass) drifting (along x). I'm running a LANGEVIN  (SGCPMD) on my 

heterostructure, consisting of a GST + TiTe2 (633 atoms), 

where I consider van der Waals interactions (with DFT-D3).  The cutoff is 320 Ry: maybe is it too small? Maybe at a temperature of 600K a gamma of 1.9e-3 is too big?<br><br>I noticed that by using COMVEL_TOL=0 I have the following behaviour: the temperature of the whole system does not correspond to the average (weighted according to the number of atoms) of the temperatures of the individual atomic species. <br><br><br>Many thanks in advance,<br>Simone<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/9ba34af5-66f0-40c9-aff5-dcfe339fbe90n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/9ba34af5-66f0-40c9-aff5-dcfe339fbe90n%40googlegroups.com</a>.<br />