<div>Well, CP2K was setup to do Born-Oppenheimer MD for bio systems - so for wide gap insulators you should for sure only need a few SCF cycles to converge.</div><div><br></div><div>Is your system insulating (band gap, let's say > 0.2 eV)? If not different approach needed.<br></div><div><br></div><div>If so, improve the OT setup - especially the preconditioner<br></div><div><br></div><div>remove<br></div><div>PRECONDITIONER FULL_KINETIC<br>ENERGY_GAP 0.01<br><br>try <br></div><div><br></div><div>PRECONDITIONER FULL_SINGLE_INVERSE</div><div><br></div><div>or (can be expensive for larger systems)<br></div><div><br></div><div>PRECONDITIONER FULL_ALL</div><div><br></div><div>and leave out the energy_gap (if you are using an up-to-date cp2k)</div><div><br></div><div>Also check that you XC setup is correct ...</div><div><br></div><div>Matt<br></div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, 16 November 2021 at 17:58:59 UTC nwi...@berkeley.edu 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;">Hello all, <div><br></div><div>When I run a relaxation for a complex system, e.g. in my case a defective transition metal oxide with hybrid dft on a modest ~200 atom supercell, each SCF still takes quite a long time to converge. </div><div><br></div><div>I originally came from vasp before switching to cp2k, and it in that code, my experience was that once the first scf loop converged, the next relaxation step would take much less time, and so on until the last few steps converged in just a few iterations. This isn't my experience with CP2K though. Running CG on my most recent relaxation I still find that after more than 20 relaxation steps, 40-50 scf steps are required to achieve convergence.</div><div><br></div><div>Is this something that I could improve by changing certain settings? Below is the relevant info for my input file. Please note that I understand FULL_KINETIC and CG might lead to slower convergence than FULL_ALL + DIIS, but my question lies in why convergence is not appreciably speeding up rather than its absolute speed.</div><div><br></div><div>-Nick</div><div><br></div><div><p>&DFT</p>
<p> BASIS_SET_FILE_NAME BASIS_MOLOPT</p>
<p> BASIS_SET_FILE_NAME BASIS_MOLOPT_UCL</p>
<p> BASIS_SET_FILE_NAME BASIS_ADMM</p>
<p> BASIS_SET_FILE_NAME BASIS_ADMM_MOLOPT</p>
<p> POTENTIAL_FILE_NAME GTH_POTENTIALS</p>
<p> UKS T</p>
<p> MULTIPLICITY 0</p>
<p> CHARGE 4</p>
<p> &SCF</p>
<p> MAX_ITER_LUMO 400</p>
<p> MAX_SCF 50</p>
<p> EPS_SCF 1e-06</p>
<p> SCF_GUESS RESTART</p>
<p> &OT</p>
<p> ALGORITHM STRICT</p>
<p> MINIMIZER CG</p>
<p> LINESEARCH 2PNT</p>
<p> PRECONDITIONER FULL_KINETIC</p>
<p> ENERGY_GAP 0.01</p>
<p> ROTATION F</p>
<p> OCCUPATION_PRECONDITIONER F</p>
<p> &END OT</p>
<p> &OUTER_SCF</p>
<p> EPS_SCF 1e-06</p>
<p> MAX_SCF 20</p>
<p> &END OUTER_SCF</p>
<p> &END SCF</p>
<p> &AUXILIARY_DENSITY_MATRIX_METHOD</p>
<p> ADMM_PURIFICATION_METHOD None</p>
<p> METHOD BASIS_PROJECTION</p>
<p> &END AUXILIARY_DENSITY_MATRIX_METHOD</p>
<p> &QS</p>
<p> EPS_DEFAULT 1e-12</p>
<p> EPS_PGF_ORB 1e-16</p>
<p> EXTRAPOLATION PS</p>
<p> METHOD GPW</p>
<p> &END QS</p>
<p> &MGRID</p>
<p> NGRIDS 5</p>
<p> CUTOFF 550.0</p>
<p> REL_CUTOFF 50.0<br></p>
<p> &END MGRID</p><p> &XC</p><p> &XC_FUNCTIONAL NO_SHORTCUT</p><p> &PBE</p><p> SCALE_X 0.0<br></p><p> SCALE_C 1.0</p><p> &END PBE</p><p> &XWPBE</p><p> SCALE_X 0.063</p><p> SCALE_X0 0.877</p><p> OMEGA 0.07</p><p> &END XWPBE</p><p> &PBE_HOLE_T_C_LR</p><p> SCALE_X 0.123</p><p> CUTOFF_RADIUS 5.930000000000001</p><p> &END PBE_HOLE_T_C_LR</p><p> &END XC_FUNCTIONAL</p><p> &HF</p><p> FRACTION 1.0</p><p> &SCREENING</p><p> EPS_SCHWARZ 1e-07</p><p> EPS_SCHWARZ_FORCES 1e-06</p><p> SCREEN_P_FORCES T</p><p> SCREEN_ON_INITIAL_P T</p><p> &END SCREENING</p><p> &INTERACTION_POTENTIAL</p><p> POTENTIAL_TYPE MIX_CL_TRUNC</p><p> OMEGA 0.07</p><p> SCALE_COULOMB 0.06</p><p> SCALE_LONGRANGE 0.063</p><p> CUTOFF_RADIUS 5.930000000000001</p><p> T_C_G_DATA t_c_g.dat</p><p> &END INTERACTION_POTENTIAL</p><p> &LOAD_BALANCE</p><p> RANDOMIZE T</p><p> &END LOAD_BALANCE</p><p> &MEMORY</p><p> EPS_STORAGE_SCALING 0.1</p><p> MAX_MEMORY 2000</p><p> &END MEMORY</p><p> &END HF</p><p> &END XC</p><p>&END DFT</p></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/0cea4d69-5edf-4c48-9ac4-e9e8c9c497b9n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/0cea4d69-5edf-4c48-9ac4-e9e8c9c497b9n%40googlegroups.com</a>.<br />