Thank you for reply! I created two Si (100) systems with alpha beta gamma 90 90 90 and 90 90 120. The first system did't show any sign of convergence within 10 hours of simulation, whereas other fully converged (with 120 gamma degrees) at DZVP level within 19 hours.  The system with 90 90 120 cell parameters performed about 66 steps  within 19 hours, but the other 90 90 90 system performed only 2 steps for 10 hours. Before doing simulation, I created slab in avogadro and put in CP2K later. I wasn't sure should I choose 90 90 90 or 90 90 120 degree, as I saw smb published a work for Si (100) with 90 90 120 degrees cell parameters.<div>Here is my and ouput file for system with 90 90 120 cell angels. <br><div><br></div><div>  Total Charge                                                            0.002</div><div> !-----------------------------------------------------------------------------!</div><div><br></div><div> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.):               -6.901440296115577</div><div><br></div><div><br></div><div> --------  Informations at step =    66 ------------</div><div>  Optimization Method        =                LBFGS</div><div>  Total Energy               =        -6.9014402961</div><div>  Real energy change         =        -0.0000072318</div><div>  Decrease in energy         =                  YES</div><div>  Used time                  =              883.680</div><div><br></div><div>  Convergence check :</div><div>  Max. step size             =         0.0025436647</div><div>  Conv. limit for step size  =         0.0030000000</div><div>  Convergence in step size   =                  YES</div><div>  RMS step size              =         0.0006204331</div><div>  Conv. limit for RMS step   =         0.0015000000</div><div>  Convergence in RMS step    =                  YES</div><div>  Max. gradient              =         0.0001874079</div><div>  Conv. limit for gradients  =         0.0004500000</div><div>  Conv. in gradients         =                  YES</div><div>  RMS gradient               =         0.0000472944</div><div>  Conv. limit for RMS grad.  =         0.0003000000</div><div>  Conv. in RMS gradients     =                  YES</div><div> ---------------------------------------------------</div><div><br></div><div> *******************************************************************************</div><div> ***                    GEOMETRY OPTIMIZATION COMPLETED                      ***</div><div> *******************************************************************************</div></div>