Hi Marcella,<div><br></div><div>Thanks for responding! There is a giant zoo of methods people seem to use for a system like this, from PBE to relatively exotic stuff and everyone states that their method is "optimal." That said, I do agree that I should have probably started with something that is known to reproduce the electronic structure of graphene rather than what works for small organics. :)</div><div><br></div><div>Again, thanks.</div><div><br></div><div>Alex<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, August 29, 2022 at 6:10:57 AM UTC-6 Marcella Iannuzzi 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><br></div><div><br></div>Well<div>Let's say that you should have a good reason to use all-electrons and hybrid-functionals with no ADMM for a BOMD simulation of a condensed matter system. </div><div>Moreover one should consider that graphene is semi-metallic and the convergence might be strongly dependent on the system size as well as other parameters/conditions. </div><div>In the literature you can find many DFT studies of water/graphene interfaces, I would start from there to figure out which are the optimal settings to correctly reproduce electronic, structural and dynamic properties. </div><div>Regards</div><div>Marcella</div><div><br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, August 29, 2022 at 9:37:05 AM UTC+2 Alex 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">Hi everyone,<div><br></div><div>I am trying set up a dftmd simulation (room-temp NVT, rectangular box periodic in XYZ) of a tiny solvated sheet of graphene (~1700 electrons total). I am using the following as my starting point for the force calculation: <a href="https://github.com/misteliy/cp2k/blob/master/tests/QS/regtest-hybrid/H2O-hybrid-b3lyp.inp" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://github.com/misteliy/cp2k/blob/master/tests/QS/regtest-hybrid/H2O-hybrid-b3lyp.inp&source=gmail&ust=1661896760194000&usg=AOvVaw12Ccjhh2ftDAaW08pY1lcJ">https://github.com/misteliy/cp2k/blob/master/tests/QS/regtest-hybrid/H2O-hybrid-b3lyp.inp</a> -- the only two major changes being POISSON set to be periodic in xyz and I'm also adding D3 corrections. </div><div>The initial structure comes from classical MD and SCF is refusing to converge with that dreaded huge oscillation. Should I abandon this idea and use something much more familiar to me (say, PBE + GTH potentials) or am I setting it up incorrectly? </div><div>Thank you,</div><div><br></div><div>Alex</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/346a9665-9286-43c1-8d45-3abc2d06ddfen%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/346a9665-9286-43c1-8d45-3abc2d06ddfen%40googlegroups.com</a>.<br />