Cholesky decomposition failed

24 views
Skip to first unread message

Dat Ho

unread,
Sep 8, 2025, 5:22:16 AM (6 days ago) Sep 8
to cp2k

Dear all,


I am trying to perform ab initio MD simulations on a plain water box containing 512 water molecules using the PW6B95 functional. The same input works flawlessly for smaller boxes containing 64, 128 or 265 molecules, but I can't seem to get it to work for 512 molecules. I always get the error message 'Cholesky decomposition failed'. Can anyone spot an error in the input? Also, what does this error message mean in detail? 


Any help would be much appreciated!

H2O_512.log
H2O-512-1.restart
H2O_512.inp

Johann Pototschnig

unread,
Sep 8, 2025, 8:04:00 AM (6 days ago) Sep 8
to cp2k
Cholesky decomposition computes a smaller representation of a tensor, which allows more efficient computations. 

It assumes a positive definite matrix, which, off course, numerically can not be guaranteed. 

It can be switched off, but certain functionalities will not work:

Most OT preconditioners use it. 

Johann Pototschnig

unread,
Sep 8, 2025, 8:06:49 AM (6 days ago) Sep 8
to cp2k
This problem is often caused by diffuse basis functions, so it might help to change:
BASIS_SET AUX_FIT cFIT3
Reply all
Reply to author
Forward
0 new messages