<div dir="ltr"><div>Dear Torstein,</div><div><br></div><div>On Test 2:</div><div>I am no expert on LS_SCF but my understanding ist that PURIFICATION_METHOD SIGN should only be used if the chemical potential is know and set using MU. This is useful for AIMD where MU is known and does not fluctuate too much. Switching PURIFICATION_METHOD to TRS4 and enabling DYNAMIC_THRESHOLD will improve convergence if you do not know MU beforehand. I assume that the reason for the long walltime/step for LS_SCF is because of poor convergence originating in the purification.</div><div><br></div><div>On Test 3:</div><div>As OT and LS_SCF are numerical methods they will not result in identical densies, similar but not the same. So you can expect the forces to be slightely different too. If I run your input with tighter settings EPS_SCF 1.0e-8 and EPS_FILTER 1.0e-8 the max. force at step 126 (starting val 123) is already below the convergence threshold:<br></div><div><br></div><div> --------  Informations at step =   126 ------------<br>  Optimization Method        =                 BFGS<br>  Total Energy               =     -3498.3787056930<br>  Real energy change         =        -0.0000126132<br>  Predicted change in energy =        -0.0000107671<br>  Scaling factor             =         0.0000000000<br>  Step size                  =         0.0038980052<br>  Trust radius               =         0.4724315332<br>  Decrease in energy         =                  YES<br>  Used time                  =               49.665<br><br>  Convergence check :<br>  Max. step size             =         0.0038980052<br>  Conv. limit for step size  =         0.0030000000<br>  Convergence in step size   =                   NO<br>  RMS step size              =         0.0004691139<br>  Conv. limit for RMS step   =         0.0015000000<br>  Convergence in RMS step    =                  YES<br>  Max. gradient              =         0.0004187196<br>  Conv. limit for gradients  =         0.0004500000<br>  Conv. in gradients         =                  YES<br>  RMS gradient               =         0.0000483787<br>  Conv. limit for RMS grad.  =         0.0003000000<br>  Conv. in RMS gradients     =                  YES<br> ---------------------------------------------------</div><div><br></div><div>Best,</div><div>Fabian<br></div></div>