Thank your very much<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2026年1月26日星期一 UTC+8 17:21:55<Jürg Hutter> 写道:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>Please note that you forgot to specify the vdW correction (D3) in the k-point
<br>calculations.
<br>Make sure that the vdW energy is also converged in both calculations (cutoff!).
<br>
<br>best
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of <a href data-email-masked rel="nofollow">yis...@163.com</a> <<a href data-email-masked rel="nofollow">yis...@163.com</a>>
<br>Sent: Monday, January 26, 2026 4:23 AM
<br>To: cp2k
<br>Subject: [CP2K:22065] Energy Difference Between OT and Diagonalization in CP2K
<br>
<br>To the CP2K Development Team,
<br>
<br>Hello!
<br>
<br> I would like to raise another question: why is there a discrepancy in the calculated energies when using the OT (Orbital Transformation) method compared to the diagonalization method?
<br> I started with the NiO₂ unit cell (5.6512 Å × 5.6512 Å × 4.3886 Å). Due to its small size, I expanded it into supercells and performed calculations using the OT method. As shown in the attached table, I constructed supercells of sizes 2×2×2, 3×3×3, 3×3×4, and 4×4×4. The energies per NiO₂ unit cell obtained were –805.2913786 Ha, –805.2931813 Ha, –805.2932437 Ha, and –805.2931291 Ha, respectively, with relative energy differences of –1.803×10⁻³ Ha, –6.232×10⁻⁵ Ha, and 1.145×10⁻⁴ Ha. The results indicate that the energy converges to –805.2932437 Ha per unit cell when the supercell is expanded to 3×3×4.
<br>
<br>[微信图片_20260126112239_276_77.png]
<br> Additionally, I tested the diagonalization method with k-point sampling. Using k-point grids of 4×4×4, 5×5×5, 6×6×6, and 7×7×7, the energies per NiO₂ unit cell were –805.2394556 Ha, –805.2394471 Ha, –805.2395249 Ha, and –805.2394426 Ha, respectively, with relative differences of 8.512×10⁻⁶ Ha, –7.782×10⁻⁵ Ha, and 8.234×10⁻⁵ Ha. The energy converges to –805.2394556 Ha per unit cell with a 4×4×4 k-point grid.
<br> Although both methods appear to converge individually, there is a significant discrepancy of about 0.0538 Ha (~1.46 eV) between the two results, which cannot be overlooked. Could you please help explain the possible reasons for this difference?
<br> For reference, I have attached the relevant test results and input/output files:
<br> cp2k_K444.inp and cp2k_K444.out: input and output files for the diagonalization calculation with a 5×5×5 k-point grid.
<br> cp2k_SP334.inp and cp2k_SP334.out: input and output files for the OT method calculation with a 3×3×4 supercell.
<br> nio2.cif: the unit cell file used in the calculations.
<br>
<br>A figure summarizing the energy convergence with respect to supercell size and k-point sampling.
<br>
<br>I welcome any feedback or suggestions from the developers. Thank you for your time and assistance!
<br>
<br>Best regards,
<br>yisichi
<br>
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%2540googlegroups.com&source=gmail&ust=1769510436023000&usg=AOvVaw2FX9Dk82XRp__pX2sQHN8Z">https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=zh-CN&q=https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1769510436023000&usg=AOvVaw1F3xnyWfbzKap5fbteD7EV">https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></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/4519118d-e2f0-41a7-8c22-c75c0b5b8339n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/4519118d-e2f0-41a7-8c22-c75c0b5b8339n%40googlegroups.com</a>.<br />