Hi all,
When I optimized 1-butene molecular at B3lYP/6-311G** level, it took much much longer than gaussian16 did. I had test optimizer CG, BDGS and LBFGD, nothing better. It seems the optimizers have some efficiency issues. All three optimizer had taken at least 150 optimization steps to get the final structure, but the initial structure had been optimized by gaussian. How could I improve that?
Best
Linfeng
Thank you for your advice, Travis.
I had tested EPS_SCF from 1e-5 to 1e-14. I posted the 1e-5 one only because it took the shortest time to get the final structure. And I will try convergence criteria with 1e-4 as you suggest.
Linfeng
Hi,I see you have EPS_SCF commented out in your input file. The default convergence threshold is 1e-5 au which is usually not tight enough for accurate forces (and certainly not for vibrational analysis, which is computed numerically - hence it will take ages over Gaussian). The value you commented out is actually better and may take less steps to converge. Your criteria for convergence are also extremely tight, 1e-4 au is typically enough. Also, just to be clear, CP2K when setup appropriately will reproduce Gaussian output for the total energies for non-periodic calculations. So you CAN use Gaussian for the gas phase reaction and CP2K for the surface reaction.-T
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/e60c3ce3-6df6-48f1-8026-460bfc1a7798%40googlegroups.com.
<opt_geo.out><opt_geo.inp>
Thomas
To unsubscribe from this group and stop receiving emails from it, send an email to cp...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/e60c3ce3-6df6-48f1-8026-460bfc1a7798%40googlegroups.com.
<opt_geo.out><opt_geo.inp>