<div dir="ltr"><div>Hi Luis,</div><div><br></div><div>I seem to remember there is no analytical stress tensor implemented for meta-GGAs in CP2K because they use laplacian of the density, see:</div><div><br></div><div><a href="https://groups.google.com/d/msg/cp2k/V5l4_iUORkY/WuQOQGrwAQAJ">https://groups.google.com/d/msg/cp2k/V5l4_iUORkY/WuQOQGrwAQAJ</a></div><div><br></div><div>In addition, the stress tensor from MD and single point energy calculations will be different, see earlier post by Prof Hutter:</div><div><br></div><div><a href="https://groups.google.com/d/msg/cp2k/rifbK2sdjJ4/O3nUipT0M8EJ">https://groups.google.com/d/msg/cp2k/rifbK2sdjJ4/O3nUipT0M8EJ</a></div><div><br></div><div>If you decide to run the numerical stress tensor calculations, please keep in mind that you will need to use tighter convergence criteria (e.g. plane-wave cutoff, energy, force).</div><div><br></div><div>I believe you will be able to find more useful information if you search in the mailing list.</div><div><br></div><div>SL</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 14 October 2015 at 21:54, Luis Ruiz Pestana <span dir="ltr"><<a href="mailto:lar...@lbl.gov" target="_blank">lar...@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks Matt.<div><br></div><div>It would definitely be useful to implement a way in cp2k to calculate stresses when using metaGGA's (either through the analytical calculation or indirectly through the REFTRAJ ensemble using the numerical stress calculation). Let me know if you get it done any time soon, it would be a life saver!</div><div><br></div><div>I really like cp2k and I highly recommend it to other users but in the current state I might be forced to change software to simulate systems with mGGA's.</div><div><br></div><div>Best,</div><div>Luis</div></div><div class="gmail_extra"><br><div class="gmail_quote"><span>On Wed, Oct 14, 2015 at 1:35 PM, Matt W <span dir="ltr"><<a href="mailto:MattWa...@gmail.com" target="_blank">MattWa...@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><span><div dir="ltr"><div>Hi Luis,</div><div><br></div>seems like this might not be doable at the moment then, maybe only energy-force. Sorry if the suggestion was bad. I'd think it would be easy to add the extra ability to calculate stress tensor too, I'll try and have a look, but not sure I'll have time for a bit.<div><span><font color="#888888"><div><br></div></font></span><div><span><font color="#888888">Matt </font></span><div><div><br><br>On Wednesday, October 14, 2015 at 5:27:04 PM UTC+1, <a href="mailto:lar...@lbl.gov" target="_blank">lar...@lbl.gov</a> wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div dir="ltr">Hi Matt,<div><br></div><div>I have been trying to use REFTRAJ to calculate the stress on some frames from a trajectory that I already ran but I don't seem to get it to work. The simulation in the REFTRAJ ensemble runs fine but the outputs are zero for the stress and other quantities in the *.ener file. I have attached here the input file as well as the outputs for energy and stress. I am not sure what I'm doing wrong so any pointer would be highly appreciated. Thanks in advance for your help!</div><div><br></div><div>Best,</div><div>Luis</div><div><br></div><div><br></div><div><br><br>On Sunday, September 13, 2015 at 5:53:08 AM UTC-7, Matt W wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div dir="ltr"><div>Hi,</div><div><br></div>the simplest way might be to post process - use the REFTRAJ ensemble in the MD to rerun over the frames you like. The overhead of having to redo the SCF from not quite as good a guess should be negligible.<div><br></div><div>Matt<br><br>On Friday, September 11, 2015 at 7:47:06 PM UTC+1, <a>lar...@lbl.gov</a> wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid"><div dir="ltr">Dear team,<div><br></div><div>I was wondering if there is any available setting in cp2k to set the frequency at which the stress is calculated in an AIMD simulation. I am currently using metaGGAs in my simulations and because the only way of calculating the pressure is numerically for this type of xc functionals the MD simulations are painfully slow. If there would be a way to not perform the numerical stress calculation every step that would make them manageable. Btw, just changing how often it's printed doesn't change anything. Thanks in advance!</div><div><br></div><div>Best,</div><div>Luis</div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></div></div></div></span><div><div> <p></p> -- <br> You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br> To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/ZHD9Y8dsH2E/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/ZHD9Y8dsH2E/unsubscribe</a>.<br> To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<span><br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br> </span></div></div></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><span style="color:rgb(136,136,136);font-size:12.8px">Luis A. Ruiz Pestana, PhD</span><div style="color:rgb(136,136,136);font-size:12.8px">Chemist Postdoctoral Fellow</div><div style="color:rgb(136,136,136);font-size:12.8px">Lawrence Berkeley National Laboratory & UC Berkeley</div><div style="color:rgb(136,136,136);font-size:12.8px">Berkeley, CA 94720</div><div style="color:rgb(136,136,136);font-size:12.8px">Site: <a style="color:rgb(17,85,204)" href="http://luisruiz2015.wix.com/personal" target="_blank">personal website</a></div></div></div> </font></span></div><div class="HOEnZb"><div class="h5"> <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+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br> Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br> </div></div></blockquote></div><br></div>