 This thread is locked.Only browsing is available.
 Top Page > Browsing geometry relaxation Date: 2007/07/04 22:08 Name: JessK Dr. Ozaki, I am trying to relax a surface. I put a force criteria 4.0d-4 Ha/bohr, and use DIIS method for that. There is no SCF convergence problem, but relaxation is very very slow - it passed about 300 steps, but forces oscillate around one value. Are there any parameters which could help? |Maximum force| (Hartree/Bohr) = 0.004734148175 |Maximum force| (Hartree/Bohr) = 0.004922389624 |Maximum force| (Hartree/Bohr) = 0.005641408054 |Maximum force| (Hartree/Bohr) = 0.006429973134 |Maximum force| (Hartree/Bohr) = 0.008364579772 |Maximum force| (Hartree/Bohr) = 0.003883459762 |Maximum force| (Hartree/Bohr) = 0.003494295620 |Maximum force| (Hartree/Bohr) = 0.003565106358 |Maximum force| (Hartree/Bohr) = 0.003808776831 |Maximum force| (Hartree/Bohr) = 0.003518174225 |Maximum force| (Hartree/Bohr) = 0.003190680359 |Maximum force| (Hartree/Bohr) = 0.003280411383 |Maximum force| (Hartree/Bohr) = 0.006639887540 |Maximum force| (Hartree/Bohr) = 0.005232858642 |Maximum force| (Hartree/Bohr) = 0.004440083714 |Maximum force| (Hartree/Bohr) = 0.004114006691 |Maximum force| (Hartree/Bohr) = 0.004327090157 |Maximum force| (Hartree/Bohr) = 0.005095942923 |Maximum force| (Hartree/Bohr) = 0.006376893247 |Maximum force| (Hartree/Bohr) = 0.008481134569 |Maximum force| (Hartree/Bohr) = 0.011454567615 |Maximum force| (Hartree/Bohr) = 0.015366560529 |Maximum force| (Hartree/Bohr) = 0.016115039728 |Maximum force| (Hartree/Bohr) = 0.016434052761 |Maximum force| (Hartree/Bohr) = 0.014963018847 |Maximum force| (Hartree/Bohr) = 0.010052082220 |Maximum force| (Hartree/Bohr) = 0.007140709704 |Maximum force| (Hartree/Bohr) = 0.006686724870 |Maximum force| (Hartree/Bohr) = 0.005499107286 |Maximum force| (Hartree/Bohr) = 0.005175443078 |Maximum force| (Hartree/Bohr) = 0.004016993191 |Maximum force| (Hartree/Bohr) = 0.002961278589 |Maximum force| (Hartree/Bohr) = 0.004001748498 |Maximum force| (Hartree/Bohr) = 0.003882077410 |Maximum force| (Hartree/Bohr) = 0.003893408575 |Maximum force| (Hartree/Bohr) = 0.004033210610 |Maximum force| (Hartree/Bohr) = 0.004114949478 |Maximum force| (Hartree/Bohr) = 0.004298120234 |Maximum force| (Hartree/Bohr) = 0.004485083905 |Maximum force| (Hartree/Bohr) = 0.004570423744 |Maximum force| (Hartree/Bohr) = 0.004824545345 |Maximum force| (Hartree/Bohr) = 0.004987265999 |Maximum force| (Hartree/Bohr) = 0.003842487255 |Maximum force| (Hartree/Bohr) = 0.004056208063 |Maximum force| (Hartree/Bohr) = 0.003679225273 |Maximum force| (Hartree/Bohr) = 0.002953778425 |Maximum force| (Hartree/Bohr) = 0.003514826512 |Maximum force| (Hartree/Bohr) = 0.003140333926 |Maximum force| (Hartree/Bohr) = 0.003137441515 |Maximum force| (Hartree/Bohr) = 0.003154661583 |Maximum force| (Hartree/Bohr) = 0.003237588895 |Maximum force| (Hartree/Bohr) = 0.003342009854 |Maximum force| (Hartree/Bohr) = 0.004031431833 |Maximum force| (Hartree/Bohr) = 0.004068472953 |Maximum force| (Hartree/Bohr) = 0.003923797293 |Maximum force| (Hartree/Bohr) = 0.003958764119 |Maximum force| (Hartree/Bohr) = 0.004035131394 |Maximum force| (Hartree/Bohr) = 0.004063828708 |Maximum force| (Hartree/Bohr) = 0.004056360143 |Maximum force| (Hartree/Bohr) = 0.004047979732 |Maximum force| (Hartree/Bohr) = 0.004345450222 |Maximum force| (Hartree/Bohr) = 0.004775394905 |Maximum force| (Hartree/Bohr) = 0.004753100400 |Maximum force| (Hartree/Bohr) = 0.004095718748 |Maximum force| (Hartree/Bohr) = 0.005284600282 |Maximum force| (Hartree/Bohr) = 0.004442768011 |Maximum force| (Hartree/Bohr) = 0.003670183641 |Maximum force| (Hartree/Bohr) = 0.003570650150 |Maximum force| (Hartree/Bohr) = 0.003445412072 |Maximum force| (Hartree/Bohr) = 0.003522990464 |Maximum force| (Hartree/Bohr) = 0.003538940523 |Maximum force| (Hartree/Bohr) = 0.004506627192 |Maximum force| (Hartree/Bohr) = 0.004551521577 |Maximum force| (Hartree/Bohr) = 0.004666136842 |Maximum force| (Hartree/Bohr) = 0.004783940535 |Maximum force| (Hartree/Bohr) = 0.004894237769 |Maximum force| (Hartree/Bohr) = 0.004874370042 |Maximum force| (Hartree/Bohr) = 0.004849225718 |Maximum force| (Hartree/Bohr) = 0.004977702338 |Maximum force| (Hartree/Bohr) = 0.004946515199 |Maximum force| (Hartree/Bohr) = 0.005077942377 |Maximum force| (Hartree/Bohr) = 0.005053019091 |Maximum force| (Hartree/Bohr) = 0.004256245025 |Maximum force| (Hartree/Bohr) = 0.004478513265 |Maximum force| (Hartree/Bohr) = 0.005547334818 |Maximum force| (Hartree/Bohr) = 0.003987146430 |Maximum force| (Hartree/Bohr) = 0.003475534544 |Maximum force| (Hartree/Bohr) = 0.003387796841 |Maximum force| (Hartree/Bohr) = 0.003519442917 |Maximum force| (Hartree/Bohr) = 0.003814156830 |Maximum force| (Hartree/Bohr) = 0.003440362284 |Maximum force| (Hartree/Bohr) = 0.003562940194 |Maximum force| (Hartree/Bohr) = 0.005014759405 |Maximum force| (Hartree/Bohr) = 0.004986902644 |Maximum force| (Hartree/Bohr) = 0.004990464445 |Maximum force| (Hartree/Bohr) = 0.004906274298 |Maximum force| (Hartree/Bohr) = 0.004944988498 |Maximum force| (Hartree/Bohr) = 0.004895685639 thanks, jk Page: Re: geometry relaxation ( No.1 ) Date: 2007/07/05 03:59 Name: Vasilii Artyukhov I think that this problem has already been discussed several times here. This maybe a feature of the algorithm... But! you could try and tell us what happens if you increase the integration grid cutoff, this parameter _can_ cause this kind of a problem, especially with GGA.  Re: geometry relaxation ( No.2 ) Date: 2007/07/06 02:43 Name: JessK Hi, I increased cutoff from 150 Ry to 320 Ry and it looks better now (forces slowly go down). But I wonder, why steepest descent (SD) method is used when DIIS methos is chosen. From manual, it said that first N steps are SD method, and next - DIIS. But from output I see that relaxation method switches from SD to DIIS, and so on. jk  Re: geometry relaxation ( No.3 ) Date: 2007/07/06 06:10 Name: Vasilii Artyukhov Perhaps you should see some literature on DIIS, it might be the normal behavior. It can still be a bug, though... Anyway, the important thing is that the problem is more or less solved - is it? Page: