Hi Chelsea<div><br></div><div>I'm not sure about metal oxides, but the manual does mention that OT converges poorly with metals. I also had to avoid OT for pure Au.<br><br></div><div>Do you or anyone else know whether CELL_OPT is appropriate for non-periodic systems (e.g., molecules, nanoparticles).</div><div><br></div><div><div>Kind Regards</div><div><br></div><div>Sam</div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, May 19, 2022 at 2:37:22 PM UTC+2 uri...@gmail.com 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;">Dear  Marcella<div><br></div><div>Thank you very much for your help and the scf-coverge problem for cel-opt has been solved. The results are listed as follows:</div><div>-DIAG would perfectly solve the problem for cel-opt yet the following geo-opt still fluctuated without converge under same parameter.</div><div>- tighter EPS_DEFAULT didn't work with OT but would benefit for converge with DIAG method.</div><div>- higher PW-cutoff and changing the U_MINUS_J value both failed to converge the scf cycle with OT method in my cel-opt case</div><div><br></div><div>Thank you again for your help!</div><div><br></div><div>Regards</div><div>Chelsea<br></div><div><br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年5月19日星期四 UTC+8 20:20:08<Faye Guo> 写道:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>Dear all<div><br></div><div>The problem has been solved, the Diagonalisation rather than OT would converge the scf cycle and CEL-OPT.</div><div><br></div><div>However, a new one occurred in the following GEO-OPT which could hardly converge with an increasing total force and DFT+U energy. I've attached the input and output and I would really appreciate the suggestions.</div><div><div><p>Thanks in advance. </p></div><div>Bests</div><div><br></div><div>Chelsea </div></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年5月10日星期二 UTC+8 10:23:59<Faye Guo> 写道:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear 

Marcella<div><br></div><div>Thank you very much for your help. I've watched some of your lectures on youtube which are of great help to me as a very beginner in cp2k and computaional method.</div><div><br></div><div>As to my case:</div><div>1. Smaller U-J value from 5-6.8 still failed to converge the scf cycle</div><div>2. I'm confused about whether DFT+U method is right for Cu2O optimization as the magnetic moment is only about 0.012uB( data from material project). And I do read some papers discussing Cu2O modeling without +U information </div><div>3. I've tried the DIIS minimizer method which did converge the scf cycle but brought about an abnormal GEO-optimization result. I guess this might related to that DIIS method is likely to give an unstable solution</div><div>I'll try the other options to see if there's anyone could make it.</div><div><br></div><div>I've uploaded the input and output files for both DIIS and CG methods and I really appreciated your help.</div><div><br></div><div>Bests</div><div>Chelsea</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2022年5月9日星期一 UTC+8 16:38:08<Marcella Iannuzzi> 写道:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Chelsea<div><br></div><div>Some things that you can try are:</div><div>-use diagonalisation instead of OT</div><div>- use tighter EPS_DEFAULT</div><div>- higher PW-cutoff</div><div>- change the U_MINUS_J value</div><div><br></div><div>Without having the output it is hard to guess what is going on. </div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Sunday, May 8, 2022 at 5:16:12 PM UTC+2 <a rel="nofollow">uri...@gmail.com</a> 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">Dear All, <div><br></div><div>I've been trying to perform a CELL_OPT for  Cu2O  but the calculation sadly cannot converge the SCF cycles. </div><div><br></div><div> I tried all possible OT preconditioners with CG as minimizer by cp2k 6.1 edition which all failed and spin contamination even occurred under full_all condition.  </div><div><br></div><div>I also tried the increase of max_scf from 100 to 400 as well as the energy gap from 1*10-4 to 1*10-2, yet neither of them worked.</div><div><br></div><div>So I'm wondering if any else had these issues and found ways to overcome them? And the input file is attached as follows.</div><div><br></div><div><p>Any help would be much appreciated. Thanks in advance. </p></div><div>Bests</div><div><br></div><div>Chelsea </div></blockquote></div></blockquote></div></blockquote></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/ef62c6f3-b728-4737-bac5-afeea55d42e9n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ef62c6f3-b728-4737-bac5-afeea55d42e9n%40googlegroups.com</a>.<br />