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>

<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/8c2d7967-c051-4157-a531-b9bdfb83ff27n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/8c2d7967-c051-4157-a531-b9bdfb83ff27n%40googlegroups.com</a>.<br />