Hi,<div><br></div><div>I'm a new user of cp2k 8.1. I'm calculating a P doped Si system (output file attached). The system contains 1000 atoms. The calculation goes well and the SCF converges fairly smoothly. However, different from my previous calculations, I found in this calculation there are many SCF cycles (SCF converged again and again) for each geometry optimization step. In my previous calculations, the basic steps are: SCF ->  OPTIMIZATION STEP 1 -> SCF -> OPTIMIZATION STEP 2 -> SCF -> OPTIMIZATION STEP 3. But in this calculation, the steps are: SCF -> OPTIMIZATION STEP 1 -> SCF (converged) -> SCF (converged) -> SCF (converged) -> SCF (converged) ...... -> OPTIMIZATION STEP 2. I suspect this maybe caused by the different OPTIMISER. In my previous calculations, I used BFGS. But in this calculation I used LBFGS. However, based on the tutorial, LBFGS should be faster in large system calculation but the results are not what I expected. Could someone please provide some suggestions on improving this issue? Thanks</div><div><br></div><div>Regards,</div><div>Hongyang</div>