[CP2K-user] [CP2K:19014] Re: TiO2 slab not converging

[Léon Luntadila Lufungula]

I forgot to add the cif file which is read with the input, so here it is if 
someone needs it.

On Wednesday, 28 June 2023 at 13:54:37 UTC+2 Léon Luntadila Lufungula wrote:

> Update:
>
> 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!
>
> Kind regards,
> Léon 
>
> 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.
> On Wednesday, 28 June 2023 at 10:04:07 UTC+2 Léon Luntadila Lufungula 
> wrote:
>
>> Dear all,
>>
>> 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...
>>
>> If anyone could share their thoughts on the situation so I can continue 
>> my work, I would greatly appreciate the help! Thanks in advance!
>>
>> Best regards,
>> Léon
>>
>

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/12b3aec6-38d8-404f-b68f-aa681ccd6266n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230628/6ed20446/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: a101-3-layer.cif
Type: chemical/x-cif
Size: 6522 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230628/6ed20446/attachment.cif>