Benzene molecule †This tutorial explains how to optimize the geometry by using different algorithms. In this tutorial, a benzene molecule is placed in a cubic box with the cell edges of 10 Angstrom. pot_C.pbe1 and pot_H.lda1 are used with the cutoff energies of 36 Ry (GMAX=6) and 400 Ry (GMAXP=20) for the wave functions and the argmentation charge, respectively.
Monitoring the maximum force by using grep as (supposing the output file is "nfout_1") > grep -A1 f_max nfout_1 we obtain NIT TotalEnergy f_max f_rms edel vdel fdel 1 -40.14215406 0.002809 0.001947 0.10D-08 0.12D-06 0.10D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 2 -40.14215395 0.005734 0.004041 0.22D-08 0.17D-06 0.22D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 3 -40.14224266 0.001880 0.001393 0.47D-08 0.13D-06 0.47D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 4 -40.14225368 0.000871 0.000552 0.11D-08 0.15D-06 0.11D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 5 -40.14225523 0.000046 0.000030 0.51D-09 0.85D-07 0.51D-09 and we can see that the geometry optimization converges within 5 steps. Then, let us use the quenched molecular dynamics (QMD) for comparison. In this case we use IMDALG=2 and the time step (DTIO) of 200 atomic unit as 200.00 2 1 0.10D-08 : DTIO IMDALG IEXPL EDELTA and the resulting maximum forces are NIT TotalEnergy f_max f_rms edel vdel fdel 1 -40.14215406 0.002809 0.001947 0.10D-08 0.12D-06 0.10D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 2 -40.14218382 0.001971 0.001448 0.22D-08 0.10D-06 0.22D-08 -- ... -- NIT TotalEnergy f_max f_rms edel vdel fdel 14 -40.14225511 0.000195 0.000135 0.15D-08 0.66D-07 0.15D-08 -- NIT TotalEnergy f_max f_rms edel vdel fdel 15 -40.14225523 0.000071 0.000033 0.21D-09 0.23D-07 0.21D-09 The convergence is achieved with 15 steps and we can see that the GDIIS method is much more efficient than QMD. |