[CP2K-user] [CP2K:22065] Energy Difference Between OT and Diagonalization in CP2K

[Jürg Hutter]

Hi

Please note that you forgot to specify the vdW correction (D3) in the k-point
calculations.
Make sure that the vdW energy is also converged in both calculations (cutoff!).

best
JH

________________________________________
From: cp2k at googlegroups.com <cp2k at googlegroups.com> on behalf of yis... at 163.com <yisichi at 163.com>
Sent: Monday, January 26, 2026 4:23 AM
To: cp2k
Subject: [CP2K:22065] Energy Difference Between OT and Diagonalization in CP2K

To the CP2K Development Team,

Hello!

      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?
        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.

[微信图片_20260126112239_276_77.png]
          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.
         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?
          For reference, I have attached the relevant test results and input/output files:
         cp2k_K444.inp and cp2k_K444.out: input and output files for the diagonalization calculation with a 5×5×5 k-point grid.
        cp2k_SP334.inp and cp2k_SP334.out: input and output files for the OT method calculation with a 3×3×4 supercell.
        nio2.cif: the unit cell file used in the calculations.

A figure summarizing the energy convergence with respect to supercell size and k-point sampling.

I welcome any feedback or suggestions from the developers. Thank you for your time and assistance!

Best regards,
yisichi

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com<mailto:cp2k+unsubscribe at googlegroups.com>.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/3e06c02d-456d-4634-aaf4-e4284e2a7d1fn%40googlegroups.com?utm_medium=email&utm_source=footer>.

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB07596F1492A953E4D15C9B269F93A%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ????_20260126112239_276_77.png
Type: image/png
Size: 22543 bytes
Desc: ????_20260126112239_276_77.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20260126/e6e70053/attachment-0001.png>