<p>Dear CP2K Developers,</p>
<p>I am writing to ask for your advice regarding some small but systematic discrepancies I observed in density matrices and Mulliken population analysis obtained from a single‐point calculation performed using CP2K.</p>
<p>In my test for CO, I compared the following files generated from the same single-point SCF calculation:</p>
<ol>
<li>
<p><strong>The MOLDEN file</strong> (CO_molopt_smearing-mos.molden-1_0.molden)</p>
</li>
<li>
<p><strong>The Mulliken file</strong> (CO_molopt_smearing-mulliken-1.mulliken)</p>
</li>
<li>
<p><strong>The largest iteration AO log file</strong> (e.g., CO_molopt_smearing-ao-1_0_12.Log)</p>
</li>
<li>
<p><strong>The 1_0 AO log file</strong> (CO_molopt_smearing-ao-1_0.Log)</p>
</li>
</ol>
<p>I found that:</p>
<ul>
<li>
<p>The density matrices read from the last AO log file and from the 1_0 AO file differ by approximately <strong>1×10⁻⁷</strong>, even though they originate from the same single-point calculation.</p>
</li>
<li>
<p>Constructing Mulliken populations by <strong>element-wise multiplication</strong> of the density matrix <span aria-hidden="true">P</span> and the overlap matrix <span aria-hidden="true">S</span> (i.e., computing <span aria-hidden="true">PμνSμν</span> for each pair of atomic orbitals) leads to slightly different values compared to the <strong>MULLIKEN NET ATOMIC ORBITAL AND OVERLAP POPULATION MATRIX</strong> reported in CP2K’s Mulliken output file.</p>
</li>
<li>
<p>Additionally, when processing the MOLDEN file using <strong>Multiwfn</strong>, the computed Mulliken bond orders are again slightly different from those obtained directly from CP2K.</p>
</li>
</ul>
<p>Before I proceed further, I would like to confirm:</p>
<ol>
<li>
<p>Are these ~1e−7 differences expected due to formatting, truncation, numerical precision, or different internal representations used when writing the log files, MOLDEN file, and Mulliken file?</p>
</li>
<li>
<p>Is there a recommended “most reliable” source of the density matrix for reproducing CP2K Mulliken results (Mulliken file vs. AO log file vs. MOLDEN file)?</p>
</li>
<li>
<p>Are the Mulliken matrices produced by CP2K post‐processed in any additional way (e.g., symmetrization, thresholding) that could explain these differences?</p>
</li>
<li>
<p>Are there known limitations or differences in the basis function ordering between MOLDEN export and the AO log files?</p>
</li>
</ol>
<p>For reference, I attach (or can provide) the following files:</p>
<div><div><div><div></div></div></div><div dir="ltr">CO_molopt_smearing.inp<br />
CO_molopt_smearing.out<br />CO_molopt_smearing-mulliken-1.mulliken
<br />CO_molopt_smearing-mos.molden-1_0.molden
<br />CO_molopt_smearing-ao-1_0.Log
<br />CO_molopt_smearing-ao-1_0_12.Log
</div></div>
<p>Any guidance on how to interpret these discrepancies would be greatly kind of you. Thank you very much for your time and for your continuous development of CP2K.</p>
<p>Best regards</p>
<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/b263b616-bf70-463b-a9e5-7ad8561fa1e7n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b263b616-bf70-463b-a9e5-7ad8561fa1e7n%40googlegroups.com</a>.<br />