<div><br /></div>Hi ..<div><br /></div><div>Maybe it simply needs to run for more iterations to converge.</div><div>With the information you provide it is hard to guess.</div><div>Is the electronic structure calculation of the bulk working fine?</div><div>Can you reproduce with your settings (BS, PP, XC etc) the known bulk properties?</div><div><br /></div><div>Regards</div><div>Marcella</div><div><br /></div><div><br /></div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, October 16, 2024 at 7:40:32 PM UTC+2 lheid...@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;"><p>Hello CP2K community,</p><p>I am running molecular dynamics (MD) simulations on a <strong>Cu2ZnSnS4 (CZTS)</strong> system using <strong>DFT</strong> in CP2K. Below is a detailed description of my system and the modifications I applied, followed by the SCF convergence issue I am facing.</p><p>Cu2ZnSnS4 (CZTS) system modeled in a periodic box.<br>Unit cell dimensions: 10.8 × 10.8 × 10.8 Å.<br>The goal is to study surface interactions and electronic properties with a vacuum layer.<br><br>A 20 Å vacuum layer was added in the Z direction to simulate surface effects: 10.8 × 10.8 × 30.8 Å<br><br>A layer of sulfur (S) atoms was added to the surface to stabilize the system and account for surface states.<br><br>I attempted hydrogen passivation by capping the dangling bonds with H atoms to further stabilize the surface.<br>I tried running the SCF loop with and without hydrogen passivation, but both cases failed to converge.</p><p>SCF Settings and Methods Tried:<br><b>Orbital Transformation (OT):</b><br>MINIMIZER: DIIS<br>PRECONDITIONER: FULL_SINGLE_INVERSE<br>ENERGY_GAP: 0.001<br>N_HISTORY_VEC: 7<br><b>Diagonalization:</b><br>I disabled the OT section and enabled diagonalization as a fallback method, but the SCF still did not converge. ( I tried different parameters setting)<br>SCF Parameters:<br>SCF_GUESS: ATOMIC<br>EPS_SCF: 1.0E-6<br>MAX_SCF: 100</p><p><b>The SCF loop exits after a few minutes, failing to converge under both OT and diagonalization methods.</b><br>Are there specific SCF settings or preconditioners that can improve convergence for systems with large vacuum gaps?<br>Are there alternative strategies for handling surfaces and vacuum layers that could make the system more stable for electronic structure calculations?<br>Has anyone successfully applied hydrogen passivation to stabilize surfaces and improve SCF convergence in CP2K?<br></p><p>Any suggestions or advice would be greatly appreciated!</p><p>Thank you for your help and support.</p><p><br></p></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/67fea74e-5fc0-4a47-8739-d8ae68961729n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/67fea74e-5fc0-4a47-8739-d8ae68961729n%40googlegroups.com</a>.<br />