<div dir="ltr">Dear cp2k group, <div><br></div><div>I am trying to perform some cell optimization using cp2k interfaced with ASE. The command lines I am using are the following:</div><div><br></div><div><font face="monospace">sys=bulk('Cu', 'fcc', a=3.7, cubic=True)<br>sys.calc = calc</font></div><div><font face="monospace">opt=BFGS(sys, trajectory='CuBulk.traj', logfile='CuBulk.log', restart='Cubulk.rst')</font></div><div><font face="monospace">opt.run(fmax=0.01)<br></font><br></div><div>This is based on this ASE tutorial: <a href="https://wiki.fysik.dtu.dk/ase/gettingstarted/surface.html" target="_blank">https://wiki.fysik.dtu.dk/ase/gettingstarted/surface.html</a>.</div><div><br></div><div>However, as I have seen, the maximum force on each atom oscillates about the same range 0.6-0.8 eV/A (for a Cu unit cell), and it also depends on the metal I am working with. To be honest, I chose this value for fmax because I thought it was the convergence criteria for the maximum force between two consecutive optimization steps. My question is: How do I choose the maximum value for the force on each atom, or any other convergence criteria? Is this the right way to perform cell optimization using ASE?</div><div><br></div><div>Yours faithfully,</div><div><br></div><div>Lucas Bandeira</div><div><br></div></div>