Just wanted to say that you were right about the dipole correction. You need to perform the computation first without the dipole correction, use that optimized geometry then with dipole correction included. Just wanted to add this in case someone has any problems with it too.<div><br>Thanks Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, March 3, 2021 at 5:37:50 PM UTC+2 Matt W 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;">Wavelet should just work in my experience.<div><br></div><div>For the dipole correction - In 'other codes' you had to (approximately) converge the electronic stucture without correction and then turn it on. Just starting a fresh calculation tended to diverge. As I say, not tried in CP2K.<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, March 3, 2021 at 5:08:55 AM UTC ASSIDUO Network 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 dir="ltr">I will look into those things, maybe I will get lucky.<br><br>In terms of dipole correction, it's been a nightmare for me to work. I have posted a few questions in this forum about it.</div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 2, 2021 at 6:37 PM Matt W <<a rel="nofollow">mat...@gmail.com</a>> wrote:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I normally use WAVELET poisson solver - then you need Y as the non-periodic direction (so PERIODIC XZ), the system must be in the centre of the cell and you need ~5 A vacuum each side to get converged results.<div><br></div><div>The ANALYTIC psolver is also usable for 2D poisson, but it requires a much bigger cell to converge (~2x), to my knowledge. </div><div><br></div><div>Do check that you are converged with cell size if you go these routes. There is also a dipole correction method available somewhere (I've not used).</div><div><br></div><div>Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, March 2, 2021 at 4:04:03 PM UTC Lucas Lodeiro wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Yes, and in &CELL too, to avoid interaction between slab replicas and their electrostatic fields.</div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 2 mar 2021 a las 12:07, Lenard Carroll (<<a rel="nofollow">len...@gmail.com</a>>) escribió:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Just to be 100% clear, are you saying that for the &POISSON section which deals with the electrostatic potential, that I select PERIODIC XY?</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 2, 2021 at 4:55 PM Lucas Lodeiro <<a rel="nofollow">el...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Just a comment, see your electrostatic potential profile. You are using a periodical slab picture, then when you add a polar molecule, there is not a flat potential into the vacuum, in that case, disable the periodicity in this direction or add an electrostatic field to correct it.<div>Also, to check the convergence of your quantities with respect to your model, you could explore bigger slabs than 6 layers.</div><div><br></div><div>Regards - Lucas Lodeiro</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 2 mar 2021 a las 6:37, ASSIDUO Network (<<a rel="nofollow">len...@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thank you for the quick response.<div><br></div><div>I understand what you are saying and it makes total sense to me. Thank you for clarifying.<br><br>The 3x3x1 was just to test how my final adsorption energy changes when I change various settings, including alpha, beta, kpoints and cutoff. I got similar results to the the 5x5x1 case, however the 5x5x1 result was more accurate compared to the ltierature.<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, March 2, 2021 at 11:33:40 AM UTC+2 <a rel="nofollow">fa...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Dear Lenard,</div><div><br></div><div>Yes, you have correctly answered your question. For metals you need a quite good sampling of the brillouin zone. Therefore a gamma-point calculation of a 5x5 angstrom unit cell is expected to yield unreasonable results. <br></div><div>The fact that you managed to get a good structure with a cutoff of 300 is probably due to a coincidental error cancelling. But as soon as you alter the configuration by adding a CO the error cancelling does not work anymore and the structure becomes unreasonable again.<br></div><div><br></div><div>You should make sure that your quantity of interest is converged w.r.t the size of the k-point grid. 3x3x1still seems small to me, but maybe its already enough for you.</div><div><br></div><div>Cheers,</div><div>Fabian<br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, 2 March 2021 at 06:42:45 UTC+1 ASSIDUO Network wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi there everyone, I am currently trying to geometry optimize a stack of Au (111) 2x2 surfaces, specifically by fixing the bottom 3 surfaces and relaxing the top 3.<div><br></div><div>Unfortunately, the top surfaces completely mess up and the final structure looks bad. When I turn KPOINTS on (3x3x1 or 5x5x1), it no longer happens.</div><div>Finally, when I lower the cutoff to 300, but I leave KPOINTS off, I get a good optimized geometry, but when I try to perform a CO adsorption with that optimized structure, it severely deforms as well.</div><div><br></div><div>Does anyone have a solution to my problem, particularly one that doesn't change the ELECTRONIC_TEMPERATURE, KPOINTS or CUTOFF. Is it the need for KPOINTS in a small structure like this that messes up the results?<br><br>I have attached all my relevant files.<br><br>Kind regards,<br>Lenard</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 rel="nofollow">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d1100c7c-514e-4fbc-82e5-528bcdad45abn%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/d1100c7c-514e-4fbc-82e5-528bcdad45abn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1615880664929000&usg=AFQjCNHjCp-hWx2PHiMOz4r9L4uestwyAg">https://groups.google.com/d/msgid/cp2k/d1100c7c-514e-4fbc-82e5-528bcdad45abn%40googlegroups.com</a>.<br> </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 rel="nofollow">cp...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PJSHNzMu9F08E0UO-rBJZJO7puF6Z9UMUxQDxzu10L_bA%40mail.gmail.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/CAOFT4PJSHNzMu9F08E0UO-rBJZJO7puF6Z9UMUxQDxzu10L_bA%2540mail.gmail.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1615880664929000&usg=AFQjCNFu7rKlOnqZHNsg2foxfVSZGZnFsQ">https://groups.google.com/d/msgid/cp2k/CAOFT4PJSHNzMu9F08E0UO-rBJZJO7puF6Z9UMUxQDxzu10L_bA%40mail.gmail.com</a>.<br> </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 rel="nofollow">cp...@googlegroups.com</a>.<br></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAKur158Q1dvprpLA5QyzBVSMy1qJd14oxfFPUaFygaKF2TcYBA%40mail.gmail.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/CAKur158Q1dvprpLA5QyzBVSMy1qJd14oxfFPUaFygaKF2TcYBA%2540mail.gmail.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1615880664929000&usg=AFQjCNHUW_M5YSpNXeHhNyhBq2vl2qZr0g">https://groups.google.com/d/msgid/cp2k/CAKur158Q1dvprpLA5QyzBVSMy1qJd14oxfFPUaFygaKF2TcYBA%40mail.gmail.com</a>.<br> </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 rel="nofollow">cp...@googlegroups.com</a>.<br></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/31940cf9-268f-443d-a5ff-17f5759e0969n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/31940cf9-268f-443d-a5ff-17f5759e0969n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1615880664929000&usg=AFQjCNF1MunBxatkxscFP2e3de0Al32rdQ">https://groups.google.com/d/msgid/cp2k/31940cf9-268f-443d-a5ff-17f5759e0969n%40googlegroups.com</a>.<br> </blockquote></div> </blockquote></div></blockquote></div>