<div dir="ltr"><span style="color: rgb(0, 0, 255);">Hi all,<br><br>I've been using the OT[CP2K(2.4)] eigensolver successfully for calculating energies for quite large systems (60-70 atoms) and now I need to find the optimized geometry using BFGS with OT (attached 3.inp). But I get the following error message:<br></span><br><br>  <br> *************************************************************************<br> *** 20:24:42 ERRORL2 in cp_dbcsr_cholesky:cp_dbcsr_cholesky_decompose ***<br> *** processor 0  :: err=-300 condition FAILED at line 123             ***<br> *************************************************************************<br><br><br> ===== Routine Calling Stack ===== <br><br>           13 cp_dbcsr_cholesky_decompose<br>           12 qs_ot_get_derivative<br>           11 ot_mini<br>           10 ot_scf_mini<br>            9 qs_scf_loop_do_ot<br>            8 scf_env_do_scf_inner_loop<br>            7 scf_env_do_scf<br>            6 qs_energies_scf<br>            5 qs_forces<br>            4 cp_eval_at<br>            3 geoopt_bfgs<br>            2 cp_geo_opt<br>            1 CP2K<br><br><br><span style="color: rgb(0, 0, 255);">However, when I constrain all atoms in the large system I do not get this error and the calculations go OK.  Changing the solver to Digonalization proceeds on successfully (attached 1.inp), but I would like to use OT throughout for its fast performance. I tested the combination OT+BFGS for the small system of H2O (attached 2.inp), with and without pinning the atoms, and that worked fine: so why doesn't work for the large system? Please help me figure out where the problem lies. <br><br>Thanks,<br>ER<br></span><br></div>