<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=Windows-1252"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"> <div class="WordSection1"> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US">Hi Léon<o:p></o:p></span></p> <p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">You can export cell and coordinates as a POSCAR file in VESTA and from that file you can easily cut&paste CP2K’s &CELL and &COORD section.<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p> <div id="mail-editor-reference-message-container"> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">cp2k@googlegroups.com <cp2k@googlegroups.com> on behalf of Léon Luntadila Lufungula <Leon.luntadilalufungula@uantwerpen.be><br> <b>Date: </b>Wednesday, 28 June 2023 at 16:10<br> <b>To: </b>cp2k <cp2k@googlegroups.com><br> <b>Subject: </b>Re: [CP2K:19020] Re: TiO2 slab not converging<o:p></o:p></span></p> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Hi Matthias,<o:p></o:p></span></p> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Thanks a lot for the tip! The output when using the pre-compiled binary shows that some atoms are indeed too close together, which is probably the cause of the bad convergence. Is there an easy way to fix this? <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">I made the slabs with ASE so I'm quite surprised that it did not make the slabs correctly. Here is the script I used to make them with ASE, perhaps you can see where I made a mistake, resulting in the faulty CIF files...<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><br> from ase.io import read, write<br> from ase.visualize import view<br> from ase.build.surfaces_with_termination import surfaces_with_termination<br> from ase.build.tools import sort<br> <br> ana = read('anatase_opt_tzvp.pdb')<br> <br> for N in range(2,13):<br> s = surfaces_with_termination(ana, (1,0,1), N, vacuum=20, termination='O')<br> s = s[1]<br> s = s.repeat((1,3,1))<br> s = sort(s, tags=s.positions[:, 2])<br> name = f'a101-{N}-layer'<br> write(f'{name}.cif', s)<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt"> <span style="font-size:11.0pt"><br> Léon<o:p></o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">On Wednesday, 28 June 2023 at 15:58:25 UTC+2 Krack Matthias wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Hi Léon </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">You are using CP2K version 8.2. Newer CP2K versions perform a distance check for CP2K/Quickstep by default (see also this <a href="https://manual.cp2k.org/cp2k-2023_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/PRINT/INTERATOMIC_DISTANCES.html#CHECK_INTERATOMIC_DISTANCES" target="_blank"> input keyword</a>) when the system is smaller than 2000 atoms and prints a warning if a small distance is found. </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">You can download a pre-compiled binary of CP2K 2023.1 <a href="https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp" target="_blank"> here</a>, which is sufficient for that initial check.</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">HTH</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:36.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> </div> <div> <div> <div> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:72.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><span class="MsoHyperlink">cp...@googlegroups.com</span> <<span class="MsoHyperlink">cp...@googlegroups.com</span>> on behalf of Léon Luntadila Lufungula <<span class="MsoHyperlink">Leon.luntad...@uantwerpen.be</span>><br> <b>Date: </b>Wednesday, 28 June 2023 at 15:46<br> <b>To: </b>cp2k <<span class="MsoHyperlink">cp...@googlegroups.com</span>><br> <b>Subject: </b>Re: [CP2K:19017] Re: TiO2 slab not converging</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p style="margin-left:72.0pt"><span lang="EN-US">Hi Matthias,</span><o:p></o:p></p> <p style="margin-left:72.0pt"><span lang="EN-US"> </span><o:p></o:p></p> <p style="margin-left:72.0pt"><span lang="EN-US">Thanks for you reply!</span><o:p></o:p></p> <p style="margin-left:72.0pt"> <o:p></o:p></p> <p style="margin-left:72.0pt"><span lang="EN-US">I looked at my CIF file in various visualization software and it looked okay there. Is there a way to check this via CP2K to be sure?</span><o:p></o:p></p> <p style="margin-left:72.0pt"><span lang="EN-US">I also have some preliminary results of some ongoing calculations with DIIS/FULL_ALL and DIIS/FULL_SINGLE_INVERSE and these seem to be a bit faster than CG but reaching convergence is still a big issue (see attachments).<br> <br> Best,<br> Léon</span><o:p></o:p></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">On Wednesday, 28 June 2023 at 15:35:32 UTC+2 Krack Matthias wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="DE-CH" style="font-size:11.0pt">Hi Léon</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="DE-CH" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt">Did you check the set of atomic coordinates read from the cif file for close contacts possibly caused by the periodic boundary conditions?</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt">Best</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><o:p></o:p></span></p> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> </div> <div> <div> <div> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:108.0pt"> <b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><span class="MsoHyperlink">cp...@googlegroups.com</span> <<span class="MsoHyperlink">cp...@googlegroups.com</span>> on behalf of Léon Luntadila Lufungula <<span class="MsoHyperlink">Leon.luntad...@uantwerpen.be</span>><br> <b>Date: </b>Wednesday, 28 June 2023 at 14:16<br> <b>To: </b>cp2k <<span class="MsoHyperlink">cp...@googlegroups.com</span>><br> <b>Subject: </b>[CP2K:19014] Re: TiO2 slab not converging</span><span style="font-size:11.0pt"><o:p></o:p></span></p> </div> <p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt;margin-left:108.0pt"> <span style="font-size:11.0pt">I forgot to add the cif file which is read with the input, so here it is if someone needs it.<o:p></o:p></span></p> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">On Wednesday, 28 June 2023 at 13:54:37 UTC+2 Léon Luntadila Lufungula wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Update:<br> <br> I tried to optimize a 3-layer slab with altered SCF cycle settings (inner loop: 15 steps, outer loop: 45 steps) and still it does not converge. The convergence goes towards 2E-7 after 22 outer SCF cycles, but then doesn't seem to get lower (see attachment a101-3-layer-jobID_188.out). I tried to restart the calculation with altered EPS_SCF (= 1E-6) and then the restarted SCF cycle converges as the convergence was already below 1E-6. Unfortunately, after that the SCF again seems to only very slowly converge, after 27 outer SCF steps it's still only at 8E-5 convergence (see attachment a101-3-layer-jobID_189.out). Can anybody please help me out? I feel like the calculation is not converging properly, as after nearly 2 hours it still hasn't converged one SCF cycle but I don't know what might be wrong orhow to fix it... Requiring 300+ steps to reach SCF convergence seems quite high to me and makes me start to think there is something fundamentally wrong with my input. If anyone could look at my inputs and tell me what's wrong with them, I would be forever thankful!<br> <br> Kind regards,<o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Léon <br> <br> P.S. Could the TZVP basis set pose a problem? I often see people use DZVP-SR for similar calculations, but my promotor and co-promotor both told me that TZVP is required for proper energetics and should be used in production runs.<o:p></o:p></span></p> </div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">On Wednesday, 28 June 2023 at 10:04:07 UTC+2 Léon Luntadila Lufungula wrote:<o:p></o:p></span></p> </div> <blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Dear all,<o:p></o:p></span></p> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">I am trying to determine the ideal slab thickness for an anatase (101) slab but my SCF cycles are not converging. I optimized the cutoffs for the anatase bulk and did a CELL_OPT calculation on the unit cell to get the converged structure. From this optimized unit cell, I then used ASE to build the (101) slabs with a 20 Å vacuum width and thicknesses varying from 2 to 12 layers. During the GEO_OPT calculation I keep the bottom layer frozen to act as the bulk. I started the 2-layer calculation yesterday but when I looked at the output today, I see that the SCF cycle fails at every step of the geometry optimization (see output in attachment). Is there something wrong with my setup (too tight convergence criteria, wrong SCF setup, ...), or is it because a 2-layer structure is perhaps too small? The changes in energy go to 1E-5 - 1E-6 but never reach the convergence criteria of 1E-7 even after 300 steps (inner cycle: 10 steps, outer cycle 30 steps). I would guess that it could converge if I would increase the maximum number of steps, but 300 steps already seems quite a lot to me, although I could be wrong...<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">If anyone could share their thoughts on the situation so I can continue my work, I would greatly appreciate the help! Thanks in advance!<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt"> <o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Best regards,<o:p></o:p></span></p> </div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">Léon<o:p></o:p></span></p> </div> </blockquote> </div> </blockquote> </div> </div> </div> </div> </div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:108.0pt"> <span style="font-size:11.0pt">-- <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 <span class="MsoHyperlink">cp2k+uns...@googlegroups.com</span>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> </div> </div> </blockquote> </div> </div> </div> </div> </div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">-- <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 <span class="MsoHyperlink">cp2k+uns...@googlegroups.com</span>.<o:p></o:p></span></p> </div> </div> </div> </div> <div> <div> <div> <div> <p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;margin-left:72.0pt"> <span style="font-size:11.0pt">To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank"> https://groups.google.com/d/msgid/cp2k/94c9ad9e-58b3-4af1-844e-34b82c2fbed6n%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> </div> </div> </blockquote> </div> <p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">-- <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/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/2fb6de53-e6de-44ce-84d4-34ebbe1fb14bn%40googlegroups.com</a>.<o:p></o:p></span></p> </div> </div> </div> </body> </html> <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/ZRAP278MB08276A36D2C6A22184667559F424A%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ZRAP278MB08276A36D2C6A22184667559F424A%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<br />