<div>There is a howto on optimizing the CUTOFF in the manual (https://manual.cp2k.org/trunk/methods/dft/cutoff.html). In your case, you should check the forces and stress tensor (especially the diagonal elements) as well (signed values and sum of atomic forces) by turning on the respective PRINT section (https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/PRINT/FORCES.html, https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/PRINT/STRESS_TENSOR.html). Beware that stress tensors are generally not very accurate because with changes of the box dimensions, the integration grids will also change.<br /></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Michela Benazzi schrieb am Freitag, 26. Juli 2024 um 17:48:41 UTC+2:<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 Frederick, thank you very much for sharing this!<br><br>The spike occurs as soon as the simulation starts, hence why I shared the first couple frames. Would you still like me to include more?<br><br><div>- CP2K is not printing any warnings.</div><div>- Do you have any suggestions on what to look for when checking for accuracy? I have changed CUTOFF in the past.</div><div>- I could probably afford a larger basis set. The current one I am using was the same used in prior research in my lab, where admittedly they did not do much benchmarking! I can read more about zeta quality.</div><div>- I will also fix smearing. <br></div><div><br></div><div>Sincerely thank you so much!! This gave me a good starting point and different opinions. Best,</div><div><br></div><div>Michela<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Friday, July 26, 2024 at 10:43:53 AM UTC-4 Frederick Stein 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 Michela,</div><div>Your output is not really helpful to me, it would be better to have a few frames around the spike. Does CP2K print any warnings? Spikes within MDs may hint to too inaccurate setups (or just bad luck). My suggestions are:</div><div>- did you check the accuracy of your setup (especially CUTOFF, EPS_DEFAULT)</div><div>- are you sure that a single zeta basis set provides a sufficient accuracy or can't you afford a larger basis set? For PBE, one should employ a basis set of at least double-zeta, maybe even triple-zeta quality.<br></div><div>- if you are already using smearing, you could use the actual temperature as the electronic temperature (seems more physical to me)</div><div>For more advanced suggestions, people more experienced in MD are open to jump in.</div><div>Best,</div><div>Frederick<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">Michela Benazzi schrieb am Freitag, 26. Juli 2024 um 16:18:24 UTC+2:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello, Can anyone please give me their opinion/advice? I am trying to run with TIMECON 50 (it was 1000 before), but no improvements yet.<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, July 23, 2024 at 7:23:56 AM UTC-4 Michela Benazzi 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>Hello everyone,</div><div><br>I am running DFT at 1200 K, 1 bar with a periodic system of 64 Al atoms + 2 CO2 atoms. However, I am seeing my temperature slowly increase to >4000 K and my pressure immediately spiking to ca. 2*10^5 bar.</div><div><br></div><div>I was wondering if you could help me diagnose this issue - I thought perhaps it was because I was not accounting for vdW interactions, but I just added that term to my newer simulations and I am still incurring in the same problem!</div><div><br></div><div>I included below my input and shortened output.<br><br></div></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/b0964eac-152a-4837-9154-08a5ac66e9cbn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b0964eac-152a-4837-9154-08a5ac66e9cbn%40googlegroups.com</a>.<br />