Hi ..<div><br /></div><div>The warning for the CO calculation is not a real problem. Since the number of orbital functions, and consequently the size of  H, is 62, it is not possible to generate more than 62 MOs, including occupied and virtual.  The advise provided by the warning message is to set ADDED_MOS -1 to request the maximum number of available MOs. On the other hand, the code already makes this choice and the calculation runs to the end.</div><div><br /></div><div>The ill conditioned S matrix in the calculation of  the slab is due to the choice of the basis set, which is not really suitable for condensed matter calculations. To successfully run such systems,  basis sets optimised to reduce linear dependencies, as the MOLOPT series (but you can find others in the literatures), need to be adopted. </div><div><br /></div><div>Regards</div><div>Marcella</div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, November 19, 2025 at 9:38:22 AM UTC+1 scarlet...@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>Dear CP2K developers,</p>
<p>I am using CP2K with PBE-D3(BJ) to calculate single-point energies and try to make an output of AO overlap, density matrices, Mulliken and Löwdin populations for CO and CO adsorbed on a Cu slab.</p>
<p>I am facing two problems:</p>
<ol>
<li>
<p><strong>MO warning for CO calculations:</strong><br>
When using the def2-TZVP basis set with smearing, CP2K reports:</p>
<div><div><div><div></div></div></div><div dir="ltr">*** WARNING in qs_environment.F:1486 :: More added MOs requested than ***
*** available. The full set of unoccupied MOs will be used. Use         ***
*** 'ADDED_MOS -1' to always use all available MOs and to get rid of this ***
*** warning. 
</div></div>
<p>This occurs even though I explicitly requested 30 MOs in the input.</p>
</li>
<li>
<p><strong>Cholesky decomposition failure for Cu slab:</strong><br>
When using def2-TZVP for Cu (with or without smearing and CO), the calculation aborts with:</p>
<div><div><div><div></div></div></div><div dir="ltr">[ABORT] Cholesky decompose failed: matrix is not positive definite or ill-conditioned
fm/cp_cfm_basic_linalg.F:877
</div></div>
</li>
</ol>
<p>I would like to ask for guidance on how to avoid these issues. Should I adjust the number of MOs, EPS_DEFAULT, smearing, or other SCF parameters for metal systems with large basis sets? I have no idea about this.</p>
<p>I have attached the input XYZ files and example output files for reference.</p>
<p>Thank you very much for your time and help.</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 visit <a href="https://groups.google.com/d/msgid/cp2k/ab4f1c83-d180-4d28-b343-892c7019e7e2n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ab4f1c83-d180-4d28-b343-892c7019e7e2n%40googlegroups.com</a>.<br />