<div>Dear Lucas,</div><div><br></div><div>(I'm not a CP2K developer, but I have worked quite a bit with ASE.)<br></div><div><br></div><div>Your code snippet looks fine, except that this will simply relax the ionic positions.</div><div></div><div>If you also want to optimize the cell vectors, you will need to look at e.g. the StrainFilter</div><div>or UnitCellFilter in ase.constraints. An "fmax" value of 0.01 eV/Angstrom should</div><div>generally give well-converged minima. Fmax=0.05 eV/Angstrom is also sufficiently</div><div>accurate in many cases (but please check it yourself to be sure).<br></div><div><br></div><div>But the fact that the forces do not readily decrease in this calculation suggests</div><div>that something may be wrong with your CP2K input. But without a minimal working</div><div>example (including all input and output files you get when running the example),</div><div>it is impossible to diagnose that further.<br></div><div><br></div><div>Best,<br></div>Maxime<br><div class="gmail_quote"><div dir="auto" class="gmail_attr">Op woensdag 12 mei 2021 om 19:13:26 UTC+2 schreef ban...@gmail.com:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><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" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=nl&q=https://wiki.fysik.dtu.dk/ase/gettingstarted/surface.html&source=gmail&ust=1621324566993000&usg=AFQjCNF_0emU0HPC4HePdDz0R_BN7lIWiQ">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>
</blockquote></div>