<div dir="ltr">Dear all, <div><br></div><div>I am trying to run a Cell Optimisation on a periodic molecular crystal (beta-resorcinol: C6H6O2, 4 molecules in unit cell) using K-points and BFGS as optimiser. I am using the vdW-DF-optPBE functional (<a href="https://github.com/cp2k/cp2k/commit/74ef05b96a46523a0daf7da04e75804420b5cd5e">https://github.com/cp2k/cp2k/commit/74ef05b96a46523a0daf7da04e75804420b5cd5e</a>). After 83 optimisation steps, my system does converge to the required MOTION thresholds, however from step 7 onwards the total energy increases at every step i.e.:</div><div><br></div><div><div>--------  Informations at step =     7 ------------</div><div>  Optimization Method        =                 BFGS</div><div>  Total Energy               =      -280.1489665596</div><div>  Internal Pressure [bar]    =     -2667.4573178334</div><div>  Real energy change         =         0.0001037473</div><div>  Predicted change in energy =        -0.0001063956</div><div>  Scaling factor             =         0.0000000000</div><div>  Step size                  =         0.0124910496</div><div>  Trust radius               =         0.4724315332</div><div><b style="background-color: rgb(255, 255, 0);">  Decrease in energy         =                   NO</b></div><div>  Used time                  =               60.978</div><div><br></div><div>  Convergence check :</div><div>  Max. step size             =         0.0124910496</div><div>  Conv. limit for step size  =         0.0010000000</div><div>  Convergence in step size   =                   NO</div><div>  RMS step size              =         0.0029192108</div><div>  Conv. limit for RMS step   =         0.0010000000</div><div>  Convergence in RMS step    =                   NO</div><div>  Max. gradient              =         0.0034738629</div><div>  Conv. limit for gradients  =         0.0001000000</div><div>  Conv. for gradients        =                   NO</div><div>  RMS gradient               =         0.0006857738</div><div>  Conv. limit for RMS grad.  =         0.0001000000</div><div>  Conv. for gradients        =                   NO</div><div>  Pressure Deviation [bar]   =     -2668.4705678334</div><div>  Pressure Tolerance [bar]   =       200.0000000000</div><div>  Conv. for  PRESSURE        =                   NO</div><div> ---------------------------------------------------</div></div><div><br></div><div>The resulting cell-optimised structure is therefore not a minimum in energy. I've attached the input and output files. Some points which may (or may not) be relevant:</div><div><ul><li>I did not encounter this problem when using a different functional (PBE + D3 instead of optPBE-vdW as used here), but otherwise identical settings.</li><li>I encountered the same issue (energy increasing at every step) using CUTOFF/REL_CUTOFF = 800/60</li><li>I initially tried to converge CUTOFF and REL_CUTOFF as in the tutorial (<a href="https://www.cp2k.org/howto:converging_cutoff">https://www.cp2k.org/howto:converging_cutoff</a>), however my single point calculations did not converge quite as nicely as in the tutorial - instead oscillating with an amplitude of order ~ e-4 Ha (I've attached a figure showing this). This is why I use such large CUTOFF/REL_CUTOFF values. I note that this was also the case when I ran PBE + D3 calculations, however in the latter case my cell optimisation did work. </li></ul><div><br></div><div>Any input would be greatly appreciated. </div></div><div><br></div><div>Best wishes, </div><div><br></div><div>Jan</div><div><br></div><div><br></div></div>